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

CTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGA‑GGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGC  >  NC_000913/2677894‑2678189
                                                                                                                                                  |                                                                                                                                                       
atgctTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGACGGGGGAGGtttt                                                                                                                                                      >  1:810684/2‑150 (MQ=255)
   ctTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACTGACCGC‑GGGCGAGGgttttctcc                                                                                                                                                 >  2:750882/1‑146 (MQ=255)
         cgcCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGA‑GGCCGAGGGTTTCCCCCCCCGTc                                                                                                                                           >  1:131461/1‑151 (MQ=255)
             gggCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCTTCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGA‑GGCAGAGGTTTTCCGCCACTGTCAGcc                                                                                                                                       <  2:680277/151‑1 (MQ=255)
              ggCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACCGACCGA‑GGCAGAGGGTTTTCGCCACCGTCAGccc                                                                                                                                      >  2:153819/1‑151 (MQ=255)
               gCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGAAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCTCCCCG‑GGCAGAGGGTTTTCGCCCCTTTTAGccc                                                                                                                                      >  1:378474/1‑150 (MQ=255)
                 ttAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACTGACCCA‑GGCAGAGGGTTTTCGCCACTGTCCGcccct                                                                                                                                    >  1:374989/1‑148 (MQ=255)
                  tAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCCCTGACCGG‑GGGAGAGGGTTTTCGCCCCTGTCAGACCTTcc                                                                                                                                  >  2:454931/1‑151 (MQ=255)
                     cGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCTACCCA‑GGCCGGGGGTTTTCCCCCCCCTCAGCCCCTccccc                                                                                                                               >  1:880833/1‑151 (MQ=255)
                          aaaCAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGA‑GGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑t                                                                                                                         <  2:464186/151‑1 (MQ=255)
                              aaGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTTTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGA‑GGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTaa                                                                                                                      <  1:994426/150‑1 (MQ=255)
                              aaGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACCGACCCA‑GGCAGAGGGTTTTCGCCCCCGTCAGCCCCTCCACCCGAC‑CTaa                                                                                                                      >  1:777949/1‑150 (MQ=255)
                                 ttgttCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCCCTGACCGG‑GGCGGAGGTTTTTCGCCACTCTCAGCCCCCccaccatccctaaatcc                                                                                                                   >  1:674661/1‑140 (MQ=255)
                                 ttgttCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCCA‑GGGAGAGGTTTTCCCCCCCCGTCAGCCCCTTCACCCGAC‑TTAACTcc                                                                                                                  >  2:754183/1‑151 (MQ=255)
                                     tCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCGA‑GGCAGAGGGTTTTCGCCCCCGTCAGCCCCTCCACCCAAC‑TTAACTCCTGAt                                                                                                              >  2:1055535/1‑151 (MQ=255)
                                             aaCGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGA‑GGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACcgcg                                                                                                      <  1:516097/151‑1 (MQ=255)
                                             aaCGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGA‑GGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACcgc                                                                                                       <  2:570371/150‑1 (MQ=255)
                                             aaCGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGA‑GGCAGAGGGTTTTCGCCACTGTCAGCCCCTCCACCCGAC‑TTAACTCCTGATAAAAcgc                                                                                                       >  2:1073979/1‑150 (MQ=255)
                                              aCGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCCCCCCC‑GGCCGCGGGTTTTCTCCCCCGTCCGCCCTTCCCCCAGAC‑TTAAATCCCGTTaaaac                                                                                                         >  1:296736/1‑145 (MQ=255)
                                               cGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAGTCAATCATGCCGTTGCGGCGGGGTCACCGACCGA‑GGCGAAGGGTTTTCGCCACTGTCAGCCCCTCCCCCAGAA‑TTAACCCCTGATAAAAcgc                                                                                                       >  1:928931/1‑148 (MQ=255)
                                               cGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGA‑GGCAGAGGGTTTTCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACcgc                                                                                                       >  2:433474/1‑148 (MQ=255)
                                                   aaCGCCCAACGCCTGGAAGCAACGTTGGGAATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGA‑GGCAGGGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCa                                                                                                  <  1:613762/149‑1 (MQ=255)
                                                    aCGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCGA‑GGGAGAGGGTTTTCGCCACTGTCCGCCCCTCCCCCCGAC‑TTAACTCCTGGTAAACCCCGCGCcccc                                                                                               >  1:247593/1‑151 (MQ=255)
                                                       cccAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGA‑GGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCa                                                                                            <  2:5395/151‑1 (MQ=255)
                                                        ccAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTTACCCA‑GGCAGAGGTTTTTCGCCCCTGTCCGCCCTTCCCCCCGAC‑CTAACTCCTGATAAAACGCGCGCcccccc                                                                                             >  2:541945/1‑149 (MQ=255)
                                                        ccAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGA‑GGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAg                                                                                           <  1:838517/151‑1 (MQ=255)
                                                         cAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCGCCCCA‑GGCCGAGGGTTTTCCCCCCCCTTAGCCCCTCCCCCCGTCTTTAACCCCT‑ATAAACCCCGCGCAcccccc                                                                                            >  1:919709/1‑148 (MQ=255)
                                                         cAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACGGA‑GGGAGGGGGTTTTCGCCACCTTCAGCCCCTCCCCCAGAC‑TTAACTCCCGATAAAACGCCCGCAccccc                                                                                             >  1:885405/1‑148 (MQ=255)
                                                         cAGCGCCTGTAAGCAACGTTGGGCATCCCGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCCA‑GGCAGAGGTTTTTCCCCACCGTCAGCCCCTCCCCCAGAC‑CTAACTCCTGATAAACCCCGCCCAccccc                                                                                             >  2:486448/1‑148 (MQ=255)
                                                              ccTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGA‑GGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAg                                                                                      <  1:1025593/150‑1 (MQ=255)
                                                                 gTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGA‑GGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCg                                                                                    <  2:731902/149‑1 (MQ=255)
                                                                     gCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTCACCGA‑GGCAGAGGGTTTTCCCCCCTGTCAGCCCCTCCACCCGAC‑TTAACTCCTGATAAACCCCGCGCACCCCCCGCTCAGccg                                                                                   >  1:516255/1‑144 (MQ=255)
                                                                        aCGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCCACCGA‑GGCAGAGGGTTTTCTCCCCTGTCCGCCCCTCCCCCCCAC‑TTAACTCCTGATAAAACCCGCGCACCCCCCGTTCAGCCGCCCGGc                                                                             >  1:921765/1‑149 (MQ=255)
                                                                         cGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGA‑GGCAGAGGGTTTTCGCCACTGTCAGCCCCTCCCCCCGAC‑TTAACTCCTGATAAACCCCGCGCCCCCCCCGTTCAGCGGCACGGc                                                                             >  2:259248/1‑148 (MQ=255)
                                                                          gTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCGA‑GGCAGAGGGTTTTCGCCACTGTCAGCCCCTCCCCCAGAC‑TTAACTCCTGATAAAACGCGCGCACCCCCCGCTCAGCGGCCCGGCGAg                                                                          >  2:661654/1‑150 (MQ=255)
                                                                                cATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCGACCCA‑GGCAGAGGGTTTTCGCCCCCGTCAGCCCTTCCACCAGAC‑TTAAATCCTGATAAAACGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCg                                                                    >  2:891191/1‑150 (MQ=255)
                                                                                cATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACCCC‑GGCAGAGGTTTTTCCCCCCTGTCCGCCCCTCCCCCAGAC‑TTAACTCCTGATAAAACGCCCCCACCCCCCCTTCCGCGGCCCGGCGAGTCAGCg                                                                    >  1:941546/1‑150 (MQ=255)
                                                                                 aTCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGA‑GGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTa                                                                  <  1:249722/151‑1 (MQ=255)
                                                                                  tCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGA‑GGCAGAGGGTTTCCCCCCCTGTCAGCCCTTCCACCCGAC‑TTAACTCCTGATAAAACCCGCGCCCCCCCAGTTCAGCGGCACCGCGAGGCAGc                                                                     >  2:386238/1‑147 (MQ=255)
                                                                                       cgcCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCAC‑GGCGCAGGGTTTTCCCCACCCCTTCCCCTTTCACCACAC‑TCCAAACCTGATAAACCGCGCGCACCCCCCGTTCACCTGCCCGCCGGGTCCGCGTa                                                                  >  1:334104/1‑145 (MQ=255)
                                                                                         ccATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGA‑GGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGt                                                           <  1:160424/150‑1 (MQ=255)
                                                                                         ccATTTGCAGGTAATCAACCCTGCCGTTGCGGCGGGGCCACTGACCGA‑GGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTc                                                          <  1:24578/151‑1 (MQ=255)
                                                                                             ttGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCGACCGC‑GGCAGAGGGTTTTCCCCCCTGTCAGCCCCTCCCCCCGAC‑TTAACTCCCGATAAAACCCGCGCCCCCCCCGGTCAGCGGCACGGCGGGGCAGCGGAGCTTCGTCAcc                                                       >  2:1031889/1‑150 (MQ=255)
                                                                                                cAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGA‑GGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTcc                                                   <  2:880878/151‑1 (MQ=255)
                                                                                                  ggTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACCGA‑GGCAGAGGGTTTTCGCCACTGTCAGCCCCTCCACCCGAC‑TTAACTCCTGATAAAACCCGCGCACCCCCCCTTCAGCGGCACCGCGGGGCAGCGTAGCCTCGTCACCCTcc                                                   >  2:1053085/1‑149 (MQ=255)
                                                                                                     aatcaatcATGCCGTTGCGGCGGGGCCACTGACCGA‑GGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTc                                               <  1:1031823/150‑1 (MQ=255)
                                                                                                        caatcaTGCCGTTGCGGCGGGGCCCCTGACCGA‑GGCAGAGGGTTTTCGCCACCGTCAGCCCTTCCCCCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCCGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCCTCCCGTc                                               >  2:182962/1‑147 (MQ=255)
                                                                                                        caatcaTGCCGTTGCGGCGGGGCCACTGACCGA‑GGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGt                                           <  1:844807/151‑1 (MQ=255)
                                                                                                          atcaTGCCGTTGCGGCGGGGCCACTGACCGA‑GGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGtct                                         <  2:828022/151‑1 (MQ=255)
                                                                                                            caTGCCGTTGCGGCGGGGCCCCCGACCCA‑GGCAGAGGGTTTTCGCCCCCCTCAGCCCCTCCCCCAGAC‑TTAACTCCTGATAAAACCCGCGCACCCCCCCCTCAGCGGGACCGCGGGGCCGCGGAGCTTCGTCACCTTCCCCTCGCGcct                                         >  1:758007/1‑149 (MQ=255)
                                                                                                             aTGCCGTTGCGGCGGGGCCACTGACCGA‑GGCAGAGGGTTTTCGCCACTGTCAGCCCTTCCACCCGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCCGTTCAGCGGCACGGCGGGGCAGCGGAGCTTCGTCACCCTCCCGTCGCGTCTCGc                                      >  2:950075/1‑151 (MQ=255)
                                                                                                                 cGTTGCGGCGGGGCCACTGACCGA‑GGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAAt                                  <  2:904492/151‑1 (MQ=255)
                                                                                                                  gTTGCGGCGGGGCCACTGACCGA‑GGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATc                                 <  2:716301/151‑1 (MQ=255)
                                                                                                                     gcggcgGGGCCACTGACCGA‑GGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCa                               <  2:663117/150‑1 (MQ=255)
                                                                                                                       ggcggGGCCACTGACCGA‑GGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGat                            <  1:290276/151‑1 (MQ=255)
                                                                                                                                 cTGACCGA‑GGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCt                   <  1:957963/150‑1 (MQ=255)
                                                                                                                                   gACCGA‑GGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTa                  <  1:1079489/149‑1 (MQ=255)
                                                                                                                                   gACCGA‑GGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATc                <  2:2036/151‑1 (MQ=255)
                                                                                                                                      cGA‑GGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCgg             >  1:246388/1‑151 (MQ=255)
                                                                                                                                            cAGAGTTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGt         <  1:636743/150‑1 (MQ=255)
                                                                                                                                            cAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGtt        <  1:693336/151‑1 (MQ=255)
                                                                                                                                             agagGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTc       <  2:375009/151‑1 (MQ=255)
                                                                                                                                              gagGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCg      <  2:988800/151‑1 (MQ=255)
                                                                                                                                                ggTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCt    <  1:757400/151‑1 (MQ=255)
                                                                                                                                                  ttttCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGc  >  1:760631/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                       
CTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACCGA‑GGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGAC‑TTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGCTGC  >  NC_000913/2677894‑2678189

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: ‑