Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,846,6740TG42.9% 90.5 / 14.3 65T247P (ACC→CCC) hycDhydrogenase 3, membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (12/24);  new base G (27/0);  total (41/24)
Fisher's exact test for biased strand distribution p-value = 8.32e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGCCAAACCCAGCCCAGGTAATGCGCGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTGAAGGTTTCCATTTGTCCCCACGGAATAAACACCCCGACGAACATCTGCAACACCACCAGCTGTTTCAGGCTGATACCCCATTTCATGACGCCAAAGCCGCTGCCGCTGTATTCAGAGAGCGGGCCTTCCTGCAGCTCCTGCTCGGCTTCCGCC  >  NC_000913/2846526‑2846822
                                                                                                                                                     |                                                                                                                                                      
aaGCCAAACCCAGCCCAGGTAATGCGCGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTGAAGGttt                                                                                                                                                      <  1:73555/151‑1 (MQ=255)
  gCCAAACCCAGCCCAGGTAATGCGCGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTGAAGGTTTcc                                                                                                                                                    <  2:533555/151‑1 (MQ=255)
    cAAACCCAGCCCAGGTAATGCGCGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTGTAAGGTTTc                                                                                                                                                     >  1:996654/1‑149 (MQ=255)
      aaCCCAGCCCAGGTAATGCGCGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTGAAGGTTTCCAttt                                                                                                                                                <  2:633242/151‑1 (MQ=255)
      aaCCCAGCCCAGGTAATGCGCGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTGAAGGTTTCCAttt                                                                                                                                                <  2:20070/151‑1 (MQ=255)
       aCCCAGCCCAGGTAATGCGCGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTGAAGGGTTCCAttt                                                                                                                                                >  1:678850/1‑150 (MQ=255)
         ccagcccagGTAATGCGCGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTGAAGGTTTTCAttttt                                                                                                                                              >  2:318483/1‑148 (MQ=255)
           agcccagGTAATGCGCGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGCGAAGGGTTCCAttttttcc                                                                                                                                           >  2:806919/1‑146 (MQ=255)
                    aaTGCGCGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTGAAGGTTTCCATTTGTCCCCACGGAATa                                                                                                                                 <  1:857623/151‑1 (MQ=255)
                     aTGCGCGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTGAAGGGTTCCCTTTTTCCCCCCGGAATa                                                                                                                                 >  2:945651/1‑150 (MQ=255)
                       gcgcGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTGAAGGGTTCCATTTTTTCCCCCGGAATaaa                                                                                                                               >  2:324334/1‑150 (MQ=255)
                          cgGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTGAAGGGTTTCATTTTTTCCCCCGGAAAAAAAAcc                                                                                                                           >  2:613438/1‑151 (MQ=255)
                          cgGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTGAAGGGTTTCATTTTTCCCCACGGAATAAAAAcc                                                                                                                           >  2:574330/1‑151 (MQ=255)
                           ggAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTGAAGGTTTCCATTTGTCCCCACGGAATAAACAccc                                                                                                                          <  1:723826/151‑1 (MQ=255)
                            gAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGCGGAAGTATCCCTTTCGCCCCACGCAATaaaacccac                                                                                                                         >  2:597280/1‑145 (MQ=255)
                              gTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTGAAGGTTTCCATTTGTCCCCACGGAATAAACACCCcga                                                                                                                       <  1:174175/151‑1 (MQ=255)
                                aaTATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTGAAAGTTTCCATTTTTCCCCCCGGCAAAAAAccccccac                                                                                                                       >  1:744312/1‑148 (MQ=255)
                                 atatCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTGAAGGTTTCCATTTGTCCCCACGGAATAAACAcccccccg                                                                                                                     >  2:579675/1‑146 (MQ=255)
                                 atatCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTGAAGGGTTTCATTTTTTCCCACGGGATAAAAAcccc                                                                                                                         >  1:176855/1‑146 (MQ=255)
                                 atatCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTGAAGGGTTCCATTTGTCCCCACGGAATAAACACCCccac                                                                                                                      >  1:191125/1‑149 (MQ=255)
                                  tatCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTGAAGGGTTTCCTTTGTCCCCCCCGAATAAAAAcccc                                                                                                                         >  2:523514/1‑145 (MQ=255)
                                  tatCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTGAAGGGTTCCATTTGTCCCCCCGGGATAAACAcccc                                                                                                                         >  1:474917/1‑145 (MQ=255)
                                   atCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGAGAAGGATTTCCTTTGTTCCCCCGGAAAAAAAACCCccacaca                                                                                                                   >  2:437356/1‑147 (MQ=255)
                                      aaGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTGAAGGGTTTCATTTTTTCCCCCCGAATAAAAAcccccccga                                                                                                                    >  2:83335/1‑141 (MQ=255)
                                            cAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTGAAGGGTTCCATTTGTCCCCACGGAATAAACACCCCGACGAACATCCGCaa                                                                                                          >  2:373101/1‑150 (MQ=255)
                                                cGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTGGAGGGTTTCATTTTTCCCCCCGGAAAAAAAAACCCCCCCCACATCTGCAAaacc                                                                                                      >  2:321864/1‑150 (MQ=255)
                                                 gggCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTGAAGGTTTCCATTTGTCCCCACGGAATAAACACCCCGACGAACATCTGCAAcaccac                                                                                                    <  2:491151/151‑1 (MQ=255)
                                                  ggCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTGAAGGTTTCCATTTGTCCCCACGGAATAAACACCCCGACGAACATCTGCAACaccacc                                                                                                   <  1:558187/151‑1 (MQ=255)
                                                   gCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTGAAGGTTTCCATTTGTCCCCACGGAATAAACACCCCGACGAACATCTGCAACAccacca                                                                                                  <  1:479487/151‑1 (MQ=255)
                                                         cTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTGAAGGGTTTCATTTTTCCCCCCCGGAAAAACACCCCCACGAACATCTGCCACACCACCAGCTgct                                                                                            >  2:1000828/1‑149 (MQ=255)
                                                           gTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGCGAAGGGTTCCATTTGTCCCCCCCGGATAAACACCCCCCCCAACATCTGCCACACCCCCACCTGtt                                                                                            >  2:636808/1‑149 (MQ=255)
                                                                  aaCAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTGAAGGTTTCCATTTGTCCCCACGGAATAAACACCCCGACGAACATCTGCAACACCACCAGCTGTTTCAGGCTGa                                                                                   <  1:893787/151‑1 (MQ=255)
                                                                  aaCAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTGAAGGGTTTCATTTTTCCCCACGGAATAAAAACCCCCACGAACAACTGCAACACCACCAGCTGttt                                                                                           >  2:348843/1‑143 (MQ=255)
                                                                      gcgcGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTGAAGGTTTCCATTTGTCCCCACGGAATAAACACCCCGACGAACATCTGCAACACCACCAGCTGTTTCAGGCTGATAcc                                                                              <  1:523668/151‑1 (MQ=255)
                                                                          gATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTAAAAGGTTTCATTTGTACCCCCGGAATAAAAACACCCACTAACACCTGCTAAAACACCAGCGGGTTTAGGTCGATAcccc                                                                            >  1:746347/1‑149 (MQ=255)
                                                                               ccAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTGAAGGTTTCCATTTGTCCCCACGGAATAAACACCCCGACGAACATCTGCAACACCACCAGCTGTTTCAGGCTGATACCCCATTTCAt                                                                     <  1:695120/151‑1 (MQ=255)
                                                                               ccAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTGAAGGGTTCCATTTGTCCCCCCGGAATAAACACCCCCACGAACATCTGCAACACCACCAGCTGTTTCAGGCTGATACCCCCttt                                                                        >  1:105336/1‑148 (MQ=255)
                                                                                cAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTGAAGGGTTTCATTTTTCCCCACGGAATAAACACCCCGACGAACATCTGCAACACCACCCGCCGTTTCAGGCTGATACCCCCTTTCAt                                                                     >  1:250016/1‑150 (MQ=255)
                                                                                   gACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTGAAGGGTTCCCTTTTTCCCCCCGGAATAAACACCCCGACGAACATCTGCAACACCACCAGCCGTTTCAGGCTGATACCCCCTTTCATGAc                                                                  >  1:405940/1‑150 (MQ=255)
                                                                                         gACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTGAAGGTTTTCCTTTTTCCCCCCGGGATAAACACCCCGCCGAACAACTGCAACACCACCAGCTGTTTTAGGCTGGTACCCCATTTCATGACGccca                                                             >  1:1042281/1‑147 (MQ=255)
                                                                                         gACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTGAAGGTTTCCATTTGTCCCCACGGAATAAACACCCCGACGAACATCTGCAACACCACCAGCTGTTTCAGGCTGATACCCCATTTCATGACGCCAAAg                                                           <  2:697159/151‑1 (MQ=255)
                                                                                           ccaccaGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTGAAGGTTTCCATTTGTCCCCACGGAATAAACACCCCGACGAACATCTGCAACACCACCAGCTGTTTCAGGCTGATACCCCATTTCATGACGCCAAAgcc                                                         <  2:474780/151‑1 (MQ=255)
                                                                                            caccaGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTGAAGGGTTCCCTTTTTCCCCCCGGAAAAAAAACCCCCACGAACATCTGCCAAACCCCCACCCGTTTCAGGCCGGTACCCCATTTTCTGACGCCCAAgcc                                                         >  2:139440/1‑150 (MQ=255)
                                                                                               caGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTGAAGGTTTCCATTTGTCCCCACGGAATAAACAACCCGACGAACATCTGCAACACCACCAACTGTTTCAGGCTGATAACCCATTTCATGACGCCAAAaccgcc                                                      >  2:422649/1‑149 (MQ=255)
                                                                                                   tttACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTGAAAGTTTCCATTTTTTCCCCCGGAAAAAAAACCCCGACGAACAACTGCAACACCACCACCTGTTTCAGGCTGATACCCCATTTTCTGACGCCCAAgccgccgccgc                                                 >  2:793242/1‑151 (MQ=255)
                                                                                                      aCGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTGAAGGTTTCCATTTGTCCCCACGGAATAAACACCCCGACGAACATCTGCAACACCACCAGCTGTTTCAGGCTGATACCCCATTTCATGACGCCAAAGccgctgccgctg                                               <  2:624099/150‑1 (MQ=255)
                                                                                                      aCGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTGAAGGTTTCCATTTGTCCCCACGGAATAAACACCCCGACGAACATCTGCAACACCACCAGCTGTTTCAGGCTGATACCCCATTTCATGACGCCAAAGCCGCTGCCGCTGt                                              <  2:457819/151‑1 (MQ=255)
                                                                                                      aCGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTGAAGGGTTCCATTTTTTCCCCCGGAATAAAAACCCCCACCAACATCTGCAACACCACCCGCTGTTTCCGGCTGATACCCCACTTCATTACCCCAAAacccccgccgc                                                 >  2:310390/1‑148 (MQ=255)
                                                                                                       cGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTGTAGGAGTTCATTTTTCCCCCCGGAATAAAAAACCCCCCGAACAACTTCCACACCCCCCGCGGGTTCACGCGGAAATCCCCATTTATAACGCCCAAACCGCTGCCGCCGt                                              >  2:57482/1‑150 (MQ=255)
                                                                                                          tGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTGAAGGGTTTCATTTTTCCCCCCGGAATAAACACCCCCACGGACATCTGCAACACCACCCGCTGTTTCAGGCTGATACCCCCTTTCATTACGCCCAAGCCGCCGCCGCTGGAtt                                           >  2:800705/1‑150 (MQ=255)
                                                                                                          tGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTGAAGGGTTCCATTTTTCCCCACCGAATAAAAACCCCCACGAACATCTTCAACACCACCAGCTGTTTCAGGCTGATACCCCATTTCATGACGCCAAAACCGCTGCCGCTGTAtt                                           >  1:555333/1‑150 (MQ=255)
                                                                                                                tCACCAGCGCCAGCAGCAGTCCACCGGCGGTGAAGGGTTCCATTTGTCCCCACGGAATAAAAACCCCGACGCACATCTGCAACACCACCAGCTGTTTCAGGCTGATACCCCATTTCATGACGCCCAAACCGCTGCCCCTGTATTCaga                                       >  2:399545/1‑148 (MQ=255)
                                                                                                                   ccAGCGCCAGCAGCAGTCCACCGGCGGTGAAGGTTTCCATTTGTCCCCACGGAATAAACACCCCGACGAACATCTGCAACACCACCAGCTGTTTCAGGCTGATACCCCATTTCATGACGCCAAAGCCGCTGCCGCTGTATTCAGAGAGCgg                                 <  2:705009/151‑1 (MQ=255)
                                                                                                                   ccAGCGCCAGCAGCAATCCAACGGCGGTGAAGGTTTCCATTTGTCCCCCCGGAATAAACACCCCGACGAACATCTGCAACACCACCAGCTGTTTCAGGCTGATACCCCATTTCATGACGCCAAAGCCGCTGCCGCTGTATTCAGAGAGc                                   <  2:1117945/149‑1 (MQ=255)
                                                                                                                         ccAGCAGCAGTCCACCGGCGGTGAAGGTTTTCATTTGTCCCCCCGGAATAAACACCCCGCCGCACATCTGCAACACCACCAGCTGTTTTCGGCTGATACCCCCTTTCATGACGCCAAAGCCGCTGCCGCTGTATTTAGAGAGCGGGCCtt                            >  2:224359/1‑150 (MQ=255)
                                                                                                                         ccAGCAGCAGTCCACCGGCGGTGAAGGGTTCCATTTGTCCCCACGGAATAAAAACCCCCACGAACATCTGCAAAACCACCAGCTGTTTCAGGCTGATACCCCATTTCATGACGCCAAAACCGCTGCCGCTGTATTCAGAGAGCGGGCCtt                            >  1:338784/1‑150 (MQ=255)
                                                                                                                           agcagcagTCCACCGGCGGTGAAGGTTTTCATTTTTCCCCCCGGAATAAAAACCCCCCCGAACCTCTGCAACACCCCCAGCTGTTTCAGGCTGGTACCCCCTTTTATTACCCCAAAACCGCTGCCGCTGTATTTAGAGAGCCGGCCTTcc                          >  1:478924/1‑150 (MQ=255)
                                                                                                                            gcagcagTCCACCGGCGGTGAAGGTTTCCATTTGTCCCCACGGAATAAACACCCCGACGAACATCTGCAACACCACCAGCTGTTTCAGGCTGATACCCCATTTCATGACGCCAAAGCCGCTGCCGCTGTATTCAGAGAGCGGGCCTTCCTg                        <  2:986616/151‑1 (MQ=255)
                                                                                                                              agcagTCCACCGGCGGTGAAGGTTTCCATTTGTCCCCACGGAATAAACACCCCGACGAACATCTGCAACACCACCAGCTGTTTCAGGCTGATACCCCATTTCATGACGCCAAAGCCGCTGCCGCTGTATTCAGAGAGCGGGCCTTCCTGCa                      <  1:214055/151‑1 (MQ=255)
                                                                                                                                      aCCGGCGGTGAAGGGTTTCATTTTTTCCCCCGGAAAAAAAACCCCCACGGACATCTGCAACACCACCAGCTGTTTTTGGCTGATAACCCCCTTCATGAAGCCAAAGCCGCTGCCGCCGTATTTAGAGAGCGGGGCCTCCTCCCGCTcccg               >  2:858308/1‑148 (MQ=255)
                                                                                                                                         ggcggTGAAGGTTTCCATTTGTCCCCACGGAATAAACACCCCGACGAACATCTGCAACACCACCAGCTGTTTCAGGCTGATACCCCATTTCATGACGCCAAAGCCGCTGCCGCTGTATTCAGAGAGCGGGCCTTCCTGCAGCTCCTGCTCg           <  1:646271/151‑1 (MQ=255)
                                                                                                                                         ggcggTGAAGGTTTCCATTTGTCCCCACGGAATAAACACCCCGACGAACATCTGCAACACCACCAGCTGTTTCAGGCTGATACCCCATTTCATGACGCCAAAGCCGCTGCCGCTGTATTCAGAGAGCGGGCCTTCCTGCAGCTCCTGCTCg           <  1:638603/151‑1 (MQ=255)
                                                                                                                                          gcggTGAAGGTTTCCATTTGTCCCCACGGAATAAACACCCCGACGAACATCTGCAACACCACCAGCTGTTTCAGGCTGATACCCCATTTCATGACGCCAAAGCCGCTGCCGCTGTATTCAGAGAGCGGGCCTTCCTGCAGCTCCTGCTCgg          <  2:1146119/151‑1 (MQ=255)
                                                                                                                                                  aGGTTTCCATTTGTCCCCACGGAATAAACACCCCGACGAACATCTGCAACACCACCAGCTGTTTCAGGCTGATACCCCATTTCATGACGCCAAAGCCGCTGCCGCTGTATTCAGAGAGCGGGCCTTCCTGCAGCTCCTGCTCGGCTTCCGc   <  1:437480/151‑1 (MQ=255)
                                                                                                                                                    gTTTCCATTTGTCCCCACGGAATAAACACCCCGACGAACATCTGCAACACCACCAGCTGTTTCAGGCTGATACCCCATTTCATGACGCCAAAGCCGCTGCCGCTGTATTCAGAGAGCGGGCCTTCCTGCAGCTCCTGCTCGGCTTCCGcc  >  2:241908/1‑150 (MQ=255)
                                                                                                                                                    gTTTCCATTTGTCCCCACGGAATAAACACCCCGACGAACATCTGCAACACCACCAGCTGTTTCAGGCTGATACCCCATTTCATGACGCCAAAGCCGCTGCCGCTGTATTCAGAGAGCGGGCCTTCCTGCAGCTCCTGCTCGGCTTCCGcc  >  2:950444/1‑150 (MQ=255)
                                                                                                                                                     tttCCATTTGTCCCCACGGAATAAACACCCCGACGAACATCTGCAACACCACCAGCTGTTTCAGGCTGATACCCCATTTCATGACGCCAAAGCCGCTGCCGCTGTATTCAGAGAGCGGGCCTTCCTGCAGCTCCTGCTCGGCTTCCGcc  >  2:161403/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
AAGCCAAACCCAGCCCAGGTAATGCGCGGAGTAATATCAAGACGCAGACGGGCCATGCTGTTTTCGAACAGCGCGATAACCAGGACGCCGACCACCAGTTTTACGATGGCAATCACCAGCGCCAGCAGCAGTCCACCGGCGGTGAAGGTTTCCATTTGTCCCCACGGAATAAACACCCCGACGAACATCTGCAACACCACCAGCTGTTTCAGGCTGATACCCCATTTCATGACGCCAAAGCCGCTGCCGCTGTATTCAGAGAGCGGGCCTTCCTGCAGCTCCTGCTCGGCTTCCGCC  >  NC_000913/2846526‑2846822

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