Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,830,064 A→C 7.7% P234P (CCT→CCG) ‡ astA ← arginine succinyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,830,0640AC 7.7% 275.1 / 10.5 100P234P (CCT→CCG) ‡astAarginine succinyltransferase
Reads supporting (aligned to +/- strand):  ref base A (46/45);  new base C (5/3);  total (51/49)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

GCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCGGCTAAAATCCA  >  NC_000913/1829916‑1830209
                                                                                                                                                    |                                                                                                                                                 
gcgcAGGCTGCCCTTCTGCCACTTCCACCCGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGGCAAAGATGTCGATAGAGTTTCGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAg                                                                                                                                                  >  2:1202098/1‑150 (MQ=255)
   cAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGgcg                                                                                                                                               >  1:748300/1‑150 (MQ=255)
       acgCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCggg                                                                                                                                           >  2:1754129/3‑149 (MQ=255)
        ggCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCCTGCCCGCGCGGtg                                                                                                                                          <  1:1864856/149‑2 (MQ=255)
         gCCCTTCTGCCACTTCCACCCGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATCGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGttt                                                                                                                                         <  2:1380485/150‑1 (MQ=255)
         gCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGttt                                                                                                                                         <  1:248613/150‑1 (MQ=255)
          cccTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCCTGCCCGCGCGGtgct                                                                                                                                        <  2:32436/150‑4 (MQ=255)
           ccTTCTGCCACTTCCACCAGCCGACTTTTACGGATTGCGCGCACGAGGTCGATGTCACACTCAAGCGTCGGCCCACCGGCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCCGCACCGCGAGGGCAGGCGCGGTTTGc                                                                                                                                       <  2:1253201/150‑1 (MQ=255)
           ccTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCCCCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGc                                                                                                                                       >  2:466025/1‑150 (MQ=255)
            cTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCCCGCGGTCGCTGTCACACTCAAGCGTCGGCCCACCGTCGAAGATGTCGCTATAGTTACGGTAGCGAAAACGTGCTTTCTCCAGCACACCGCGGGCAGGCGCGGGTTGCg                                                                                                                                      >  1:1799409/1‑150 (MQ=255)
            cTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCCTGCCCGCGCGGtgctgg                                                                                                                                      <  2:1879688/150‑6 (MQ=255)
             ttCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCgg                                                                                                                                     <  2:1854648/150‑1 (MQ=255)
             ttCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCgg                                                                                                                                     <  1:1324828/150‑1 (MQ=255)
              tCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGa                                                                                                                                    <  2:737946/150‑1 (MQ=255)
                tGCCACTTCCAGCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTAGGGCCACCTTCAAAGATGTCGATATAGTTAAGGTAGCGAATACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATg                                                                                                                                  <  2:123745/150‑1 (MQ=255)
                ggCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCAGGCGGTCGATGTCACACTCAAGCGGCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCTAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATg                                                                                                                                  <  2:942676/149‑1 (MQ=255)
                 gCCACTTCCCCCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGt                                                                                                                                 <  1:860062/150‑1 (MQ=255)
                 gCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTAGGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGt                                                                                                                                 >  1:956782/1‑150 (MQ=255)
                 gCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCGCCGCGCGGGCAGGCGCGGTTTGCGGATGt                                                                                                                                 >  2:362338/1‑150 (MQ=255)
                 gCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGt                                                                                                                                 >  1:61871/1‑150 (MQ=255)
                  ccACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGAGCAGGCGCGGTTTGCGGATGTa                                                                                                                                >  2:1630812/1‑150 (MQ=255)
                    aCTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTAcc                                                                                                                              <  2:1741135/150‑1 (MQ=255)
                    aCTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGGCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTAcc                                                                                                                              <  2:1526032/150‑1 (MQ=255)
                       tCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACAGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTGTCTCCAGCACCGCGCGGGCAGGCGCGGGTTGCGGATGTACCTGa                                                                                                                           >  1:1540383/1‑150 (MQ=255)
                        ccaccaGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGAc                                                                                                                          <  2:65992/150‑1 (MQ=255)
                        ccaccaGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTGGCGGGTGTACCGGAc                                                                                                                          >  1:154599/1‑150 (MQ=255)
                                     tttACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTg                                                                                                             <  2:1116036/150‑1 (MQ=255)
                                      ttACGGATGGCGAGCAAGCGGTCGTTGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGTTGACGTCCTgg                                                                                                            <  2:1670751/150‑1 (MQ=255)
                                         cGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCATTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCt                                                                                                         >  2:1888646/1‑150 (MQ=255)
                                            aTGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGAGGTACCTGACCGATGACGTCCTGGGCTTcc                                                                                                      >  2:133380/1‑150 (MQ=255)
                                             tGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTTCCTGCCCGATGACGTCCTGGGCTTCCt                                                                                                     >  2:1033860/1‑150 (MQ=255)
                                               gcgcgcACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTgg                                                                                                   <  2:1461256/150‑1 (MQ=255)
                                                    cACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGCCGTCCTGGGCTTCCTGGGATaa                                                                                              >  1:709126/1‑150 (MQ=255)
                                                        cgGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGAGGTACCTGACCGATGACGGCCTGGGCTTCCTGGGATTAAAAg                                                                                          >  1:245425/1‑150 (MQ=255)
                                                           tCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTgg                                                                                       <  2:863848/150‑1 (MQ=255)
                                                           gcGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTAGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTgg                                                                                       <  2:1460885/149‑1 (MQ=255)
                                                            cgctGTCACACTCAAGCGTCGGCCCACCGTCAAAGAAGTCGCTATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTggg                                                                                      <  2:754958/147‑1 (MQ=255)
                                                            cgctGTCACACGCAAGCGTCGGCCCACCGTCAAAGATGGCGATATAGGTTCGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGAGGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTggg                                                                                      <  1:791986/147‑1 (MQ=255)
                                                            cGATGTCACACTCAAGCGTCGGCCCACCTTCAAAGATGTCGATATAGATACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTggg                                                                                      >  2:1888749/1‑150 (MQ=255)
                                                            cGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTggg                                                                                      <  2:449534/150‑1 (MQ=255)
                                                             gaagTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGt                                                                                     <  1:1798681/147‑1 (MQ=255)
                                                             gATCTCACACTCAAGCGTCGGCCCACCGTCAAAGAAGTCGATATAGTAACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGt                                                                                     <  1:154045/150‑1 (MQ=255)
                                                               tGTCACACTCATGCGTCGGCCCACCGTCAATGATGTCGATATAGTTAAGGTAGCGAAAACCTTCTTTCTCCAGCTCCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTTGGCTTCCTGGGATAAAAAGTGGGtat                                                                                   <  2:1204267/150‑1 (MQ=255)
                                                                gTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGtata                                                                                  <  1:915615/150‑1 (MQ=255)
                                                                 tCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAACCCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGGTTGCGGATGTACCTGAACGATGCAGTCCTCGGCTTCCTGGGATAAAAAGGGGGTATAg                                                                                 >  2:361085/1‑150 (MQ=255)
                                                                 tCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATTAAAAGTGGGTATAg                                                                                 >  2:251738/1‑150 (MQ=255)
                                                                 tCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAg                                                                                 >  2:1721682/1‑150 (MQ=255)
                                                                 tCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAg                                                                                 >  1:21859/1‑150 (MQ=255)
                                                                      cTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCgg                                                                            <  2:713415/150‑1 (MQ=255)
                                                                      cTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCgg                                                                            >  1:438720/1‑150 (MQ=255)
                                                                      cTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCgg                                                                            <  1:235119/150‑1 (MQ=255)
                                                                       tCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGa                                                                           <  1:1355107/150‑1 (MQ=255)
                                                                        cAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTGTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGAt                                                                          >  1:978995/1‑150 (MQ=255)
                                                                          aGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGt                                                                        <  1:1490792/150‑1 (MQ=255)
                                                                             gTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTc                                                                     >  1:725612/1‑150 (MQ=255)
                                                                               cGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCgg                                                                   <  1:4905/150‑1 (MQ=255)
                                                                                ggCCCACCGTCAAAGAGGTCGATTTAGATACGGGAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCTGATGTACCTGACCAATGACGGCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGc                                                                  <  2:1047368/150‑1 (MQ=255)
                                                                                 gCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCa                                                                 <  1:1474559/150‑1 (MQ=255)
                                                                                  cccACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCAt                                                                >  1:1528046/1‑150 (MQ=255)
                                                                                   ccACCGTCATAGAGGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTTGGCTTCATGGGATAAAAAGTGGGTATATATCGGATGTTTCGGCATc                                                               <  1:322201/150‑1 (MQ=255)
                                                                                    cACCGTCAAAGAGGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCAACGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGTTATAGATCGGATGTATCGGCATCa                                                              >  2:1723487/1‑150 (MQ=255)
                                                                                         tCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCt                                                         >  1:1322615/1‑150 (MQ=255)
                                                                                           aaaGCTGTCGAATTAGTTACTGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAACGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGc                                                       <  2:1752031/150‑1 (MQ=255)
                                                                                           aaaGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGCTGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGc                                                       <  2:334528/150‑1 (MQ=255)
                                                                                           aaaGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGc                                                       <  2:475122/150‑1 (MQ=255)
                                                                                           aaaGATGTCGATATAGTGACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGGAGGCGCGGTTTGCGGATTTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAATAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGc                                                       <  2:681541/150‑1 (MQ=255)
                                                                                                        taGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCttt                                          <  1:71771/150‑1 (MQ=255)
                                                                                                          gTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTg                                        >  1:667963/1‑150 (MQ=255)
                                                                                                               ggTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGAACGGATGTTTCGGCATCAGTTCTGCAATAAATGACTTTTGCCCgg                                   >  1:594678/1‑150 (MQ=255)
                                                                                                                   gCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCCATAAATGCCTTTTGCCCGGTGcc                               >  1:1563483/1‑150 (MQ=255)
                                                                                                                   gCGAAAACCTTCTTTCTCCAGCACCGCGAGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTTGGATAAAAAGTGGGTATAGATAGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCAGGTGcc                               >  1:591074/1‑150 (MQ=255)
                                                                                                                    cGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGcag                              >  1:1434423/1‑148 (MQ=255)
                                                                                                                    cGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCg                              >  1:1442471/1‑150 (MQ=255)
                                                                                                                     gAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAGAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGc                             >  1:407889/1‑150 (MQ=255)
                                                                                                                          ccTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAgaga                        >  1:1399775/1‑150 (MQ=255)
                                                                                                                             tctttctCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCGTGGGCTTCGTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGTCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGaaaa                     <  1:1327673/150‑1 (MQ=255)
                                                                                                                             tctttctCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGaaaa                     >  2:478522/1‑150 (MQ=255)
                                                                                                                             tctttctCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCAGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGaaaa                     >  2:1668327/1‑150 (MQ=255)
                                                                                                                             tcatgctcCCGCCCCGGGCGTGCTTGTGCGGTTTTCGTTTTTAGCTGACGGTTGACGTCCTGGTCTTCGTGTGATAAAACGTGGGTATTGTTGGGATGTTTCGGCATCAGTTCTGCAATAACTGCCTTTTGCCCGGTGCCGCAGAGaaaa                     <  2:1040434/145‑1 (MQ=255)
                                                                                                                                ttctCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCg                  >  1:168302/1‑150 (MQ=255)
                                                                                                                                 tctcCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCgg                 <  1:703442/150‑1 (MQ=255)
                                                                                                                                   tcCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGgcg               >  2:23171/1‑150 (MQ=255)
                                                                                                                                   tcCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGgcg               >  1:1653209/1‑150 (MQ=255)
                                                                                                                                   tcCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGgcg               >  2:438562/1‑150 (MQ=255)
                                                                                                                                   ccgcaaaCCGCGCCTGCCCGCGCGGTTTGCGGATGTACCTGACCGCTGACGTCCTGGGCGTCCTGGGATAAAAAGGGGGTTTAGAGCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTGTTGCCCGGTGCCGCAGAGAAAATCGgcg               >  1:283132/7‑150 (MQ=255)
                                                                                                                                    cgcaaaCCGCGCCTGCCCGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGgcgc              >  1:402921/6‑150 (MQ=255)
                                                                                                                                    ccAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGgcgc              >  2:1084510/1‑150 (MQ=255)
                                                                                                                                    ccAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGgcgc              >  1:692191/1‑150 (MQ=255)
                                                                                                                                     gcaaaCCGCGCCTGCCCGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGcgcg             >  2:1786093/5‑150 (MQ=255)
                                                                                                                                       gcgccGCGCGGGCAGGCGCGGGTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCAGGGATAAAAAGTGGGGATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCGGc           <  1:362338/147‑1 (MQ=255)
                                                                                                                                       gCACCGCGCGGGCAGGCGCGGTTTGCGGATGTTCCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGAGGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGATAATCGGCGCGGc           >  1:1679010/1‑150 (MQ=255)
                                                                                                                                       gCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGTCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCGGc           <  1:1196441/150‑1 (MQ=255)
                                                                                                                                       gCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGGCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCGGc           <  1:1557990/150‑1 (MQ=255)
                                                                                                                                       ccACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCGGc           <  2:774332/149‑1 (MQ=255)
                                                                                                                                       aaaCCGCGCCTGCCCGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCGGc           >  2:750126/3‑150 (MQ=255)
                                                                                                                                        aaCCGCGCCTGCCCGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCGGCt          >  2:1841439/2‑150 (MQ=255)
                                                                                                                                           cgcgcgGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGACCGGTGCCGCAGAGAAAATCGGCGCGGCTaaa       >  2:477483/1‑150 (MQ=255)
                                                                                                                                                gggCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATTAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCGGCTAAAATCCa  >  1:544033/1‑150 (MQ=255)
                                                                                                                                                gggCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCGGCTAAAATCCa  <  1:1022340/150‑1 (MQ=255)
                                                                                                                                                gggCAGGCGCGGATTGCGGATGTACCTGACCGATGACGTCATGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAACTGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCGGCTAAAATCCa  <  2:301477/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCGGCTAAAATCCA  >  NC_000913/1829916‑1830209

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: