Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,221,149 G→A 7.4% I33I (ATC→ATT) ‡ 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,1490GA 7.4% 240.8 / 7.0 81I33I (ATC→ATT) ‡ygjHputative tRNA‑binding protein YgjH
Reads supporting (aligned to +/- strand):  ref base G (32/42);  new base A (4/2);  total (37/44)
Fisher's exact test for biased strand distribution p-value = 4.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.11e-01

GCATTCCGACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACTATAGCAAAAGTTTTTCT  >  NC_000913/3221005‑3221298
                                                                                                                                                |                                                                                                                                                     
gcatTCCGACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTa                                                                                                                                                  >  2:1109529/1‑150 (MQ=255)
 catTCCGACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCCCCAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTAc                                                                                                                                                 >  2:1256673/1‑150 (MQ=255)
 catTCCGACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTAc                                                                                                                                                 <  1:432355/150‑1 (MQ=255)
  atTCCGACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACa                                                                                                                                                <  1:738272/150‑1 (MQ=255)
       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:496840/150‑1 (MQ=255)
                                                                                                                         ttttGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTAc                         <  1:1568039/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:1603345/1‑150 (MQ=255)
                                                                                                                                   caTCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACTATAg               >  2:1481317/1‑150 (MQ=255)
                                                                                                                                    aTCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACTATAGc              <  1:944929/150‑1 (MQ=255)
                                                                                                                                     tCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACTATAGCa             <  1:951770/150‑1 (MQ=255)
                                                                                                                                     tCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACTATAGCa             >  1:288350/1‑150 (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:1223764/150‑1 (MQ=255)
                                                                                                                                               cGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACTATAGCAAAAGTTTTtc   >  1:539638/1‑150 (MQ=255)
                                                                                                                                               cGATGTACAGCTTGGCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACTATAGCAAAAGTTTTtc   <  1:798127/150‑1 (MQ=255)
                                                                                                                                                gATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACTATAGCAAAAGTTTTtct  <  2:378416/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
GCATTCCGACGTTTCACCACGCATTTTCGCTTTTTGCAGATTGCATAAGACCACCACCGTTTTCCCCATCAGCTCTTCTTCGCTGTAGTACGGCACCAGGCTGGTCACGGTTTGCAGCGTTTTTTGCCCCACATCAACCTGTACGATGTACAGCTTGTCGGCGTTTTCATGGCGTTTCACTTCCACAATCTTTCCGACGCGCATTTCCAGACGTGCAAAATCAGCGTAAGCCACGGTTTCCATTGCTTCCTCCCTTAGGGTAAAATTTTACTAAACTATAGCAAAAGTTTTTCT  >  NC_000913/3221005‑3221298

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: