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

CCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCT  >  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)
                cGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAG‑c                                                                                                                                     <  2:339636/151‑1 (MQ=255)
                  ggCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCCCGGCCGAGGGTTTTCGCCCCCGTCCG‑ccc                                                                                                                                   >  1:1579736/1‑151 (MQ=255)
                       aGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTcc                                                                                                                               <  2:57328/150‑1 (MQ=255)
                         cGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACTGACCGAGGCAGAGGGTTTTCGCCACTGTCAG‑CCCTTccacc                                                                                                                            >  2:178646/1‑151 (MQ=255)
                         cGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTTACCGAGGCAGAGGGTTTTCGCCACTGTCCG‑CCCCTccccc                                                                                                                            >  2:221020/1‑151 (MQ=255)
                           ttGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGGTTTTCCCCCCCGTCAG‑CCCTTccccc                                                                                                                            >  2:818019/1‑149 (MQ=255)
                            tGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAg                                                                                                                          <  2:251603/150‑1 (MQ=255)
                              aaaCAAGTTGTTCAGGCCTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAg                                                                                                                          <  1:907814/148‑1 (MQ=255)
                                aCAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGCGGCAGAGGGTTTTCGCCCCTGTCAG‑CCCCTCCACCCGACtt                                                                                                                      >  2:1033902/1‑150 (MQ=255)
                                aCAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTGCCGCCACTGTCAG‑CCCTTCCACCAGACTTa                                                                                                                     <  1:300584/151‑1 (MQ=255)
                                 cAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTaa                                                                                                                    <  2:1594258/151‑1 (MQ=255)
                                  aaGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTaa                                                                                                                    <  2:1231217/150‑1 (MQ=255)
                                   aGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTAACt                                                                                                                  <  1:1127095/151‑1 (MQ=255)
                                   aGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGGTTTTCGCCACTGTCAG‑CCCTTCCCCCAGACTTaa                                                                                                                    >  2:1070877/1‑149 (MQ=255)
                                   aGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGGTTTTCGCCACCGTCAG‑CCCTTCCCCCAGACttta                                                                                                                    >  2:129421/1‑147 (MQ=255)
                                           aGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCAATGACCGAGGGAGAGGGTTTTCGCCCCTGTCAG‑CCCTTCCACCAGACCTAACTCCTGa                                                                                                             >  2:834861/1‑148 (MQ=255)
                                           aGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGCAATCAATCATGCCGTTGCGGCGGGGCCACCGACCGCGGCAGAGGGTTTCCGCCCCTGTCAG‑CCCCTCCCCCCCACTTAACTCCTGAt                                                                                                            >  1:1205727/1‑149 (MQ=255)
                                            ggACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTAACTCCTGATaaa                                                                                                         <  1:1045572/151‑1 (MQ=255)
                                                tAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGGAGAGGGTTTTCGCCACTGTCAG‑CCCCTCCCCCACACTTAACTCCCGATaaa                                                                                                         >  2:490458/1‑147 (MQ=255)
                                                 aaCGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTAACTCCTGATAAACcgcg                                                                                                    <  2:1312540/151‑1 (MQ=255)
                                                 aaCGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCGCGGCGGAGGGTTTTCGCCCCTGCTAG‑CCCCTCCCCCCGCCTTAACACCTGATaaa                                                                                                         >  1:313073/1‑146 (MQ=255)
                                                    gTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCCAGGCAGAGGGTTTTCGCCCCCGTCAG‑CCCCTCCACCCGACTTAAATCCCGATAAACcgcgc                                                                                                   >  1:1672602/1‑149 (MQ=255)
                                                       aaCGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCCCCGAACGCGGCAGAGGGTTTTCGCCCCCGTCAA‑CCCTTCCACCAACCTTAAATCCTCATAAAACCCCCGCAcc                                                                                              >  1:1103327/1‑151 (MQ=255)
                                                            ccAGCGCCTGTAAGCAGCGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCAGGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAg                                                                                         <  1:1676357/151‑1 (MQ=255)
                                                            ccAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCGAGCCAGGCAGAGGGTTTCCCGCCCTGTCAA‑CCCCTCCACCCGAATTAAATCCCGATAAAACGCGCGCAccccc                                                                                           >  2:914769/1‑149 (MQ=255)
                                                              aGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGGTTTTCGCCCCCGTCAG‑CCCTTCCACCAGACTTAAATCCTGATAAACCGCGCGCAccccccggt                                                                                       >  2:614091/1‑147 (MQ=255)
                                                                   cTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCCGAGGGTTTTCGCCACTGTCAG‑CCCTTCCCCCAGACTTAACTCCTGATAAACCCCGCGCCCCCCCAGTTCAGCg                                                                                  >  2:1490226/1‑151 (MQ=255)
                                                                    tGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCGACCGGGGCAGGGGGTTTCCGCCCCCGTCAG‑CCCTTCCCCCCCACCTTACTCCTGGTAAAACGCGCGCccccccccttccac                                                                                   >  1:444914/1‑141 (MQ=255)
                                                                       aaGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCAc                                                                              <  1:26246/151‑1 (MQ=255)
                                                                       aaGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCAc                                                                              <  1:1590550/151‑1 (MQ=255)
                                                                       aaGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCTACCCCGGGCGAGGGTTTTCGCCCCTGTCAG‑CCCCTCCACCCCACCTAAATCCTGCTAAAACCCGCccccccccccgt                                                                                       >  1:1672542/1‑138 (MQ=255)
                                                                        aGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGGCCCAGGCAGAGGGTTTTCGCCACTGTCAG‑CCCCTCCCCCACACTTAACTCCTGATAAACCGCGCGCCCCCCCCGCTCAGCGGccc                                                                              >  1:556490/1‑148 (MQ=255)
                                                                          cAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACgg                                                                            <  1:602762/150‑1 (MQ=255)
                                                                          cAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAAGCAGAGGGGTTTCCCCCCGCACCG‑CCCTTCCCCCCGACTTAAATCCTGATAAAACGCGCGCCCCCCCCGGTCacccggccg                                                                             >  2:1647416/1‑141 (MQ=255)
                                                                           aaCGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGGTTTTCGCCACTGTCAG‑CCCTTCCACCAGACTTAACTCCTGATAAAACGCGCGCACCCCCCGTTCAGCgg                                                                                 >  2:984719/1‑144 (MQ=255)
                                                                                tGGGCATCCTGCGCCATTTGCAGGAAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTc                                                                      <  2:1277380/150‑1 (MQ=255)
                                                                                    cATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCGCGGCAGAGGGTTTTCCCCACCGTCCG‑CCCCTCCCCCCCACTTAAATCCTTATAAACCCCCCCCACCCCCCGTTCAACGGCACCGCGAGTCacc                                                                   >  2:822968/1‑147 (MQ=255)
                                                                                      tCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCGCCCCAGGCCGAGGTTTTTCCCCACCGTCAG‑CCCCTCCACCCGACTTAACTCCTGGTAAACCCCGCGCACCCCCCCTTCAGCGGCCCCGCGAGTCAGCg                                                                  >  1:302499/1‑148 (MQ=255)
                                                                                         tGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGGAGAGGGTTTTCGCCCCTGTCAG‑CCCTTCCACCCGACTTAAATCCTGATAAAACGCGCGCACCCCCCGTTCAGCGGCACGGCGAGTCagcttagctt                                                            >  2:131281/1‑151 (MQ=255)
                                                                                         tGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGGTTTTCGCCCCCGTCAG‑CCCTTCCACCCCACTTAACTCCCGATAAAACGCCCGCACCCCCCCTTCAGCGGCACCGCGGGTCAGCGCAGc                                                              >  1:1096144/1‑149 (MQ=255)
                                                                                                 ttGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACCGACCGAGGCCGAGGGTTTTCGCCCCTGTCCGTCCCTTCC‑CCACACTTAACACCTGTTAAACCGCCCCCACCCCCCCGTCAGCGGCACGGCGGGGCACGGGAGCTTCGTCAcc                                                     >  1:176941/1‑150 (MQ=255)
                                                                                                   gCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCtt                                                   <  1:284158/150‑1 (MQ=255)
                                                                                                     aGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCTAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCa                                                <  1:659934/151‑1 (MQ=255)
                                                                                                         aatcaatcATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTcg                                            <  1:1007672/151‑1 (MQ=255)
                                                                                                         aatcaatcATGCCGTTGCGGCGGGGCCACTGACCCAGGCAGAGGGTTTTCCCCCCCCTCAG‑ACCTTCCCCCCGACTTAAATCCTGATAAAACCCGCGCACCCCCCCTTCAACGGGACGGCGAGGCAACGTAGCCTCGTCACCCTcccgcc                                             >  1:1233782/1‑146 (MQ=255)
                                                                                                              atcaTGCCGTTGCGGCGGGGCCACTGACCCAGGGAGAGGGTTTTCGCCCCTGTCAG‑CCCCTCCCCCCCACTTAACTCCTGATAAAACCCGCGCACCCCCCGTTCAGCGGCACCGGGGGGCAGCGGAGCTTCGTCACCCTcccctc                                             >  1:1148811/1‑141 (MQ=255)
                                                                                                               tcaTGCCGTTGCGGCGGGGCCCCAGACCGAGGGAGAGGGTTTTCCCCCCTGTCAG‑CCCCTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCCCTTCAGCGGCACGGCGAGGCAGCGGAGCTTCGTCACCCTCCCGTCGCGgctc                                      >  1:887460/1‑151 (MQ=255)
                                                                                                                caTGCCGTTGCGGCGGGGCCCCCGACCCGGGCAGAGGGTTTTCGCCCCCGTCCG‑CCCCTCCCCCCGACTTAACTCCTGATAAAACGCGCGCACCCCCCGGTCAGCGGGACGGCGAGTCAGCGGAGCTTCGTCACCCTCCCGt                                              >  1:421638/1‑142 (MQ=255)
                                                                                                                caTGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAGGCTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCg                                     <  1:948509/151‑1 (MQ=255)
                                                                                                                caTGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGtctc                                      <  1:1329503/150‑1 (MQ=255)
                                                                                                                 aTGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGtctc                                      <  1:59480/149‑1 (MQ=255)
                                                                                                                  tGCCGTTGCGGCGGGGCCACTCACCGAGGGAGAGGGTTTCCGCCACTGTCCG‑CCCCTCCCCCAGACTTAAATCCTGATTAACCACCCGCACCCCCAGCTCTGCGGCACGGCGGGGCAGCCGAGCTTCCTCACCTTCCAGTCCCGTCCCGcc                                   >  1:23328/1‑151 (MQ=255)
                                                                                                                          cggcggGGCCCCTGACCGAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTAACTCCTGATAAAACGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCCCCTTCCAGTCGCGTCTCGCCAATcc                              >  2:530247/1‑148 (MQ=255)
                                                                                                                             cggGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGatat                        >  1:1418507/1‑151 (MQ=255)
                                                                                                                                     cTGACCGAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCt                 <  1:656533/150‑1 (MQ=255)
                                                                                                                                       gACCGAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATc              <  1:796730/151‑1 (MQ=255)
                                                                                                                                         ccGAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCg            <  1:921527/151‑1 (MQ=255)
                                                                                                                                          cGAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCgg           <  1:723496/151‑1 (MQ=255)
                                                                                                                                              gCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGt       >  2:453465/1‑151 (MQ=255)
                                                                                                                                                 gggattttCCGCCACTGTCAG‑CCCCGCCGCCAGACTTAACTCCGGATAAGCCGCGCACACCCCCAGTTCAGCGGCACGGGGAGTTAGCGTAGCTTCGTCACCTTCCAGTGGCGTCTCGCCAATCCAGACATCACCGCTATCCGGGCGtt      <  1:1185999/145‑1 (MQ=255)
                                                                                                                                                 gagTTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTc     <  2:542621/150‑1 (MQ=255)
                                                                                                                                                   ggTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCt  <  2:4346/151‑1 (MQ=255)
                                                                                                                                                   ggTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCt  <  1:746049/151‑1 (MQ=255)
                                                                                                                                                   ggTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCt  <  2:928715/151‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGAGGCAGAGGTTTTCCGCCACTGTCAG‑CCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCT  >  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: ‑