Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,477,460 G→C 7.4% G406G (GGC→GGG PP_3087 ← excinuclease ABC subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,477,4600GC 7.4% 278.8 / 10.0 92G406G (GGC→GGGPP_3087excinuclease ABC subunit A
Reads supporting (aligned to +/- strand):  ref base G (64/20);  new base C (5/2);  total (69/23)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

TGTTCGCCGGCAGCGGGCCCTACATCCACCAGCCAGTCGGCGCGGCGCATGGTGTCCAGGTCATGCTCCACCACGTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCGGGGTGCTGCGCTCCAGGGCCAGGTAGCCGAGGCCGAGGTCAAGCAGCGT  >  NC_002947/3477314‑3477609
                                                                                                                                                  |                                                                                                                                                     
tgtTCGCCGGCAGCGGGCCCTACATCCACCAGCCAGTCGGCGCGGCGCATGGTGTCCAGGTCATGCTCCACCACGTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAgccg                                                                                                                                                    <  1:4862687/150‑1 (MQ=255)
 gcccGCCGGCAGCGGGCCCTACATCCACCAGCCAGTCGGCGCGGCGCATGGTGTCCAGGTCATGCTCCACCACGTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAgccgg                                                                                                                                                   <  1:861109/147‑1 (MQ=255)
  ttCGCCGGCAGCGGGCCCTACATCCACCAGCCAGTCGGCGCGGCGCATGGTGTCCAGGTCATGCTCCACCACGTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAgccggc                                                                                                                                                  >  1:4686134/1‑150 (MQ=255)
  ttCGCCGGCAGCGGGCCCTACATCCACCAGCCAGTCGGCGCGGCGCATGGTGTCCAGGTCATGCTCCACCACGTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAgccggc                                                                                                                                                  <  2:3052366/150‑1 (MQ=255)
   tCGCCGGCAGCGGGCCCTACATCCACCAGCCAGTCGGCGCGGCGCATGGTGTCCAGGTCATGCTCCACCACGTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAgccggcc                                                                                                                                                 <  2:1400206/150‑1 (MQ=255)
      ccGGCAGCGGGCCCTACATCCACCAGCCAGTCGGCGCGGCGCATGGTGTCCAGGTCATGCTCCACCACGTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGAc                                                                                                                                              <  1:5432275/150‑1 (MQ=255)
             cGGGCCCTACATCCACCAGCCAGTCGGCGCGGCGCATGGTGTCCAGGTCATGCTCCACCACGTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCtcgt                                                                                                                                       >  2:468035/1‑150 (MQ=255)
                gCCCTACATCCACCAGCCAGTCGGCGCGGCGCATGGTGTCCAGGTCATGCTCCACCACGTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGa                                                                                                                                    >  1:2585925/1‑150 (MQ=255)
                 cccTACATCCACCAGCCAGTCGGCGCGGCGCATGGTGTCCAGGTCATGCTCCACCACGTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAg                                                                                                                                   >  1:2520328/1‑150 (MQ=255)
                  ccTACATCCACCAGCCAGTCGGCGCGGCGCATGGTGTCCAGGTCATGCTCCACCACGTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGt                                                                                                                                  >  2:724111/1‑150 (MQ=255)
                  ccTACATCCACCAGCCAGTCGGCGCGGCGCATGGTGTCCAGGTCATGCTCCACCACGTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGt                                                                                                                                  >  2:3457366/1‑150 (MQ=255)
                   cTACATCCACCAGCCAGTCGGCGCGGCGCATGGTGTCCAGGTCATGCTCCACCAGGTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCCGGTGCAAGCCGGCCGACGGCTCGTCGAGTa                                                                                                                                 >  2:1462020/1‑150 (MQ=255)
                      cATCCACCAGCCAGTCGGCGCGGCGCATGGTGTCCAGGTCATGCTCCACCACGTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGt                                                                                                                              >  2:226850/1‑150 (MQ=255)
                       aTCCACCAGCCAGTCGGCGCGGCGCATGGTGTCCAGGTCATGCTCCACCACGTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTa                                                                                                                             >  2:4674750/1‑150 (MQ=255)
                        tCCACCAGCCAGTCGGCGCGGCGCATGGTGTCCAGGTCATGCTCCACCACGTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTaa                                                                                                                            >  2:3122483/1‑150 (MQ=255)
                          caccaGCCAGTCGGCGCGGCGCATGGTGTCCAGGTCATGCTCCACCACGTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAac                                                                                                                          >  1:3059086/1‑150 (MQ=255)
                          caccaGCCAGTCGGCGCGGCGCATGGTGTCCAGGTCATGCTCCACCACGTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAac                                                                                                                          >  2:3469560/1‑150 (MQ=255)
                           accaGCCAGTCGGCGCGGCGCATGGTGTCCAGGTCATGCTCCACCACGTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAgccggcc                                                                                                                                                 <  2:413509/126‑1 (MQ=255)
                           accaGCCAGTCGGCGCGGCGCATGGTGTCCAGGTCATGCTCCACCACGTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAgccggcc                                                                                                                                                 >  1:413509/1‑126 (MQ=255)
                           accaGCCAGTCGGCGCGGCGCATGGTGTCCAGGTCATGCTCCACCACGTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAacc                                                                                                                         >  1:3063342/1‑150 (MQ=255)
                           accaGCCAGTCGGCGCGGCGCATGGTGTCCAGGTCATGCTCCACCACGTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAacc                                                                                                                         >  1:4978361/1‑150 (MQ=255)
                           accaGCCAGTCGGCGCGGCGCATGGTGTCCAGGTCATGCTCCACCACGTACACCGAGTTGCCAGCCTGTTGCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAacc                                                                                                                         >  2:3559449/1‑150 (MQ=255)
                                      ggcgcggcgcATGGTGTCCAGGTCATGCTCCACCACGTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAg                                                                                                              >  2:424997/1‑150 (MQ=255)
                                       gcgcggcgcATGGTGTCCAGGTCATGCTCCACCACGTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGc                                                                                                             >  1:437657/1‑150 (MQ=255)
                                          cggcgcTGGGTGTCCTGGTCTTGCACCACCACGTACACCGCGTTGCCTGCTTTTTTCTGGGGGTGCTGGGCTTCTAACGGCGCTTCACTGTCGGCGGGGTGAACGCGGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGc                                                                                                          <  2:4187388/150‑1 (MQ=255)
                                                  ggTGTCCAGGTCATGCTCCACCACGTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAg                                                                                                  >  1:2804535/1‑150 (MQ=255)
                                                  ggTGTCCAGGTCATGCTCCACCACGTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAg                                                                                                  >  2:2556162/1‑150 (MQ=255)
                                                    tttCCTGGTCATGCACCACCACTTACACCGTGTTGCCAGCCTTTTTCTGGCGCTGCTGGGCTTCTAACAGCTCTACACTGTCGGCGGGGTGAACGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCt                                                                                                <  2:1210538/148‑1 (MQ=255)
                                                    tgtCCAGGTCATGCTCCACCACGTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACTGCGCTTCACTGTCAGCCGGGTGCAATCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCt                                                                                                <  2:5125100/150‑1 (MQ=255)
                                                    tgtCCAGGTCATGCTCCACCACGTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCCAACAGCTGCGAATTCAGcg                                                                                                >  2:1083418/1‑149 (MQ=255)
                                                    tgtCCAGGTCATGCTCCACCACGTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCt                                                                                                >  2:3744824/1‑150 (MQ=255)
                                                    ggtacAGGTCTTGCTCCACCACGTACACCGTGTGGCCAGCCTTTTTCTGGCGGTGCTGGGTTTGTAACTGGTCTTCTCTGTGGGCGGGGTGCAAGCGGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCt                                                                                                <  2:2208083/146‑1 (MQ=255)
                                                        cAGGTCATGCTCCACCACGTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGt                                                                                            >  2:769833/1‑150 (MQ=255)
                                                         aGGTCATGCTCCACCACGTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTg                                                                                           >  2:4732367/1‑150 (MQ=255)
                                                               tccTCCCCCACTTACACCGTGTTGCCAGCCTGTTTCGGGCGCTGCGGGGCTACGAACCGGGTTTCTCTGTCGGCGGGGTGCACGCCGCCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCcagg                                                                                     <  1:2255757/148‑1 (MQ=255)
                                                               tGCTCCACCACGTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCcagg                                                                                     >  1:460104/1‑150 (MQ=255)
                                                                gCTCCACCCCGTACACGGTTTTGCCAGCTTTTTTCCGGCGCTGCTGGGATTCTAACAGCGCTTCACTGTCGGGCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCcaggc                                                                                    <  2:2650322/150‑1 (MQ=255)
                                                                gCTCCACCACGTACCCGGTGTGGCCTGCGTGTTTCAGGCGCTGCCGGGCTTCGAACCGCGCTTCACTGTGGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCcaggc                                                                                    <  2:5375476/150‑1 (MQ=255)
                                                                  tCCACCACGTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCcaggcgc                                                                                  >  1:1066689/1‑150 (MQ=255)
                                                                   ccaccaCGTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCcaggcgca                                                                                 >  1:308406/1‑150 (MQ=255)
                                                                    caccacGTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCACGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCcaggcgcag                                                                                >  1:947298/1‑150 (MQ=255)
                                                                    caccacGTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCcaggcgcag                                                                                >  1:554812/1‑150 (MQ=255)
                                                                    caccacGTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCcaggcgcag                                                                                >  2:4644497/1‑150 (MQ=255)
                                                                    caccacGTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCcaggcgcag                                                                                >  1:3961785/1‑150 (MQ=255)
                                                                    caccacGTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCcaggcgcag                                                                                >  2:4644389/1‑150 (MQ=255)
                                                                     accacGTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCcaggcgcagg                                                                               >  1:5011900/1‑150 (MQ=255)
                                                                       cacGTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCcaggcgcaggcg                                                                             >  2:1017392/1‑150 (MQ=255)
                                                                          gtccacCGCGTGCCCAGCCTGTTCCCGGCGCGGCCGGGCTCCGACCCGCGCTTCACGGGCGGCCGGGTCCCGCCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTg                                                                          <  2:1637395/147‑1 (MQ=255)
                                                                          gTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTg                                                                          >  1:3636413/1‑150 (MQ=255)
                                                                              acCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGt                                                                      >  2:4426741/1‑150 (MQ=255)
                                                                                    ttGCCAGCCTTTTTCCGGCGCTGCTGGGCTCCCACCTGCGCTACACTGTGGGCGGGGGGCACGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCa                                                                <  1:2710432/150‑1 (MQ=255)
                                                                                    ttGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCa                                                                >  2:915129/1‑150 (MQ=255)
                                                                                    gggCCAGCTTTTTTCTGGGGCTGCTGGGCTTCTACCCGCGCTTCACGGGCGGCGGGGTGCCCGCGGCCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCa                                                                <  2:2337453/148‑1 (MQ=255)
                                                                                      gCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCaga                                                              >  2:533219/1‑150 (MQ=255)
                                                                                       ccAGCGTTTTTCCGCCGCTGCTGGTCTACGAACAGGGCTTCCCTGTCGGCGGGGTGCACGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCagaa                                                             <  2:2102900/150‑1 (MQ=255)
                                                                                         aGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCagaaga                                                           >  2:2050125/1‑150 (MQ=255)
                                                                                             tGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAg                                                       >  1:1227825/1‑150 (MQ=255)
                                                                                               tttCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGt                                                     >  1:4645936/1‑150 (MQ=255)
                                                                                                   tggCGCTGGTGGGTTTCTAACTGGGCTTCAGTGGGGGCGGGTTGCATGCGGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCggg                                                 <  2:4610345/149‑1 (MQ=255)
                                                                                                    ggCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGTCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCgggg                                                >  1:1792402/1‑150 (MQ=255)
                                                                                                       gctgctGGGCTTCTAACTGGGGTTCACGGTCGGCGGGGTGCACGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCGGGGtgc                                             <  2:1730757/150‑1 (MQ=255)
                                                                                                         tGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCGGGGtgctg                                           >  2:2079003/1‑150 (MQ=255)
                                                                                                            aGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCGGGGTGCTgcgc                                        >  1:242420/1‑150 (MQ=255)
                                                                                                              ggCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCGGGGTGCTGCGCTc                                      >  1:1990412/1‑150 (MQ=255)
                                                                                                              ggCTACGAACTGGGCTACCCGGGCGGCCGGTTCAACGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCGGGGTGCTGCGCTc                                      <  1:4987875/150‑1 (MQ=255)
                                                                                                               gCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCGGGGTGCTGCGCTcc                                     >  1:620052/1‑150 (MQ=255)
                                                                                                                     aaCAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCGGGGTGCTGCGCTCCAGGGcc                               >  1:477259/1‑150 (MQ=255)
                                                                                                                     aaCAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCGGGGTGCTGCGCTCCAGGGcc                               >  1:477245/1‑150 (MQ=255)
                                                                                                                     aaCAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCGGGGTGCTGCGCTCCAGGGcc                               >  1:477248/1‑150 (MQ=255)
                                                                                                                        aGCGCTTCCCTGTCGGCGGGTTTCAAGCGGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGAGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCGGGGTGCTGCGCTCCAGGGCCAgg                            <  2:419720/150‑1 (MQ=255)
                                                                                                                        aGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCGGGGTGCTGCGCTCCAGGGCCAgg                            >  1:273701/1‑150 (MQ=255)
                                                                                                                        aGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCGGGGTGCTGCGCTCCAGGGCCAgg                            >  1:2958279/1‑150 (MQ=255)
                                                                                                                        aGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCGGGGTGCTGCGCTCCAGGGCCAgg                            >  2:4563039/1‑150 (MQ=255)
                                                                                                                         gcgcTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCGGGGTGCTGCGCTCCAGGGCCAGGt                           >  2:4038294/1‑150 (MQ=255)
                                                                                                                            cTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCGGGGTGCTGCGCTCCAGGGCCAGGTCgc                        >  1:4508179/1‑150 (MQ=255)
                                                                                                                             ttCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCGGGGTGCTGCGCTCCAGGGCCAGGTAgcc                       >  1:5253942/1‑150 (MQ=255)
                                                                                                                               ctcTGTCGGCGGGGTGCATGCGGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCGGGGTGCTGCGCTCCAGGGCCAGGTAgccga                     <  2:3636413/148‑1 (MQ=255)
                                                                                                                               ccgggcgggCGGGGTGCACGCCGCCCGACGGCTCGTCGAGTACGTAAACCACCCCAACCAGCTGCGCATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCGGGGTGCTGCGCTCCAGGGCCAGGTAgccga                     <  2:437657/143‑1 (MQ=255)
                                                                                                                                aCTGTCGGCCGGGTGCAAGCGGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCGGGGTGCTGCGCTCCAGGGCCAGGTAgccgag                    <  2:4844125/150‑1 (MQ=255)
                                                                                                                                 ccgTCGGCCGGCTTGCACCCGGCCGATAGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCGGGGtgctg                                           >  2:4344528/3‑126 (MQ=255)
                                                                                                                                 ccgTCGGCCGGCTTGCACCCGGCCGATAGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCGGGGtgctg                                           <  1:4344528/124‑1 (MQ=255)
                                                                                                                                   gTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCGGGGTGCTGCGCTCCAGGGCCAGGTAgccgaggcc                 <  2:4978361/150‑1 (MQ=255)
                                                                                                                                   gTCGGCCGGCTTGCACCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCGGGGTGCTGCGCTCCAGGGCCAGGTAgccgaggcc                 >  2:3722415/1‑150 (MQ=255)
                                                                                                                                   gTCGGCCGGCTTGCACCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCGGGGTGCTGCGCTCCAGGGCCAGGTAgccgaggcc                 >  2:4314557/1‑150 (MQ=255)
                                                                                                                                   gTCGGCCGGCTTGCACCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCGGGGTGCTGCGCTCCAGGGCCAGGTAgccgaggcc                 >  2:4314685/1‑150 (MQ=255)
                                                                                                                                   gTCGGCCGGCTTGCACCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCGGGGTGCTGCGCTCCAGGGCCAGGTAgccgaggcc                 >  2:4314822/1‑150 (MQ=255)
                                                                                                                                    tCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCGGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCGGGGTGCTGCGCTCCAGGGCCAGCTAgccgaggccg                >  2:797668/1‑150 (MQ=255)
                                                                                                                                     cggcggGTTGAAAGCCGGCCGACGTCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCGGGGTGCTGCGCTCCAGGGCCAGGTAgccgaggccga               <  2:1792402/150‑1 (MQ=255)
                                                                                                                                      ggCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCGGGGTGCTGCGCTCCAGGGCCAGGTAgccgaggccgag              >  1:1647756/1‑150 (MQ=255)
                                                                                                                                       gCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCGGGGTGCTGCGCTCCAGGGCCAGGTAGccgaggccgagg             >  2:4439648/1‑150 (MQ=255)
                                                                                                                                       gCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCGGGGTGCTGCGCTCCAGGGCCAGGTAGccgaggccgagg             >  2:3611722/1‑150 (MQ=255)
                                                                                                                                            gTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCGGGGTGCTGCGCTCCAGGGCCAGGTAGCCGAGGCCGAGGTCAag        <  2:4878866/150‑1 (MQ=255)
                                                                                                                                            gTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGCTGCAGTTCACCAGAAGACAAGGTCGGGGTGCTGCGCTCCAGGGCCAGGTAGCCGAGGCCGAGGTCAag        >  1:949026/1‑150 (MQ=255)
                                                                                                                                             tGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCGGGGTGCTGCGCTCCAGGGCCAGGTAGCCGAGGCCGAGGTCAagc       >  1:516684/1‑150 (MQ=255)
                                                                                                                                             tGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCGGGGTGCTGCGCTCCAGGGCCAGGTAGCCGAGGCCGAGGTCAagc       >  1:516688/1‑150 (MQ=255)
                                                                                                                                              gCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCGGGGTGCTGCGCTCCAGGGCCAGGTAGCCGAGGCCGAGGTCAagca      >  1:4854105/1‑150 (MQ=255)
                                                                                                                                              gCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCGGGGTGCTGCGCTCCAGGGCCAGGTAGCCGAGGCCGAGGTCAagca      >  1:162043/1‑150 (MQ=255)
                                                                                                                                              gCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCGGGGTGCTGCGCTCCAGGGCCAGGTAGCCGAGGCCGAGGTCAagca      >  2:1548976/1‑150 (MQ=255)
                                                                                                                                                aaGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCGGGGTGCTGCGCTCCAGGGCCAGGTAGCCGAGGCCGAGGTCAagcagc    >  2:2262214/1‑150 (MQ=255)
                                                                                                                                                  gccggccgACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCGGGGTGCTGCGCTCCAGGGCCAGGTAGCCGAGGCCGAGGTCAAGCAGCGt  >  1:2879208/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TGTTCGCCGGCAGCGGGCCCTACATCCACCAGCCAGTCGGCGCGGCGCATGGTGTCCAGGTCATGCTCCACCACGTACACCGAGTTGCCAGCCTGTTTCAGGCGCTGCAGGGCTTCGAACAGCGCTTCACTGTCGGCCGGGTGCAAGCCGGCCGACGGCTCGTCGAGTACGTAAACCACCCCAAACAGCTGCGAATTCAGCTGGGTGGCCAGGCGCAGGCGTTGCAGTTCACCAGAAGACAAGGTCGGGGTGCTGCGCTCCAGGGCCAGGTAGCCGAGGCCGAGGTCAAGCAGCGT  >  NC_002947/3477314‑3477609

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: