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

GCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCG  >  NC_000913/3508507‑3508797
                                                                                                                                                       |                                                                                                                                             
gcAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGa                                                                                                                                               <  2:571029/151‑1 (MQ=255)
     ccACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGgaaaa                                                                                                                                           >  1:852935/1‑150 (MQ=255)
     ccACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGgaaa                                                                                                                                            >  2:1099652/1‑149 (MQ=255)
     ccACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGAGGCGCGCGGaaaaa                                                                                                                                           >  1:119460/1‑150 (MQ=255)
      cACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAtt                                                                                                                                         <  1:463283/151‑1 (MQ=255)
            gttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTcacaca                                                                                                                                   <  1:131325/151‑1 (MQ=255)
              ttttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCacacac                                                                                                                                 >  1:582263/1‑150 (MQ=255)
                ttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGAGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTCGGGCGCGGGAAAAATcac                                                                                                                                      >  2:904189/1‑144 (MQ=255)
                 tAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATAaaacaccg                                                                                                                               >  1:641757/1‑147 (MQ=255)
                  aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAg                                                                                                                            <  1:148558/151‑1 (MQ=255)
                        ggAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCAAACACCGTAGccccc                                                                                                                       >  1:704745/1‑150 (MQ=255)
                          agagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCACACACCGTAGcccccc                                                                                                                      >  1:1041819/1‑148 (MQ=255)
                            agaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGCGCGGGAAAAATCAAAAACGGTAGgccccc                                                                                                                      >  2:723175/1‑146 (MQ=255)
                             gaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCa                                                                                                                 <  1:30825/151‑1 (MQ=255)
                             gaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCa                                                                                                                 <  1:825496/151‑1 (MQ=255)
                             gaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCa                                                                                                                 <  1:191658/151‑1 (MQ=255)
                             gaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCa                                                                                                                 <  1:580522/151‑1 (MQ=255)
                             gaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAACCCCcac                                                                                                                     >  1:1044960/1‑147 (MQ=255)
                                 aaTATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTAAAACCGGTAGCCCCCCCCTCaaaa                                                                                                              >  1:591649/1‑150 (MQ=255)
                                 aaTATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGCGCGGGAAAAATCACAAACCGTAGCCCCCCCATCAGaa                                                                                                              >  1:465659/1‑150 (MQ=255)
                                   tatGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCACACACGGTAGCCCCCCCATCAGaaa                                                                                                             >  1:298394/1‑149 (MQ=255)
                                   tatGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGaaa                                                                                                             >  1:232327/1‑149 (MQ=255)
                                    atGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGc                                                                                                           <  1:452186/150‑1 (MQ=255)
                                    atGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGCGGGGGGCGGGAAAAATTACAAACCGTGGCCCCCCCAACacgaa                                                                                                             >  2:951163/1‑144 (MQ=255)
                                     tGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCACAAACGGGAGCCCCCCCATCAGAAAACgg                                                                                                         >  2:531217/1‑150 (MQ=255)
                                          aTGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACCGTAGCCCCCCCATCAGAAAACggggtg                                                                                                     >  1:212416/1‑150 (MQ=255)
                                           tGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGCGCGCGGGAAAAATCACAAACGGGAGCCCCCACATCACGAAGCGCGGGgcg                                                                                                   >  2:530770/1‑151 (MQ=255)
                                            gCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTgcg                                                                                                   <  2:197305/150‑1 (MQ=255)
                                              tCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACAAGGTAGCCCCCACATCAGAAAACGTGGGGCGcac                                                                                                >  1:828104/1‑151 (MQ=255)
                                                 tGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGAGCGCGCGGGAAAATTCAAACACGGTGGACCCCCCAGCAGGAAGGGAGGGGCGGACAcc                                                                                             >  1:436197/1‑151 (MQ=255)
                                                   ggCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGa                                                                                           <  2:648232/151‑1 (MQ=255)
                                                        gggTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACGGGAGCCCCCCCACCAGAAAACGGGGGGGGCACACCGAGCTGc                                                                                      >  2:746577/1‑151 (MQ=255)
                                                         ggTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGCGCGCGGGAAAAATTACACACGGAAGGCCCCCCAACAGAAAAAGCGGGGCGCACAcc                                                                                             >  2:134035/1‑143 (MQ=255)
                                                           tgtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATCACACACGGTTGCCCCCCCCTCAAAAAGGGTGGGGCGGACACCGCGCCGCGt                                                                                    >  1:670580/1‑150 (MQ=255)
                                                            gtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAg                                                                                  <  2:992339/151‑1 (MQ=255)
                                                               aCCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGa                                                                               <  1:833847/151‑1 (MQ=255)
                                                                       gATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGTAGCCCCCACATCAGAAAACGTGGGGCGCACACCGAGCTGCGTAGCGATCATccc                                                                        >  1:638797/1‑150 (MQ=255)
                                                                            aaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCACACACGGTAGGCCCCACATCAGAAAGGGTGGGGGGCACACCGAGCCGCGGAGCGATCATCCCCTGc                                                                    >  2:17913/1‑149 (MQ=255)
                                                                                   tAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGGGCGGGAAAAATCACACACCGTAAACCCCCCAACAGAAAACGTGGTGCGCGCACCGAGCCGCGTAGGGATGATCCCCTGCTGATaaa                                                             >  2:1105251/1‑149 (MQ=255)
                                                                                      ggTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGATGGGGGCGGGAAAAATCAAACACCGTAATCCCCACACCACCAAGTGCGGTGCGGACACCGACCCGCATAGCGATCATCCCTTTCTCCTaaaac                                                           >  2:217735/1‑146 (MQ=255)
                                                                                         gtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGGAGCCCCCCCATCAGAAAACGTGGGGCGCCCACCGAGCTGCGTAGCGATCATCCCCTGCTGATaaaac                                                           >  1:216864/1‑143 (MQ=255)
                                                                                            tcCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGTAGCCCCCACATCAAAAAACGTGGTGCGCACACCCAGCTGCGGAGCGATCATCCCCTGCTGATAAACCAGCACTg                                                    >  2:21056/1‑149 (MQ=255)
                                                                                                gCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGcacca                                              <  2:19512/151‑1 (MQ=255)
                                                                                                  aaaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAATTCACACACCGGTACCCCCCCAACAGAAAGAGTGGGGCGGCCACCCCGCCGCGGCACGGTCCTCCCTTGCTGATAAAACaacccccgcacc                                               >  2:502749/1‑137 (MQ=255)
                                                                                                    aCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGGAGCCCCCACATCAGAAAGCGTGGTGGGCACACCCAGCTGCGTAGCGATCATCCCTTTCTGGTAAAACAGCACTGGcaccac                                             >  1:556759/1‑148 (MQ=255)
                                                                                                      cGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTcgc                                        <  1:10161/151‑1 (MQ=255)
                                                                                                        gagggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCgg                                      <  2:440679/151‑1 (MQ=255)
                                                                                                         agggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGTAGCCCCCCCATCAGAAAACGGGGGGCGCACACCGAGCTGCGGAGCGATCATCCCCTGCTGATAAAACAGCACTGGCACCACCTCGCggg                                     >  1:647319/1‑151 (MQ=255)
                                                                                                              ggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTAtgg                                <  1:1092047/151‑1 (MQ=255)
                                                                                                                 aagaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATCACACACGGTAGCCCCCACATCAGAAAGCGTGGGGCGCACACCCAGCTGCGTAGCGATCATCCCCTGCTGATAAAACAACACTGGCACCACTTCGCGGGGAtggtgg                             >  1:825580/1‑151 (MQ=255)
                                                                                                                 aagaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATCACACACGGGAGCCCCCCCATCAGAAAAAGTGGGGCGCACACCCAGCTGCGTAGCGATCATCCCTTGCTGATAAAAAAGCACTGGCACCACCTCGCGGGTAtggtgg                             >  2:1057212/1‑151 (MQ=255)
                                                                                                                   gaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGGAGCCCCCACATCAGAAAGCGTGGTGCGCACACCCAGCTGCGGAGCGATCATCCCTTGCTGATAAAACAGCACTGGCACCCCCTCGCGGGGAtggtgg                             >  2:561237/1‑149 (MQ=255)
                                                                                                                   gaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTTACACACGGGAGCCCCCCCATCAGAAAACGGGGGGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGGAtggtgg                             >  1:564025/1‑149 (MQ=255)
                                                                                                                   gaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATTACACACGGGAGCCCCCACATCAGAAAGCGTGGGGCGCACACCCAGCTGCGTAGCGATCATCCCCTGCTGATAAAACAGCACCGGCACCCCCTCGCGGGGAtggtgg                             >  1:825681/1‑149 (MQ=255)
                                                                                                                     aaagaGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTAAAAAAGGTAGCCCCCCCATCAAAAAACGGGGGGCGCACACCCAGGCGCGGAGCGATCATCCCCTGCTGGTAAAACACCACTGGCACCACCTCCCCGGGAtggggg                             >  2:441276/1‑147 (MQ=255)
                                                                                                                      aagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCtgtg                        <  2:1001932/151‑1 (MQ=255)
                                                                                                                         agCGACCATTTTGTGGCGGTGGCGCGCGGGAAAATTTACACCCGGGAGCCCCCACCTCAGAAAGCGGGGTGCGGACACCCCGCTGCGGAGCGGGCATTCCCTTCTGATAAACAAACACTGGCAACACTTCGCGGGTCTGGTGGCTGTGAcc                     >  2:1098973/1‑151 (MQ=255)
                                                                                                                         agCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGAcc                     <  2:135704/151‑1 (MQ=255)
                                                                                                                             aCCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGTAGCCCCCACATCAGAAAACGTGGGGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAAACAGCACTGGCACCACTTCGCGGGGATGGGGGCTGTGAcc                     >  1:1054960/1‑147 (MQ=255)
                                                                                                                               cATTTTGTGGCGGTGGCGCGCGGAAAAAATCACACACGGTAGCCCCCCCATCAGAAAGAGTGGGGCGCACACCCAGCTGCGTAGCGATCATCCCTTGCTGATAAAACAGCACTGGCACCACTTCGCGGGGATGGTGGGTGTGACCAAAgg                >  1:1145247/1‑150 (MQ=255)
                                                                                                                                aTTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTc              <  1:531306/151‑1 (MQ=255)
                                                                                                                                     gtggcggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATc         <  2:1083702/151‑1 (MQ=255)
                                                                                                                                           gtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCTCTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGcc   <  2:334426/151‑1 (MQ=255)
                                                                                                                                           gtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGcc   >  2:7207/1‑151 (MQ=255)
                                                                                                                                            tggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCg  >  1:502017/1‑151 (MQ=255)
                                                                                                                                            tggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCg  >  1:978024/1‑151 (MQ=255)
                                                                                                                                             ggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCg  >  2:149995/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                             
GCAGTCCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCG  >  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: