Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913514,6240AC40.6% 174.2 / 11.6 96V79G (GTG→GGG) ybbJinner membrane protein; stimulator of the QmcA suppressor of ftsH‑htpX
Reads supporting (aligned to +/- strand):  ref base A (40/17);  new base C (0/39);  total (40/56)
Fisher's exact test for biased strand distribution p-value = 9.41e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAATAACCAT  >  NC_000913/514475‑514774
                                                                                                                                                        |                                                                                                                                                      
gaTCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCAc                                                                                                                                                       <  2:1140352/151‑1 (MQ=255)
 aTCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCAcc                                                                                                                                                      >  4:1653321/1‑151 (MQ=255)
   ccTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCAccc                                                                                                                                                     <  3:423483/150‑1 (MQ=255)
   ccTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCAccc                                                                                                                                                     <  4:1510135/150‑1 (MQ=255)
   ccTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCAccc                                                                                                                                                     <  4:912210/150‑1 (MQ=255)
   ccTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACccg                                                                                                                                                    >  2:841081/1‑151 (MQ=255)
      cGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGAACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACccgc                                                                                                                                                   <  4:167584/149‑1 (MQ=255)
       gCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACCcgccgc                                                                                                                                                >  1:186085/1‑151 (MQ=255)
               tGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAAcc                                                                                                                                         >  2:1511608/1‑150 (MQ=255)
                    ggCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTcc                                                                                                                                   >  4:588885/1‑151 (MQ=255)
                           aaCTGCCCCCGACGCGCATATGCCCGCCCCCGTTGCCCAGCGGAGATTCCAGCAAAAAACGCCGCCAAATCAGCTGCCCCCCGCGCTGTTTAAAATGACTATCGCTGTGTTTTTGTCCGCGCCCCCGCCGCGACAACCATTTCCACcaca                                                                                                                             <  4:186365/150‑1 (MQ=255)
                            aCTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGc                                                                                                                           >  2:1534869/1‑151 (MQ=255)
                              tGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGcc                                                                                                                          >  2:648069/1‑150 (MQ=255)
                               gtcgcggcctcgtatatGACCGCCCCCTTTGACCTGGGAAGATTCCAGCACAAAACGGCGACAAACCGGCTGTCGCCGGCTTTGTTTAAATTAATTAGCGCTGTGTTTTGCTCGGCTCACCCGCCGCGACAACCATTTCCACCACAGCCa                                                                                                                         <  4:1247212/137‑1 (MQ=255)
                               gccACCCACGCGCATATGCCCGCCCCCGTTGACCGGCGGAGATTCCAGCACAAACCGTCGGCCAATCAGCTGCTGCCCGCGTTGGTTAAAATGACTATCGCTGTGTTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCa                                                                                                                         <  2:1516848/148‑1 (MQ=255)
                                   ccgccgCGCATAGGCCCGCGCCCGTGGACCAGGGGAGTTCCCACCACAAACCGTCGCCCAATCAGCTGCCCCCCGCGTTGTTTTAAATGACTATCGCTGTGTTTTTGTTCGCGCCCCCGCCGCGACAACCATTTCCACCACAGcc                                                                                                                          <  2:1051914/145‑1 (MQ=255)
                                     gACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGgcggcg                                                                                                                  >  2:1053564/1‑151 (MQ=255)
                                        gcgcATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGAGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCTGCGagc                                                                                                               >  2:544755/1‑151 (MQ=255)
                                        gcgcATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCGGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCACGACAACCCTTTCCACCACAGCCAGGCGGCGagc                                                                                                               >  2:416511/1‑151 (MQ=255)
                                         cgcATATGCCCGCCCCCGTTGACCAGGGGAGATCCCAGCAAAAACCGCCGGCCAATCGCCTGCTGCCCGCGCTGGTTAAAATGACTACCGCTGTGTTTTTGTTCGCGCCCCCGCCGCGACAACCATTTCCACCACAGcc                                                                                                                          <  1:308993/139‑1 (MQ=255)
                                          gcATATGCCCGCGCCCGTTGACCAGCGGAGATTCCAGCAAAAAACGTCGGCCAATCAGCTGCTCCCCGCGTTGGTTAAAATGTCTATCGCTGTGTTTTTGTTCGCGCCCCCGCCGCGACAACCATTTCTACCACAGCCAGGCGGCGagc                                                                                                               <  4:1320528/149‑1 (MQ=255)
                                           cgatatgccGCGCCCGGTTGACCGGGGGAGATTCCGGCAAAAACCGTGGCCAAAGCTCCTGCTCCCGCGGCTGTTTTAAATGCCTCTCGTTGTGCTTTTTTTCGCCCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGagca                                                                                                              <  2:1016377/143‑1 (MQ=255)
                                            atatGCCCGCGCCGGTGGACCGGCGGAGTTCCCAGCAAAAAACGCCGCCAAATCAGCTGCTGCCCGCGCTGTTTTAAATAACTACCGCTGTGTTTTTGTCCGCGCACCCGCCGCGACAACCATTTccacc                                                                                                                                <  2:1331805/130‑1 (MQ=255)
                                              atGCCCGCCCCCGTTGCCCGGGGGAGATCCCACCACAAACCGCCGGCAAATCAGCTGCCCCCCGCGCTGTTTAAAATGATTACCGCTTTTTTTTTGTCCGCGCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGagca                                                                                                              <  2:1339697/146‑1 (MQ=255)
                                                gcccGCGCCCGTTGCCCAGCGGAGTTTCCAGCAAAAAACGTCGGCCAATCAGCTGCTCCCCGCGCTGGTTAAAATGACTATCGCTGTGTTTTTGTCCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGt                                                                                                          <  4:970897/146‑1 (MQ=255)
                                                  ccGCGCCGTTGGCCCAGCGGGAATCCCAGCACAAAACGGCGCCCAATCGGCTGCCGCCCGCGCGGGTTTAAAGGCCTAGCCGGGTGTTTTTGTTCGCACCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGa                                                                                                     <  2:1080572/151‑1 (MQ=255)
                                                  ccGCGCCCGTTGACGGGGGGGGTTCCCAGCAAAAACCGGCGGCAAATCGGCTCCTGCCGGCGGTGGTTTAATTGCTTATCGCTTTGTTTTGGTCCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGt                                                                                                          <  2:1213689/146‑1 (MQ=255)
                                                   cgcgCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGAt                                                                                                    >  3:38380/1‑151 (MQ=255)
                                                       cccTTTGACCAGGGGAGATTCCAGCAAAAAACGTGGGCAAATCAGCTGCTCCCCGCGCTGGTTTAAATGACTATCGTTGTCTTTTTGTCCGCGCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATgg                                                                                                  <  4:746770/149‑1 (MQ=255)
                                                       cccGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGc                                                                                                 >  1:174467/1‑150 (MQ=255)
                                                            ttGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAAc                                                                                             >  4:1041895/1‑150 (MQ=255)
                                                            ttGACCAGCGGAGATCCCAGCACAAACCGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTAAAATGACTTCCGCTGTTTTTTTGTCCGCCCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAAc                                                                                             <  3:109498/150‑1 (MQ=255)
                                                            ttGACCAGCGGAGATCCCAGCAAAAAACGCCGGCCAATCGGCTGCTGCCCGCGCTGTTTTAAATGCTTATCGCTGTGTTTTTGTCCGCGCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCaa                                                                                               <  2:1623407/148‑1 (MQ=255)
                                                              gACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACAtt                                                                                          >  4:606919/1‑151 (MQ=255)
                                                                ccAGCGGAGTTTCCAGCACAAACCGCCGGCCAATCAGCTGCTGCCCGCGCTGTTTTAAATGACTATCGCTGTGTTTTTGTCCGCGCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACAtt                                                                                          <  4:301965/149‑1 (MQ=255)
                                                                ccAGCGGAGATCCCAGCAAAAACCGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGCCTATCGCTGTTTTTTTGTTCGCGCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTa                                                                                         <  4:1085263/150‑1 (MQ=255)
                                                                   ggggAGATCCCACAAAAAAGCGGCGGCCATCCAGTTGCGCCCCCCGCTGTTTTAAATGACTATCGTTGTGTTTTTGTTCGCCCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACa                                                                                            <  4:700317/142‑1 (MQ=255)
                                                                   gCGGAGTTCCCAGCAAAAAACGTCGGCCAATCAGCGGCTCCCCGCGGTGGTTAAAATGACTACCGCTGTGTTTTTGTCCGCGCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTAccc                                                                                      <  4:39706/150‑1 (MQ=255)
                                                                     ggAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCtt                                                                                   >  4:1476534/1‑151 (MQ=255)
                                                                     ggAGATCCCAGCACAAACCGGCGGCAAATCGGCTGCGGCCCGCGCTGTTTTAAATGACTACCGCTGTTTTTTTGTCCGCCCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCtt                                                                                   <  2:1685393/151‑1 (MQ=255)
                                                                     ggAAATCCCACCAAAAAGCGTGGGCCAATCAGCTCCGGCCCGCGTTGTTTAAAATGACTATCGCTGTGTTTTTGTTCGCGCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCt                                                                                    <  1:476572/150‑1 (MQ=255)
                                                                     gaagaTCCCACCACAAACCGCCGCCCAACCGCCTGCCGCCCGCGGTGTTTTAAAGACCTACCGCTGGTTTTTTTTCCGCCCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCt                                                                                    <  1:30055/148‑1 (MQ=255)
                                                                      gagaTTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTg                                                                                  >  4:1602509/1‑151 (MQ=255)
                                                                       agaTTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGc                                                                                 >  4:1364277/1‑151 (MQ=255)
                                                                       agaTTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGc                                                                                 >  1:102289/1‑151 (MQ=255)
                                                                       agaTTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGc                                                                                 >  2:459068/1‑151 (MQ=255)
                                                                       agaTTCCAGCAAAAACCGTCGCCCAACCCGCTGCTCCCCGCGTTGTTTAAAATGACTATCGCTGTGTTTTTGTCCGCCCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTg                                                                                  <  2:399133/150‑1 (MQ=255)
                                                                       agaTCCCAGCAAAAAACGCCGCCCAATCAGCTGCTGCCCGCGCTGTTTAAAATGACTATCGTTGTGTTTTTGTTCGCGCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTg                                                                                  <  3:174481/150‑1 (MQ=255)
                                                                        gttcccAGCAAAAACCGCGGCCCAACCAGCTGCTGCCCGCGTTGGTTTAAATGATTATCGCTGTTTTTTTTTTCGCCCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGcc                                                                                <  4:726947/147‑1 (MQ=255)
                                                                        gaTTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGcc                                                                                >  2:1641276/1‑151 (MQ=255)
                                                                            ccAGCAAAAACCGCCGGCCAATCAGCTGCTCCCCGCGTTGTTTTAAATGACTATCGCTGTGTTTTTGTCCGCGCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCAc                                                                              <  4:717540/149‑1 (MQ=255)
                                                                              aGCACAAAGCGTCGCCCAACCGCCTGCTCCCCGCGCTGTTTAAAATGATTACCGCTGTTTTTTTGTGCGCGCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTccc                                                                          <  1:238916/151‑1 (MQ=255)
                                                                              aGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTccc                                                                          >  2:272632/1‑151 (MQ=255)
                                                                                caaaaacggccgcccAACCAGCTGCTCCCCGCGCGGGTTTAAATGCCTATCGCTTGTTTTTTGTTCGCGCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTccc                                                                          <  2:1107079/136‑1 (MQ=255)
                                                                                   aaacCGCCGGCAAATCCGCTGCTCCCCGCGCTGTTTAAAATCACTATCGCTGTGTTTTTGTCCGCCCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTcca                                                                          <  4:268870/146‑2 (MQ=255)
                                                                                   aaacCGCCGGCAAATCAGCTGCTCCCCGCGTTGTTTTAAATGCCTCTCGCTGTGTTTTTGTCCGCCCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAAcc                                                                      <  4:1466703/150‑1 (MQ=255)
                                                                                   aaaaCGTCGGCCAATCAGCTGCTGCCCGCGCTGTTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAAcc                                                                      <  3:403010/150‑1 (MQ=255)
                                                                                     aaCGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAAccc                                                                     >  1:247727/1‑149 (MQ=255)
                                                                                      aCGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGc                                                                  >  3:389348/1‑151 (MQ=255)
                                                                                       cGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCg                                                                 >  3:69250/1‑151 (MQ=255)
                                                                                        gTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCgg                                                                >  1:148229/1‑151 (MQ=255)
                                                                                         tCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGc                                                               >  1:27455/1‑151 (MQ=255)
                                                                                           ggCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCAc                                                             >  3:331183/1‑151 (MQ=255)
                                                                                           ggCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCAc                                                             >  2:1273619/1‑151 (MQ=255)
                                                                                            gCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCAcc                                                            >  4:140358/1‑151 (MQ=255)
                                                                                             ccAATCCGCTCCCCCCCGCGTTGTTTAAAAGGCCGTCCGCTGTCTTTTGGTCCGCCCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCAcca                                                           <  2:191305/151‑1 (MQ=255)
                                                                                              cAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCAcc                                                            >  3:231768/1‑149 (MQ=255)
                                                                                              cAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCAcc                                                            >  4:1359559/1‑149 (MQ=255)
                                                                                                        ctcccggcgcTGTTTTAATTGACTCTCGCTGTGTTTTTGTCCGCCCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGacc                                                   <  4:265419/142‑1 (MQ=255)
                                                                                                           cccGCGCTGTTTTAAATGCCTACCGTTGTGTTTTTTTCCCCCCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGaccac                                                 <  4:979428/147‑1 (MQ=255)
                                                                                                           cccGCGCTGGTTTAAATGACTATCGCTTTGTTTTTGTCCGCGCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAg                                              <  3:216081/150‑1 (MQ=255)
                                                                                                            cccgcctttttttAAATGACTACCGCTGTCTTTTTTTGCGCCCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCAccagcca                                                       <  4:345611/131‑1 (MQ=255)
                                                                                                            ccGCGTTGTTTAAAATTCCTACCGCTGTGCTTTTGCCCCCCCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGc                                            <  2:1002865/151‑1 (MQ=255)
                                                                                                             ccgcgtttttttAAATGTTTTTCGCTGTTTTTTTGTTCGCCCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGc                                            <  4:1633034/142‑1 (MQ=255)
                                                                                                               cgcTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCgg                                         >  2:32149/1‑151 (MQ=255)
                                                                                                                  cTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGt                                        >  2:1604343/1‑150 (MQ=255)
                                                                                                                   tGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTa                                       >  2:1140785/1‑150 (MQ=255)
                                                                                                                        tAAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCAcc                                 >  2:1416109/1‑151 (MQ=255)
                                                                                                                         aaaTGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCg                                >  2:1412835/1‑151 (MQ=255)
                                                                                                                          aatcccTATCGCTTTGTTTTTGTCCGCCCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGc                               <  4:356768/146‑1 (MQ=255)
                                                                                                                           atcaatATCGCTTTCGTTTCGTGCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGGGCAGCGTCAGGATGGCAAACATTACACCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCt                              <  2:514618/146‑1 (MQ=255)
                                                                                                                           aTGACTGTCGCGGTGTTTTTGTCCGCGCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCt                              <  2:1283778/151‑1 (MQ=255)
                                                                                                                                  tCGCTGTGTTTTTTTCCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCAcg                        <  3:244542/150‑1 (MQ=255)
                                                                                                                                   cgttttgtttttGTCCCCGCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCAcg                        <  4:920849/141‑1 (MQ=255)
                                                                                                                                   cggtgtttttttGTCCGCCCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCAc                         <  4:1056789/140‑1 (MQ=255)
                                                                                                                                    gCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCAcgccg                     >  2:924438/1‑151 (MQ=255)
                                                                                                                                        tgtttttGTCCGCGCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCAcg                        <  3:472109/141‑1 (MQ=255)
                                                                                                                                          tttttGTGCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCAcgc                       <  1:294564/142‑1 (MQ=255)
                                                                                                                                           ttttttCCGCCCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCcac               <  3:246826/150‑1 (MQ=255)
                                                                                                                                           ttttGTCCGCCCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCAcgccg                     <  1:377358/144‑1 (MQ=255)
                                                                                                                                           ttttGTCCGCCCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCAcgcc                      <  4:423211/143‑1 (MQ=255)
                                                                                                                                            tttGTCCGCCCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCcaca              <  2:1326918/150‑1 (MQ=255)
                                                                                                                                              tGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCacaac            >  3:328551/1‑150 (MQ=255)
                                                                                                                                                  cgcgCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAAt        >  4:1357045/1‑150 (MQ=255)
                                                                                                                                                    gcgcACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAATa       >  4:1661439/1‑150 (MQ=255)
                                                                                                                                                     cgcCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACaaa         <  4:459713/147‑1 (MQ=255)
                                                                                                                                                        aCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAATAACCAt  >  4:1044730/1‑151 (MQ=255)
                                                                                                                                                        aCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAATAACCAt  >  2:1011877/1‑151 (MQ=255)
                                                                                                                                                        aCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAATAACCAt  >  1:423416/1‑151 (MQ=255)
                                                                                                                                                        aCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAATAACCAt  >  1:404117/1‑151 (MQ=255)
                                                                                                                                                        |                                                                                                                                                      
GATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCCGTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTTAAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAATAACCAT  >  NC_000913/514475‑514774

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