Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,508,6570AG44.4% 75.2 / 18.4 63F386S (TTC→TCC) yhfWphosphopentomutase‑related metalloenzyme superfamily protein
Reads supporting (aligned to +/- strand):  ref base A (14/21);  new base G (28/0);  total (42/21)
Fisher's exact test for biased strand distribution p-value = 8.80e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCG  >  NC_000913/3508507‑3508797
                                                                                                                                                       |                                                                                                                                             
gcAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGa                                                                                                                                               <  2:571029/151‑1 (MQ=255)
     ccACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGgaaaa                                                                                                                                           >  1:852935/1‑150 (MQ=255)
     ccACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGgaaa                                                                                                                                            >  2:1099652/1‑149 (MQ=255)
     ccACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGAGGCGCGC‑GGaaaaa                                                                                                                                           >  1:119460/1‑150 (MQ=255)
      cACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAAtt                                                                                                                                         <  1:463283/151‑1 (MQ=255)
            gttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTcacaca                                                                                                                                   <  1:131325/151‑1 (MQ=255)
              ttttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATCacaca‑c                                                                                                                                 >  1:582263/1‑150 (MQ=255)
                ttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGAGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTCGGGCGC‑GGGAAAAATcac                                                                                                                                      >  2:904189/1‑144 (MQ=255)
                 tAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATAaaaca‑ccg                                                                                                                               >  1:641757/1‑147 (MQ=255)
                  aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACA‑CGGTAg                                                                                                                            <  1:148558/151‑1 (MQ=255)
                        ggAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATCAAACA‑CCGTAGccccc                                                                                                                       >  1:704745/1‑150 (MQ=255)
                          agagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATCACACA‑CCGTAGcccccc                                                                                                                      >  1:1041819/1‑148 (MQ=255)
                            agaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGCGC‑GGGAAAAATCAAAAA‑CGGTAGgccccc                                                                                                                      >  2:723175/1‑146 (MQ=255)
                             gaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACA‑CGGTAGCCCCCACATCa                                                                                                                 <  1:30825/151‑1 (MQ=255)
                             gaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACA‑CGGTAGCCCCCACATCa                                                                                                                 <  1:825496/151‑1 (MQ=255)
                             gaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACA‑CGGTAGCCCCCACATCa                                                                                                                 <  1:191658/151‑1 (MQ=255)
                             gaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACA‑CGGTAGCCCCCACATCa                                                                                                                 <  1:580522/151‑1 (MQ=255)
                             gaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACA‑CGGTAACCCCcac                                                                                                                     >  1:1044960/1‑147 (MQ=255)
                                 aaTATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTAAAAC‑‑CGGTAGCCCCCCCCTCaaaa                                                                                                              >  1:591649/1‑150 (MQ=255)
                                 aaTATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGCGC‑GGGAAAAATCACAAA‑CCGTAGCCCCCCCATCAGaa                                                                                                              >  1:465659/1‑150 (MQ=255)
                                   tatGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATCACACA‑CGGTAGCCCCCCCATCAGaaa                                                                                                             >  1:298394/1‑149 (MQ=255)
                                   tatGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACA‑CGGTAGCCCCCACATCAGaaa                                                                                                             >  1:232327/1‑149 (MQ=255)
                                    atGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGc                                                                                                           <  1:452186/150‑1 (MQ=255)
                                    atGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGCGGGGGGC‑GGGAAAAATTACAAA‑CCGTGGCCCCCCCAACacgaa                                                                                                             >  2:951163/1‑144 (MQ=255)
                                     tGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATCACAAA‑CGGGAGCCCCCCCATCAGAAAACgg                                                                                                         >  2:531217/1‑150 (MQ=255)
                                          aTGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACA‑CCGTAGCCCCCCCATCAGAAAACggggtg                                                                                                     >  1:212416/1‑150 (MQ=255)
                                           tGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGCGCGC‑GGGAAAAATCACAAA‑CGGGAGCCCCCACATCACGAAGCGCGGGgcg                                                                                                   >  2:530770/1‑151 (MQ=255)
                                            gCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTgcg                                                                                                   <  2:197305/150‑1 (MQ=255)
                                              tCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATCACACA‑AGGTAGCCCCCACATCAGAAAACGTGGGGCGcac                                                                                                >  1:828104/1‑151 (MQ=255)
                                                 tGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGAGCGCGC‑GGGAAAATTCAAACA‑CGGTGGACCCCCCAGCAGGAAGGGAGGGGCGGACAcc                                                                                             >  1:436197/1‑151 (MQ=255)
                                                   ggCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGa                                                                                           <  2:648232/151‑1 (MQ=255)
                                                        gggTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATTACACA‑CGGGAGCCCCCCCACCAGAAAACGGGGGGGGCACACCGAGCTGc                                                                                      >  2:746577/1‑151 (MQ=255)
                                                         ggTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGCGCGC‑GGGAAAAATTACACA‑CGGAAGGCCCCCCAACAGAAAAAGCGGGGCGCACAcc                                                                                             >  2:134035/1‑143 (MQ=255)
                                                           tgtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAAATCACACA‑CGGTTGCCCCCCCCTCAAAAAGGGTGGGGCGGACACCGCGCCGCGt                                                                                    >  1:670580/1‑150 (MQ=255)
                                                            gtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAg                                                                                  <  2:992339/151‑1 (MQ=255)
                                                               aCCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGa                                                                               <  1:833847/151‑1 (MQ=255)
                                                                       gATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATCACACA‑CGGTAGCCCCCACATCAGAAAACGTGGGGCGCACACCGAGCTGCGTAGCGATCATccc                                                                        >  1:638797/1‑150 (MQ=255)
                                                                            aaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATCACACA‑CGGTAGGCCCCACATCAGAAAGGGTGGGGGGCACACCGAGCCGCGGAGCGATCATCCCCTGc                                                                    >  2:17913/1‑149 (MQ=255)
                                                                                   tAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGGGC‑GGGAAAAATCACACA‑CCGTAAACCCCCCAACAGAAAACGTGGTGCGCGCACCGAGCCGCGTAGGGATGATCCCCTGCTGATaaa                                                             >  2:1105251/1‑149 (MQ=255)
                                                                                      ggTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGATGGGGGC‑GGGAAAAATCAAACA‑CCGTAATCCCCACACCACCAAGTGCGGTGCGGACACCGACCCGCATAGCGATCATCCCTTTCTCCTaaaac                                                           >  2:217735/1‑146 (MQ=255)
                                                                                         gtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATCACACA‑CGGGAGCCCCCCCATCAGAAAACGTGGGGCGCCCACCGAGCTGCGTAGCGATCATCCCCTGCTGATaaaac                                                           >  1:216864/1‑143 (MQ=255)
                                                                                            tcCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATCACACA‑CGGTAGCCCCCACATCAAAAAACGTGGTGCGCACACCCAGCTGCGGAGCGATCATCCCCTGCTGATAAACCAGCACTg                                                    >  2:21056/1‑149 (MQ=255)
                                                                                                gCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGcacca                                              <  2:19512/151‑1 (MQ=255)
                                                                                                  aaaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAATTCACACACCGGTA‑CCCCCCCAACAGAAAGAGTGGGGCGGCCACCCCGCCGCGGCACGGTCCTCCCTTGCTGATAAAACaacccccgcacc                                               >  2:502749/1‑137 (MQ=255)
                                                                                                    aCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATCACACA‑CGGGAGCCCCCACATCAGAAAGCGTGGTGGGCACACCCAGCTGCGTAGCGATCATCCCTTTCTGGTAAAACAGCACTGGcaccac                                             >  1:556759/1‑148 (MQ=255)
                                                                                                      cGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTcgc                                        <  1:10161/151‑1 (MQ=255)
                                                                                                        gagggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCgg                                      <  2:440679/151‑1 (MQ=255)
                                                                                                         agggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATCACACA‑CGGTAGCCCCCCCATCAGAAAACGGGGGGCGCACACCGAGCTGCGGAGCGATCATCCCCTGCTGATAAAACAGCACTGGCACCACCTCGCggg                                     >  1:647319/1‑151 (MQ=255)
                                                                                                              ggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTAtgg                                <  1:1092047/151‑1 (MQ=255)
                                                                                                                 aagaaaagaGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAAATCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGGGCGCACACCCAGCTGCGTAGCGATCATCCCCTGCTGATAAAACAACACTGGCACCACTTCGCGGGGAtggtgg                             >  1:825580/1‑151 (MQ=255)
                                                                                                                 aagaaaagaGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAAATCACACA‑CGGGAGCCCCCCCATCAGAAAAAGTGGGGCGCACACCCAGCTGCGTAGCGATCATCCCTTGCTGATAAAAAAGCACTGGCACCACCTCGCGGGTAtggtgg                             >  2:1057212/1‑151 (MQ=255)
                                                                                                                   gaaaagaGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATCACACA‑CGGGAGCCCCCACATCAGAAAGCGTGGTGCGCACACCCAGCTGCGGAGCGATCATCCCTTGCTGATAAAACAGCACTGGCACCCCCTCGCGGGGAtggtgg                             >  2:561237/1‑149 (MQ=255)
                                                                                                                   gaaaagaGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTTACACA‑CGGGAGCCCCCCCATCAGAAAACGGGGGGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGGAtggtgg                             >  1:564025/1‑149 (MQ=255)
                                                                                                                   gaaaagaGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAAATTACACA‑CGGGAGCCCCCACATCAGAAAGCGTGGGGCGCACACCCAGCTGCGTAGCGATCATCCCCTGCTGATAAAACAGCACCGGCACCCCCTCGCGGGGAtggtgg                             >  1:825681/1‑149 (MQ=255)
                                                                                                                     aaagaGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATTAAAAA‑AGGTAGCCCCCCCATCAAAAAACGGGGGGCGCACACCCAGGCGCGGAGCGATCATCCCCTGCTGGTAAAACACCACTGGCACCACCTCCCCGGGAtggggg                             >  2:441276/1‑147 (MQ=255)
                                                                                                                      aagaGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCtgtg                        <  2:1001932/151‑1 (MQ=255)
                                                                                                                         agCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAATTTACACC‑CGGGAGCCCCCACCTCAGAAAGCGGGGTGCGGACACCCCGCTGCGGAGCGGGCATTCCCTTCTGATAAACAAACACTGGCAACACTTCGCGGGTCTGGTGGCTGTGAcc                     >  2:1098973/1‑151 (MQ=255)
                                                                                                                         agCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGAcc                     <  2:135704/151‑1 (MQ=255)
                                                                                                                             aCCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATCACACA‑CGGTAGCCCCCACATCAGAAAACGTGGGGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAAACAGCACTGGCACCACTTCGCGGGGATGGGGGCTGTGAcc                     >  1:1054960/1‑147 (MQ=255)
                                                                                                                               cATTTTGTGGCGGTGGCGCGC‑GGAAAAAATCACACA‑CGGTAGCCCCCCCATCAGAAAGAGTGGGGCGCACACCCAGCTGCGTAGCGATCATCCCTTGCTGATAAAACAGCACTGGCACCACTTCGCGGGGATGGTGGGTGTGACCAAAgg                >  1:1145247/1‑150 (MQ=255)
                                                                                                                                aTTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTc              <  1:531306/151‑1 (MQ=255)
                                                                                                                                     gtggcggtggcgCGC‑GGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATc         <  2:1083702/151‑1 (MQ=255)
                                                                                                                                           gtggcgCGC‑GGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCTCTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGcc   <  2:334426/151‑1 (MQ=255)
                                                                                                                                           gtggcgCGC‑GGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGcc   >  2:7207/1‑151 (MQ=255)
                                                                                                                                            tggcgCGC‑GGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCg  >  1:502017/1‑151 (MQ=255)
                                                                                                                                            tggcgCGC‑GGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCg  >  1:978024/1‑151 (MQ=255)
                                                                                                                                             ggcgCGC‑GGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCg  >  2:149995/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                             
GCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACA‑CGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCG  >  NC_000913/3508507‑3508797

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