Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,000,850 A→C Q169P (CAG→CCG)  xdhA → xanthine dehydrogenase, molybdenum binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,000,8500AC100.0% 266.4 / NA 68Q169P (CAG→CCG) xdhAxanthine dehydrogenase, molybdenum binding subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (35/33);  total (35/33)

AGTGGCAAGAATTACCCGTTATCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAAAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTCGCCGCGTGGTTGGTCAGGC  >  NC_000913/3000706‑3000997
                                                                                                                                                |                                                                                                                                                   
aGTGGCAAGAATTACCCGTTATCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCCGGTACa                                                                                                                                               <  1:1077137/151‑1 (MQ=255)
   ggCAAGAATTACCCGTTATCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCCGGTACAggg                                                                                                                                            >  2:694784/1‑151 (MQ=255)
      aaGAATTACCCGTTATCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCCGGTACAGGGGCa                                                                                                                                         <  2:397487/151‑1 (MQ=255)
       aGAATTACCCGTTATCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCCGGTACAGGGGCAc                                                                                                                                        >  1:605681/1‑151 (MQ=255)
        gAATTACCCGTTATCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCCGGTACAGGGGCACt                                                                                                                                       <  1:261868/151‑1 (MQ=255)
         aaTTACCCGTTATCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCCGGTACAGGGGCACt                                                                                                                                       <  2:755353/150‑1 (MQ=255)
         aaTTACCCGTTATCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCCGGTACAGGGGCACt                                                                                                                                       >  2:277360/1‑150 (MQ=255)
          aTTACCCGTTATCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCCGGTACAGGGGCACTAt                                                                                                                                     >  1:658496/1‑151 (MQ=255)
            tACCCGTTATCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCCGGTACAGGGGCACTATc                                                                                                                                    <  2:265827/150‑1 (MQ=255)
             aCCCGTTATCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCCGGTACAGGGGCACTATc                                                                                                                                    >  2:63878/1‑149 (MQ=255)
             aCCCGTTATCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCCGGTACAGGGGCACTATCa                                                                                                                                   <  2:95963/150‑1 (MQ=255)
                cGTTATCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCCGGTACAGGGGCACTATCAGAcc                                                                                                                               >  2:743916/1‑151 (MQ=255)
                   tATCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCCGGTACAGGGGCACTATCAGAccccc                                                                                                                            >  1:843843/1‑151 (MQ=255)
                   tATCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCCGGTACAGGGGCACTATCAGAccccc                                                                                                                            >  2:1019285/1‑151 (MQ=255)
                    aTCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCCGGTACAGGGGCACTATCAGAccccc                                                                                                                            <  2:302539/150‑1 (MQ=255)
                      caccacGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCCGGTACAGGGGCACTATCAGACCCCCGtt                                                                                                                         <  2:342094/151‑1 (MQ=255)
                      caccacGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCCGGTACAGGGGCACTATCAGACCCCCGt                                                                                                                          <  1:907264/150‑1 (MQ=255)
                         cacGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCCGGTACAGGGGCACTATCAGACCCCCGTTAtt                                                                                                                      <  1:106213/151‑1 (MQ=255)
                          acGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCCGGTACAGGGGCACTATCAGACCCCCGTTATTc                                                                                                                     >  2:592324/1‑151 (MQ=255)
                             ccAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCCGGTACAGGGGCACTATCAGACCCCCGTTATTCAAc                                                                                                                  <  1:816457/151‑1 (MQ=255)
                                gAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCCGGTACAGGGGCACTATCAGACCCCCGTTATTCAACAtt                                                                                                               >  2:116957/1‑151 (MQ=255)
                                 aaGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCCGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTg                                                                                                              >  2:702438/1‑151 (MQ=255)
                                  aGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCCGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGt                                                                                                             <  1:823472/151‑1 (MQ=255)
                                  aGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCCGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGt                                                                                                             >  1:412854/1‑151 (MQ=255)
                                   gcggcgCTGGCAGAAGACGCTGCACCAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCCGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTc                                                                                                            <  2:638257/151‑1 (MQ=255)
                                    cggcgCTGGCAGAAGACGCTGCACCAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCCGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTca                                                                                                           <  2:1036784/151‑1 (MQ=255)
                                     ggcgCTGGCAGAAGACGCTGCACCAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCCGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTcac                                                                                                          >  1:951513/1‑151 (MQ=255)
                                      gcgcTGGCAGAAGACGCTGCACCAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCCGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTcac                                                                                                          >  1:574147/1‑150 (MQ=255)
                                        gcTGGCAGAAGACGCTGCACCAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCCGGTACAGGGGCACTATCAGACCCGCGTTATTCAACATTGGCACATg                                                                                                       >  2:184176/1‑151 (MQ=255)
                                         cTGGCAGAAGACGCTGCACCAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCCGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATgg                                                                                                      >  1:53697/1‑151 (MQ=255)
                                            gCAGAATACGCTGCACCAAGCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCCGGTACAGAGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGaaa                                                                                                   <  1:111955/151‑1 (MQ=255)
                                            gCAGAAGACGCTGCACCAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCCGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGaaa                                                                                                   <  2:632087/151‑1 (MQ=255)
                                                  gaCGCTGCACCAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCCGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAAAGCGTa                                                                                              >  1:977997/1‑150 (MQ=255)
                                                     gCTGCACCAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCCGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAAAGCGTAACAt                                                                                          >  1:1104870/1‑151 (MQ=255)
                                                          aCCAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCCGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAAAGCGTAACATCGCt                                                                                      <  2:1117713/150‑1 (MQ=255)
                                                           ccAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCCGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAAAGCGTAACATCGCTgg                                                                                    >  1:816015/1‑151 (MQ=255)
                                                             aaTCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCCGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAAAGCGTAACATCGCTGGCg                                                                                  <  2:1061716/151‑1 (MQ=255)
                                                                   tAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCCGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAAAGCGTAACATCGCTGGCGtgga                                                                              >  1:659009/1‑149 (MQ=255)
                                                                    aaCGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCCGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAAAGCGTAACATCGCTGGCGtggatgg                                                                           >  1:724323/1‑151 (MQ=255)
                                                                    aaCGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATATCCAACAAACAATCGATGCCGCCGACTACCCGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAAAGCGTAACATCGCTGGCGtggatgg                                                                           <  1:32682/151‑1 (MQ=255)
                                                                        gTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCCGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAAAGCGTAACATCGCTGGCGTGGATggagg                                                                        >  1:19610/1‑150 (MQ=255)
                                                                          ggCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCCGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAAAGCGTAACATCGCTGGCGTGGATGGAGGAt                                                                      <  2:688755/150‑1 (MQ=255)
                                                                           gCAATTGACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGACGGCGCCGGCTACCCGGTACAGGGGCACTACCAGACCCCCGATATTCAACATTGTCACATGGAAAGCGTAACATCGCTGGCGTGGATGGGCGATg                                                                     >  2:211989/1‑150 (MQ=255)
                                                                                 tACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCCGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAAAGCGTAACATCGCTGGCGTGGATGGAGGATGACTcgcg                                                              >  2:369891/1‑151 (MQ=255)
                                                                                  aCTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCCGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAAAGCGTAACATCGCTGGCGTGGATGGAGGATGACTcgcg                                                              >  1:86658/1‑150 (MQ=255)
                                                                                    tGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCCGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAAAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAAt                                                           >  2:797770/1‑151 (MQ=255)
                                                                                     gAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCCGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAAAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAAt                                                           <  1:965710/150‑1 (MQ=255)
                                                                                               aCGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCCGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAAAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGttt                                                >  1:631727/1‑151 (MQ=255)
                                                                                               aCGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCCGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAAAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGttt                                                >  2:1101256/1‑151 (MQ=255)
                                                                                                cGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCCGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAAAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTc                                               <  1:997298/151‑1 (MQ=255)
                                                                                                          gggTAATGTCCAACAAACAATCGATGCCGCCGACTACCCGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAAAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAg                                     <  1:252079/151‑1 (MQ=255)
                                                                                                            gTAATGTCCAACAAACAATCGATGCCGCCGACTACCCGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAAAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGAt                                   <  2:1070030/151‑1 (MQ=255)
                                                                                                                 gTCCAACAAACAATCGATGCCGCCGACTACCCGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAAAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGc                              >  1:991557/1‑151 (MQ=255)
                                                                                                                    caacaaACAATCGATGCCGCCGACTACCCGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAAAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGcaca                           <  2:944632/151‑1 (MQ=255)
                                                                                                                    caacaaACAATCGATGCCGCCGACTACCCGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAAAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGcaca                           <  1:985596/151‑1 (MQ=255)
                                                                                                                     aacaaacaATCGATGCCGCCGACTTCCCGGTACAGGGGCGCTATCAGACCCCCGTTATTCAACATTGTCACATGGAAAGCGTAACATCGCTGGCGTGGATAGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACAt                          <  2:447049/151‑1 (MQ=255)
                                                                                                                        aaacaaTCGATGCCGCCGACTACCCGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAAAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGt                       >  1:995799/1‑151 (MQ=255)
                                                                                                                          acaaTCGATGCCGCCGACTACCCGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAAAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTc                     >  1:1066695/1‑151 (MQ=255)
                                                                                                                           caaTCGATGCCGCCGACTACCCGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAAAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTcg                    <  1:853770/151‑1 (MQ=255)
                                                                                                                            aaTCGATGCCGCCGACTACCCGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAAAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTcgc                   >  1:953477/1‑151 (MQ=255)
                                                                                                                             aTCGATGCCGCCGACTACCCGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAAAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTcgcc                  >  1:51581/1‑151 (MQ=255)
                                                                                                                               cGATGCCGCCGACTACCCGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAAAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTcgccg                 <  1:959428/150‑1 (MQ=255)
                                                                                                                                gATGCCGCCGACTACCCGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAAAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTCGCcgcg               >  1:72237/1‑151 (MQ=255)
                                                                                                                                gATGCCGCCGACTACCCGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAAAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTCGCcgcg               <  1:806083/151‑1 (MQ=255)
                                                                                                                                gATGCCGCCGACTACCCGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAAAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTCGCcgcg               <  1:923811/151‑1 (MQ=255)
                                                                                                                                          aCTACCCGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAAAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTCGCCGCGTGGTTGGTCa     <  1:888494/151‑1 (MQ=255)
                                                                                                                                             aCCCGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAAAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTCGCCGCGTGGTTGGTCAGgc  >  1:814666/1‑151 (MQ=255)
                                                                                                                                             aCCCGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAAAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTCGCCGCGTGGTTGGTCAGgc  <  2:225944/151‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
AGTGGCAAGAATTACCCGTTATCACCACGCCAGAAGCGGCGCTGGCAGAAGACGCTGCACCAATCCATAACGGTGGCAATTTACTGAAACAAAGCACGATGTCGACGGGTAATGTCCAACAAACAATCGATGCCGCCGACTACCAGGTACAGGGGCACTATCAGACCCCCGTTATTCAACATTGTCACATGGAAAGCGTAACATCGCTGGCGTGGATGGAGGATGACTCGCGAATTACCATCGTTTCCAGCACCCAGATCCCGCACATTGTTCGCCGCGTGGTTGGTCAGGC  >  NC_000913/3000706‑3000997

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: