Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,221,158 C→G 6.2% K30N (AAG→AAC ygjH ← putative tRNA‑binding protein YgjH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,221,1580CG 6.2% 272.6 / 6.1 81K30N (AAG→AACygjHputative tRNA‑binding protein YgjH
Reads supporting (aligned to +/- strand):  ref base C (34/42);  new base G (3/2);  total (37/44)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

GACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACTATAGCAAAAGTTTTTCTCAATCCTGT  >  NC_000913/3221012‑3221307
                                                                                                                                                  |                                                                                                                                                     
gACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACATCGTACAGGTTg                                                                                                                                                    <  2:1473326/150‑1 (MQ=255)
gACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACATCGTACAGGTTg                                                                                                                                                    <  1:650264/150‑1 (MQ=255)
    tttCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCgg                                                                                                                                                <  1:783911/150‑1 (MQ=255)
     ttCACCACGCATTTTCGCTTTTTGCAGATTGTATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGc                                                                                                                                               >  1:62311/1‑150 (MQ=255)
         ccacGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGttt                                                                                                                                           >  1:1570369/1‑150 (MQ=255)
          cacGCATTTTCGCTTTTTGCAGATTGCATAAGCCCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCACGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACCATGTACAGCTTGTCGGCGtttt                                                                                                                                          >  2:1112583/1‑150 (MQ=255)
            cGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCa                                                                                                                                        <  1:836955/150‑1 (MQ=255)
             gCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCAt                                                                                                                                       <  1:1602805/150‑1 (MQ=255)
               aTTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATgg                                                                                                                                     <  2:421922/150‑1 (MQ=255)
               aTTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATgg                                                                                                                                     <  1:144755/150‑1 (MQ=255)
                      cTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCa                                                                                                                              >  2:1458451/1‑150 (MQ=255)
                       tttttGGAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGGCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCAc                                                                                                                             <  1:1048229/150‑1 (MQ=255)
                       tttttGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCAc                                                                                                                             <  2:1599625/150‑1 (MQ=255)
                         tttGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACtt                                                                                                                           <  1:1216894/150‑1 (MQ=255)
                         tttGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACtt                                                                                                                           <  2:813203/150‑1 (MQ=255)
                           tGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTcc                                                                                                                         <  2:605062/150‑1 (MQ=255)
                            gCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCca                                                                                                                        <  2:67200/150‑1 (MQ=255)
                            gCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCca                                                                                                                        >  2:754586/1‑150 (MQ=255)
                                  tGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCt                                                                                                                  <  1:1000006/150‑1 (MQ=255)
                                   gCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCtt                                                                                                                 >  1:378416/1‑150 (MQ=255)
                                      tAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTcc                                                                                                              <  1:929196/150‑1 (MQ=255)
                                           ccaccaccGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACGAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGAcgc                                                                                                         <  1:1061927/150‑1 (MQ=255)
                                               caccGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCAt                                                                                                     >  2:1188995/1‑150 (MQ=255)
                                                 ccGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCAttt                                                                                                   <  1:1041014/150‑1 (MQ=255)
                                                  cGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTc                                                                                                  >  2:951770/1‑150 (MQ=255)
                                                         cccATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGt                                                                                           >  2:568067/1‑150 (MQ=255)
                                                          ccATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGGACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTg                                                                                          >  1:347297/1‑150 (MQ=255)
                                                            aTCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCa                                                                                        <  1:1370667/150‑1 (MQ=255)
                                                                 gtcttctTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATc                                                                                   <  2:411704/149‑1 (MQ=255)
                                                                 ctctTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATc                                                                                   >  2:607859/1‑150 (MQ=255)
                                                                  tcttctTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCa                                                                                  >  1:1555966/1‑150 (MQ=255)
                                                                    ttcttcGCTGTAGTACGGCACCAGGCTGGTCCCGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGc                                                                                >  1:820837/1‑150 (MQ=255)
                                                                     tcttcGCTGTAGTACGGGACCAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCACCCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCg                                                                               <  2:683325/150‑1 (MQ=255)
                                                                         cGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAg                                                                           >  1:712572/1‑150 (MQ=255)
                                                                            tGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCa                                                                        >  2:438702/1‑150 (MQ=255)
                                                                              tagtaCGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACg                                                                      <  1:568067/150‑1 (MQ=255)
                                                                               agtaCGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACgg                                                                     <  1:574628/150‑1 (MQ=255)
                                                                                  aCGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGttt                                                                  >  2:466532/1‑150 (MQ=255)
                                                                                  aCGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGCTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGttt                                                                  >  2:968174/1‑150 (MQ=255)
                                                                                   cGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTc                                                                 <  1:1188995/150‑1 (MQ=255)
                                                                                     gCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCa                                                               >  1:3877/1‑150 (MQ=255)
                                                                                      cACCAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCAt                                                              >  2:1300993/1‑150 (MQ=255)
                                                                                       cccAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCAtt                                                             <  1:1301276/149‑1 (MQ=255)
                                                                                        ccAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTg                                                            <  1:754586/150‑1 (MQ=255)
                                                                                           ggCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTTTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCtt                                                         <  2:1091199/150‑1 (MQ=255)
                                                                                             cTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTtcc                                                       <  2:1414886/150‑1 (MQ=255)
                                                                                               ggTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTtcctc                                                     <  2:1570369/150‑1 (MQ=255)
                                                                                                gTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTtcctcc                                                    <  2:1555966/150‑1 (MQ=255)
                                                                                                  cACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTGCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCt                                                  >  2:186709/1‑150 (MQ=255)
                                                                                                     ggTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAg                                               >  2:36907/1‑150 (MQ=255)
                                                                                                      gTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAgg                                              >  2:617803/1‑150 (MQ=255)
                                                                                                           cAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTaaa                                         <  1:1256673/150‑1 (MQ=255)
                                                                                                           cAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTaaa                                         >  1:205956/1‑150 (MQ=255)
                                                                                                                ttttttGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAAtttt                                    <  2:50359/150‑1 (MQ=255)
                                                                                                                  ttttGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTAc                                  <  1:1568039/150‑1 (MQ=255)
                                                                                                                  ttttGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTAc                                  <  1:496840/150‑1 (MQ=255)
                                                                                                                       ccccACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAAc                             <  1:1109529/150‑1 (MQ=255)
                                                                                                                       ccccACATCAACCTGTACGATGTACACCTTGGCGGCGGTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCACACGTGGAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAAc                             <  1:1112583/150‑1 (MQ=255)
                                                                                                                         ccACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACta                           >  2:298876/1‑150 (MQ=255)
                                                                                                                          cacaTCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACtat                          <  2:347297/150‑1 (MQ=255)
                                                                                                                           ccgacAAGCTGTACATCGTACAGGTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACtata                         >  1:1265826/5‑150 (MQ=255)
                                                                                                                           ccgacAAGCTGTACATCGTACAGGTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACtata                         >  2:589561/5‑150 (MQ=255)
                                                                                                                           ccgacAAGCTGTACATCGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACtata                         >  2:758791/5‑150 (MQ=255)
                                                                                                                           acaTCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACtata                         <  2:1462117/150‑1 (MQ=255)
                                                                                                                            caTCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACTATAg                        >  2:1481317/1‑150 (MQ=255)
                                                                                                                            caTCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACTATAg                        >  2:1603345/1‑150 (MQ=255)
                                                                                                                             aTCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACTATAGc                       <  1:944929/150‑1 (MQ=255)
                                                                                                                              tCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACTATAGCa                      >  1:288350/1‑150 (MQ=255)
                                                                                                                              tCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACTATAGCa                      <  1:951770/150‑1 (MQ=255)
                                                                                                                               cAAGCTGTACATCGTACAGGTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACTATAGCaa                     >  1:1474466/1‑150 (MQ=255)
                                                                                                                               cAACCTGTACGATCTACAGCTTGTCCGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGCCGCGCATTTCCAGCCGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTAATAAACTATAGCaa                     >  2:1193047/1‑150 (MQ=255)
                                                                                                                                aaCCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACTATAGCaaa                    >  2:1307649/1‑150 (MQ=255)
                                                                                                                                     gTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACTATAGCAAAAGttt               >  2:41405/1‑150 (MQ=255)
                                                                                                                                      tatgATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACTATAGCAAAAGtttt              >  1:507229/4‑150 (MQ=255)
                                                                                                                                      tACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACTATAGCAAAAGtttt              <  2:3877/150‑1 (MQ=255)
                                                                                                                                        cGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACTATAGCAAAAGTTTTtc            >  1:539638/1‑150 (MQ=255)
                                                                                                                                        cGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACTATAGCAAAAGTTTTtc            <  1:1223764/150‑1 (MQ=255)
                                                                                                                                        cGATGTACAGCTTGGCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACTATAGCAAAAGTTTTtc            <  1:798127/150‑1 (MQ=255)
                                                                                                                                         gATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACTATAGCAAAAGTTTTtct           <  2:378416/150‑1 (MQ=255)
                                                                                                                                          aTGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACTATAGCAAAAGTTTTtctc          <  1:438702/150‑1 (MQ=255)
                                                                                                                                               cAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACTATAGCAAAAGTTTTTCTCAATcc     >  1:875028/1‑150 (MQ=255)
                                                                                                                                                  cTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACTATAGCAAAAGTTTTTCTCAATCCTGt  <  2:288350/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACTATAGCAAAAGTTTTTCTCAATCCTGT  >  NC_000913/3221012‑3221307

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: