Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,238,0490AG42.9% 119.8 / 21.7 77E139G (GAA→GGA) yjbHDUF940 family extracellular polysaccharide protein
Reads supporting (aligned to +/- strand):  ref base A (7/37);  new base G (33/0);  total (40/37)
Fisher's exact test for biased strand distribution p-value = 5.75e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTACAGCAGCGTCGAAGCGTTCTCTGGCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGG  >  NC_000913/4237900‑4238199
                                                                                                                                                        |                                                                                                                                                          
gTACAGCAGCATCGAAGCGTTCTCTGGCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTACGGGGCTGTTTGATGCGGaa                                                                                                                                                           <  2:947792/151‑1 (MQ=255)
     gcagcGTCGAAGCGTTCTCTGGCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGGACGGGGGTGTTTTGTGCGGGATATCt                                                                                                                                                      >  1:769062/1‑151 (MQ=255)
       agcGTCGAAGCGTTCTCTGGCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTACGGGGCTGTTTGATGCGGAATATCttg                                                                                                                                                    <  1:803048/151‑1 (MQ=255)
       agcGTCGAAGCGTTCTCTGGCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGGACGGGGGTGTTTGTTGGGGGATATCtt                                                                                                                                                     >  2:776782/1‑150 (MQ=255)
          gTCGAAGCGTTCTCTGGCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGGTATCGGCGGGACGGGGGGGTTTGATGCGGGATATCtt                                                                                                                                                     >  2:337453/1‑147 (MQ=255)
                   ttCTCTGGCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTACCGGGGTGTTTGATGCGGGATATATTTTTTCCAGGaaa                                                                                                                                      >  2:203855/1‑150 (MQ=255)
                      tctGGCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGGACGGGGGTGTTTTATGCGGGATATATTGTTTGCAGCAAAgac                                                                                                                                  >  1:296386/1‑148 (MQ=255)
                       ctGGCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGcc                                                                                                                                  <  1:788446/149‑1 (MQ=255)
                       ctGGCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGGACGGGGGTGTTTGGTGGGGGATATTTTTTTGGCCGCAAAAcc                                                                                                                                  >  1:443607/1‑149 (MQ=255)
                            gATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGc                                                                                                                           <  1:891457/151‑1 (MQ=255)
                             aTCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGGACGGGGGTGTTTTATGCGGGATATCTTTTTGCCCGCAAAACCCggggg                                                                                                                            >  1:672244/1‑149 (MQ=255)
                               cAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCg                                                                                                                         <  2:850581/150‑1 (MQ=255)
                                 aaCGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGGACGGGGGTGTTTGGTGCGGGATATCTTGTTGCCCGCAAAACCTGGGGGcc                                                                                                                          >  1:465635/1‑147 (MQ=255)
                                  aCGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGGACGGGGGTGTTTGGTGGGGGATATTTTTTTTGCAGCCAAAACCTGGGGGCGGt                                                                                                                        >  2:888531/1‑149 (MQ=255)
                                  aCGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGGACGGGGGTGGTTTATGGGGGATATCTTTTTTCCCGCAAAGCCTGGGGGGCg                                                                                                                         >  2:12967/1‑147 (MQ=255)
                                    gTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGGTATCGGCGGGAAGGGGGCGTTTTGTGGGGGATAACTTTTTGGCCGCAAAACCCGGGGGGCggttgg                                                                                                                    >  2:560446/1‑145 (MQ=255)
                                     tataAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGGACGGGGGTGTTTTATTCGGGATATCTTTTTTGCCGCCAAACCCGGGGGGCGtttggtt                                                                                                                  >  1:400294/1‑146 (MQ=255)
                                      ataaagataaagCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGGACGGGGGTGGTTTATGCGGGATATCTTTTTTGCCGCAAAACCCGGGGGGCGGTCGAttt                                                                                                                 >  2:890477/1‑151 (MQ=255)
                                          agataaagCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGGACGGGGGCGTTTTGTGGGGGATATTTTGTTGCCCGCCAAACCCGGGGGGCGGTCGGtttt                                                                                                                >  2:785330/1‑148 (MQ=255)
                                            ataaagCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGGACGGGGGTGTTTTTTGCGGGATATCTTGTTGGCCGCAAAACCTGGGGGCCGGTCGGTTTTAcac                                                                                                            >  2:347935/1‑148 (MQ=255)
                                            ataaagCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGGAAGGGGGTGGTTTGTGGGGGATATTTTTTTGGCCGCAAAACCTGGGGGGCGTTCGGttttttgc                                                                                                            >  2:425078/1‑146 (MQ=255)
                                             taaagCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGGACGGGGGTGTTTTGTGCGGGATATCTTTTTGGCCGCAAAACCTGGGGGGCGTTCGGTTTTACg                                                                                                             >  1:195095/1‑148 (MQ=255)
                                              aaagCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGGTCCGGGGCGGTTTGTGGGGGATATCCTGTTGGCAGCAAAGCCTGGGGGGCGTTCGGTTTTTCGCTc                                                                                                          >  2:282151/1‑150 (MQ=255)
                                               aagCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGGTATCGGCGGGACGGGGGTGTTTGGTGGGGGATATTTTTTTTGCCGGAAAACCTGGGGGGCGGTTGGTTTTTTGCTCgg                                                                                                        >  2:853495/1‑151 (MQ=255)
                                                agCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGGACGGGGGCGGTTTGTGGGGGATATCTTGTTGGCCGCAAAAGCCGGGGGGCGGTTGGTTTTTAGCTCgg                                                                                                        >  2:213061/1‑150 (MQ=255)
                                                agCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGGACGGGGGTGTTTTGTGGGGGATATTTTTTTGGCAGCAAAACCTGGGGGGCGGTCGGTTTTTCGCTccg                                                                                                        >  2:687063/1‑148 (MQ=255)
                                                agCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGGACGGGGGGGGTTTGTGCGGGGTATTTTGTTGGCAGGCAAGACTGGGGGGCGGTTGGtttttagacc                                                                                                          >  2:500023/1‑142 (MQ=255)
                                                  ccTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGGACGGGGGTGTTTGATGCGGGATATCTTGTTGCCCGCAAAACCCGGGGGGCGGTCGATTTTACGCCCGGcc                                                                                                      >  1:889400/1‑150 (MQ=255)
                                                            cAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGGTATCGGCGGGACGGGGGCGTTTTGTGGGGGAAATCTTGTTGCCCGCCAAAGCTGGGGGGCGTTTGATTTTTTGCTCGGGCTGGGGTtgggg                                                                                           >  2:740206/1‑151 (MQ=255)
                                                            cAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCctgggctggg                                                                                             <  2:402162/149‑1 (MQ=255)
                                                            cAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCctgggctggg                                                                                             <  1:967481/149‑1 (MQ=255)
                                                             aaaCTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGt                                                                                          <  1:875122/151‑1 (MQ=255)
                                                             aaaCTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGt                                                                                          <  1:801402/151‑1 (MQ=255)
                                                              aaCTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTa                                                                                         <  1:407770/151‑1 (MQ=255)
                                                                   cGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTACCGGGCTGTTTTGTGCGGGATATCTTGTTGCCAGCAAAACCTGGGGGCCGTTTGATTTTTCGCTCGGCCTTGGGTgggggggttt                                                                                      >  1:869678/1‑144 (MQ=255)
                                                                        gtgGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGGGCGGGGTATCGGCGGTTCCGGGGCTGTTGTATCGGGATAACCTTTTGTCAGCAAACGCCGGGGGCCCGTTGATTTTTCGCTTCGCCCGGGGGGGGGGGGCTTGGGGAcc                                                                               >  2:505187/1‑151 (MQ=255)
                                                                          gggAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTACGGGGGCGTTTTATGCGGGATATTTTGTTGCCCGCAAAACCCGGGGGGCGGTCGGTTTTTCCCCCTGCCTGGGCTGGGGGGATTTTGGCAcc                                                                               >  2:140360/1‑149 (MQ=255)
                                                                            gAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCg                                                                           <  2:732786/151‑1 (MQ=255)
                                                                             aaGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCgg                                                                          <  1:754558/151‑1 (MQ=255)
                                                                               gagagTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTa                                                                        <  2:967009/151‑1 (MQ=255)
                                                                                  agTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAAtg                                                                     <  1:967389/151‑1 (MQ=255)
                                                                                  agTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGGACCGGGGTGTTTTGTGGGGGAAATCTTGTTTGCAGGCAAACCCGGGGGGCGGTCGGTTTTTCGCCCCGCCTGGGGTGGGGGGATTTTGGGACCAACGGGaa                                                                       >  2:948338/1‑149 (MQ=255)
                                                                                   gTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGGACGGGGGTGTTTGGTGCGGGATATTTTGTTGGCCGCAAAGCCTGGGGGGCGGTTGGTTTTACGCTCGGCCTGGGCTGGGGGGATTTGGGCACCAGCGGGaa                                                                       >  1:604031/1‑148 (MQ=255)
                                                                                       cTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGGACGGGGCTGGTTTATGCGGGATATCTTGTTGCCAGCAAAGCCTGGGGGGCGGTCGATTTTTAGCTCGGCCTGGGCTGGGGGGATTTTGGCACCAGCGGGAAtgtgga                                                                 >  1:77549/1‑148 (MQ=255)
                                                                                        tGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGGACGGGGGTGTTTTGTGGGGGAAATCTTGTTGCCAGCAAAACCTGGGGGGCGTTCGGTTTTTTGCTCGGCCTGGGGGGGGGGGATTTGGGCACCAGCGGGAATGTGgaa                                                                >  2:649703/1‑150 (MQ=255)
                                                                                         ggCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGaaaa                                                               <  1:62898/150‑1 (MQ=255)
                                                                                          gCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAAt                                                             <  1:697895/151‑1 (MQ=255)
                                                                                          gCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAAt                                                             <  1:12969/151‑1 (MQ=255)
                                                                                            tGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGGACGGGGGTGTTTTGTGGGGGAAATCTTGTTGGCAGCCAAACCTGGGGGGCGTTCGATTTTTCGCTCGGGCCGGGGTGGGGGGGTTTTGGCCCCCGCGGGAAAGGGAAAAATcc                                                           >  1:579178/1‑151 (MQ=255)
                                                                                               cGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGct                                                        <  1:200041/151‑1 (MQ=255)
                                                                                                   aGTGGCGGTTGGTGCGCGGGATATCGGCGGTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGt                                                    <  2:244018/151‑1 (MQ=255)
                                                                                                    gTGGCGGTTGGTGCGCGGGATATCGGCGGGACCGGGGGGTTTTGTGGGGGAAATTTTTTTGGCCGCCAAACCCGGGGGGCGGTTGGTTTTTAGCCCGGCCCTGGGTGGGGGGATTTTGGGACCCGCGGGAATGTGGAAAAACCGCCCTGtt                                                   >  1:156560/1‑151 (MQ=255)
                                                                                                      ggCGGTTGGTGCGCGGGATATCGGCGGTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTc                                                  <  1:425120/150‑1 (MQ=255)
                                                                                                      ggCGGTTGGTGCGCGGGATATCGGCGGTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCa                                                 <  2:672179/151‑1 (MQ=255)
                                                                                                        cGGTTGGTGCGCGGGATATCGGCGGTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAg                                                >  2:160789/1‑150 (MQ=255)
                                                                                                          gTTGGTGCGCGGGATATCGGCGGTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCa                                             <  1:495092/151‑1 (MQ=255)
                                                                                                            tGGTGCGCGGGATATCGGCGGTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGGATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGcc                                              >  2:435014/1‑148 (MQ=255)
                                                                                                                gcgcGGGATATCGGCGGTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGAt                                        <  1:689483/150‑1 (MQ=255)
                                                                                                                gcgcGGGATATCGGCGGTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATa                                       <  2:420359/151‑1 (MQ=255)
                                                                                                                 cgcgGGATATCGGCGGTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATaa                                      <  1:160801/151‑1 (MQ=255)
                                                                                                                    gggATATCGGCGGTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAata                                   <  1:953392/151‑1 (MQ=255)
                                                                                                                      gATATCGGCGGTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATAtt                                 <  1:776857/151‑1 (MQ=255)
                                                                                                                      gATATCGGCGGTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATAtt                                 <  1:775800/151‑1 (MQ=255)
                                                                                                                       atatCGGCGGTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAat                                    <  2:91449/147‑1 (MQ=255)
                                                                                                                        tatCGGCGGTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGGATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATAtt                                 >  1:812912/1‑149 (MQ=255)
                                                                                                                            ggcggTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATc                           <  1:147736/151‑1 (MQ=255)
                                                                                                                                 gTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGa                       <  2:400255/151‑1 (MQ=255)
                                                                                                                                    aCGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGAt                      <  1:453912/150‑1 (MQ=255)
                                                                                                                                        gggCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACa                  <  1:538377/151‑1 (MQ=255)
                                                                                                                                           cTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCt               >  1:762239/1‑151 (MQ=255)
                                                                                                                                             gTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTa              >  1:810093/1‑150 (MQ=255)
                                                                                                                                               ttGCTGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACa            >  1:555303/1‑150 (MQ=255)
                                                                                                                                               ttGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACa            <  1:524478/150‑1 (MQ=255)
                                                                                                                                                tGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGGTACaa           <  2:310399/150‑1 (MQ=255)
                                                                                                                                                tGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACaa           <  2:719714/150‑1 (MQ=255)
                                                                                                                                                tGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACaa           <  1:835051/150‑1 (MQ=255)
                                                                                                                                                  aTGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAAc         >  1:68200/1‑150 (MQ=255)
                                                                                                                                                      ggAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCg    <  1:560501/151‑1 (MQ=255)
                                                                                                                                                        aaTATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCggg  >  1:259124/1‑151 (MQ=255)
                                                                                                                                                        |                                                                                                                                                          
GTACAGCAGCGTCGAAGCGTTCTCTGGCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTACGGGGCTGTTTGATGCGGAATATCTTGTTGCCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGG  >  NC_000913/4237900‑4238199

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