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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,830,0650GC 8.7% 222.0 / 10.5 69P234R (CCT→CGT) ‡astAarginine N‑succinyltransferase
Reads supporting (aligned to +/- strand):  ref base G (33/30);  new base C (3/3);  total (36/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.72e-01

CTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCGGCTAAAATCCATCGA  >  NC_000913/1829923‑1830213
                                                                                                                                              |                                                                                                                                                    
cTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCCTGCCCGCGCGGt                                                                                                                                               <  1:1259360/150‑1 (MQ=255)
    ccTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCCTGCCCGCGCGGtgctg                                                                                                                                           <  1:987195/150‑5 (MQ=255)
    ccTTCTGCCACTTCCACCAGCCGACTTTTACGGAGGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGc                                                                                                                                           <  2:1527772/150‑1 (MQ=255)
     cTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCg                                                                                                                                          <  2:247459/150‑1 (MQ=255)
     cTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCg                                                                                                                                          >  2:1172630/1‑150 (MQ=255)
      ttCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCgg                                                                                                                                         <  1:1090624/150‑1 (MQ=255)
        cTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGAt                                                                                                                                       >  2:932462/1‑150 (MQ=255)
        cTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGAt                                                                                                                                       >  2:398725/1‑150 (MQ=255)
          gCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGt                                                                                                                                     >  2:1268850/1‑150 (MQ=255)
                 ccaccaGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGAc                                                                                                                              <  1:795097/150‑1 (MQ=255)
                  caccaGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGAcc                                                                                                                             <  1:86701/150‑1 (MQ=255)
                  caccaGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGAcc                                                                                                                             <  1:616367/150‑1 (MQ=255)
                     caGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTAcc                                                                                                                                  >  2:564609/1‑142 (MQ=255)
                     caGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTAcc                                                                                                                                  <  1:564609/142‑1 (MQ=255)
                            cTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTcc                                                                                                                   <  1:337728/150‑1 (MQ=255)
                            cTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTcc                                                                                                                   >  1:742083/1‑150 (MQ=255)
                                tACGGATGGCGCGCACGCGGTCGATGTCACCCTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGCCCGATGCCGTCCTggg                                                                                                               >  1:922124/1‑150 (MQ=255)
                                tACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTggg                                                                                                               <  1:1200860/150‑1 (MQ=255)
                                  cGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCt                                                                                                             >  2:1431899/1‑150 (MQ=255)
                                      tGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCt                                                                                                         >  2:909879/1‑150 (MQ=255)
                                         cgcgcACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTggg                                                                                                      <  1:358212/150‑1 (MQ=255)
                                           cgcACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGAt                                                                                                    >  1:1072014/1‑150 (MQ=255)
                                            gcACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATa                                                                                                   <  2:1500532/150‑1 (MQ=255)
                                               cgcgGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATaaaa                                                                                                >  2:1069807/1‑150 (MQ=255)
                                                    tCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTgg                                                                                           <  1:932462/150‑1 (MQ=255)
                                                     cGATGTCAGACTCAAGCTTCGGCCCACCGTGAAAGATGTCGATATTGTTTCGGTTGTTAAAACCTTCTTTATCCAGCAGCGCGCGTGCAGGGGCGGTTTGTGGATGTACCTGACCGATGACGTCCTGGTCTTCCTGGGATAAAAATTggg                                                                                          <  2:786771/150‑1 (MQ=255)
                                                     cGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTggg                                                                                          <  1:1268850/150‑1 (MQ=255)
                                                     cGATGTCACACTCAAGCGGCGGCCCACCGTCCAAGAGGTCGATATAGGTACGGTAGCGAAAACCTGCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTggg                                                                                          <  1:644954/150‑1 (MQ=255)
                                                         gTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGtata                                                                                      >  2:308064/1‑150 (MQ=255)
                                                          tCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGtatat                                                                                     >  2:776213/1‑149 (MQ=255)
                                                          tCACACTAAAGCGTCGGCCCACAGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTTACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAg                                                                                     >  2:752533/1‑150 (MQ=255)
                                                           cacaCTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGCCGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGa                                                                                    >  1:1399145/1‑150 (MQ=255)
                                                              acTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCg                                                                                 <  2:1491118/150‑1 (MQ=255)
                                                                 cAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGAt                                                                              >  1:227095/1‑150 (MQ=255)
                                                                     cGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAGCCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGttt                                                                          >  2:350195/1‑150 (MQ=255)
                                                                      gTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCCGTTTGCGGATGTCCCTGACCGATGCCGTCCTGGGCTTCCTGGGATAAAAAGAGGGTATAGATAGGATGTTTc                                                                         >  2:356120/1‑150 (MQ=255)
                                                                       gcGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCg                                                                        <  2:1483866/149‑1 (MQ=255)
                                                                         ggCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGc                                                                      >  2:338916/1‑150 (MQ=255)
                                                                          gCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCa                                                                     <  1:1210583/150‑1 (MQ=255)
                                                                           cccACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCAt                                                                    >  1:470905/1‑150 (MQ=255)
                                                                             cACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCa                                                                  <  1:186362/150‑1 (MQ=255)
                                                                              aCCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAg                                                                 <  2:11365/150‑1 (MQ=255)
                                                                                 gTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTc                                                              <  1:752533/150‑1 (MQ=255)
                                                                                 gTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTc                                                              >  1:680129/1‑150 (MQ=255)
                                                                                   cAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTg                                                            >  2:535806/1‑150 (MQ=255)
                                                                                     aaGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCa                                                          >  1:580046/1‑150 (MQ=255)
                                                                                      aGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCaa                                                         >  1:1207222/1‑150 (MQ=255)
                                                                                             gATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGc                                                  >  1:728465/1‑150 (MQ=255)
                                                                                                 taGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCttt                                              >  1:1348268/1‑150 (MQ=255)
                                                                                                             cGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCg                                  >  2:919207/1‑150 (MQ=255)
                                                                                                                    cttcttTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGaa                           >  2:321273/1‑150 (MQ=255)
                                                                                                                       ctttctCCAGCACCGCGCGGGCAGGCGCGCTTTCCGGATGTACCTGACCGATGACGTCCTGGGCTTCATGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAAt                        >  2:1275787/1‑150 (MQ=255)
                                                                                                                         ttctCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCgg                                       >  1:1538841/1‑133 (MQ=255)
                                                                                                                         ttctCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCgg                                       <  2:1538841/133‑1 (MQ=255)
                                                                                                                         ttctCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCg                      >  1:858122/1‑150 (MQ=255)
                                                                                                                           ctcCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGgc                    >  2:1304873/1‑150 (MQ=255)
                                                                                                                           ctcCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGgc                    <  1:308064/150‑1 (MQ=255)
                                                                                                                             cgcaaaCCGCGCCTGCCCGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGgcgc                  >  1:503369/6‑150 (MQ=255)
                                                                                                                             cgcaaaCCGCGCCTGCCCGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGgcgc                  >  2:136925/6‑150 (MQ=255)
                                                                                                                             ccAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGgcgc                  >  2:264897/1‑150 (MQ=255)
                                                                                                                             ccAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGgcgc                  <  1:398725/150‑1 (MQ=255)
                                                                                                                              cAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGcgcg                 <  2:478861/150‑1 (MQ=255)
                                                                                                                              cAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGcgcg                 <  1:798149/150‑1 (MQ=255)
                                                                                                                               caaaCCGCGCCTGCCCGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCgg                >  1:518799/4‑150 (MQ=255)
                                                                                                                               caaaCCGCGCCTGCCCGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCgg                <  1:136925/147‑1 (MQ=255)
                                                                                                                                      cgcgGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCGGCTAAAAt         <  2:728465/150‑1 (MQ=255)
                                                                                                                                           gCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCGGCTAAAATCCATc    <  2:1072014/150‑1 (MQ=255)
                                                                                                                                            cAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCGGCTAAAATCCATCg   <  2:1207222/150‑1 (MQ=255)
                                                                                                                                             aGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCGGCTAAAATCCATCGa  <  1:1172630/150‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
CTGCCCTTCTGCCACTTCCACCAGCCGACTTTTACGGATGGCGCGCACGCGGTCGATGTCACACTCAAGCGTCGGCCCACCGTCAAAGATGTCGATATAGTTACGGTAGCGAAAACCTTCTTTCTCCAGCACCGCGCGGGCAGGCGCGGTTTGCGGATGTACCTGACCGATGACGTCCTGGGCTTCCTGGGATAAAAAGTGGGTATAGATCGGATGTTTCGGCATCAGTTCTGCAATAAATGCCTTTTGCCCGGTGCCGCAGAGAAAATCGGCGCGGCTAAAATCCATCGA  >  NC_000913/1829923‑1830213

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