Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,678,0390TG42.6% 100.5 / 15.3 61N108T (AAC→ACC) yphEputative sugar ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (6/29);  new base G (26/0);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 1.37e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTC  >  NC_000913/2677892‑2678184
                                                                                                                                                   |                                                                                                                                                  
aGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGtttt                                                                                                                                                 <  1:503876/151‑1 (MQ=255)
  ctgctTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGgttt                                                                                                                                                 >  2:615253/1‑149 (MQ=255)
    gctTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGc                                                                                                                                             <  2:771096/151‑1 (MQ=255)
    gctTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGc                                                                                                                                             <  1:254978/151‑1 (MQ=255)
    gctTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCCAGGCAGGGGgttt                                                                                                                                                 >  1:1121359/1‑147 (MQ=255)
      ttttGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGGCCGAGGCAGAGGGTTTTCGcc                                                                                                                                            >  1:322626/1‑150 (MQ=255)
       tttGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCa                                                                                                                                           <  2:300493/150‑1 (MQ=255)
        ttGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACTGACCGAGGGAGAGGGTTTTCGCCAc                                                                                                                                          >  2:756549/1‑150 (MQ=255)
               gggCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCGTCCGAGGGAGAGGGTTTCTTCCCCTGTc                                                                                                                                      >  2:577283/1‑147 (MQ=255)
                 gCTTAGCGTTGAAACAAGTTTTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCt                                                                                                                                <  2:735664/151‑1 (MQ=255)
                 gCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCt                                                                                                                                <  2:288148/151‑1 (MQ=255)
                         ttGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTccacca                                                                                                                         <  2:1020437/150‑1 (MQ=255)
                          tGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGCGGCCGAGGGTTTCCGCCACTGTCAGCCCTTccccc                                                                                                                          >  2:932607/1‑148 (MQ=255)
                          tGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAg                                                                                                                        <  2:999719/150‑1 (MQ=255)
                           gAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGa                                                                                                                       <  1:615488/150‑1 (MQ=255)
                           gAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTCACCCAGGGAGAGGGTTTCCCCCACCCTCCGCCCCTCCCCCCGa                                                                                                                       >  2:5447/1‑150 (MQ=255)
                            aaccAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCAGCGCCTTTTGCAGGTAATCACCCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAc                                                                                                                      <  1:299211/147‑1 (MQ=255)
                                aaGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCCAACCAGGCGGAGGGGTTCCGCCCCCATCAGCCCCTCCACCCGACTTAAc                                                                                                                 >  2:859383/1‑151 (MQ=255)
                                  gttgttCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGTGCCCCCCCCCCCGGGAGGGGGTTTTTCCCCCTGTCAGTCCTTCCCCCAGaagcaa                                                                                                                  >  1:732918/1‑143 (MQ=255)
                                     gttCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCCCCGACCGGGGCAGAGGGTTTTCGCCCCCGTCCGCCCCTCCCCCCGACTTAACTcc                                                                                                              >  2:190418/1‑149 (MQ=255)
                                     gttCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCCCCGACCGAGGGAGAGGGTTTTCGCCACCGTCAGCCCCTCCCCCCCCCCTAACTcc                                                                                                              >  1:902073/1‑149 (MQ=255)
                                      ttCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACCGACCCAGGCCGAGGGTTTTCGCCCCTGTCAGCCCCTCCCCCAGACCTAACTcc                                                                                                              >  1:1223034/1‑148 (MQ=255)
                                       tCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGa                                                                                                           <  1:655198/150‑1 (MQ=255)
                                        cAGGACTAACGCCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCGGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTcc                                                                                                              <  1:831653/146‑1 (MQ=255)
                                          ggACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATaaa                                                                                                       <  1:1029466/151‑1 (MQ=255)
                                           gACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACCGACCGAGGCAGAGGGTTTTCGCCACTGTCAGCCCCTCCACCCCACTTAAATCCTGATaaa                                                                                                       >  2:213740/1‑150 (MQ=255)
                                             cTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCGAGGGAGAGGGTTTTCGCCCCTGTCAACCCCTCCCCCCGACCTAACTCCTGATaaa                                                                                                       >  2:1246047/1‑148 (MQ=255)
                                              tAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCGACCGGGGCAGAGGGTTTTTCCCCCCGTTAGGCCCTCCACCCGACTTAACACCTCATaaaac                                                                                                     >  1:1087210/1‑147 (MQ=255)
                                               aaCGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACTGACCGAGGCAGAGGGTTTTCGCCCCTGTCAGCCCTTCCACCAGACTTAAATCCTGATaaaaccc                                                                                                   >  2:706446/1‑146 (MQ=255)
                                               aaCGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACcgcg                                                                                                  <  2:1041228/151‑1 (MQ=255)
                                                   tCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCCACCCGGGCAGAGGGTTTCTGCCCCCGTTAGCCCCTCCACCCGACATTAATCCTCGTAAAAcgcgcg                                                                                                >  1:248606/1‑149 (MQ=255)
                                                     aaCGCCCAGCGCTTGTAAGCAACGTTGAGCATCCGGCGCCATTGGCAGGTAATCAATCAGGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCgcgcgc                                                                                               <  2:675548/148‑1 (MQ=255)
                                                         cccAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGGTTTTCCCCCCCCTCAGCCCTTCCACCCGACATAAATCCTGATAAAACGCGCGCcccccc                                                                                         >  2:775370/1‑150 (MQ=255)
                                                           cAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCCCGGCGGCGGGTTTTCGCCACTGTCAGCCCTTCCACCAGACTTAACTTCTGATAAAACACGCcccccccc                                                                                         >  2:153801/1‑148 (MQ=255)
                                                               gcCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCCAGGCAGAGGGTTTTCCCCCCCGTCAGCCCTTCCACCCGACTTAACTCCTGATAAACCGCGCGCCCCCCCAGCTc                                                                                    >  2:457331/1‑149 (MQ=255)
                                                                ccTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAg                                                                                  <  2:672590/150‑1 (MQ=255)
                                                                   gTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCGAGGCAGAGGGTTTTCGCCCCTGTCCGCCCTTCCCCCAGACCTAACTCCTGATAAAACGCGCGCCCCCCCCCTTTAGCgg                                                                               >  1:960939/1‑150 (MQ=255)
                                                                     aggCAACGTTGGGCATCCTGCGCCATTTGCAGGGAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGAGAAACCGCGCGCACCCCCAGTTCAGCGGCAc                                                                            <  1:203265/149‑1 (MQ=255)
                                                                     aaGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTCACCGAGGCAGAGGGTTTCCCCCACCGTTAGCCCTTCCACCCCACTTAACACCTGATAAAACGCGCGCCCCCCCCGCTCAGCgg                                                                               >  2:256800/1‑148 (MQ=255)
                                                                              tGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCCAGGCAGAGGGTTTTCGCCACTGTCAGCCCCTCCACCCGACTTAACTCCTGATAAAACCCCCGCCCCCCCCGTTCAGCGGCAccg                                                                          >  2:955463/1‑142 (MQ=255)
                                                                               gggCATCCTGCGCCATTTGCAGTTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCa                                                                   <  1:437010/150‑1 (MQ=255)
                                                                               gggCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACCCACCGAGGGGGAGGGTTTTCGCCCCCGTCATCCCCTCCCCCCGACTTAACTCCTGATAAAACGCGCGCACCCCCCCCTCAGCGGCGCGGCGAgcc                                                                    >  1:278365/1‑147 (MQ=255)
                                                                               gggCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCa                                                                   <  2:900345/150‑1 (MQ=255)
                                                                                   aTCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCCCTGACCGAGGGAGAGGGTTTTCCCCCCTCTCCGACCCTCCCCCCCACCTAAACCCCCATAAAACCCGCGCCCCCCCCCTTCAGCGGCCCGGCggggcc                                                                   >  1:761970/1‑141 (MQ=255)
                                                                                        gcgcCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACCGAACGAGGGCGGGGGTTTTCGCCCCCGTCCGCCCTTCCACCCGACTTAACTCCTGATAAAACGCCCCCCCCCCCCCGTCAGCGGCACCGCGGGGCAGCGGAGcct                                                          >  1:98353/1‑148 (MQ=255)
                                                                                         cgcCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCCACCGAGGCCGAGGGTTTTCGCCCCCCTCAGGCCCTCCACCACACTTCACTCCTGTTAAAACGCGCGCCCCCCCCGTTTATCCGCACGGCGAGGCAGCGGAGcgt                                                          >  1:1124676/1‑147 (MQ=255)
                                                                                           ccATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTc                                                     <  2:1138761/151‑1 (MQ=255)
                                                                                            cATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCa                                                    <  2:906032/151‑1 (MQ=255)
                                                                                                       aatcaatcATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTc                                          <  2:1134096/150‑1 (MQ=255)
                                                                                                        atcaatcaTGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTcgc                                        <  2:108451/151‑1 (MQ=255)
                                                                                                           attcaTGCCGTTGCGGCGGGGCCACTGCCCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCGCCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTcgcg                                       <  1:970231/147‑1 (MQ=255)
                                                                                                           aatcaTGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGtc                                     <  2:6548/151‑1 (MQ=255)
                                                                                                             tcaTGCCGTTGCGGCGGGGCCGCTGCCCCAGGGAGAGGTTTTTCGCCCCCGTCAGGCCTTCCCCCCCACTTCTCCCCTGCTTAAACCCGCGCACCCCCCGTTCTGCGCCACGGCGCGGCCGCGTAGCTTTGTCCACCTTccctcccgtgtc                                    >  1:579857/1‑141 (MQ=255)
                                                                                                                 gCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGcc                                <  2:902121/150‑1 (MQ=255)
                                                                                                                  ccGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCaa                              <  1:924621/151‑1 (MQ=255)
                                                                                                                      tGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATcc                           <  2:79008/150‑1 (MQ=255)
                                                                                                                       gcggcgGGGCCACTGACCGAGGCAGAGGGTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAAACGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGGCTCGCCAATc                            >  2:246447/1‑148 (MQ=255)
                                                                                                                               gCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCAcc                 >  1:1229726/1‑151 (MQ=255)
                                                                                                                                  aCTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGc               >  1:948570/1‑150 (MQ=255)
                                                                                                                                       ccGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCg         >  1:692207/1‑151 (MQ=255)
                                                                                                                                         gAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCggg       <  2:547907/151‑1 (MQ=255)
                                                                                                                                          aGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGc      >  2:496487/1‑151 (MQ=255)
                                                                                                                                           ggCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGc      <  1:538993/150‑1 (MQ=255)
                                                                                                                                              agagGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTc  >  1:403610/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTC  >  NC_000913/2677892‑2678184

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