Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 161,733 A→C 15.6% intergenic (‑228/+18) PP_0152 ← / ← PP_0153 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947161,7330AC15.6% 934.9 / 49.5 382intergenic (‑228/+18)PP_0152/PP_0153hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (213/96);  new base C (6/64);  total (219/163)
Fisher's exact test for biased strand distribution p-value = 1.88e-21
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 383 total aligned reads displayed.

CTGTGCCAGACAACTGTTAAGGCGAATTTTAATTTTTGTAAAACTGTACTGGTCGCAATCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCC‑GCTGCAGTCACCTGCAACGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTCCAGGCGAGCCAGTTGACGACGGCTGGTCATACGGCGCTGCCACAGCAGCAGGGT  >  NC_002947/161586‑161882
                                                                                                                                                     |                                                                                                                                                     
 tgtgCCAGACAACTGTTAAGGCGAATTTTAATTTTTGTAAAACTGTACTGGTCGCAATCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCC‑GCTGCAGTCACCTGCAACGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGcc                                                                                                                                                    <  2:3919089/150‑1 (MQ=255)
   tgCCAGACAACTGTTAAGGCGACTTTTAATTTTTGTAAAACTGTACTGGTCGCAATCTGAAGCACGCGACTGCGTGTCTCTCGGGCAAGAAAAAACCC‑GCTGCAGTCACCTGCAACGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGCCgg                                                                                                                                                  <  2:1636567/150‑1 (MQ=255)
   tgCCAGACAACTGTTAAGGCGAATTTTAATTTTTGTAAAACTGTACTGGTCGCAATCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCC‑GCTGCAGTCACCTGCAACGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGCCgg                                                                                                                                                  <  2:558006/150‑1 (MQ=255)
    gCCAGACAACTGTTAAGGCGAATTTTAATTTTTGTAAAACTGTACTGGTCGCAATCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAATAACCC‑GCTGCAGTCACCTGCAACGGG‑TTTGAGAAATTCACGGTGTCGGCGAAGCCGGc                                                                                                                                                 >  1:5199470/1‑150 (MQ=255)
     ccAGACAACTGTTAAGGCGAATTTTAATTTTTGTAAAACTGTACTGGTCGCAATCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCC‑GCTGCAGTCACCTGCAACGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGCCGGCa                                                                                                                                                >  1:2208324/1‑150 (MQ=255)
       aGACAACTGTTAAGGCGAATTTTAATTTTTGTAAAACTGTACTGGTCGCAATCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAATACCC‑GCTGCAGTCACCTGCAACGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGCCGGCAAg                                                                                                                                              >  2:5348945/1‑150 (MQ=255)
       aGACAACTGTTAAGGCGAATTTTAATTTTTGTAAAACTGTACTGGTCGCAATCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCC‑GCTGCAGTCACCTGCAACGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGCCGGCAAg                                                                                                                                              >  2:6028411/1‑150 (MQ=255)
           aaCTGTTAAGGCGAATTTTAATTTTTGTAAAACTGTACTGGTCGCAATCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCC‑GCTGCAGTCACCTGCAACGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGCCGGCAAGCCgg                                                                                                                                          >  2:8562982/1‑150 (MQ=255)
              tGTTAAGGCGAATTTTAATTTTTGTAAAACTGTACTGGTCGCAATCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCC‑GCTGCAGTCACCTGCAACGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAg                                                                                                                                       >  2:7935952/1‑150 (MQ=255)
                ttAAGGCGAATTTTAATTTTTGTAAAACTGTACTGGTCGCAATCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCC‑GCTGCAGTCACCTGCAACGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCt                                                                                                                                     >  2:5294766/1‑150 (MQ=255)
                    ggCTAATTTTATTTTTTGTAAAACTGTACTGGTCGAAATCTGAAGCACGCGCCTGCGCGTCTCTCGGGCAATAAAAAACCC‑GCTGCAGTCACCTGCAACGGG‑TTTTAGAAATTCACGTTGCCGGCTAAGCCGGCAAGCCGGGAGCTTAgc                                                                                                                                 <  2:4328011/150‑1 (MQ=255)
                     gCGAATTTTAATTTTTGTAAAACTGTACTGGTCGCAATCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCC‑GCTGCAGTCACCTGCAACGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAgcg                                                                                                                                >  2:8743479/1‑150 (MQ=255)
                     gCGAATTTTAATTTTTGTAAAACTGTACTGGTCGCAATCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCC‑GCTGCAGTCACCTGCAACGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAgcg                                                                                                                                >  1:4708917/1‑150 (MQ=255)
                     gCGAATTTTAATTTTTGTAAAACTGTACTGGTCGCAATCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCC‑GCTGCAGTCACCTGCAACGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGCCGGCAAGCCCGGAGCTTAgcg                                                                                                                                >  2:7403880/1‑150 (MQ=255)
                      cGAATTTTAATTTTTGTAAAACTGTACTGGTCGCAATCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCC‑GCTGCAGTCACCTGCAACGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAgcgc                                                                                                                               >  1:1095592/1‑150 (MQ=255)
                      cGAATTTTAATTTTTGTAAAACTGTACTGGTCGCAATCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCC‑GCTGCAGTCACCTGCAACGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAgcgc                                                                                                                               >  2:6585085/1‑150 (MQ=255)
                       taattttaattttTGAAAAACTTTACTGCTCTAAATCTTAACCCCGCCCCGGCGCTTCTCTCGGCAAAAAAAAAACCC‑CCTGCCGTCACCTGCAACGGG‑TTTTACAAATACCCGGTCCCGGCGACGCCGGCAAGCCGGGAGCTTAGCGcc                                                                                                                              <  2:279948/149‑1 (MQ=255)
                       gAATTTTAATTTTTGTAAAACTGTACTGGTCGCAATCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCC‑GCTGCAGTCACCTGCAACGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGcc                                                                                                                              >  2:5659664/1‑150 (MQ=255)
                        aattttaattttTGTAAAACTGTACTGGTCGCAATCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCC‑GCTGCAGTCACCTGCAACGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCa                                                                                                                             >  1:8040217/1‑150 (MQ=255)
                            ttatttttTGTAAAACTGTACTGGTCGCAATCTGACGCCCGCGACTGCGCGTCTCTGGGGAAATAAAAACCCG‑GCTGCCGTCACCTGCAAGGGT‑TTTGACAAATCCACGGGGCCGGCGACGCCGGCAAGCCGGGAGCTTAGCGCCAGaaa                                                                                                                         <  2:7508849/150‑1 (MQ=255)
                               attttTTTAAAACTGTCCGGGTCGAAATCTGAACCACGCGACTGCGCGTCTCTCGGGAAATAAAAAACCC‑GCTGCCGACACCTGCACCGGG‑TTTGACAAATACACGGTGCCGGCGACGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGc                                                                                                                      <  1:5079879/150‑1 (MQ=255)
                               attttTGTAAAACTGTACTGGTCGCAATCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCC‑GCTGCAGTCACCTGCAACGGG‑TTTGAGAAATTCACGGTACCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGc                                                                                                                      <  1:5560484/150‑1 (MQ=255)
                                    tGTAAAACTGTACTGGTCGCAATCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCC‑GCTGCAGTCACCTGCAACGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCt                                                                                                                 >  2:4130813/1‑150 (MQ=255)
                                       aaaaCTGTCCTGGGCGCAATCTGACGCACGCGACTGCGCGTCTCTCGGGCAACAAAAAACCC‑GCGGCCGTCACCTG‑AACGGGTTTTGAGAAATTCCCGGTGCCGGCGACGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAg                                                                                                              <  2:991110/150‑1 (MQ=255)
                                        aaaCTGTAGTGGTCGCAATCTGAAGCCCGGGACGGCGCGTCGCTGGGGAAATAAAAAACCG‑GCTGCCGCCACCGGCAAGGGG‑TTTGCGAAATGCCCGGTGCCGGCGATGCGGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGc                                                                                                             <  2:429857/150‑1 (MQ=255)
                                         aaCTGTACTGGTCGCAATCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCC‑GCTGCAGTCACCTGCAACGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGct                                                                                                            >  1:4295947/1‑150 (MQ=255)
                                         aaCTGTACTGGTCGCAATCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCC‑GCTGCAGTCACCTGCAACGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGct                                                                                                            >  1:7393164/1‑150 (MQ=255)
                                            tGTACTGGTCGCAATCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCC‑GCTGCAGTCCCCTGCACCGGG‑TTTGAGAAATCCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCtt                                                                                                         <  2:1976722/150‑1 (MQ=255)
                                             gTACTGGTCGCAATCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCC‑GCTGCAGTCACCTGCAACGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTc                                                                                                        <  2:7302863/150‑1 (MQ=255)
                                             gTACTGGTCGAAATCTGAAGCACGCGACTGCGCGTCTCTGGGGAAATAAAAAACCC‑GCTGCAGTCACCTGCAAGGGG‑TTTGAGAAATACACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTc                                                                                                        <  1:6587359/150‑1 (MQ=255)
                                                 tGGTCGCAATCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCC‑GCTGCAGTCACCTGCAACGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTaa                                                                                                    >  1:124659/1‑150 (MQ=255)
                                                   gTCGCAATCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCC‑GCTGCAGTCACCTGCAACGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACg                                                                                                  >  2:7133160/1‑150 (MQ=255)
                                                   gTCGCAATCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCC‑GCTGCAGTCACCTGCAACGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACg                                                                                                  >  1:3554967/1‑150 (MQ=255)
                                                      gCAATCTGAAGCACGCCACTGCGCGTCTCTCGGGAAATAAAAAACCC‑GCTGAAGTCACCTGAAAGGGT‑TTTGATAAATTCACGGTGCCGGCCAAGCCGCCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTg                                                                                               <  2:5936074/150‑1 (MQ=255)
                                                          tCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCC‑GCTGCAGTCACCTGCAACGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGAc                                                                                           >  1:322636/1‑150 (MQ=255)
                                                            tGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCC‑GCTGCAGTCACCTGCAACGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTc                                                                                         >  1:1518521/1‑150 (MQ=255)
                                                             gAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCC‑GCTGCAGTCACCTGCAACGGG‑TTTGACAAATTCAGGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCg                                                                                        <  2:503001/150‑1 (MQ=255)
                                                                gCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCC‑GCTGCAGTCACCTGCAACGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTg                                                                                     <  2:8619011/150‑1 (MQ=255)
                                                                    gcgACTGCGCGTCTCTCGGGCAATAAAAAACCC‑GCTGCAGTCACCTGCAACGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATAc                                                                                 >  1:6166690/1‑150 (MQ=255)
                                                                    gcgACTGCGCGTCTCTCGGGCAATAAAAAACCC‑GCTGCAGTCACCTGCAACGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATAc                                                                                 >  1:5399468/1‑150 (MQ=255)
                                                                       acggcgcGACTCTCGGGAAATAAAAAACCC‑GCGGCCGTCACGTGAACCGGT‑TTTTATAATTTCAGGGTGCCGGCGACGCGGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCgg                                                                              <  2:6096061/147‑1 (MQ=255)
                                                                              gTCTCTCGGGCAATAAAAAACCC‑GCTGCAGTCACCTGCAAGGGG‑TTTGAGAATTTCACGGTGCGGGCGACGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGc                                                                       <  2:8430467/150‑1 (MQ=255)
                                                                                 tctcGGGCAATAAAAAACCC‑GCTGCAGTCACCTGCAACGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCcag                                                                    >  2:3542524/1‑150 (MQ=255)
                                                                                   tcGGGCAATAAAAAACCC‑GCTGCAGTCACCTGCAACGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCcagca                                                                  >  2:7983621/1‑150 (MQ=255)
                                                                                   gcggcaaaTAAAAACCCC‑GCGGCCGTCACCTCCACGGGG‑TTTGAGAACTCCCGGGGGCCGGCCACGCCGCCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCcagca                                                                  <  1:418827/144‑1 (MQ=255)
                                                                                    cGGGCAATAAAAAACCC‑GCTGCAGTCACCTGCAACGGG‑TTTGAGAAATTCACGGTGCCGGCGACGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCcagcag                                                                 <  2:1346389/150‑1 (MQ=255)
                                                                                     gggCAATAAAAAACCC‑GCTGCAGTCACCTGCAACGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAgg                                                                >  1:4416800/1‑150 (MQ=255)
                                                                                     gggCAATAAAAAACCC‑GCTGCAGTCACCTGCAACGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCGTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAgg                                                                >  2:2653572/1‑150 (MQ=255)
                                                                                      ggaaaTAAAAAACCC‑GCTGCCGGCACCTGCAACGGG‑TTTGAGAAATTCCCGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGc                                                               <  2:8207551/147‑1 (MQ=255)
                                                                                      ggCAATAAAAAACCC‑GCTGCAGTCACCTGCAACGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGc                                                               >  1:6107032/1‑150 (MQ=255)
                                                                                       gCAATAAAAAACCC‑GCTGCAGTCACCTGCAACGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCg                                                              >  1:5667092/1‑150 (MQ=255)
                                                                                       caaaTAAAAACCCC‑GCTGCTGTCACCTGAAAGGGT‑TTTGAGAAATTCACGGGGCCGGCGAAGCCGCCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCg                                                              <  1:8514479/148‑1 (MQ=255)
                                                                                        gaaaaaaaaaCCCGGCTGCAGTCACCTGCACGGGT‑TTTTAGAATTTCCCGGGGCGGGCGACGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCg                                                              <  1:579205/146‑1 (MQ=255)
                                                                                        caaacaaaaaCCC‑GCGGCCGGCCCCGGCACGGGG‑TTTCCGAACTCCCCGGGGCGGGCGCCGCGGCCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCgg                                                             <  2:18499/145‑1 (MQ=255)
                                                                                        cAATAAAAACCCC‑GCGGCCGGCACCTGCACCGGG‑TTTGAGAAATTCCCGGTGCCGGCGACGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCgg                                                             <  2:6225174/150‑1 (MQ=255)
                                                                                        cAATAAAAAACCC‑GCTGCAGTCACCTGCAACGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCgg                                                             >  1:6887970/1‑150 (MQ=255)
                                                                                        cAAAAAAAAACCC‑GCTGCAGTCACCTGCACGGGG‑TTTGCTAAATTCCCGGTGCCGGCTAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCgg                                                             <  2:4836956/150‑1 (MQ=255)
                                                                                          aTAAAAAACCC‑GCTGCAGTCACCTGCAACGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGa                                                           >  2:3127445/1‑150 (MQ=255)
                                                                                          aTAAAAAACCC‑GCTGCAGTCACCTGCAACGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGa                                                           >  2:3802857/1‑150 (MQ=255)
                                                                                           tctcaaaCCC‑GTTGCAGGTGACTGCAGCGGG‑TTTTAGAAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAg                                                          >  1:2683517/5‑150 (MQ=255)
                                                                                           tAAAAAACCC‑GCTGCAGTCACCTGCAACGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAg                                                          >  1:4894197/1‑150 (MQ=255)
                                                                                           tAAAAAACCC‑GCTGCAGTCACCTGCAACGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAg                                                          >  2:6751258/1‑150 (MQ=255)
                                                                                           tAAAAAACCC‑GCTGCAGTCACCTGCAACGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAg                                                          >  1:1982835/1‑150 (MQ=255)
                                                                                           tAAAAAACCC‑GCTGCAGTCACCTGCAACGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAg                                                          >  1:901307/1‑150 (MQ=255)
                                                                                           tAAAAAACCC‑GCTGCAGTCACCTGCAACGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAg                                                          >  1:76307/1‑150 (MQ=255)
                                                                                           tAAAAAACCC‑GCTGCAGTCACCTGCAACGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAg                                                          >  1:1645997/1‑150 (MQ=255)
                                                                                            aaaaaaCCC‑GCTGCAGTCACCTGCAACGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGt                                                         >  1:1443604/1‑150 (MQ=255)
                                                                                             tcaaaCCC‑GTTGCAGGTGACTGCAGCGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTc                                                        >  2:7189317/3‑150 (MQ=255)
                                                                                             aaaaaCCC‑GCTGCAGACACCTGCAACGGT‑TTTGAGAAATACACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTc                                                        <  1:2868596/150‑1 (MQ=255)
                                                                                              aaaaCCC‑GCTGCAGTCACCTGCAACGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTcc                                                       >  1:6737124/1‑150 (MQ=255)
                                                                                              aaaaCCC‑GCTGCAGTCACCTGCAACGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTcc                                                       >  2:229861/1‑150 (MQ=255)
                                                                                               aaaCCC‑GCTGCTGTCACCTGCATCGGG‑TTTGAGAAATTCAGGGTGCGGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTCCa                                                      <  2:7796073/150‑1 (MQ=255)
                                                                                                aaCCC‑GTTGCAGGTGACTGCAGCGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTCCAg                                                     >  2:5786431/1‑150 (MQ=255)
                                                                                                    c‑GCTGCAGTCACCTGCAACGGT‑TTTGAGAAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTCCAGGCGa                                                 <  2:4683446/150‑1 (MQ=255)
                                                                                                       cTGCAGCCACCTGCACGGGT‑TTTGAGAAATTCACGGTGCCGGCGACGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTCCAGGCGAGc                                               <  2:2430395/150‑1 (MQ=255)
                                                                                                            gTCACCTGCACCGGG‑TTTGAGAAATTCCCGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTCCAGGCGAGCCAGtt                                          <  2:6892870/150‑1 (MQ=255)
                                                                                                               aCCTGCAACGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTCCAGGCGAGCCAGTTgac                                       >  2:1072096/1‑150 (MQ=255)
                                                                                                                 cTGCAACGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTCCAGGCGAGCCAGTTgacga                                     >  2:548228/1‑150 (MQ=255)
                                                                                                                    aaaCGGG‑TTTTATAATTTCACGGTGCCGGCTATGCCGGCAAGCCGTGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTCCAGGCGAGCCAGTTGacgacgt                                  <  2:93449/149‑2 (MQ=255)
                                                                                                                     aaCGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTCCAGGCGAGCCAGTTGACGACGGc                                 >  1:279497/1‑150 (MQ=255)
                                                                                                                       cGGG‑TTTGACAAATTCAGGGTGCCGGCGACGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTCCAGGCGAGCCAGTTGACGACGGCTg                               <  2:2208324/150‑1 (MQ=255)
                                                                                                                          gttttGATAAATTCTGGGTGCCGGCGCAGCGGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTCCAGGCGAGCCAGTTGACGACGGCTGGt                             <  1:3377216/148‑1 (MQ=255)
                                                                                                                          g‑TTTGAGAAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTCCAGGCGAGCCAGTTGACGACGGCTGGTc                            >  1:7123804/1‑150 (MQ=255)
                                                                                                                            tttGAGAAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTCCAGGCGAGCCAGTTGACGACGGCTGGTCa                           >  2:4740736/1‑150 (MQ=255)
                                                                                                                              tGAGAAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGAAGTCCAGGCGAGCCAGTTGACGACGGCTGGTCATa                         <  2:653481/150‑1 (MQ=255)
                                                                                                                                agaAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTCCAGGCGAGCCAGTTGACGACGGCTGGTCATACg                       >  2:7003350/1‑150 (MQ=255)
                                                                                                                                 gaAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTCCAGGCGAGCCAGTTGACGACGGCTGGTCATACgg                      >  1:2017588/1‑150 (MQ=255)
                                                                                                                                 gaAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTCCAGGCGAGCCAGTTGACGACGGCTGGTCATACgg                      >  1:6482936/1‑150 (MQ=255)
                                                                                                                                  aaaTTCAGGGGGCGGGCGAAGCCGCCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTCCAGGCGAGCCAGTTGACGACGGCTGGTCATACGgc                     <  1:4001057/150‑1 (MQ=255)
                                                                                                                                   aaTTCTCGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTCCAGGCGAGCCAGTTGACGACGGCTGGTCATACGgcg                    <  1:3121497/150‑1 (MQ=255)
                                                                                                                                    aTTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTCCAGGCGAGCCAGTTGACGACGGCTGGTCATACGgcgc                   >  2:1117165/1‑150 (MQ=255)
                                                                                                                                     tccACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTCCAGGCGAGCCAGTTGACGACGGCTGGTCATACGGCGCt                  <  2:3365696/148‑1 (MQ=255)
                                                                                                                                      tCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTCCAGGCGAGCCAGTTGACGACGGCTGGTCATACGGCGCTg                 >  2:5612389/1‑150 (MQ=255)
                                                                                                                                       cACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTCCAGGCGAGCCAGTTGACGACGGCTGGTCATACGGCGCTGc                >  1:1867421/1‑150 (MQ=255)
                                                                                                                                        aCGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGTGTCTGCCAGCAGGCGGGAGTCCAGGCGAGCCAGTTGACGACGGCTGGTCATACGGCGCTGcc               >  2:3066495/1‑150 (MQ=255)
                                                                                                                                         cGGTGCGGGCGACGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTCCAGGCGAGCCAGTTGACGACGGCTGGTCATACGGCGCTGCca              <  2:3175609/150‑1 (MQ=255)
                                                                                                                                           ttGCCGGCTTCGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTCCAGGCGAGCCAGTTGACGACGGCTGGTCATACGGCGCTGCcaca            >  1:2429420/2‑150 (MQ=255)
                                                                                                                                           ttGCCGGCTTCGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTCCAGGCGAGCCAGTTGACGACGGCTGGTCATACGGCGCTGCcaca            >  1:6922084/2‑150 (MQ=255)
                                                                                                                                                ggCGAAGCGGGCAAGCCGGGAGCTTAGCGCCAG‑AAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTCCAGGCGAGCCAGTTGACGACGGCTGGTCATACGGCGCTGCCAcagcagca      <  2:1525904/150‑1 (MQ=255)
                                                                                                                                                   gAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTCCAGGCGAGCCAGTTGACGACGGCTGGTCATACGGCGCTGCCACAGCAGCAgg    <  2:8065765/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTGTGCCAGACAACTGTTAAGGCGAATTTTAATTTTTGTAAAACTGTACTGGTCGCAATCTGAAGCACGCGACTGCGCGTCTCTCGGGCAATAAAAAACCC‑GCTGCAGTCACCTGCAACGGG‑TTTGAGAAATTCACGGTGCCGGCGAAGCCGGCAAGCCGGGAGCTTAGCGCCAGAAAGGCTTGCTCAGCTCTTCGTAACGCTGCGACTCGCTGATACCGGCGTCTGCCAGCAGGCGGGAGTCCAGGCGAGCCAGTTGACGACGGCTGGTCATACGGCGCTGCCACAGCAGCAGGGT  >  NC_002947/161586‑161882

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑