Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913993,9130GC39.7% 137.0 / 20.8 74G44G (GGC→GGGssuBaliphatic sulfonate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (41/3);  new base C (1/28);  total (43/31)
Fisher's exact test for biased strand distribution p-value = 9.95e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTT  >  NC_000913/993770‑994058
                                                                                                                                                 |                                                                                                                                                 
aGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCccacca                                                                                                                                            <  2:1074802/151‑1 (MQ=255)
aGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCccacca                                                                                                                                            <  2:1139942/151‑1 (MQ=255)
  cAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCaccacc                                                                                                                                          >  2:973478/1‑151 (MQ=255)
   aGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCg                                                                                                                                         >  1:759794/1‑151 (MQ=255)
      cgcgCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGcc                                                                                                                                       >  1:1488769/1‑150 (MQ=255)
          cATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACaa                                                                                                                                   >  1:328750/1‑150 (MQ=255)
                gacaCATCATTCGCGTCTTTTCCTGAATTTCCGCCAGCGGTGTGGTGCCCGCTCACACATTGCCTGCGCTTGGCGCTTCCCAACCTGCCCGCCCGCCCAGCAGCGTACTTTTGCCACCACCGCTGCGCTTCACCACCGCCACAAACTGAcc                                                                                                                            >  2:1422874/4‑151 (MQ=255)
                aaaaCATCATTCGGGTTTCTCCTGAAATTCCGGCCGGGGGTGTGGGCCCCGCAACCACACCGCGGGGGGTTGGTTTTCCCACCCCCGCCAGCAGGCGCGGGGGGTTACTTTTCCCCCCCCCGCGGGGCCCCACCACCGCCACAAACTGAcc                                                                                                                            <  2:543608/150‑1 (MQ=255)
                 aaaCATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCt                                                                                                                           >  2:324668/1‑151 (MQ=255)
                 aaaCATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCt                                                                                                                           >  2:921544/1‑151 (MQ=255)
                 aaaCATCATTCGCACGCCTTCCTGAATTACAGCCACCAGTGTGGGGCCCGCTACCACATCGCCTGCGTTTAGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGCACTTTTGCGACCACCCATGCGTCCCACCACACCcacgcaccgac                                                                                                                             >  2:318919/1‑141 (MQ=255)
                  aacaccttccgcgTGTCTCCCTAATTTCCAGCCAGGGGTGTGGGCCCCGCAAACACACCCCCGGGGGTTGGCTTTCCCAGCCCCGCCAGCGGGCGCAGCGGGGTTCTTTTCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTg                                                                                                                          <  1:704553/142‑1 (MQ=255)
                   aCATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGc                                                                                                                         >  2:1337118/1‑151 (MQ=255)
                     atcatTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCg                                                                                                                       >  2:691892/1‑151 (MQ=255)
                               tgtCTTCCTAAATTCCAGCCAGCGGTGTGGTGCCCGCTAACACACCGCCTGGGTTGGGCTTTCCCAGCCCGGCCAGCGGGCGCAGCGGGTTTTTTTTCCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATg                                                                                                                <  1:872725/148‑1 (MQ=255)
                                 tCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAAt                                                                                                           >  1:335261/1‑151 (MQ=255)
                                  cTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATc                                                                                                          >  2:1350854/1‑151 (MQ=255)
                                   ttCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAAGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGGATATGTAAATcc                                                                                                         >  2:1424771/1‑151 (MQ=255)
                                     ccTGATTTTCAGCCGGGGGTGTGGGCCCCGCAAACCCATCGCCGGCGGTGGGCTTTCCCAGCCCTGCCAGCGGGCGCAGCGGGTTTCTTTTCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCa                                                                                                        <  1:1321814/150‑1 (MQ=255)
                                     ccTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATcc                                                                                                         >  1:470470/1‑149 (MQ=255)
                                        gAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGCACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTg                                                                                                    >  2:604537/1‑151 (MQ=255)
                                          attttCACCCAGGGGTGGGGTCCCCGCAAACCCACCGCCGGCGGTGGGCGTTCCCAGCCCTGCCAGAAGGCGCGGCGGGTTCTTTTTCCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATg                                                                                                                <  1:843902/136‑1 (MQ=255)
                                          atttccacccgggggTGGGGTCCCCGCAAACACCCCCCCTGGGTTTGGCTTTCCCAGCCCTCCCGGCGGGCGCAGCGGGTTCTTTTGCCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTg                                                                                                    <  1:1565428/137‑1 (MQ=255)
                                             ttCAGCCCGCCGTGTGGTGCCCGCTAACACATCTCCTGCGCTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTCCCAGCACCGCTCCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATTTAAATCCAGTTGGtt                                                                                                 >  2:1092021/1‑150 (MQ=255)
                                              tCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTc                                                                                                >  1:668486/1‑150 (MQ=255)
                                               cagccagcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTc                                                                                                >  1:559599/1‑149 (MQ=255)
                                               cagccagcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTc                                                                                                >  1:520790/1‑149 (MQ=255)
                                                  ccgggggggtggtcCCCCCAAACCCTCCGCCTGGGGTGGGCGTTCCCACCCCTGCCAGCAGGCGCAGCGGGGTATTTTTCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTc                                                                                                <  1:93316/138‑1 (MQ=255)
                                                  cccggggTGTGGCCCCCCATAACACCCCGCCTGGGTTGGGCTTTCCCACCCCCGCCAGCGGGCCCAGCGGGGTCCTTTCCCCACCCCCGCTCCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTc                                                                                                <  2:1313777/141‑1 (MQ=255)
                                                  ccagcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGAc                                                                                           >  2:775897/1‑151 (MQ=255)
                                                   cggggggggggtcCCCCCAAACCCACCGCCGGGGGTTGGCTTTCCCACCCCGCCCGGCGGGCGCAGAGGGGTTTTTTCGCCCCCACCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGa                                                                                            <  2:696532/140‑1 (MQ=255)
                                                   cagcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACg                                                                                          >  1:575039/1‑151 (MQ=255)
                                                   cagcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACg                                                                                          >  1:1513982/1‑151 (MQ=255)
                                                   cagcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACg                                                                                          >  1:1135512/1‑151 (MQ=255)
                                                    agcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACg                                                                                          >  2:141038/1‑150 (MQ=255)
                                                       gggggggtcCCCGCAACAACACCGCCGGGGTTTGGCTTTCCCACCCCTGCCAGCAGCCGCAGCGGGTTTCTTTTCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGatat                                                                                      <  1:1369953/146‑1 (MQ=255)
                                                                   gCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAAt                                                                          >  2:1325996/1‑151 (MQ=255)
                                                                      aaCACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCAGTTGGTTCAGGACGATATTTTCCGCGTAATGttt                                                                      >  1:182590/1‑151 (MQ=255)
                                                                          caTCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCt                                                                   >  2:1211086/1‑151 (MQ=255)
                                                                           attgCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCtt                                                                  >  2:110224/4‑151 (MQ=255)
                                                                           aTCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCtt                                                                  >  2:68833/1‑151 (MQ=255)
                                                                             cGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTAc                                                                >  1:1094058/1‑151 (MQ=255)
                                                                              gtcTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCAAAAACTGACCTGACGGAATATGTAAATCCAGTTGGCTCAGGACGACATTTTCCGCGTAATGTTTGCTTACt                                                               >  2:1397838/3‑151 (MQ=255)
                                                                              gCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACt                                                               >  2:68749/1‑151 (MQ=255)
                                                                                  gcgttgggCTTTCCCAGCCCGGCCAGCAGGCGCAGCAGGTTTTTTTTGCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTg                                                              <  1:115427/142‑1 (MQ=255)
                                                                                    gtttGGCTTTCCCAGCCCTCCCAGCGGGCGCAGCGGGTTCTTTTTCCCCCCCCCGCTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTg                                                              <  1:896156/144‑1 (MQ=255)
                                                                                    gtttGGCTTTCCCACCCCGGCCAGCGGGCGCGGCGGGGTTTTTTTCCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTg                                                         <  1:1376251/149‑1 (MQ=255)
                                                                                      tgggCTTTCCCAGCCCTGCCGGCGGGCGCAGCGGGGTTTTTTTCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGc                                                                    <  1:416275/136‑1 (MQ=255)
                                                                                       tGGAGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGc                                                      >  2:242831/1‑151 (MQ=255)
                                                                                            ttcccAGCCCTGCCGGCGGGCGGAGCGGGGTACTTTGCCCCCCCCCGGTGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaac                                                   <  1:657755/146‑1 (MQ=255)
                                                                                                              ggCGCGGCGGGGTTTTTTGCCCCCCCCCGCCGCGCCCCCCCTCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCa                                <  1:1211850/150‑1 (MQ=255)
                                                                                                               gcgcAGCGGGGTACTTTTCCCCCCCCCGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAg                               <  1:640884/150‑1 (MQ=255)
                                                                                                                cgcAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGAc                             >  1:189277/1‑151 (MQ=255)
                                                                                                                  cagcagGTTTCTTCCGCCCCCCCCGCGGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGa                              <  2:987686/148‑1 (MQ=255)
                                                                                                                   agcagGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACg                            >  1:343896/1‑149 (MQ=255)
                                                                                                                   agcagGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACg                            >  2:388676/1‑149 (MQ=255)
                                                                                                                     gggggTATTTTTGCCCCCCCCCGCGCCGCCCCGCCGCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACg                            <  2:759781/146‑1 (MQ=255)
                                                                                                                     cagGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGc                         >  1:1107473/1‑150 (MQ=255)
                                                                                                                      agGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGc                         >  1:777452/1‑149 (MQ=255)
                                                                                                                       ggttttttttCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGa                              <  1:820682/137‑1 (MQ=255)
                                                                                                                       gggtctttttCCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGc                         <  1:793733/142‑1 (MQ=255)
                                                                                                                       gggtctttgcgCCCCCCCCGCGGCCCCCCACCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGAc                             <  2:1483730/134‑1 (MQ=255)
                                                                                                                          tatttttCCCCCCCCGGCTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTa                     <  2:17225/146‑1 (MQ=255)
                                                                                                                          tACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTAtt                   >  1:502254/1‑151 (MQ=255)
                                                                                                                              tttgcccccccccGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTAtt                   <  1:489728/137‑1 (MQ=255)
                                                                                                                              tttGCCACCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTc                  <  1:111528/148‑1 (MQ=255)
                                                                                                                              tgtacccccccccgccgCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCa                 <  1:763864/139‑1 (MQ=255)
                                                                                                                               ttgcccccccccGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTc                  <  1:1145170/138‑1 (MQ=255)
                                                                                                                                 gCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCg            >  2:1210199/1‑151 (MQ=255)
                                                                                                                                   ccccccccGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCg            <  1:1096202/144‑1 (MQ=255)
                                                                                                                                       ccccGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcc       <  1:163290/149‑1 (MQ=255)
                                                                                                                                       caccGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcc       >  1:7930/1‑151 (MQ=255)
                                                                                                                                           gCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcctcct   >  2:916311/1‑151 (MQ=255)
                                                                                                                                            cTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCtt  >  1:1521411/1‑151 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
AGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTT  >  NC_000913/993770‑994058

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: