Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,061,025 T→C 5.2% intergenic (‑106/+43) PP_5312 ← / ← hupA hypothetical protein/DNA‑binding protein HU‑alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,061,0250TC 5.2% 294.5 / 3.5 96intergenic (‑106/+43)PP_5312/hupAhypothetical protein/DNA‑binding protein HU‑alpha
Reads supporting (aligned to +/- strand):  ref base T (43/48);  new base C (2/3);  total (45/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATCAACAGCCCCATGGCCCAGAGAAGAAAGCGGAACTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAGGTTCTTGCCAGGCTTGAAGGCCACGGTATTGCTGGCCTTGATCTTGACCGGTTC  >  NC_002947/6060878‑6061172
                                                                                                                                                    |                                                                                                                                                   
cATCAACAGCCCCATGGCCCAGAGAAGAAAGCGGAACTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCtt                                                                                                                                                    >  4:141409/1‑149 (MQ=255)
cATCAACAGCCCCATGGCCCAGAGAAGAAAGCGGAACTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCtt                                                                                                                                                    <  1:2881/149‑1 (MQ=255)
cATCAACAGCCCCATGGCCCAGAGAAGAAAGCGGAACTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCtt                                                                                                                                                    >  6:230434/1‑149 (MQ=255)
 aTCAACAGCCCCATGGCCCAGAGAAGAAAGCGGAACTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTg                                                                                                                                                   >  5:23333/1‑149 (MQ=255)
  tCAACAGCCCCATGGCCCAGAGAAGAAAGCGGAACTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTgg                                                                                                                                                  <  5:230434/149‑1 (MQ=255)
    aaCAGCCCCATGGCCCAGAGAAGAAAGCGGAACTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTgggg                                                                                                                                                >  4:147438/1‑149 (MQ=255)
    aaCAGCCCCATGGCCCAGAGAAGAAAGCGGAACTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTggg                                                                                                                                                 <  6:23333/148‑1 (MQ=255)
     aCAGCCCCATTGCCCAGAGAACAAAGCGTAACTTCATCAACTCGACTCGAATGAAAAGTGACTAATGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGgc                                                                                                                                               <  6:173446/149‑1 (MQ=255)
      cAGCCCCATGGCCCAGAGAAGAAAGCGGAACTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGgcg                                                                                                                                              >  6:31876/1‑149 (MQ=255)
         ccccATGGCCCAGAGAAGAAACCGGAACTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGccc                                                                                                                                           <  8:54816/149‑1 (MQ=255)
           ccATGGCCCAGAGAAGAAAGCGGAACTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGTTGACTTGGGGCGCCCg                                                                                                                                          >  1:109079/1‑149 (MQ=255)
           ccATGGCCCAGAGAAGAAAGCGGAACTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGc                                                                                                                                         >  1:317198/1‑149 (MQ=255)
            cATGGCCCAGAGAAGAAAGCGGAACTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCt                                                                                                                                        >  3:36097/1‑149 (MQ=255)
              tGGCCCATAGAAGAATGCGGAACTTCATAAACTGGCCTCGAATTAAAAGTGAATAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCTAAACATTTTGCCTGCAGAAAGGACAACGGGCGCCCCAAGGCGCCCgttgt                                                                                                                                      <  5:177730/149‑5 (MQ=255)
                gCCCAGAGAAGAAAGCGGAACTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTgg                                                                                                                                     <  2:16457/148‑1 (MQ=255)
                 cccATACAACAAAGCGCAACTTCATCACCGCGCCACTAATCAAACGTTACTACGCCCCGCATTTTATCCTTTACCATAACTATTTATCCAAGTTCTAAACATTTTGCCTGCACAAACCACAACGGGCGCCTTGGGGCGCCCGCTTGGcc                                                                                                                                   <  2:295739/149‑1 (MQ=255)
                 cccAGAGAAGAAAGCGGAACTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGcc                                                                                                                                   <  2:85549/149‑1 (MQ=255)
                 cccAGAGAAGAAAGCGGAACTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTAATGCAAGTTAGAAACATTTTGCCTGCAGAAAGCACAACCGGCGCCTTGGGGCGCCCGCTTGGcc                                                                                                                                   >  3:233853/1‑149 (MQ=255)
                   cAGAGAAGAAAGCGGAACTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGa                                                                                                                                 >  8:167640/1‑149 (MQ=255)
                   cAGAGAAGAAAGCGGAACTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGa                                                                                                                                 >  6:124404/1‑149 (MQ=255)
                   cAGAGAAGAAAGCGGAACTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGa                                                                                                                                 <  5:144280/149‑1 (MQ=255)
                        aagaaAGCGGACCTTCATCACCCCGCCTCGAATGAAAAGTGACTACGGCGCGCAGTTTATCCTTTTCCAGAACTATTTTTGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTg                                                                                                                            <  7:182955/149‑1 (MQ=255)
                          gaaAGCGGAACTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTggg                                                                                                                          >  7:54310/1‑149 (MQ=255)
                          gaaAGCGGAACTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTggg                                                                                                                          <  5:105202/149‑1 (MQ=255)
                           aaaGCGGAACTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCCCAAGGCGCCCgttgt                                                                                                                                      >  5:45626/1‑132 (MQ=255)
                           aaaGCGGAACTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCCCAAGGCGCCCgttgt                                                                                                                                      <  6:45626/136‑5 (MQ=255)
                             aGCGGAACTTCATCACCGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACACCGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGtt                                                                                                                       >  4:246022/1‑149 (MQ=255)
                               cGGAACTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTgggttgt                                                                                                                      >  4:124273/1‑144 (MQ=255)
                               cGGAACTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCATTTTCAAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTg                                                                                                                      >  2:283696/1‑148 (MQ=255)
                                ggAACTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGt                                                                                                                     >  8:150500/1‑148 (MQ=255)
                                ggAACTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGt                                                                                                                     >  5:310762/1‑148 (MQ=255)
                                ggAACTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCCCAAGGCGCCCg                                                                                                                                          >  5:293686/1‑127 (MQ=255)
                                ggAACTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCCCAAGGCGCCCg                                                                                                                                          <  6:293686/127‑1 (MQ=255)
                                 gAACTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGt                                                                                                                     >  1:290960/1‑147 (MQ=255)
                                      tcatcaACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTgggg                                                                                                                         >  1:298704/1‑138 (MQ=255)
                                       catcaACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCagaga                                                                                                               <  6:273806/147‑1 (MQ=255)
                                        atcaACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGc                                                                                                             >  7:240062/1‑148 (MQ=255)
                                        atcaACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGc                                                                                                             >  7:249219/1‑148 (MQ=255)
                                        atcaACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGc                                                                                                             >  2:271574/1‑148 (MQ=255)
                                          caACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGAt                                                                                                          <  1:271574/149‑1 (MQ=255)
                                             cgcgCCGCGAATGAAAAGTGACTAAGCCGCGCAGTTAATCCTTTTCCAGAACTATTTATGCAAGTCCTAAACATTTTGCCGGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCAc                                                                                                       <  4:48290/149‑1 (MQ=255)
                                               cgcCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCATAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACtt                                                                                                     <  1:289693/149‑1 (MQ=255)
                                                   tCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCt                                                                                                 >  7:70177/1‑149 (MQ=255)
                                                   tCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCt                                                                                                 >  5:271334/1‑149 (MQ=255)
                                                     gAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGTTGTCAGAGAGCGATCACTTGCCTTTc                                                                                              <  3:124273/149‑1 (MQ=255)
                                                       aTGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCt                                                                                             <  7:150500/149‑1 (MQ=255)
                                                        tGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCtt                                                                                            <  8:54310/149‑1 (MQ=255)
                                                        tGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCtt                                                                                            <  8:249219/149‑1 (MQ=255)
                                                           aaaGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTtggc                                                                                         <  3:147438/149‑1 (MQ=255)
                                                            aaGTGTCTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTtggct                                                                                        >  6:186359/1‑149 (MQ=255)
                                                               tGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTtggctggc                                                                                     >  1:167895/1‑149 (MQ=255)
                                                                  cTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTg                                                                                  <  3:141409/149‑1 (MQ=255)
                                                                   tAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTgg                                                                                 >  7:172606/1‑149 (MQ=255)
                                                                     aGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCTAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGc                                                                                <  5:124404/148‑1 (MQ=255)
                                                                     aGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGc                                                                                <  6:271334/148‑1 (MQ=255)
                                                                       gcgcgcAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCggg                                                                             >  4:244947/1‑149 (MQ=255)
                                                                        ctcgctttttATCCTTTTCCAGAAGTATTTATGCATGTTCGAAACATTTTGCCTGCAGAATGCACAACGGGGGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCggg                                                                             <  8:257171/141‑1 (MQ=255)
                                                                           gcAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCggg                                                                             >  7:257171/1‑145 (MQ=255)
                                                                            cAGTTTATCCTTTTCAAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGt                                                                        <  1:283696/149‑1 (MQ=255)
                                                                               tttATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGGCAGAGAGCGATCACTTGCCTTTCTtg                                                                                           <  7:253559/127‑1 (MQ=255)
                                                                               tttATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGACAGAGAGCGATCACTTGCCTTTCTtg                                                                                           >  8:253559/1‑127 (MQ=255)
                                                                                 tATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAAc                                                                    <  7:167640/148‑1 (MQ=255)
                                                                                  aTCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGc                                                                  >  1:169773/1‑149 (MQ=255)
                                                                                  aTCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGc                                                                  >  6:168042/1‑149 (MQ=255)
                                                                                   tCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCt                                                                 <  2:169773/149‑1 (MQ=255)
                                                                                      ttttCCAGAACTATTTATGCAAGTTCGAAACATTTTGTCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGt                                                               >  8:103144/1‑148 (MQ=255)
                                                                                      ttttCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTggg                                                                                                                          <  6:112178/89‑1 (MQ=255)
                                                                                      ttttCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTggg                                                                                                                          >  5:112178/1‑89 (MQ=255)
                                                                                           cAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACg                                                           >  1:127169/1‑147 (MQ=255)
                                                                                           cAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGc                                                          >  1:272500/1‑148 (MQ=255)
                                                                                           cAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGc                                                          <  2:272500/148‑1 (MQ=255)
                                                                                            aGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAg                                                        >  7:118975/1‑149 (MQ=255)
                                                                                              aaCTATTTATGCAAGTTCGAAACATATTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAGGt                                                      >  5:24362/1‑149 (MQ=255)
                                                                                                  aTTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTg                                                                                                                                      <  6:108713/65‑1 (MQ=255)
                                                                                                  aTTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTg                                                                                                                                      >  5:108713/1‑65 (MQ=255)
                                                                                                       tGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAGGTTCTTGCCAg                                             <  2:127169/149‑1 (MQ=255)
                                                                                                         cacgTTCGAAACATTTTGCCTGCAGAAACCACACCGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAGGTTCTTGCCAgg                                            <  5:31876/145‑1 (MQ=255)
                                                                                                          aaGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAGGTTCTTGCCAGGCt                                          <  8:240062/149‑1 (MQ=255)
                                                                                                            gTTCTACACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAGGTTCTTGCCACGCTTg                                        <  7:286159/149‑1 (MQ=255)
                                                                                                             ttCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGAt                                                                                                          >  1:200708/1‑82 (MQ=255)
                                                                                                             ttCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGAt                                                                                                          <  2:200708/82‑1 (MQ=255)
                                                                                                              tCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAGGTTCTTGCCAGGCTTGa                                       <  2:167895/148‑1 (MQ=255)
                                                                                                               cGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAGGTTCTTGCCAGGCTTGaa                                      <  8:70177/148‑1 (MQ=255)
                                                                                                               cGAAACATATTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAGGTTCTTGCCAGGCTTGaa                                      <  6:24362/148‑1 (MQ=255)
                                                                                                                gAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAGGTTCtt                                                  >  4:9351/1‑135 (MQ=255)
                                                                                                                gAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAGGTTCtt                                                  <  3:9351/135‑1 (MQ=255)
                                                                                                                gAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAGGTTCTTGCCAGGCTTGAAgg                                    >  5:227291/1‑149 (MQ=255)
                                                                                                                gAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAGGTTCTTGCCAGGCTTGAAg                                     <  5:168042/148‑1 (MQ=255)
                                                                                                                 aaaCATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAGGTTCtt                                                  <  2:16741/134‑1 (MQ=255)
                                                                                                                 aaaCATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAGGTTCtt                                                  >  1:16741/1‑134 (MQ=255)
                                                                                                                     tttttGCCGCCATAAAGCACACGGGGCGCGTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAGGTTCTTGCCAGGCTTGAAGGCCACg                               <  2:290960/148‑1 (MQ=255)
                                                                                                                            cTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCATTAACGCTGTCACGCAGGTTCTTGCCAGGCTTGAAGGCCACGGTATTGc                        >  6:307877/1‑149 (MQ=255)
                                                                                                                                  aaaGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAGGTTCTTGCCAGGCTTGAAGGCCACGGTATTGCTGGCCt                  <  3:246022/149‑1 (MQ=255)
                                                                                                                                        caACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAGGTTCTTGCCAGGCTTGAAGGCCACGGTATTGCTGGCCTTGATCt            <  5:114073/149‑1 (MQ=255)
                                                                                                                                          aCGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAGGTTCTTGCCAGGCTTGAAGGCCACGGTATTGCTGGCCTTGATCTTg          <  7:198387/149‑1 (MQ=255)
                                                                                                                                            gggCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAGGTTCTTGCCAGGCTTGAAGGCCACGGTATTGCTGGCCTTGATCTTGAc        <  2:109079/149‑1 (MQ=255)
                                                                                                                                                 gcCTTGGGGCGCCCGCTTGGCCGACACTGGTGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCTGGTTCTTGCCAGGATTGAAGGCGACGGTATTGCTGGCCTTg                <  5:192843/137‑1 (MQ=255)
                                                                                                                                                 gcCTTGGGGCGCCCGCTTGGCCGACACTCGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTGGGCGGGGGCGTTAACGCTGTCACGCAGGTTCTAGCCAGGCTTGAAGGCAACTGTATTCCTGGCCTTg                >  6:192843/1‑137 (MQ=255)
                                                                                                                                                   cTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAGGTTCTTGCCAGGCTTGAAGGCCACGGTATTGCTGGCCTTGATCTTGACCGGTTc  <  8:172606/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CATCAACAGCCCCATGGCCCAGAGAAGAAAGCGGAACTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAGGTTCTTGCCAGGCTTGAAGGCCACGGTATTGCTGGCCTTGATCTTGACCGGTTC  >  NC_002947/6060878‑6061172

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