Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,678,0390TG47.0% 95.4 / 18.2 66N108T (AAC→ACC) yphEputative sugar ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (4/31);  new base G (31/0);  total (35/31)
Fisher's exact test for biased strand distribution p-value = 8.34e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCT  >  NC_000913/2677890‑2678187
                                                                                                                                                     |                                                                                                                                                     
ccaGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGt                                                                                                                                                       <  2:57558/150‑1 (MQ=255)
 caGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTCAACGAGGCAGAGGgtt                                                                                                                                                     >  1:126000/1‑151 (MQ=255)
 caGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCGAGGGAGAGGgtt                                                                                                                                                     >  1:604175/1‑151 (MQ=255)
 caGCTGCTTTTGCGCCGGGCTTAGCGCTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACCGAACCAGGCAGAGGgtt                                                                                                                                                     >  2:1672415/1‑151 (MQ=255)
  aGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGtttt                                                                                                                                                    <  1:194511/151‑1 (MQ=255)
   gctgctTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCCCCGACCGAGGCCGAGGgtttt                                                                                                                                                   >  1:88861/1‑150 (MQ=255)
   gctgctTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCCAGGGCGAGGgttt                                                                                                                                                    >  1:647405/1‑150 (MQ=255)
            gcgcCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCGACCCAGGGAGAGGgtttttccccc                                                                                                                                             >  1:1668369/1‑141 (MQ=255)
             cgcCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTc                                                                                                                                         <  2:242136/151‑1 (MQ=255)
              gcCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCCCCGGCCCAGGCAGAGGGGTTTCCCCACTGt                                                                                                                                          >  1:1149027/1‑149 (MQ=255)
               ccGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAg                                                                                                                                       <  2:402192/151‑1 (MQ=255)
                cGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGc                                                                                                                                     <  2:339636/151‑1 (MQ=255)
                  ggCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCCCGGCCGAGGGTTTTCGCCCCCGTCCGccc                                                                                                                                   >  1:1579736/1‑151 (MQ=255)
                       aGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTcc                                                                                                                               <  2:57328/150‑1 (MQ=255)
                         cGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACTGACCGAGGCAGAGGGTTTTCGCCACTGTCAGCCCTTccacc                                                                                                                            >  2:178646/1‑151 (MQ=255)
                         cGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTTACCGAGGCAGAGGGTTTTCGCCACTGTCCGCCCCTccccc                                                                                                                            >  2:221020/1‑151 (MQ=255)
                           ttGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGGTTTTCCCCCCCGTCAGCCCTTccccc                                                                                                                            >  2:818019/1‑149 (MQ=255)
                            tGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAg                                                                                                                          <  2:251603/150‑1 (MQ=255)
                              aaaCAAGTTGTTCAGGCCTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAg                                                                                                                          <  1:907814/148‑1 (MQ=255)
                                aCAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGCGGCAGAGGGTTTTCGCCCCTGTCAGCCCCTCCACCCGACtt                                                                                                                      >  2:1033902/1‑150 (MQ=255)
                                aCAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTGCCGCCACTGTCAGCCCTTCCACCAGACTTa                                                                                                                     <  1:300584/151‑1 (MQ=255)
                                 cAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTaa                                                                                                                    <  2:1594258/151‑1 (MQ=255)
                                  aaGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTaa                                                                                                                    <  2:1231217/150‑1 (MQ=255)
                                   aGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACt                                                                                                                  <  1:1127095/151‑1 (MQ=255)
                                   aGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGGTTTTCGCCACTGTCAGCCCTTCCCCCAGACTTaa                                                                                                                    >  2:1070877/1‑149 (MQ=255)
                                   aGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGGTTTTCGCCACCGTCAGCCCTTCCCCCAGACttta                                                                                                                    >  2:129421/1‑147 (MQ=255)
                                           aGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCAATGACCGAGGGAGAGGGTTTTCGCCCCTGTCAGCCCTTCCACCAGACCTAACTCCTGa                                                                                                             >  2:834861/1‑148 (MQ=255)
                                           aGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGCAATCAATCATGCCGTTGCGGCGGGGCCACCGACCGCGGCAGAGGGTTTCCGCCCCTGTCAGCCCCTCCCCCCCACTTAACTCCTGAt                                                                                                            >  1:1205727/1‑149 (MQ=255)
                                            ggACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATaaa                                                                                                         <  1:1045572/151‑1 (MQ=255)
                                                tAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGGAGAGGGTTTTCGCCACTGTCAGCCCCTCCCCCACACTTAACTCCCGATaaa                                                                                                         >  2:490458/1‑147 (MQ=255)
                                                 aaCGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACcgcg                                                                                                    <  2:1312540/151‑1 (MQ=255)
                                                 aaCGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCGCGGCGGAGGGTTTTCGCCCCTGCTAGCCCCTCCCCCCGCCTTAACACCTGATaaa                                                                                                         >  1:313073/1‑146 (MQ=255)
                                                    gTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCCAGGCAGAGGGTTTTCGCCCCCGTCAGCCCCTCCACCCGACTTAAATCCCGATAAACcgcgc                                                                                                   >  1:1672602/1‑149 (MQ=255)
                                                       aaCGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCCCCGAACGCGGCAGAGGGTTTTCGCCCCCGTCAACCCTTCCACCAACCTTAAATCCTCATAAAACCCCCGCAcc                                                                                              >  1:1103327/1‑151 (MQ=255)
                                                            ccAGCGCCTGTAAGCAGCGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCAGGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAg                                                                                         <  1:1676357/151‑1 (MQ=255)
                                                            ccAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCGAGCCAGGCAGAGGGTTTCCCGCCCTGTCAACCCCTCCACCCGAATTAAATCCCGATAAAACGCGCGCAccccc                                                                                           >  2:914769/1‑149 (MQ=255)
                                                              aGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGGTTTTCGCCCCCGTCAGCCCTTCCACCAGACTTAAATCCTGATAAACCGCGCGCAccccccggt                                                                                       >  2:614091/1‑147 (MQ=255)
                                                                   cTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCCGAGGGTTTTCGCCACTGTCAGCCCTTCCCCCAGACTTAACTCCTGATAAACCCCGCGCCCCCCCAGTTCAGCg                                                                                  >  2:1490226/1‑151 (MQ=255)
                                                                    tGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCGACCGGGGCAGGGGGTTTCCGCCCCCGTCAGCCCTTCCCCCCCACCTTACTCCTGGTAAAACGCGCGCccccccccttccac                                                                                   >  1:444914/1‑141 (MQ=255)
                                                                       aaGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCAc                                                                              <  1:26246/151‑1 (MQ=255)
                                                                       aaGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCAc                                                                              <  1:1590550/151‑1 (MQ=255)
                                                                       aaGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCTACCCCGGGCGAGGGTTTTCGCCCCTGTCAGCCCCTCCACCCCACCTAAATCCTGCTAAAACCCGCccccccccccgt                                                                                       >  1:1672542/1‑138 (MQ=255)
                                                                        aGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGGCCCAGGCAGAGGGTTTTCGCCACTGTCAGCCCCTCCCCCACACTTAACTCCTGATAAACCGCGCGCCCCCCCCGCTCAGCGGccc                                                                              >  1:556490/1‑148 (MQ=255)
                                                                          cAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACgg                                                                            <  1:602762/150‑1 (MQ=255)
                                                                          cAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAAGCAGAGGGGTTTCCCCCCGCACCGCCCTTCCCCCCGACTTAAATCCTGATAAAACGCGCGCCCCCCCCGGTCacccggccg                                                                             >  2:1647416/1‑141 (MQ=255)
                                                                           aaCGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGGTTTTCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAAACGCGCGCACCCCCCGTTCAGCgg                                                                                 >  2:984719/1‑144 (MQ=255)
                                                                                tGGGCATCCTGCGCCATTTGCAGGAAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTc                                                                      <  2:1277380/150‑1 (MQ=255)
                                                                                    cATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCGCGGCAGAGGGTTTTCCCCACCGTCCGCCCCTCCCCCCCACTTAAATCCTTATAAACCCCCCCCACCCCCCGTTCAACGGCACCGCGAGTCacc                                                                   >  2:822968/1‑147 (MQ=255)
                                                                                      tCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCGCCCCAGGCCGAGGTTTTTCCCCACCGTCAGCCCCTCCACCCGACTTAACTCCTGGTAAACCCCGCGCACCCCCCCTTCAGCGGCCCCGCGAGTCAGCg                                                                  >  1:302499/1‑148 (MQ=255)
                                                                                         tGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGGAGAGGGTTTTCGCCCCTGTCAGCCCTTCCACCCGACTTAAATCCTGATAAAACGCGCGCACCCCCCGTTCAGCGGCACGGCGAGTCagcttagctt                                                            >  2:131281/1‑151 (MQ=255)
                                                                                         tGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGGTTTTCGCCCCCGTCAGCCCTTCCACCCCACTTAACTCCCGATAAAACGCCCGCACCCCCCCTTCAGCGGCACCGCGGGTCAGCGCAGc                                                              >  1:1096144/1‑149 (MQ=255)
                                                                                                 ttGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACCGACCGAGGCCGAGGGTTTTCGCCCCTGTCCGTCCCTTCCCCACACTTAACACCTGTTAAACCGCCCCCACCCCCCCGTCAGCGGCACGGCGGGGCACGGGAGCTTCGTCAcc                                                     >  1:176941/1‑150 (MQ=255)
                                                                                                   gCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCtt                                                   <  1:284158/150‑1 (MQ=255)
                                                                                                     aGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCTAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCa                                                <  1:659934/151‑1 (MQ=255)
                                                                                                         aatcaatcATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTcg                                            <  1:1007672/151‑1 (MQ=255)
                                                                                                         aatcaatcATGCCGTTGCGGCGGGGCCACTGACCCAGGCAGAGGGTTTTCCCCCCCCTCAGACCTTCCCCCCGACTTAAATCCTGATAAAACCCGCGCACCCCCCCTTCAACGGGACGGCGAGGCAACGTAGCCTCGTCACCCTcccgcc                                             >  1:1233782/1‑146 (MQ=255)
                                                                                                              atcaTGCCGTTGCGGCGGGGCCACTGACCCAGGGAGAGGGTTTTCGCCCCTGTCAGCCCCTCCCCCCCACTTAACTCCTGATAAAACCCGCGCACCCCCCGTTCAGCGGCACCGGGGGGCAGCGGAGCTTCGTCACCCTcccctc                                             >  1:1148811/1‑141 (MQ=255)
                                                                                                               tcaTGCCGTTGCGGCGGGGCCCCAGACCGAGGGAGAGGGTTTTCCCCCCTGTCAGCCCCTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCCCTTCAGCGGCACGGCGAGGCAGCGGAGCTTCGTCACCCTCCCGTCGCGgctc                                      >  1:887460/1‑151 (MQ=255)
                                                                                                                caTGCCGTTGCGGCGGGGCCCCCGACCCGGGCAGAGGGTTTTCGCCCCCGTCCGCCCCTCCCCCCGACTTAACTCCTGATAAAACGCGCGCACCCCCCGGTCAGCGGGACGGCGAGTCAGCGGAGCTTCGTCACCCTCCCGt                                              >  1:421638/1‑142 (MQ=255)
                                                                                                                caTGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGGCTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCg                                     <  1:948509/151‑1 (MQ=255)
                                                                                                                caTGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGtctc                                      <  1:1329503/150‑1 (MQ=255)
                                                                                                                 aTGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGtctc                                      <  1:59480/149‑1 (MQ=255)
                                                                                                                  tGCCGTTGCGGCGGGGCCACTCACCGAGGGAGAGGGTTTCCGCCACTGTCCGCCCCTCCCCCAGACTTAAATCCTGATTAACCACCCGCACCCCCAGCTCTGCGGCACGGCGGGGCAGCCGAGCTTCCTCACCTTCCAGTCCCGTCCCGcc                                   >  1:23328/1‑151 (MQ=255)
                                                                                                                          cggcggGGCCCCTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAAACGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCCCCTTCCAGTCGCGTCTCGCCAATcc                              >  2:530247/1‑148 (MQ=255)
                                                                                                                             cggGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGatat                        >  1:1418507/1‑151 (MQ=255)
                                                                                                                                     cTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCt                 <  1:656533/150‑1 (MQ=255)
                                                                                                                                       gACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATc              <  1:796730/151‑1 (MQ=255)
                                                                                                                                         ccGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCg            <  1:921527/151‑1 (MQ=255)
                                                                                                                                          cGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCgg           <  1:723496/151‑1 (MQ=255)
                                                                                                                                              gCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGt       >  2:453465/1‑151 (MQ=255)
                                                                                                                                                 gggattttCCGCCACTGTCAGCCCCGCCGCCAGACTTAACTCCGGATAAGCCGCGCACACCCCCAGTTCAGCGGCACGGGGAGTTAGCGTAGCTTCGTCACCTTCCAGTGGCGTCTCGCCAATCCAGACATCACCGCTATCCGGGCGtt      <  1:1185999/145‑1 (MQ=255)
                                                                                                                                                 gagTTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTc     <  2:542621/150‑1 (MQ=255)
                                                                                                                                                   ggTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCt  <  2:4346/151‑1 (MQ=255)
                                                                                                                                                   ggTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCt  <  1:746049/151‑1 (MQ=255)
                                                                                                                                                   ggTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCt  <  2:928715/151‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCT  >  NC_000913/2677890‑2678187

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