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

GATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCC‑GTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGC‑GCTGGTTTAAATGACTATCGCTGTGCTTTTGTTC‑GCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAATAACCAT  >  NC_000913/514475‑514774
                                                                                                                                                        |                                                                                                                                                      
gaTCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCC‑GTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGC‑GCTGGTTTAAATGACTATCGCTGTGCTTTTGTTC‑GCGCAc                                                                                                                                                       <  2:1140352/151‑1 (MQ=255)
 aTCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCC‑GTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGC‑GCTGGTTTAAATGACTATCGCTGTGCTTTTGTTC‑GCGCAcc                                                                                                                                                      >  4:1653321/1‑151 (MQ=255)
   ccTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCC‑GTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGC‑GCTGGTTTAAATGACTATCGCTGTGCTTTTGTTC‑GCGCAccc                                                                                                                                                     <  3:423483/150‑1 (MQ=255)
   ccTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCC‑GTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGC‑GCTGGTTTAAATGACTATCGCTGTGCTTTTGTTC‑GCGCAccc                                                                                                                                                     <  4:1510135/150‑1 (MQ=255)
   ccTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCC‑GTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGC‑GCTGGTTTAAATGACTATCGCTGTGCTTTTGTTC‑GCGCAccc                                                                                                                                                     <  4:912210/150‑1 (MQ=255)
   ccTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCC‑GTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGC‑GCTGGTTTAAATGACTATCGCTGTGCTTTTGTTC‑GCGCACccg                                                                                                                                                    >  2:841081/1‑151 (MQ=255)
      cGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCC‑GTTGACCAGCGGAGATTCCAGAACAAAACGTCGGCCAATCAGCTGCTGCCCGC‑GCTGGTTTAAATGACTATCGCTGTGCTTTTGTTC‑GCGCACccgc                                                                                                                                                   <  4:167584/149‑1 (MQ=255)
       gCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCC‑GTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGC‑GCTGGTTTAAATGACTATCGCTGTGCTTTTGTTC‑GCGCACCcgccgc                                                                                                                                                >  1:186085/1‑151 (MQ=255)
               tGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCC‑GTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGC‑GCTGGTTTAAATGACTATCGCTGTGCTTTTGTTC‑GCGCACCCGCCGCGACAAcc                                                                                                                                         >  2:1511608/1‑150 (MQ=255)
                    ggCCATGAACTGTCACCGACGCGCATATGACCGCGCCC‑GTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGC‑GCTGGTTTAAATGACTATCGCTGTGCTTTTGTTC‑GCGCACCCGCCGCGACAACCATTTcc                                                                                                                                   >  4:588885/1‑151 (MQ=255)
                           aaCTGCCCCCGACGCGCATATGCCCGCCCCC‑GTTGCCCAGCGGAGATTCCAGCAAAAAACGCCGCCAAATCAGCTGCCCCCCGC‑GCTGTTTAAAATGACTATCGCTGTGTTTTTGTCC‑GCGCCCCCGCCGCGACAACCATTTCCACcaca                                                                                                                             <  4:186365/150‑1 (MQ=255)
                            aCTGTCACCGACGCGCATATGACCGCGCCC‑GTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGC‑GCTGGTTTAAATGACTATCGCTGTGCTTTTGTTC‑GCGCACCCGCCGCGACAACCATTTCCACCACAGc                                                                                                                           >  2:1534869/1‑151 (MQ=255)
                              tGTCACCGACGCGCATATGACCGCGCCC‑GTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGC‑GCTGGTTTAAATGACTATCGCTGTGCTTTTGTTC‑GCGCACCCGCCGCGACAACCATTTCCACCACAGcc                                                                                                                          >  2:648069/1‑150 (MQ=255)
                               gtcgcggcctcgtatatGACCGCCCCC‑TTTGACCTGGGAAGATTCCAGCACAAAACGGCGACAAACCGGCTGTCGCCGGC‑TTTGTTTAAATTAATTAGCGCTGTG‑TTTTGCTCGGCTCACCCGCCGCGACAACCATTTCCACCACAGCCa                                                                                                                         <  4:1247212/137‑1 (MQ=255)
                               gccACCCACGCGCATATGCCCGCCCCC‑GTTGACCGGCGGAGATTCCAGCACAAACCGTCGGCCAATCAGCTGCTGCCCGC‑GTTGGTTAAAATGACTATCGCTGTGTTTTTGTTC‑GCGCACCCGCCGCGACAACCATTTCCACCACAGCCa                                                                                                                         <  2:1516848/148‑1 (MQ=255)
                                   ccgccgCGCATAGGCCCGCGCCC‑GTGGACCAGGGGAGTTCCCACCACAAACCGTCGCCCAATCAGCTGCCCCCCGC‑GTTGTTTTAAATGACTATCGCTGTGTTTTTGTTC‑GCGCCCCCGCCGCGACAACCATTTCCACCACAGcc                                                                                                                          <  2:1051914/145‑1 (MQ=255)
                                     gACGCGCATATGACCGCGCCC‑GTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGC‑GCTGGTTTAAATGACTATCGCTGTGCTTTTGTTC‑GCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGgcggcg                                                                                                                  >  2:1053564/1‑151 (MQ=255)
                                        gcgcATATGACCGCGCCC‑GTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGC‑GCTGGTTTAAATGACTATCGCTGAGCTTTTGTTC‑GCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCTGCGagc                                                                                                               >  2:544755/1‑151 (MQ=255)
                                        gcgcATATGACCGCGCCC‑GTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCGGCCCGC‑GCTGGTTTAAATGACTATCGCTGTGCTTTTGTTC‑GCGCACCCGCCACGACAACCCTTTCCACCACAGCCAGGCGGCGagc                                                                                                               >  2:416511/1‑151 (MQ=255)
                                         cgcATATGCCCGCCCCC‑GTTGACCAGGGGAGATCCCAGCAAAAACCGCCGGCCAATCGCCTGCTGCCCGC‑GCTGGTTAAAATGACTACCGCTGTGTTTTTGTTC‑GCGCCCCCGCCGCGACAACCATTTCCACCACAGcc                                                                                                                          <  1:308993/139‑1 (MQ=255)
                                          gcATATGCCCGCGCCC‑GTTGACCAGCGGAGATTCCAGCAAAAAACGTCGGCCAATCAGCTGCTCCCCGC‑GTTGGTTAAAATGTCTATCGCTGTGTTTTTGTTC‑GCGCCCCCGCCGCGACAACCATTTCTACCACAGCCAGGCGGCGagc                                                                                                               <  4:1320528/149‑1 (MQ=255)
                                           cgatatgccGCGCCCGGTTGACCGGGGGAGATTCCGGCAAAAACCGTGGCCAAAGCTCCTGCT‑CCCGCGGCTGTTTTAAATGCCTCTCGTTGTGCTTTTTTTC‑GCCCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGagca                                                                                                              <  2:1016377/143‑1 (MQ=255)
                                            atatGCCCGCGCCG‑GTGGACCGGCGGAGTTCCCAGCAAAAAACGCCGCCAAATCAGCTGCTGCCCGC‑GCTGTTTTAAATAACTACCGCTGTGTTTTTGTCC‑GCGCACCCGCCGCGACAACCATTTccacc                                                                                                                                <  2:1331805/130‑1 (MQ=255)
                                              atGCCCGCCCCC‑GTTGCCCGGGGGAGATCCCACCACAAACCGCCGGCAAATCAGCTGCCCCCCGC‑GCTGTTTAAAATGATTACCGCTTTTTTTTTGTCC‑GCGCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGagca                                                                                                              <  2:1339697/146‑1 (MQ=255)
                                                gcccGCGCCC‑GTTGCCCAGCGGAGTTTCCAGCAAAAAACGTCGGCCAATCAGCTGCTCCCCGC‑GCTGGTTAAAATGACTATCGCTGTGTTTTTGTCC‑GCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGt                                                                                                          <  4:970897/146‑1 (MQ=255)
                                                  ccGCGCCG‑TTGGCCCAGCGGGAATCCCAGCACAAAACGGCGCCCAATCGGCTGCCGCCCGC‑GCGGGTTTAAAGGCCTAGCCGGGTGTTTTTGTTC‑GCACCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGa                                                                                                     <  2:1080572/151‑1 (MQ=255)
                                                  ccGCGCCC‑GTTGACGGGGGGGGTTCCCAGCAAAAACCGGCGGCAAATCGGCTCCTGCCGGC‑GGTGGTTTAATTGCTTATCGCTTTGTTTTGGTCC‑GCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGt                                                                                                          <  2:1213689/146‑1 (MQ=255)
                                                   cgcgCCC‑GTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGC‑GCTGGTTTAAATGACTATCGCTGTGCTTTTGTTC‑GCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGAt                                                                                                    >  3:38380/1‑151 (MQ=255)
                                                       ccc‑TTTGACCAGGGGAGATTCCAGCAAAAAACGTGGGCAAATCAGCTGCTCCCCGC‑GCTGGTTTAAATGACTATCGTTGTCTTTTTGTCC‑GCGCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATgg                                                                                                  <  4:746770/149‑1 (MQ=255)
                                                       ccc‑GTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGC‑GCTGGTTTAAATGACTATCGCTGTGCTTTTGTTC‑GCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGc                                                                                                 >  1:174467/1‑150 (MQ=255)
                                                            ttGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGC‑GCTGGTTTAAATGACTATCGCTGTGCTTTTGTTC‑GCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAAc                                                                                             >  4:1041895/1‑150 (MQ=255)
                                                            ttGACCAGCGGAGATCCCAGCACAAACCGTCGGCCAATCAGCTGCTGCCCGC‑GCTGGTTAAAATGACTTCCGCTGTTTTTTTGTCC‑GCCCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAAc                                                                                             <  3:109498/150‑1 (MQ=255)
                                                            ttGACCAGCGGAGATCCCAGCAAAAAACGCCGGCCAATCGGCTGCTGCCCGC‑GCTGTTTTAAATGCTTATCGCTGTGTTTTTGTCC‑GCGCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCaa                                                                                               <  2:1623407/148‑1 (MQ=255)
                                                              gACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGC‑GCTGGTTTAAATGACTATCGCTGTGCTTTTGTTC‑GCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACAtt                                                                                          >  4:606919/1‑151 (MQ=255)
                                                                ccAGCGGAGTTTCCAGCACAAACCGCCGGCCAATCAGCTGCTGCCCGC‑GCTGTTTTAAATGACTATCGCTGTGTTTTTGTCC‑GCGCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACAtt                                                                                          <  4:301965/149‑1 (MQ=255)
                                                                ccAGCGGAGATCCCAGCAAAAACCGTCGGCCAATCAGCTGCTGCCCGC‑GCTGGTTTAAATGCCTATCGCTGTTTTTTTGTTC‑GCGCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTa                                                                                         <  4:1085263/150‑1 (MQ=255)
                                                                   ggggAGATCCCACAAAAAAGCGGCGGCCATCCAGTTGCGCCCCCC‑GCTGTTTTAAATGACTATCGTTGTGTTTTTGTTC‑GCCCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACa                                                                                            <  4:700317/142‑1 (MQ=255)
                                                                   gCGGAGTTCCCAGCAAAAAACGTCGGCCAATCAGCGGCTCCCCGC‑GGTGGTTAAAATGACTACCGCTGTGTTTTTGTCC‑GCGCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTAccc                                                                                      <  4:39706/150‑1 (MQ=255)
                                                                     ggAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGC‑GCTGGTTTAAATGACTATCGCTGTGCTTTTGTTC‑GCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCtt                                                                                   >  4:1476534/1‑151 (MQ=255)
                                                                     ggAGATCCCAGCACAAACCGGCGGCAAATCGGCTGCGGCCCGC‑GCTGTTTTAAATGACTACCGCTGTTTTTTTGTCC‑GCCCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCtt                                                                                   <  2:1685393/151‑1 (MQ=255)
                                                                     ggAAATCCCACCAAAAAGCGTGGGCCAATCAGCTCCGGCCCGC‑GTTGTTTAAAATGACTATCGCTGTGTTTTTGTTC‑GCGCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCt                                                                                    <  1:476572/150‑1 (MQ=255)
                                                                     gaagaTCCCACCACAAACCGCCGCCCAACCGCCTGCCGCCCGC‑GGTGTTTTAAAGACCTACCGCTGGTTTTTTTTCC‑GCCCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCt                                                                                    <  1:30055/148‑1 (MQ=255)
                                                                      gagaTTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGC‑GCTGGTTTAAATGACTATCGCTGTGCTTTTGTTC‑GCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTg                                                                                  >  4:1602509/1‑151 (MQ=255)
                                                                       agaTTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGC‑GCTGGTTTAAATGACTATCGCTGTGCTTTTGTTC‑GCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGc                                                                                 >  4:1364277/1‑151 (MQ=255)
                                                                       agaTTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGC‑GCTGGTTTAAATGACTATCGCTGTGCTTTTGTTC‑GCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGc                                                                                 >  1:102289/1‑151 (MQ=255)
                                                                       agaTTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGC‑GCTGGTTTAAATGACTATCGCTGTGCTTTTGTTC‑GCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGc                                                                                 >  2:459068/1‑151 (MQ=255)
                                                                       agaTTCCAGCAAAAACCGTCGCCCAACCCGCTGCTCCCCGC‑GTTGTTTAAAATGACTATCGCTGTGTTTTTGTCC‑GCCCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTg                                                                                  <  2:399133/150‑1 (MQ=255)
                                                                       agaTCCCAGCAAAAAACGCCGCCCAATCAGCTGCTGCCCGC‑GCTGTTTAAAATGACTATCGTTGTGTTTTTGTTC‑GCGCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTg                                                                                  <  3:174481/150‑1 (MQ=255)
                                                                        gttcccAGCAAAAACCGCGGCCCAACCAGCTGCTGCCCGC‑GTTGGTTTAAATGATTATCGCTGTTTTTTTTTTC‑GCCCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGcc                                                                                <  4:726947/147‑1 (MQ=255)
                                                                        gaTTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGC‑GCTGGTTTAAATGACTATCGCTGTGCTTTTGTTC‑GCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGcc                                                                                >  2:1641276/1‑151 (MQ=255)
                                                                            ccAGCAAAAACCGCCGGCCAATCAGCTGCTCCCCGC‑GTTGTTTTAAATGACTATCGCTGTGTTTTTGTCC‑GCGCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCAc                                                                              <  4:717540/149‑1 (MQ=255)
                                                                              aGCACAAAGCGTCGCCCAACCGCCTGCTCCCCGC‑GCTGTTTAAAATGATTACCGCTGTTTTTTTGTGC‑GCGCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTccc                                                                          <  1:238916/151‑1 (MQ=255)
                                                                              aGCACAAAACGTCGGCCAATCAGCTGCTGCCCGC‑GCTGGTTTAAATGACTATCGCTGTGCTTTTGTTC‑GCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTccc                                                                          >  2:272632/1‑151 (MQ=255)
                                                                                caaaaacggccgcccAACCAGCTGCTCCCCGC‑GCGGGTTTAAATGCCTATCGCTTGTTTTTTGTTC‑GCGCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTccc                                                                          <  2:1107079/136‑1 (MQ=255)
                                                                                   aaacCGCCGGCAAATCCGCTGCTCCCCGC‑GCTGTTTAAAATCACTATCGCTGTGTTTTTGTCC‑GCCCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTcca                                                                          <  4:268870/146‑2 (MQ=255)
                                                                                   aaacCGCCGGCAAATCAGCTGCTCCCCGC‑GTTGTTTTAAATGCCTCTCGCTGTGTTTTTGTCC‑GCCCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAAcc                                                                      <  4:1466703/150‑1 (MQ=255)
                                                                                   aaaaCGTCGGCCAATCAGCTGCTGCCCGC‑GCTGTTTTAAATGACTATCGCTGTGCTTTTGTTC‑GCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAAcc                                                                      <  3:403010/150‑1 (MQ=255)
                                                                                     aaCGTCGGCCAATCAGCTGCTGCCCGC‑GCTGGTTTAAATGACTATCGCTGTGCTTTTGTTC‑GCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAAccc                                                                     >  1:247727/1‑149 (MQ=255)
                                                                                      aCGTCGGCCAATCAGCTGCTGCCCGC‑GCTGGTTTAAATGACTATCGCTGTGCTTTTGTTC‑GCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGc                                                                  >  3:389348/1‑151 (MQ=255)
                                                                                       cGTCGGCCAATCAGCTGCTGCCCGC‑GCTGGTTTAAATGACTATCGCTGTGCTTTTGTTC‑GCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCg                                                                 >  3:69250/1‑151 (MQ=255)
                                                                                        gTCGGCCAATCAGCTGCTGCCCGC‑GCTGGTTTAAATGACTATCGCTGTGCTTTTGTTC‑GCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCgg                                                                >  1:148229/1‑151 (MQ=255)
                                                                                         tCGGCCAATCAGCTGCTGCCCGC‑GCTGGTTTAAATGACTATCGCTGTGCTTTTGTTC‑GCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGc                                                               >  1:27455/1‑151 (MQ=255)
                                                                                           ggCCAATCAGCTGCTGCCCGC‑GCTGGTTTAAATGACTATCGCTGTGCTTTTGTTC‑GCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCAc                                                             >  3:331183/1‑151 (MQ=255)
                                                                                           ggCCAATCAGCTGCTGCCCGC‑GCTGGTTTAAATGACTATCGCTGTGCTTTTGTTC‑GCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCAc                                                             >  2:1273619/1‑151 (MQ=255)
                                                                                            gCCAATCAGCTGCTGCCCGC‑GCTGGTTTAAATGACTATCGCTGTGCTTTTGTTC‑GCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCAcc                                                            >  4:140358/1‑151 (MQ=255)
                                                                                             ccAATCCGCTCCCCCCCGC‑GTTGTTTAAAAGGCCGTCCGCTGTCTTTTGGTCC‑GCCCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCAcca                                                           <  2:191305/151‑1 (MQ=255)
                                                                                              cAATCAGCTGCTGCCCGC‑GCTGGTTTAAATGACTATCGCTGTGCTTTTGTTC‑GCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCAcc                                                            >  3:231768/1‑149 (MQ=255)
                                                                                              cAATCAGCTGCTGCCCGC‑GCTGGTTTAAATGACTATCGCTGTGCTTTTGTTC‑GCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCAcc                                                            >  4:1359559/1‑149 (MQ=255)
                                                                                                        ctcccggc‑gcTGTTTTAATTGACTCTCGCTGTGTTTTTGTCC‑GCCCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGacc                                                   <  4:265419/142‑1 (MQ=255)
                                                                                                           cccGC‑GCTGTTTTAAATGCCTACCGTTGTGTTTTTTTCC‑CCCCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGaccac                                                 <  4:979428/147‑1 (MQ=255)
                                                                                                           cccGC‑GCTGGTTTAAATGACTATCGCTTTGTTTTTGTCC‑GCGCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAg                                              <  3:216081/150‑1 (MQ=255)
                                                                                                            cccgcctttttttAAATGACTACCGCTGTCTTTTTTTGC‑GCCCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCAccagcca                                                       <  4:345611/131‑1 (MQ=255)
                                                                                                            ccGC‑GTTGTTTAAAATTCCTACCGCTGTGCTTTTGCCC‑CCCCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGc                                            <  2:1002865/151‑1 (MQ=255)
                                                                                                             ccgcgtttttttAAATGTTTTTCGCTGTTTTTTTGTTC‑GCCCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGc                                            <  4:1633034/142‑1 (MQ=255)
                                                                                                               c‑gcTGGTTTAAATGACTATCGCTGTGCTTTTGTTC‑GCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCgg                                         >  2:32149/1‑151 (MQ=255)
                                                                                                                  cTGGTTTAAATGACTATCGCTGTGCTTTTGTTC‑GCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGt                                        >  2:1604343/1‑150 (MQ=255)
                                                                                                                   tGGTTTAAATGACTATCGCTGTGCTTTTGTTC‑GCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTa                                       >  2:1140785/1‑150 (MQ=255)
                                                                                                                        tAAATGACTATCGCTGTGCTTTTGTTC‑GCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCAcc                                 >  2:1416109/1‑151 (MQ=255)
                                                                                                                         aaaTGACTATCGCTGTGCTTTTGTTC‑GCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCg                                >  2:1412835/1‑151 (MQ=255)
                                                                                                                          aatcccTATCGCTTTGTTTTTGTCC‑GCCCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGc                               <  4:356768/146‑1 (MQ=255)
                                                                                                                           atcaatATCGCTTTCGTTTCGTGC‑GCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGGGCAGCGTCAGGATGGCAAACATTACACCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCt                              <  2:514618/146‑1 (MQ=255)
                                                                                                                           aTGACTGTCGCGGTGTTTTTGTCC‑GCGCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCt                              <  2:1283778/151‑1 (MQ=255)
                                                                                                                                  tCGCTGTGTTTTTTTCC‑GCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCAcg                        <  3:244542/150‑1 (MQ=255)
                                                                                                                                   cgttttgtttttGTCC‑CCGCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCAcg                        <  4:920849/141‑1 (MQ=255)
                                                                                                                                   cggtgtttttttGTCC‑GCCCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCAc                         <  4:1056789/140‑1 (MQ=255)
                                                                                                                                    gCTGTGCTTTTGTTC‑GCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCAcgccg                     >  2:924438/1‑151 (MQ=255)
                                                                                                                                        tgtttttGTCC‑GCGCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCAcg                        <  3:472109/141‑1 (MQ=255)
                                                                                                                                          tttttGTGC‑GCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCAcgc                       <  1:294564/142‑1 (MQ=255)
                                                                                                                                           ttttttCC‑GCCCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCcac               <  3:246826/150‑1 (MQ=255)
                                                                                                                                           ttttGTCC‑GCCCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCAcgccg                     <  1:377358/144‑1 (MQ=255)
                                                                                                                                           ttttGTCC‑GCCCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCAcgcc                      <  4:423211/143‑1 (MQ=255)
                                                                                                                                            tttGTCC‑GCCCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCcaca              <  2:1326918/150‑1 (MQ=255)
                                                                                                                                              tGTTC‑GCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCacaac            >  3:328551/1‑150 (MQ=255)
                                                                                                                                                  c‑gcgCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAAt        >  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)
                                                                                                                                                        |                                                                                                                                                      
GATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCGCATATGACCGCGCCC‑GTTGACCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGC‑GCTGGTTTAAATGACTATCGCTGTGCTTTTGTTC‑GCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAATAACCAT  >  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: ‑