Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,748,244 G→C 10.2% T245T (ACC→ACG etfA ← electron transfer flavoprotein subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,748,2440GC10.2% 1027.6 / 83.5 322T245T (ACC→ACGetfAelectron transfer flavoprotein subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (160/129);  new base C (15/18);  total (175/147)
Fisher's exact test for biased strand distribution p-value = 3.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 331 total aligned reads displayed.

ACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCACCGTTACCCATGCCACGGCCGCCGGAAACGACGA  >  NC_002947/4748095‑4748392
                                                                                                                                                     |                                                                                                                                                    
acgaGGCCATAATCAGCCACCTGGAAGATCGGTGCTTCTTCGTCCTTGTTGATCGCAACAATCACGCGGCTGTCTTTCATGCCCGCCAGATGCTGGATCGCACCCGAGATGCCGACTGCGATGTAAAGTTCAGGTGCAACCACCTTGCCg                                                                                                                                                      >  2:5391304/4‑150 (MQ=255)
aCCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCg                                                                                                                                                      <  1:1035785/150‑1 (MQ=255)
aCCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCg                                                                                                                                                      <  2:7314180/150‑1 (MQ=255)
aCCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGAGCTTGCCg                                                                                                                                                      >  1:4089718/1‑150 (MQ=255)
  gaGGCCATAATCAGCCACCTGGAAGATCGGTGCTTCTTCGTCCTTGTTGATCGCAACAATCACGCGGCTGTCTTTCATGCCCGCCAGATGCTGGATCGCACCCGAGATGCCGACTGCGATGTAAAGTTCAGGTGCAACCACCTTGCCGGt                                                                                                                                                    <  2:3871653/149‑1 (MQ=255)
  cAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGt                                                                                                                                                    <  2:8351380/150‑1 (MQ=255)
  cAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGt                                                                                                                                                    >  1:1566970/1‑150 (MQ=255)
     gCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTg                                                                                                                                                 >  1:5022136/1‑150 (MQ=255)
     gCCATAATCGGCAACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATGGCGACGATGACGCGTGAGTCCTTCATGCCAGCCAGATGCTGGATCGCGCCGGAAATACCGACAGCGATGTAAAGCTCCGGTGCAACCACCTTGCCCGTCTg                                                                                                                                                 >  2:6121635/1‑150 (MQ=255)
      ccataatCGGCAACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATGGCGACGATGACGCGTGAGTCCTTCATGCCAGCCAGATGCTGGATCGCGCCGGAAATACCGACAGCGATGTAAAGCTCCGGTGCAACCACCTTGCCCGTCTgg                                                                                                                                                <  1:8443915/144‑1 (MQ=255)
      ccGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTgg                                                                                                                                                <  1:6933461/150‑1 (MQ=255)
      ccGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTgg                                                                                                                                                <  2:476579/150‑1 (MQ=255)
         taatCGGCAACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATGGCGACGATGACGCGTGAGTCCTTCATGCCAGCCAGATGCTGGATCGCGCCGGAAATACCGACAGCGATGTAAAGCTCCGGTGCAACCACCTTGCCCGTCTGGCCg                                                                                                                                             >  1:6417996/4‑150 (MQ=255)
          aGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGa                                                                                                                                            >  1:5067761/1‑150 (MQ=255)
            tCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGAcc                                                                                                                                          >  2:6687003/1‑150 (MQ=255)
            tCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGAcc                                                                                                                                          >  2:4718839/1‑150 (MQ=255)
             cagCCACCTGGAAGATCGGTGCTTCTTCGTCCTTGTTGATCGCAACAACCACGCGGCTGTCTTTCATGCCCGCCAGATGCTGGATCGCACCCGAGATGCCGACTGCGATGTAAAGTTCAGGTGCAACCACCTTGCCGGTCTGGCCAACCt                                                                                                                                         >  2:1243467/3‑150 (MQ=255)
             cGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCt                                                                                                                                         <  1:4889152/150‑1 (MQ=255)
               gcaaCCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATGGCGACGATGACGCGTGAGTCCTTCATGCCAGCCAGATGCTGGATCGCGCCGGAAATACCGACAGCGATGTAAAGCTCCGGTGCAACCACCTTGCCCGTCTGGCCGACCTGc                                                                                                                                       <  2:4830818/147‑1 (MQ=255)
                 cACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCAt                                                                                                                                     <  2:7762363/150‑1 (MQ=255)
                 cACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCAt                                                                                                                                     <  1:1265332/150‑1 (MQ=255)
                  aCCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATg                                                                                                                                    >  2:8176011/1‑150 (MQ=255)
                   ccTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGt                                                                                                                                   <  2:4949114/150‑1 (MQ=255)
                       gAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGtt                                                                                                                               >  2:8918206/1‑150 (MQ=255)
                        aaGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTc                                                                                                                              <  2:1566970/150‑1 (MQ=255)
                        aaGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTc                                                                                                                              <  1:13820/150‑1 (MQ=255)
                        aaGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTc                                                                                                                              <  1:8986134/150‑1 (MQ=255)
                           aTCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGc                                                                                                                           >  1:4562974/1‑150 (MQ=255)
                            tCGGTGCTTCTTCGTCCTTGTTGATCGCAACAATCACGCGGCTGTCTTTCATGCCCGCCAGATGCTGGATCGCACCCGAGATGCCGACTGCGATGTAAAGTTCAGGTGCAACCACCTTGCCGGTCTGGCCAACCTGCCAGTCGTTTGGcg                                                                                                                          >  1:5741612/1‑149 (MQ=255)
                             cGGCGCTTCTTCGTCCTTGTTGATGGCGACGATGACGCGTGAGTCCTTCATGCCAGCCAGATGCTGGATCGCGCCGGAAATACCGACAGCGATGTAAAGCTCCGGTGCAACCACCTTGCCCGTCTGGCCGACCTGCCAGTCGTTCGGcgc                                                                                                                         >  1:8715553/1‑148 (MQ=255)
                              ggCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACg                                                                                                                        >  1:942712/1‑150 (MQ=255)
                               gtgcTTCTTCGTCCTTGTTGATCGCAACAATCACGCGGCTGTCTTTCATGCCCGCCAGATGCTGGATCGCACCCGAGATGCCGACTGCGATGTAAAGTTCAGGTGCAACCACCTTGCCGGTCTGGCCAACCTGCCAGTCGTTTGGcgcat                                                                                                                       <  1:7685617/148‑5 (MQ=255)
                               gcgcTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGa                                                                                                                       >  2:1331567/1‑150 (MQ=255)
                                cgcTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGaa                                                                                                                      <  2:8607773/150‑1 (MQ=255)
                                 gcTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAg                                                                                                                     >  1:3887468/1‑150 (MQ=255)
                                  cttcttCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGc                                                                                                                    <  1:5866893/150‑1 (MQ=255)
                                         gTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTc                                                                                                             >  2:257906/1‑150 (MQ=255)
                                         gTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTc                                                                                                             >  2:5055723/1‑150 (MQ=255)
                                           ccTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGa                                                                                                           <  2:5913057/150‑1 (MQ=255)
                                             ttgttgATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACt                                                                                                         <  1:3412413/150‑1 (MQ=255)
                                              tgttgATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTg                                                                                                        <  1:1331567/150‑1 (MQ=255)
                                                ttgATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTgcg                                                                                                      >  1:497762/1‑150 (MQ=255)
                                                ttgATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTgcg                                                                                                      <  2:3243414/150‑1 (MQ=255)
                                                gtgATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTgcg                                                                                                      <  2:5343071/149‑1 (MQ=255)
                                                 tgATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTTCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTgcgg                                                                                                     >  2:6554568/1‑150 (MQ=255)
                                                    tCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGgcgc                                                                                                  >  1:4506188/1‑150 (MQ=255)
                                                     cGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGgcgcg                                                                                                 >  2:7061027/1‑150 (MQ=255)
                                                      gCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGgcgcgc                                                                                                >  2:241092/1‑150 (MQ=255)
                                                        cACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGa                                                                                              >  2:8558609/1‑150 (MQ=255)
                                                         aCGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGaa                                                                                             <  1:532705/150‑1 (MQ=255)
                                                          cGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAg                                                                                            <  1:8025051/150‑1 (MQ=255)
                                                           gATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGc                                                                                           >  1:1017471/1‑150 (MQ=255)
                                                            aTCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCa                                                                                          <  2:3406168/150‑1 (MQ=255)
                                                            aTCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCa                                                                                          <  2:3073728/150‑1 (MQ=255)
                                                                 tttGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAAc                                                                                     >  1:7363378/1‑150 (MQ=255)
                                                                 tttGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAAc                                                                                     >  1:5153705/1‑150 (MQ=255)
                                                                 tttGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAAc                                                                                     <  1:7464041/150‑1 (MQ=255)
                                                                 tttGGAGTCTCTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAAc                                                                                     >  1:6214165/1‑150 (MQ=255)
                                                                    ggAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGc                                                                                  >  1:4677803/1‑150 (MQ=255)
                                                                    ggAGTCTTTCATGCCGGCCAGGTGCTGGATCACGCCGGAGATACTGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGc                                                                                  <  2:7538463/150‑1 (MQ=255)
                                                                     gAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCa                                                                                 <  1:6011943/150‑1 (MQ=255)
                                                                     gAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCa                                                                                 >  2:4281217/1‑150 (MQ=255)
                                                                          tttCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGcc                                                                            >  2:5809142/1‑150 (MQ=255)
                                                                          tttCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGcc                                                                            >  2:1100785/1‑150 (MQ=255)
                                                                            tCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGa                                                                          <  1:4286212/150‑1 (MQ=255)
                                                                             cATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAg                                                                         <  1:2457935/150‑1 (MQ=255)
                                                                               tGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCt                                                                       >  2:1212300/1‑150 (MQ=255)
                                                                                gCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCtt                                                                      >  1:4176278/1‑150 (MQ=255)
                                                                                gCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCtt                                                                      <  2:2260391/150‑1 (MQ=255)
                                                                                 ccGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTg                                                                     <  2:7412498/150‑1 (MQ=255)
                                                                                   ggCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTc                                                                   <  1:2781860/150‑1 (MQ=255)
                                                                                      cAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGc                                                                <  2:334681/150‑1 (MQ=255)
                                                                                         gTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAg                                                             <  1:3086624/150‑1 (MQ=255)
                                                                                           gCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGc                                                           >  1:7881262/1‑150 (MQ=255)
                                                                                                aTCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTAc                                                      >  1:737341/1‑150 (MQ=255)
                                                                                                    cgcCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGt                                                  <  1:2178907/150‑1 (MQ=255)
                                                                                                     gcCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTg                                                 >  1:2821189/1‑150 (MQ=255)
                                                                                                     gcCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTg                                                 <  1:6698644/150‑1 (MQ=255)
                                                                                                     gcCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTg                                                 >  1:8593534/1‑150 (MQ=255)
                                                                                                           gaTACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGaa                                           <  2:8456829/150‑1 (MQ=255)
                                                                                                           gaTACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGaa                                           <  2:4438233/150‑1 (MQ=255)
                                                                                                                cGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGt                                      >  1:2694472/1‑150 (MQ=255)
                                                                                                                   cGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCAc                                   <  2:379820/150‑1 (MQ=255)
                                                                                                                     gCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCACCg                                 <  2:2783489/150‑1 (MQ=255)
                                                                                                                      cGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCACCGt                                >  2:5627884/1‑150 (MQ=255)
                                                                                                                      cGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCACCGt                                >  1:7171336/1‑150 (MQ=255)
                                                                                                                           tACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCACCGTTAccc                           >  2:5852125/1‑150 (MQ=255)
                                                                                                                             cAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCACCGTTACCCAt                         <  2:7088311/150‑1 (MQ=255)
                                                                                                                             cAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCACCGTTACCCAt                         <  1:3521957/150‑1 (MQ=255)
                                                                                                                               gCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCACCGTTACCCATGc                       <  2:587701/150‑1 (MQ=255)
                                                                                                                                cTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGTAGTTGTCACCGTTACCCATGcc                      <  2:8801529/150‑1 (MQ=255)
                                                                                                                                 tgtgGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCgg                                                                 >  1:1492099/1‑106 (MQ=255)
                                                                                                                                       gcAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCACCGTTACCCATGCCACGgccg               >  2:6253202/1‑150 (MQ=255)
                                                                                                                                       gcAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCACCGTTACCCATGCCACGgccg               >  2:279792/1‑150 (MQ=255)
                                                                                                                                         aaCGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTAc                                                      <  1:7108651/109‑1 (MQ=255)
                                                                                                                                              tCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCACCGTTACCCATGCCACGGCCGCCGGaaa        <  1:9087239/150‑1 (MQ=255)
                                                                                                                                                ttGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCACCGTTACCCATGCCACGGCCGCCGGAAacg      >  2:2072897/1‑150 (MQ=255)
                                                                                                                                                ttGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCACCGTTACCCATGCCACGGCCGCCGGAAacg      >  1:3633543/1‑150 (MQ=255)
                                                                                                                                                atGCAGGTCGGCCAGACCGGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCACCGTTACCCATGCCACGGCCGCCGGAAacg      >  2:7844979/2‑150 (MQ=255)
                                                                                                                                                atGCAGGTCGGCCAGACCGGCATGTCGATCCGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGAACAGGTGCTAGAAGTTGTCACCGTTACCCATGCCACGGCCGCCGGAAacg      >  1:1900531/2‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ACCAGGCCGTAGTCGGCCACCTGGAAGATCGGCGCTTCTTCGTCCTTGTTGATCGCCACGATCACTTTGGAGTCTTTCATGCCGGCCAGGTGCTGGATCGCGCCGGAGATACCGACGGCGATGTACAGCTGTGGCGCAACGATCTTGCCGGTCTGGCCGACCTGCATGTCGTTCGGCACGAAGCCTGCGTCGACTGCGGCGCGCGAAGCACCAACCGCAGCGCCGAGCTTGTCGGCCAGGCTGTACAGGTGCTTGAAGTTGTCACCGTTACCCATGCCACGGCCGCCGGAAACGACGA  >  NC_002947/4748095‑4748392

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