Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,830,060 G→C 8.2% R236G (CGC→GGC)  astA ← arginine succinyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,830,0600GC 8.2% 288.1 / 16.8 98R236G (CGC→GGC) astAarginine succinyltransferase
Reads supporting (aligned to +/- strand):  ref base G (45/45);  new base C (5/3);  total (50/48)
Fisher's exact test for biased strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.24e-01

Only 100 of 104 total aligned reads displayed.

GCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCGGCTAAAATCCA  >  NC_000913/1829911‑1830209
                                                                                                                                                     |                                                                                                                                                     
gCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACcgcgcg                                                                                                                                                       >  1:1835413/1‑150 (MQ=255)
gCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACcgcgcg                                                                                                                                                       >  1:4521/1‑150 (MQ=255)
   cTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGTGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGc                                                                                                                                                    <  1:1705738/150‑1 (MQ=255)
   cTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGc                                                                                                                                                    >  1:696575/1‑150 (MQ=255)
     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)
                      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)
                                                                      tCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAACCCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGGTTGCGGATGTACCTGAACGATGCAGTCCTCGGCTTCCTGGGATAAAAAGGGGGTATAg                                                                                 >  2:361085/1‑150 (MQ=255)
                                                                      tCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATTAAAAGTGGGTATAg                                                                                 >  2:251738/1‑150 (MQ=255)
                                                                      tCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAg                                                                                 >  2:1721682/1‑150 (MQ=255)
                                                                           cTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCgg                                                                            >  1:438720/1‑150 (MQ=255)
                                                                           cTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCgg                                                                            <  2:713415/150‑1 (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               >  2:438562/1‑150 (MQ=255)
                                                                                                                                        tcCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGgcg               >  1:1653209/1‑150 (MQ=255)
                                                                                                                                        ccgcaaaCCGCGCCTGCCCGCGCGGTTTGCGGATGTACCTGACCGCTGACGTCCTGGGCGTCCTGGGATAAAAAGGGGGTTTAGAGCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTGTTGCCCGGTGCCGCAGAGAAAATCGgcg               >  1:283132/7‑150 (MQ=255)
                                                                                                                                         cgcaaaCCGCGCCTGCCCGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGgcgc              >  1:402921/6‑150 (MQ=255)
                                                                                                                                         ccAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGgcgc              >  1:692191/1‑150 (MQ=255)
                                                                                                                                         ccAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGgcgc              >  2:1084510/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)
                                                                                                                                                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)
                                                                                                                                                     |                                                                                                                                                     
GCCCTGCGCAGGCTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCGGCTAAAATCCA  >  NC_000913/1829911‑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: