Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,327,320 T→C 5.7% intergenic (+27/+47) atoB → / ← yfaP acetyl‑CoA acetyltransferase/DUF2135 domain‑containing protein YfaP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,327,3200TC 5.7% 270.4 / 4.7 89intergenic (+27/+47)atoB/yfaPacetyl‑CoA acetyltransferase/DUF2135 domain‑containing protein YfaP
Reads supporting (aligned to +/- strand):  ref base T (43/40);  new base C (2/3);  total (46/43)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GGTGCTCGTATTCTGGTCACACTATTACATGCCATGCAGGCACGCGATAAAACGCTGGGGCTGGCAACACTGTGCATTGGCGGCGGTCAGGGAATTGCGATGGTGATTGAACGGTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAACTTTTCACCAGCGTCAGCTCGCCAGCATTACGCATCGGTACAATAAATGTTTCCTGTTTCTCATTGACCGATCCTTCATCGGTGAT  >  NC_000913/2327171‑2327465
                                                                                                                                                     |                                                                                                                                                  
ggtgCTCGTATTCTGGTCACACTATTACATGCCATGCAGGCACGCGATAAAACGCTGGGGCTGGCAACACTGTGCATTGGCGGCGGTCAGGGAATTGCGATGGTGATTGAACGGTTGAATTAATCAATAAAAACACCCGATAGCGAAAGt                                                                                                                                                    <  1:1618496/150‑1 (MQ=255)
   gCTCGTATTCTGGTCACACTATTACATGCCATGCAGGCACGCGATAAAACGCTGGGGCTGGCAACACTGTGCATTGGCGGCGGTCAGGGAATTGCGATGGTGATTGAACGGTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTAt                                                                                                                                                 >  2:1067675/1‑150 (MQ=255)
       gTATTCTGGTCACACTATTACATGCCATGCAGGCACGCGATAAAACGCTGGGGCTGGCAACACTGTGCATTGGCGGCGGTCAGGGAATTGCGATGGTGATTGAACGGTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCggg                                                                                                                                             >  2:572955/1‑150 (MQ=255)
       gTATTCTGGTCACACTATTACATGCCATGCAGGCACGCGATAAAACGCTGGGGCTGGCAACACTGTGCATTGGCGGCGGTCAGGGAATTGCGATGGTGATTGAACGGTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCggg                                                                                                                                             <  2:430621/150‑1 (MQ=255)
       gTATTCTGGTCACACTATTACATGCCATGCAGGCACGCGATAAAACGCTGGGGCTGGCAACACTGTGCATTGGCGGCGGTCAGGGAATTGCGATGGTGATTGAACGGTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCggg                                                                                                                                             <  1:644278/150‑1 (MQ=255)
         aTTCTGGTCACACTATTACATGCCATGCAGGCACGCGATAAAACGCTGGGGCTGGCAACACTGTGCATTCGCGGCTGTCAGGGAATTGCGATTGTTATTGAACGGTTGAATTAATCAATAAAAACACCCGATAACTATCGCTATCGGgtg                                                                                                                                           <  2:440204/150‑1 (MQ=255)
           tCTGGTCACACTATTACATGCCATGCAGGCACGCGATAAAACGCTGGGGCTGGCAACACTGTGCATTGGCGGCGGTCAGGGAATTGCGATGGTGATTGAACGGTTGAATTAATCAATAAAAACACCCGATAACTTTCGCTATCGGGTGtt                                                                                                                                        <  2:1045457/150‑1 (MQ=255)
            cTGGTCACACTATTACATGCCATGCAGGCACGCGATAAAACGCTGGGGCTGGCAACACTGTGCATTGGCGGCGGTCAGGGAATTGCGATGGTGATTGAACGGTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGttt                                                                                                                                       >  2:1330512/1‑150 (MQ=255)
              ggTCACACTATTACATGCCATGCAGGCACGCGATAAAACGCTGGGGCTGGCAACACTGTGCATTGGCGGCGGTCAGGGAATTGCGATGGTGATTGAACGGTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTc                                                                                                                                     <  1:512790/150‑1 (MQ=255)
              ggTCACACTATTACATGCCATGCAGGCACGCGATAAAACGCTGGGGCTGGCAACACTGTGCATTGGCGGCGGTCAGGGAATTGCGATGGTGATTGAACGGTTGAATTAATCAATAAAAACACCCGATAACTTTCGCTATCGGGTGttttt                                                                                                                                     >  1:288733/1‑149 (MQ=255)
               gTCACACTATTACATGCCATGCAGGCACGCGATAAAACGCTGGGGCTGGCAACACTGTGCATTGGCGGCGGTCAGGGAATTGCGATGGTGATTGAACGGTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCt                                                                                                                                    <  2:913786/150‑1 (MQ=255)
               gTCACACTATTACATGCCATGCAGGCACGCGATAAAACGCTGGGGCTGGCAACACTGTGCATTGGCGGCGGTCAGGGAATTGCGATGGTGATTGAACGGTTGAATTAATCAATAAAAACACCCGATAACTTTCGCTATCGGGTGttttta                                                                                                                                    <  2:288733/150‑3 (MQ=39)
               gTCACAATATTACATGGGATGGAGGCACTCGATTAAACGCTGGTGCTGGCAACACTGTGAATTGGCGGCGGTCAGGGAATTGCGATGTTGATTGTACGGTTGAATTAATGAATAAAATCACCCGATAGCGAATGTTATGGGGTGTTTTCt                                                                                                                                    <  2:144266/150‑1 (MQ=255)
                 cacaCTATTACATGCCATGCAGGCACGCGATAAAACGCTGGGGCTGGCAACACTGTGCATTGGCGGCGGTCAGGGAATTGCGATGGTGATTGAACGGTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTg                                                                                                                                  <  1:338643/150‑1 (MQ=255)
                  acacTATTACATGCCATGCAGGCACGCGATAAAACGCTGGGGCTGGCAACACTGTGCATTGGCGGCGGTCAGGGAATTGCGATGGTGATTGAACGGTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGa                                                                                                                                 >  1:422482/1‑150 (MQ=255)
                     cTATTACATGCCATGCAGGCACGCGATAAAACGCTGGGGCTGGCAACACTGTGCATTGGCGGCGGTCAGGGAATTGCGATGGTGATTGAACGGTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACa                                                                                                                              <  1:1265647/150‑1 (MQ=255)
                     cTATTACATGCCATGCACGCACGCGATAAAACGCTGGGGCTGGCAACACTGTGCATTGGCGGCGGTCAGGGAATTGCGATCGTGAATGAACGGTTGAATTAATCAATAAAAACACACGATAGCGAAAGTTATCGGGTGTTTTCTTGAACa                                                                                                                              >  1:1052678/1‑150 (MQ=255)
                         tACATGCCATGCAGGCACGCGATAAAACGCTGGGGCTGGCAACACTGTGCATTGGCGGCGGTCAGGGAATTGCGATGGTGATTGAACGGTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGa                                                                                                                          >  1:268061/1‑150 (MQ=255)
                          aCATGCCATGCAGGCGCGCGATAAAACGCTGGGGCTGGCAACACTGTGCATTGGCGGCGGTCAGGGAATTGCGATGGTGATTGAACGGTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGAc                                                                                                                         >  1:143317/1‑150 (MQ=255)
                          aCATGCCATGCAGGCACGCGATAAAACGCTGGGGCTGGCAACACTGTGCATTGGCGGCGGTCAGGGAATTGCGATGGTGATTGAACGATTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTg                                                                                                                                  <  2:734605/141‑1 (MQ=255)
                          aCATGCCATGCAGGCACGCGATAAAACGCTGGGGCTGGCAACACTGTGCATTGGCGGCGGTCAGGGAATTGCGATGGTGATTGAACGATTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTg                                                                                                                                  >  1:734605/1‑141 (MQ=255)
                             tgccatgcAGGCACGCGATAAAACGCTGGGGCTGGCAACACTGTGCATTGGCGGCGGTCAGGGAATTGCGATGGTGATTGAACGGTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGc                                                                                                                      >  1:144273/1‑150 (MQ=255)
                                catgcAGGCACGCGATAAAACGCTGGGGCTGGCAACACTGTGCATTGGCGGCGGTCAGGGAATTGCGATGGTGATTGAACGGTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGaa                                                                                                                   >  1:1359003/1‑150 (MQ=255)
                                  tgcAGGCACGCGATAAAACGCTGGGGCTGGCAACACTGTGCATTGGCGGCGGTCAGGGAATTGCGATGGTGATTGAACGGTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAgg                                                                                                                 <  2:1282929/150‑1 (MQ=255)
                                    cAGGCACGCGATAACACGCTGGGGCTGGCAACACTGTGCATTGGCGGCGGTCAGGGAATTGCGATGGTGATTGAACGGTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCGGGGGTTTTCTTGAACATCGACGGCGAAGGTa                                                                                                               >  2:1353590/1‑150 (MQ=255)
                                         aCGCGATAAAACGCTGGGGCTGGCAACACTGTGCATTGGCGGCGGTCAGGGAATTGCGATGGTGATTGAACGGTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAAcccc                                                                                                          >  1:1046984/1‑150 (MQ=255)
                                           gcgATAAAACGCTGGGTCTGGCAACACTGTGCATTGGCGGCGGTCAGGTAATTGCGATGGTGATTGAACGGTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCAt                                                                                                        <  1:1101418/150‑1 (MQ=255)
                                           gcgATAAAACGCTGGGGCTGGCAACACTGTGCATTGGCGGCGGTCAGGGAATTGCGATGGTGATTGAACGGTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCAt                                                                                                        <  2:421760/150‑1 (MQ=255)
                                            cgATAAAACGCTGGGGCTGGCAACACTGTGCATTGGCGGCGGTCAGGGAATTGCGATGGTGATTGAACGGTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCAtt                                                                                                       >  1:642462/1‑150 (MQ=255)
                                               tAAAACGCTGGGGCTGGCAACACTGTGCATTGGCGGCGGTCAGGGAATTGCGATGGTGATTGAACGGTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCg                                                                                                                                               <  1:158601/108‑1 (MQ=255)
                                               tAAAACGCTGGGGCTGGCAACACTGTGCATTGGCGGCGGTCAGGGAATTGCGATGGTGATTGAACGGTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCg                                                                                                                                               >  2:158601/1‑108 (MQ=255)
                                                 aaaCGCTGGGGCTGGCAACACTGTGCATTGGCGGCGGTCAGGGAATTGCGATGGTGATTGAACGGTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCa                                                                                                  >  2:696319/1‑150 (MQ=255)
                                                 aaaCGCTGGGGCTGGCAACACTGTGCATTGGCGGCGGTCAGGGAATTGCGATGGTGATTGAACGGTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCa                                                                                                  >  1:1554961/1‑150 (MQ=255)
                                                  aaCGCTGGGGCTGGCAACACTGTGCATTGGCGGCGGTCAGGGAATTGCGATGGTGATTGAACGGTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCAc                                                                                                 <  1:351451/150‑1 (MQ=255)
                                                          ggCTGGCAACACTGTGCATTGGCGGCGGTCAGGGAATTGCGATGGTGATTGAACGGTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATc                                                                                                   >  2:1418959/1‑140 (MQ=255)
                                                          ggCTGGCAACACTGTGCATTGGCGGCGGTCAGGGAATTGCGATGGTGATTGAACGGTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATc                                                                                                   <  1:1418959/140‑1 (MQ=255)
                                                            cTGGCAACACTGTGCATTGGCGGCGGTCAGGGAATTGCGATGGTGATTGAACGGTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAAc                                                                                       <  1:547052/150‑1 (MQ=255)
                                                            cTGGCAACACTGTGCATTGGCGGCGGTCAGGGAATTGCGATGGTGATTGAACGGTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAAc                                                                                       >  2:866404/1‑150 (MQ=255)
                                                              ggCAACACTGTGCATTGGCGGCGGTCAGGGAATTGCGATGGTGATTGAACGGTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAACtt                                                                                     >  1:430239/1‑150 (MQ=255)
                                                                 aaCACTGTGCATTGGCGGCGGTCAGGGAATTGCGATGGTGATTGAACGGTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAACTTTTc                                                                                  <  1:521546/150‑1 (MQ=255)
                                                                        tgCATTGGCGGCGGTCAGGGAATTGCGATGGTGATTGAACGGTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAACTTTTCACCAGCg                                                                           <  1:1587411/150‑1 (MQ=255)
                                                                        tgCATTGGCGGCGGTCAGGGAATTGCGATGGTGATTGAACGGTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGACAAAACTTTTCACCAGCg                                                                           <  2:576053/150‑1 (MQ=255)
                                                                         gCATTGGCGGCGGTCAGGGAATTGCGATGGTGATTGAACGGTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAACTTTTCACCAGCGt                                                                          >  2:335207/1‑150 (MQ=255)
                                                                          cATTGGCGGCGGTCAGGGAATTGCGATGGTGATTGAACGGTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAACTTTTCACCAGCGTc                                                                         <  2:878392/150‑1 (MQ=255)
                                                                              ggcggcGGTCAGGGAATTGCGATGGTGATTGAACGGTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAACTTTTCACCAGCGTCAGCt                                                                     <  1:1067675/150‑1 (MQ=255)
                                                                                cggcggTCAGGGAATTGCGATGGTGATTGAACGGTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAACTTTTCACCAGCGTCAGCTCg                                                                   >  2:993539/1‑150 (MQ=255)
                                                                                  gcggTCAGGGAATTGCGATGGTGATTGAACGGTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAACTTTTCACCAGCGTCAGCTCGcc                                                                 <  2:430239/150‑1 (MQ=255)
                                                                                     gTCAGGGAATTGCGATGGTGATTGAACGGTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAACTTTTCACCAGCGTCAGCTCGCCAGc                                                              <  2:642462/150‑1 (MQ=255)
                                                                                       cAGGGAATTGCGATGGTGATTGAACGGTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAACTTTTCACCAGCGTCAGCTCGCCAGCAt                                                            >  2:89031/1‑150 (MQ=255)
                                                                                             aTTGCGATGGTGATTGAACGGTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAACTTTTCACCAGCGTCAGCTCGCCAGCATTACGCa                                                      >  1:749069/1‑150 (MQ=255)
                                                                                                    tGGTGATTGAACGGTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAACTTTTCACCAGCGTCAGCTCGCCAGCATTACGCATCGGTAc                                               <  1:1330512/150‑1 (MQ=255)
                                                                                                      gttaTTGAACGGTAGAAATAATCAATAAAAACACCCGATAGCGAAAGATAACCGGAGTTTTCTTGAACATCGACGGCGAAGGTAACCCCAAAAATAACCAGTCAAAACATTTCACCAGCGTCAGCTCGCCAGCATTACGCATCGGTACaa                                             >  2:478076/4‑150 (MQ=255)
                                                                                                        gATTGAACGGTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAACTTTTCACCAGCGTCAGCTCGCCAGCATTACGCATCGGTACAATa                                           >  2:17392/1‑150 (MQ=255)
                                                                                                         aTTGAACGGTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAACTTTTCACCAGCGTCAGCTCGCCAGCATTACGCATCGGTACAATaa                                          >  2:231431/1‑150 (MQ=255)
                                                                                                         aTTGAACGGTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAACTTTTCACCAGCGTCAGCTCGCCAGCATTACGCATCGGTACAATaa                                          >  1:1048268/1‑150 (MQ=255)
                                                                                                          ttGAACGGTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAACTTTTCACCAGCGTCAGCTCGCCAGCATTACGCATCGGTACAATaaa                                         >  2:1590622/1‑150 (MQ=255)
                                                                                                           tGAACGGTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAACTTTTCACCCGCGTCAGCTCGCCAGCATTACGCATCGGTACATTAAAt                                        >  2:336262/1‑150 (MQ=255)
                                                                                                                ggTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAACTTTTCACCAGCGTCAGCTCGCCAGCATTACGCATCGGTACAATAAATGTTTc                                   <  1:321610/150‑1 (MQ=255)
                                                                                                                ggTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAACTTTTCACCAGCGTCAGCTCGCCAGCATTACGCATCGGTACAATAAATGTTTc                                   <  1:696319/150‑1 (MQ=255)
                                                                                                                ggTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAACTTTTCACCAGCGTCAGCTCGCCAGCATTACGCATCGGTACAATAAATGTTTc                                   <  1:866404/150‑1 (MQ=255)
                                                                                                                 gTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAACTTTTCACCAGCGTCAGCTCGCCAGCATTACGCATCGGTACAATAAATGTTTcc                                  >  1:157843/1‑150 (MQ=255)
                                                                                                                   tGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAACTTTTCACCAGCGTCAGCTCGCCAGCATTACGCATCGGTACAATAAATGTTTCCTg                                >  1:1073737/1‑150 (MQ=255)
                                                                                                                       ttAATCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAACTTTTCACCAGCGTCAGCTCGCCAGCATTACGCATCGGTACAATAAATGTTTCCTGTTtc                            <  2:422482/150‑1 (MQ=255)
                                                                                                                        tAATCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAACTTTTCACCAGCGTCAGCTCGCCAGCATTACGCATCGGTACAATAAATGTTTCCTGTTtct                           <  1:1279373/150‑1 (MQ=255)
                                                                                                                         aaTCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAACTTTTCACCAGCGTCAGCTCGCCAGCATTACGCATCGGTACAATAAATGTTTCCTGTTtctc                          >  2:1439920/1‑150 (MQ=255)
                                                                                                                          aTCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAACTTTTCACCAGCGTCCGCTCGCCAGCATTACGCCTCGGTACAATAAATGTTTCCTGTTTCTCa                         >  1:593565/1‑150 (MQ=255)
                                                                                                                          aTCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAACTTTTCACCAGCGTCAGCTCGCCAGCATTACGCATCGGTACAATAAATGTTTCCTGTTTCTCa                         >  2:412182/1‑150 (MQ=255)
                                                                                                                           ttcaagaaaaCACCCGATAACTTTCGCTATCGGGTGTTTTTATTGAACATCGACGGCGAAGGTACCCCCATTAATCACCAGTCAAAACTTTTCACCAGCGTCAGCTCGCCAGCATTACGCATCGGTACAATAAATGTTTCCTGTTTCTCa                         >  1:1218092/7‑150 (MQ=25)
                                                                                                                            cAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAACTTTTCACCAGCGTCAGCTCGCCCGCATTACGCATCGGTACAATAAATGTTTCCTGTTTCTCAtt                       >  2:127958/1‑150 (MQ=255)
                                                                                                                            cAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAACTTTTCACCAGCGTCAGCTCGCCAGCATTACGCATCGGTACAATAAATGTTTCCTGTTTCTCAtt                       >  2:133882/1‑150 (MQ=255)
                                                                                                                             actAAACACACCCGATAGCGAAAGTTATCGGGTGGTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAACTTTTCACCAGCGGCAGCTCGCCAGCATTACGCCTCGGTACAATAAATGTTTCCTGTTTCTCATTg                      <  2:424238/148‑1 (MQ=255)
                                                                                                                             aaTAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAACTTTTCACCAGCGTCAGCTCGCCAGCATTACGCATCGGTACAATAAATGTTTCCTGTTTCTCATTg                      <  2:426901/150‑1 (MQ=255)
                                                                                                                             aaTAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAACTTTTCACCAGCGTCAGCTCGCCAGCATTACGCATCGGTACAATAAATGTTTCCTGTTTCTCATTg                      >  1:1431267/1‑150 (MQ=255)
                                                                                                                             aaTAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAACTTTTCACCAGCGTCAGCTCGCCAGCATTACGCATCGGTACAATAAATGTTTCCTGTTTCTCATTg                      <  2:1149711/150‑1 (MQ=255)
                                                                                                                              aTAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAACTTTTCACCAGCGTCAGCTCGCCAGCATTACGCATCGGTACAATAAATGTTTCCTGTTTCTCATTGa                     <  2:268061/150‑1 (MQ=255)
                                                                                                                               tAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAACTTTTCACCAGCGTCAGCTCGCCAGCATTACGCATCGGTACAATAAATGTTTCCTGTTTCTCATTGAc                    >  2:187159/1‑150 (MQ=255)
                                                                                                                               tAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAACTTTTCACCAGCGTCAGCTCGCCAGCATTACGCATCGGTACAATAAATGTTTCCTGTTTCTCATTGAc                    >  1:13255/1‑150 (MQ=255)
                                                                                                                                aaaaaCACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAACTTTTCACCAGCGTCAGCTCGCCAGCATTACGCATCGGTACAATAAATGTTTCCTGTTTCTCATTGAcc                   <  1:572955/150‑1 (MQ=255)
                                                                                                                                aaaaaCACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAACTTTTCACCAGCGTCAGCTCGCCAGCATTACGCATCGGTACAATAAATGTTTCCTGTTTCTCATTGAcc                   <  1:89031/150‑1 (MQ=255)
                                                                                                                                    acacCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAACTTTTCACCAGCGTCAGCTCGCCAGCATTACGCATCGGTACAATAAATGTTTCCTGTTTCTCATTGACCGATc               <  1:60254/150‑1 (MQ=255)
                                                                                                                                       cccGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAACTTTTCACCAGCGTCAGCTCGCCAGCATTACGCATCGGTACAATAAATGTTTCCTGTTTCTCATTGACCGATCCtt            >  1:914829/1‑150 (MQ=255)
                                                                                                                                        ccGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAACTTTTCACCAGCGTCAGCTCGCCAGCATTACGCATCGGTACAATAAATGTTTCCTGTTTCTCATTGACCGATCCTTc           >  1:265260/1‑150 (MQ=255)
                                                                                                                                        ccGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAACTTTTCACCAGCGTCAGCTCGCCAGCATTACGCATCGGTACAATAAATGTTTCCTGTTTCTCATTGACCGATCCTTc           <  2:1046984/150‑1 (MQ=255)
                                                                                                                                         cGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAACTTTTCACCAGCGTCAGCTCGCCAGCATTACGCATCGGTACAATAAATGTTTCCTGTTTCTCATTGACCGATCCTTCa          >  2:1157411/1‑150 (MQ=255)
                                                                                                                                         cGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAACTTTTCACCAGCGTCAGCTCGCCAGCATTACGCATCGGTACAATAAATGATTCCTGTTTCTCATTGACCGATCCTTCa          >  2:454867/1‑150 (MQ=255)
                                                                                                                                           aTAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAACTTTTCACCAGCGTCAGCTCGCCAGCATTACGCATCGGTACAATAAATGTTTCCTGTTTCTCATTGACCGATCCTTCATc        <  1:335207/150‑1 (MQ=255)
                                                                                                                                             agagAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAACTTTTCACCAGCGTCAGCTCGCCAGCATTACGCATCTGTACAATAAATGTTTCCTGTTTCTCATTGACCGATCCTTCATCgg      <  1:561300/147‑1 (MQ=255)
                                                                                                                                             aGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAACTTTTCACCAGCGTCAGCTCGCCAGCATTACGCATCGGTACAATAAATGTTTCCTGTTTCTCATTGACCGATCCTTCATCgg      <  1:17392/150‑1 (MQ=255)
                                                                                                                                               cGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAACTTTTCACCAGCGTCAGCTCGCCAGCATTACGCATCGGTACAATAAATGTTTCCTGTTTCTCATTGACCGATCCTTCATCGGTg    >  1:1632248/1‑150 (MQ=255)
                                                                                                                                                 aaaGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAACTTTTCACCAGCGTCAGCTCGCCAGCATTACGCATCGGTACAATAAATGTTTCCTGTTTCTCATTGACCGATCCTTCATCGGTGAt  <  1:1353392/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GGTGCTCGTATTCTGGTCACACTATTACATGCCATGCAGGCACGCGATAAAACGCTGGGGCTGGCAACACTGTGCATTGGCGGCGGTCAGGGAATTGCGATGGTGATTGAACGGTTGAATTAATCAATAAAAACACCCGATAGCGAAAGTTATCGGGTGTTTTCTTGAACATCGACGGCGAAGGTAACCCCATTAATCACCAGTCAAAACTTTTCACCAGCGTCAGCTCGCCAGCATTACGCATCGGTACAATAAATGTTTCCTGTTTCTCATTGACCGATCCTTCATCGGTGAT  >  NC_000913/2327171‑2327465

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