Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913993,9130GC41.6% 161.6 / 37.5 89G44G (GGC→GGG) ssuBaliphatic sulfonate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (48/4);  new base C (0/37);  total (48/41)
Fisher's exact test for biased strand distribution p-value = 6.81e-21
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAAT  >  NC_000913/993763‑994063
                                                                                                                                                       |                                                                                                                                                      
ccACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCgg                                                                                                                                                        >  3:177485/1‑151 (MQ=255)
  aCGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGcc                                                                                                                                                      >  3:630295/1‑151 (MQ=255)
     gTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCccac                                                                                                                                                   >  3:151464/1‑151 (MQ=255)
       aGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCccacca                                                                                                                                                 <  4:329611/151‑1 (MQ=255)
        gCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCcaccac                                                                                                                                                <  3:950994/151‑1 (MQ=255)
         cAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGCGGTGCCCGCTAACACATCGCCTGCAGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCaccacc                                                                                                                                               >  4:429362/1‑151 (MQ=255)
          aGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCg                                                                                                                                              >  1:437387/1‑151 (MQ=255)
          aGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCg                                                                                                                                              <  2:1045668/151‑1 (MQ=255)
            aCGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGcc                                                                                                                                            <  2:447208/151‑1 (MQ=255)
              gcgcATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCcac                                                                                                                                          >  2:908469/1‑151 (MQ=255)
                   tCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACt                                                                                                                                     >  3:931362/1‑151 (MQ=255)
                    ccTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTg                                                                                                                                    >  1:49441/1‑151 (MQ=255)
                     cTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTg                                                                                                                                    >  4:472093/1‑150 (MQ=255)
                      tGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGAc                                                                                                                                  >  2:770416/1‑151 (MQ=255)
                      tGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGAc                                                                                                                                  >  3:279577/1‑151 (MQ=255)
                       aaaaCATCTTTGGCGTGTCTCCCTGATTTTCAGCGAGGGGGGGGGTCCCCTCAACCACACCGCTGGGGTTTGGCTTTCCCAGCCCTGCCAGCAGGCGCAGCAGGGTTTTTTTG‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGAc                                                                                                                                  <  2:437912/149‑1 (MQ=255)
                        aaccatcatCCGCGTGTCTTCCTAAATTCCAGCCGGCGGTGGGGTGCCCGTTACCACACCGCCTGCGTTGGGCGTTTCCAGCCCTGCCAGCAGGCGCAGCGGGGTCCTTTTG‑CCACCCCCGCTGCGCCCCCCCACCGCCACAAACTGAcc                                                                                                                                 <  1:105897/147‑1 (MQ=255)
                        aaaCATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCt                                                                                                                                >  3:561358/1‑151 (MQ=255)
                            atcatTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCg                                                                                                                            >  2:788881/1‑151 (MQ=255)
                              catTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCgg                                                                                                                           >  2:1329788/1‑150 (MQ=255)
                                  cgcgTTTCTCCCTAATTTCCACCCAGGGGTGGGGTCCCCGCTAACCCACCCCCGGGGGTTGGCTTTCCCAGCCCTGCCAGCAGGCGCAGCGGGTTATTTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGa                                                                                                                          <  1:349693/147‑1 (MQ=255)
                                       gtctccgtaaaTTCCACCCGGGGGTGTGGCCCCCGCAACCCCCCCGCCTGCGGTGGGCGTTCCCAGCCCTGCCAGCGGGCGCGGGGGGTTTTTTTTC‑CCCCCCCCGCGCCGCCCCCCCACCGCCACAAACTGACCTGCCGGAata                                                                                                                       <  2:929294/136‑1 (MQ=255)
                                       gtCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTaa                                                                                                                  >  2:33009/1‑150 (MQ=255)
                                            cttaattttCACCCGGGGGTGGGGTGCCCTCTAACACATCGCCTGGGTTGGGCTTTCCCAGCCCTGCCAGCGGGCGCAGCGGGGTTTTTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAAt                                                                                                                <  2:848481/141‑1 (MQ=255)
                                            ccTGAATTTCGCCCAGGGGTGGGGTGCCCCAAACCACAGCGCCGGGGGTTGGCTTTCCCGCCCCCCCCAGAGGGCGCAGCGGGGTTTTTTTC‑CCCCCCCCCCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCa                                                                                                             <  1:307549/150‑1 (MQ=255)
                                              tGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGtt                                                                                                          >  3:445798/1‑151 (MQ=255)
                                                agtttCACCCGGCGGTGTGGGCCCCGAAAACACACCGCCGGGGGTGGGCTTTCCCACCCCCCCCAGGAGGCGAGGCGGTTTTTTTTGG‑CCCCCCCCGCGGGGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCa                                                                                                             <  2:586891/144‑1 (MQ=255)
                                                aaTTCCGGCCGGGGGTGGGGTGCCCGCAAACCCACCGCCTGGGGTGGGCTTTCCCAGCCCTGCCGGCAGGCGCAGAGGGGTCTTTTCC‑CCACCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTgg                                                                                                        <  4:15493/151‑1 (MQ=255)
                                                 ttttCACCCGGGGGTGTGGCGCCCTCAACCCCGCCGCGGGCGGTTGGGTTTCCCAGCCCCGCCAGCGGGCGCAGTGGGTTACTTTTG‑CCCCCACCGCTGCGCCCCCCCTCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGt                                                                                                           <  3:1280459/146‑1 (MQ=255)
                                                  tttCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGtt                                                                                                      >  2:348535/1‑151 (MQ=255)
                                                  tttCACCCGGGGGTGGGGTGCCCGCTAACCCATCGCCTGGGTTTGGCTTTTCCAGCCCTGCCAGCAGGCGCAGCAGGTTTTTTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGt                                                                                                       <  2:612653/150‑1 (MQ=255)
                                                     cagccagcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTc                                                                                                     >  3:1117383/1‑149 (MQ=255)
                                                      agccagcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAgg                                                                                                  >  1:132154/1‑151 (MQ=255)
                                                       gccagcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAgg                                                                                                  >  3:180589/1‑150 (MQ=255)
                                                        ccgggggtgggggccccGCTACCCCATCGCCGGGGGTTGGCTTTTCAAGCCCTGCCAGCGGGCGCAGCGGGGTTCTTTTC‑CCCCCCCCGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTc                                                                                                     <  2:919943/132‑1 (MQ=255)
                                                        ccgggggTGTGGGCCCCCCAAACACAGCGCCGGGGTTTGGCTTTCCCAGCCCCGCCAGCAGGCGAAGCGGGTTATTTTCC‑CCCCCACCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAgg                                                                                                  <  4:271106/144‑1 (MQ=255)
                                                             ggTGTGGCCCCCGCAAACACACCGCCTGCGTTTGGCTTTCCCAGACCTGCCAGCAGGCGCAGAAGGTTTTTTTTC‑CCCCCCCCGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGatat                                                                                           <  1:160072/151‑1 (MQ=255)
                                                               tgtgGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAttt                                                                                         >  2:1064777/1‑151 (MQ=255)
                                                                 tggtgCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAttt                                                                                         >  2:391672/1‑149 (MQ=255)
                                                                  gggccccGCAAACCCACCGCCGGGGTTGGGCGTTCCCAGCCCGGCCAGCAGGCGCAGCAGGGTCTTTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAtttt                                                                                        <  1:352447/145‑1 (MQ=255)
                                                                  gggccccGCAAACACTCCGCCGGGGGTGGGCTTTCCAAGCCCTCCCGGCGGGCGCAGCGGGTTACTTTGC‑CCCCCCCCGCGGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTcc                                                                                      <  2:851272/147‑1 (MQ=255)
                                                                      cccGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGt                                                                                  >  3:629682/1‑151 (MQ=255)
                                                                       ccacTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGt                                                                                  >  3:948326/4‑150 (MQ=255)
                                                                        cGCTACCCCACCCCCTGCGGTTGGCTTTCCCAGCCCTCCCAGCGGGCGCAGCAGGTTCCTTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTaa                                                                                <  1:438472/151‑1 (MQ=255)
                                                                        cGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTa                                                                                 >  2:927493/1‑150 (MQ=255)
                                                                              cacaTCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTg                                                                          >  2:59185/1‑151 (MQ=255)
                                                                               acaTCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGc                                                                         >  2:1391506/1‑151 (MQ=255)
                                                                                   cGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTAc                                                                     >  3:1271755/1‑151 (MQ=255)
                                                                                       tGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGc                                                                  >  3:763848/1‑150 (MQ=255)
                                                                                        gggtttGGCTTTCCCACCCCCGCCAGCGGGCGCGGCAGGGTTCTTTGC‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCAt                                                                <  3:348564/147‑1 (MQ=255)
                                                                                        gCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCAt                                                                >  2:1351095/1‑151 (MQ=255)
                                                                                         cGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCAt                                                                >  2:1037489/1‑150 (MQ=255)
                                                                                          gttgggCTTTCCCAGCCCTGCCAGCGGGCGCAGCGGGGTTCTTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCa                                                                 <  1:169020/144‑1 (MQ=255)
                                                                                          ggTGGGCGTTCCCAGCCCGGCCGGCGGGCGCAGCGGTTTTTTTTTC‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCAt                                                                <  2:1160567/149‑1 (MQ=255)
                                                                                           gTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAg                                                                                                            >  2:282453/1‑104 (MQ=255)
                                                                                           gTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTg                                                              >  2:488806/1‑150 (MQ=255)
                                                                                            ttgtcttctCCAGCCCTGCAAGCAGCCGCAGAAGGTTTCTTTTC‑CCACCCCGCCTCCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGa                                                             <  3:770451/142‑1 (MQ=255)
                                                                                            ttGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTg                                                              >  2:438615/1‑149 (MQ=255)
                                                                                            ttGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTg                                                              >  3:39803/1‑149 (MQ=255)
                                                                                             tGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGc                                                           >  3:1380793/1‑151 (MQ=255)
                                                                                              gccttccccAGCCCTGCCAGCGGGCGCAGCGGGGTCTTTTCC‑CCCCCCCCGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGa                                                             <  2:151446/141‑1 (MQ=255)
                                                                                                  ttcccagcccgcccAGCGGGCCCAGCGGGTTTTTTTCC‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaaca                                                       <  3:216950/138‑1 (MQ=255)
                                                                                                  ttcccAGCCCTGCCAGCAGGCGCAGCGGGTTATTTTCC‑CCCCCACCGCGGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTg                                                                   <  3:790512/135‑1 (MQ=255)
                                                                                                  tccccgccccTCCCCGCGGGCGCAGCGGGTTTCTTTCC‑CCACCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCAt                                                                <  3:831284/133‑1 (MQ=255)
                                                                                                   tcccAGCCCTCCCAGCGGGCGCAGCGGGTTATTTTTC‑CCCCCACCGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaacaa                                                      <  3:256768/148‑1 (MQ=255)
                                                                                                     ccagcccccccAGCGGGCGCAGCAGGGTCTTTTCC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATg                                                    <  1:241123/141‑1 (MQ=255)
                                                                                                         ccttcccggcagGCGCAGCGGGTTACTTTTC‑CCCCCCCCGCTGCCCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGc                                                  <  2:1129597/140‑1 (MQ=255)
                                                                                                           cTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGcc                                             >  3:995333/1‑151 (MQ=255)
                                                                                                           cTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGcc                                             >  2:447985/1‑151 (MQ=255)
                                                                                                              ccAGCAGGCCCAGCAGGGTACTTTTC‑CCCCCACCGCGGCCCCCGCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGc                                              <  4:164294/147‑1 (MQ=255)
                                                                                                                agcagGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGt                                        >  3:827366/1‑151 (MQ=255)
                                                                                                                    ggCGCAGCGGGGTACTTTTC‑CCCCCCCCCCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTc                                      <  1:240515/149‑1 (MQ=255)
                                                                                                                     gccccagagGGTTTTTTTGCCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTgg                                         <  3:253117/139‑1 (MQ=255)
                                                                                                                          gcagGTTACTTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGAc                                  <  2:1334934/147‑1 (MQ=255)
                                                                                                                            agGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGc                              >  1:259552/1‑149 (MQ=255)
                                                                                                                             ggttttttttC‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACg                                 <  2:141665/139‑1 (MQ=255)
                                                                                                                             gggTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAg                            >  3:1247269/1‑150 (MQ=255)
                                                                                                                              ggTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGt                           >  3:259141/1‑150 (MQ=255)
                                                                                                                               gTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTAt                         >  3:1270641/1‑151 (MQ=255)
                                                                                                                                 aCTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCAGGAATATGTAAATCCAGTTGGTTCAGGCCGATATTTTCCGCGTAATGTTTGCTTACTGCATTGACCAACAATGTCGTGCCCTGGTTCAGACGAGCAGTAtt                        >  4:287506/1‑150 (MQ=255)
                                                                                                                                    tttccccccccccGCGGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCa                      <  2:1255360/139‑1 (MQ=255)
                                                                                                                                    tttG‑CCCCCCCGGCGGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCAt                     <  2:445764/150‑1 (MQ=255)
                                                                                                                                     ttcgccccccccggcgcgcCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCAt                     <  2:715965/135‑1 (MQ=255)
                                                                                                                                      tcccccccaccGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCAt                     <  4:378120/143‑1 (MQ=255)
                                                                                                                                         ccccccccGCTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCg                 <  3:927527/144‑1 (MQ=255)
                                                                                                                                          cccccccGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGt                <  2:172060/145‑1 (MQ=255)
                                                                                                                                          caccacCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTg               >  3:268501/1‑151 (MQ=255)
                                                                                                                                               ccGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcctc          >  2:1309619/1‑151 (MQ=255)
                                                                                                                                                 gCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcctcct        >  2:938346/1‑151 (MQ=255)
                                                                                                                                                 gCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcctcct        >  3:1178661/1‑151 (MQ=255)
                                                                                                                                                   tGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTc      >  3:643167/1‑151 (MQ=255)
                                                                                                                                                      ggCCCCCCTCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcct           <  2:676673/143‑1 (MQ=255)
                                                                                                                                                      ggCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCaa    >  2:359034/1‑150 (MQ=255)
                                                                                                                                                      ggCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCATTCaa    >  3:578844/1‑150 (MQ=255)
                                                                                                                                                       gCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAAt  >  2:812712/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                                      
CCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAAT  >  NC_000913/993763‑994063

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