Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,830,065 G→C 10.8% P234R (CCT→CGT) ‡ astA ← arginine succinyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,830,0650GC10.8% 243.9 / 15.5 88P234R (CCT→CGT) ‡astAarginine succinyltransferase
Reads supporting (aligned to +/- strand):  ref base G (38/38);  new base C (4/5);  total (44/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCGGCTAAAATCCATCG  >  NC_000913/1829923‑1830212
                                                                                                                                              |                                                                                                                                                   
cTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGt                                                                                                                                              <  1:803099/150‑1 (MQ=255)
 tGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAACCGTCGGCCCACCGTCAAACATGTCGAAATAGTTACGGTAGCCAAAACCTTCTTACACCAGCACCGAGCGCGCACGCGCGGtt                                                                                                                                             >  2:611606/1‑150 (MQ=255)
  gCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCCTGCCCGCGCGGtgc                                                                                                                                            <  1:1081485/150‑3 (MQ=255)
  gCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCCTGCCCGCGCGGtgc                                                                                                                                            <  1:930094/150‑3 (MQ=255)
  gCCCTTCTGCCACTTCCACCAGCCGACTTATACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGttt                                                                                                                                            <  2:481501/150‑1 (MQ=255)
     cTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCCTGCCCGCGCGGtgctgg                                                                                                                                         <  1:183548/150‑6 (MQ=255)
     cTTCTGCCACTTCCACCAGCCGACTTTTACGGAAGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCCTGCCCGCGCGGtgctgg                                                                                                                                         <  1:898454/150‑6 (MQ=255)
     cTTCTGCCACATCCACCAGCCGACTTTTACGGAGGGCGCGCACGCGGTCGATGTCACACGCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCg                                                                                                                                         <  1:1001715/150‑1 (MQ=255)
      ttCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCgg                                                                                                                                        <  1:208951/150‑1 (MQ=255)
      ttCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCgg                                                                                                                                        <  1:168623/150‑1 (MQ=255)
        cTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACAGTCAAAGATGTCGATATTGTTTCGGTAGCGAAACCCTTCTTTCTCAAGCAACGCTCGGGCATGCGGGGGTTGCGGAt                                                                                                                                      >  2:103509/1‑150 (MQ=255)
         tGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATg                                                                                                                                     <  1:303354/150‑1 (MQ=255)
         tGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATg                                                                                                                                     <  2:21783/150‑1 (MQ=255)
          gCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCCAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTACAGCACCGCGCGGGCAGGCGCGGTTTGCgctcaa                                                                                                                                    >  1:1542439/1‑145 (MQ=255)
          gCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGt                                                                                                                                    >  2:1647079/1‑150 (MQ=255)
          gCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGt                                                                                                                                    >  1:1448599/1‑150 (MQ=255)
            cACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTAc                                                                                                                                  >  1:973032/1‑150 (MQ=255)
            cACTACCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGGAGGCGCGGTGTGCGGATGTAc                                                                                                                                  >  1:716009/1‑150 (MQ=255)
             aCTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACACAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTAcc                                                                                                                                 <  2:953264/150‑1 (MQ=255)
                tCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAAGCGCGGTTTGCGGATGTACCTGa                                                                                                                              >  1:1458180/1‑150 (MQ=255)
                 ccaccaGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGAc                                                                                                                             >  2:1718863/1‑150 (MQ=255)
                  caccaGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAATACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGAcc                                                                                                                            >  2:861956/1‑150 (MQ=255)
                  caccaGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGAcc                                                                                                                            >  1:998327/1‑150 (MQ=255)
                      aGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATg                                                                                                                        <  2:379164/150‑1 (MQ=255)
                        ccGACTTTTACTGATGGCGCGCACGCGGTCGATGTCACACTCAACCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGAc                                                                                                                      <  2:81166/150‑1 (MQ=255)
                           aCTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTc                                                                                                                   >  1:271849/1‑150 (MQ=255)
                               ttACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTgg                                                                                                               <  2:1361608/150‑1 (MQ=255)
                                    gATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTc                                                                                                          >  1:934105/1‑150 (MQ=255)
                                    gATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTc                                                                                                          <  2:1435910/150‑1 (MQ=255)
                                       ggCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTg                                                                                                       <  1:1536671/150‑1 (MQ=255)
                                       ggCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTg                                                                                                       <  2:1579472/150‑1 (MQ=255)
                                        gcgcgcACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTgg                                                                                                      <  1:193205/150‑1 (MQ=255)
                                           cgcACGCGGGCGATGTCACACTCACGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGAt                                                                                                   <  2:677092/150‑1 (MQ=255)
                                             cACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGCCCGATGACGTCCTGGGCTTCCTGGGATaa                                                                                                 >  2:1071552/1‑150 (MQ=255)
                                               cgcgGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATaaaa                                                                                               >  2:1365235/1‑150 (MQ=255)
                                                     cGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTggg                                                                                         <  2:1000269/150‑1 (MQ=255)
                                                     cGAGGTCACACCCCCGCGTCGGCCCACCGTCACAGATGTCGATATAGTGACGGTAGCGAACACCTTCTTTCTCCAGCACCGCGCGGGCACGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAACAGTggg                                                                                         <  2:1052835/150‑1 (MQ=255)
                                                      gATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGt                                                                                        <  1:579764/150‑1 (MQ=255)
                                                      gATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGt                                                                                        <  2:1242931/150‑1 (MQ=255)
                                                       aTGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCAATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTggtat                                                                                       <  2:1300527/150‑4 (MQ=255)
                                                           cacaCTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGa                                                                                   >  1:421490/1‑150 (MQ=255)
                                                           cacaCTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGa                                                                                   >  2:708889/1‑150 (MQ=255)
                                                           cacaCTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGa                                                                                   >  2:544652/1‑150 (MQ=255)
                                                                 aaaGCGTCGGCCCACCGTCAATGATGTCGATTTAGTTACGGTAGCGAAAACGTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGGTTCCTGGGATAAAAAGTTGGTATAGATCGGAt                                                                             <  1:611606/149‑1 (MQ=255)
                                                                    gCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGCCGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGtt                                                                          >  2:407965/1‑150 (MQ=255)
                                                                     cGTCGGCCCACCGGCAAAGATGGCGATATAGTTACGGTAGGGTCACCCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGttt                                                                         <  1:449673/150‑1 (MQ=255)
                                                                        cGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCgg                                                                      <  2:1417761/150‑1 (MQ=255)
                                                                         ggCCCACAGTCAAAGATGTCGATATAGTGTCGGTAGCGAAAACCTTCTTTCTCCAGCACCCCGCGGGCAGGCGCGGTTTGAGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGTTATAGATCGGATTTTTCGGc                                                                     <  2:76895/150‑1 (MQ=255)
                                                                           cccACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCCCCGCGCGGGCAGGCGCGGTTGGCGGATGTAGCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGGGGGGAGAGATCGGATGTTTCGGcct                                                                   >  1:215385/1‑148 (MQ=255)
                                                                            ccACCGTCAAAGATGTCGATATAGTTACGTTAGCGAAAACCTTCTTGATCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATc                                                                  <  2:1256056/150‑1 (MQ=255)
                                                                              aCCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGGCCTGGGCTTCCTGGGTTTAATTGTGGGTATAGATCGGATGTTTCGGCATCAg                                                                >  2:722539/1‑150 (MQ=255)
                                                                              aCCGTCAAAGATCACGAAATATTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGCCAGGCGCGGTTTGCGGATGTCCCTTCCCGATGACGGCCTGGGCTTCCTTGGATAAAAAGTGCGTAGAGCTCGGATGATACGGCATCac                                                                >  2:728918/1‑149 (MQ=255)
                                                                                cGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGtt                                                              >  1:1066717/1‑150 (MQ=255)
                                                                                   cAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTg                                                           <  1:1474995/150‑1 (MQ=255)
                                                                                         tGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATaa                                                     <  2:1525785/150‑1 (MQ=255)
                                                                                          gTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATaaa                                                    <  2:367126/150‑1 (MQ=255)
                                                                                            cGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCTGTTTGCGGATGTACCTGAAAGATGACGTCATGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCTGAATCAGTTCTGCAATAAATg                                                  >  1:1528454/1‑150 (MQ=255)
                                                                                            cGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATg                                                  >  1:833527/1‑150 (MQ=255)
                                                                                            cGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATg                                                  >  1:342694/1‑150 (MQ=255)
                                                                                             gATATAGTTACGGTAGCGAAACCCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGCTCGGATGTTTCGGCATCAGTTCTGCAATAAATGc                                                 >  1:887294/1‑150 (MQ=255)
                                                                                             gATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGc                                                 >  1:1403186/1‑150 (MQ=255)
                                                                                             gATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGc                                                 >  1:223457/1‑150 (MQ=255)
                                                                                                   gTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTg                                           >  2:724634/1‑150 (MQ=255)
                                                                                                      aCGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGGCCTGGGCTTCCTGGGATAAAACGTGGGAATAGAACGGATGTTTCTGCATCAGTTCTACAATAAATGCCTTTTGccc                                        <  1:102910/150‑1 (MQ=255)
                                                                                                        ggTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCgg                                      >  2:63030/1‑150 (MQ=255)
                                                                                                           aGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGGAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGc                                   <  1:79697/150‑1 (MQ=255)
                                                                                                             cGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCg                                 <  2:598942/150‑1 (MQ=255)
                                                                                                               aaaaCCTTCTTTCTCCAGCACCGCGCGGGCATGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCGTCCTGGGATAAAAAGTGGGTCTAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCa                               <  2:868910/150‑1 (MQ=255)
                                                                                                               aaaaCCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCCTCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCa                               >  1:859439/1‑150 (MQ=255)
                                                                                                               aaaaCCTTCTTTCACCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGGGGGTATAGAACGGATGTTACGGCATCAGTTATGCAATAAATGCATTTTGCCCGGTCCCGCa                               >  2:951305/1‑150 (MQ=255)
                                                                                                                  aCCTTCTTTCTCCAGCAGCGCGCGGGCAGGCGCGGTATGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAACAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCagag                            <  1:288953/150‑1 (MQ=255)
                                                                                                                  aCCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCagag                            >  2:57636/1‑150 (MQ=255)
                                                                                                                   ccTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGCTCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGACCGGTGCCGCAgaga                           >  1:678043/1‑150 (MQ=255)
                                                                                                                   ccTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAgaga                           >  2:1231584/1‑150 (MQ=255)
                                                                                                                     ttctttctCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGaaa                         <  1:133826/150‑1 (MQ=255)
                                                                                                                           ctcCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGgc                   <  1:743599/150‑1 (MQ=255)
                                                                                                                           ctcCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATAGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGAAGAGAAAATCGgc                   >  2:304310/1‑150 (MQ=255)
                                                                                                                           ctcCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGGCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGAGCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGgc                   <  1:798825/150‑1 (MQ=255)
                                                                                                                             cgcaaaCCGCGCCTGCCCGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGgcgc                 >  2:79661/6‑150 (MQ=255)
                                                                                                                             ccAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGgcgc                 >  1:536941/1‑150 (MQ=255)
                                                                                                                             ccAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGgcgc                 >  2:351051/1‑150 (MQ=255)
                                                                                                                              gcaaaCCGCGCCTGCCCGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGcgcg                >  1:1657536/5‑150 (MQ=255)
                                                                                                                                aaaCCGCGCCTGCCCGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGACATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCTCAGAGAAAATCGGCGCGGc              >  2:575327/3‑150 (MQ=255)
                                                                                                                                 cACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCGGCt             <  2:342694/150‑1 (MQ=255)
                                                                                                                                   ccGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCGGCTaa           >  1:30847/1‑150 (MQ=255)
                                                                                                                                   ccGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGAGGCCGTCCTGGGCTTCCTGGGATAAAAAGTGGGTAGAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGATAATCGGCGCGCCTaa           >  1:1623054/1‑150 (MQ=255)
                                                                                                                                           gCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCGGCTAAAATCCATc   <  1:544652/150‑1 (MQ=255)
                                                                                                                                            cAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCGGCTAAAATCCATCg  <  2:1403186/150‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
CTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCGGCTAAAATCCATCG  >  NC_000913/1829923‑1830212

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: