Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,753,062 A→G 8.1% Q11R (CAG→CGG)  PP_1555 → putative Phage integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,753,0620AG 8.1% 212.1 / 6.9 74Q11R (CAG→CGG) PP_1555putative Phage integrase
Reads supporting (aligned to +/- strand):  ref base A (37/31);  new base G (3/3);  total (40/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.41e-01

AACGACCAAGGCTACCTGGCCGTTGTGTGTCAGGGGCATTTCGACGCCATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCGGCTCATACCCTGAGCTGAGCGCTGCC  >  NC_002947/1752915‑1753210
                                                                                                                                                   |                                                                                                                                                    
aaCGACCAAGGCTACCTGGCCGTTGTGTGTCAGGGGCATTTCGACGCCATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCa                                                                                                                                                      <  6:67718/148‑1 (MQ=255)
       aaGGCTACCTGGCCGTTGTGTGTCAGGGGCATTTCGACGCCATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGccgccg                                                                                                                                              <  2:236748/149‑1 (MQ=255)
       aaGGCTACCTGGCCGTTGTGTGTCAGGGGCATTTCGACGCCATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGAAGAGATTCGCCGGCAGGccgccg                                                                                                                                              >  6:331273/1‑149 (MQ=255)
        aGGCTACCTGGCCGTTGTGTGTCAGGGGCATTTCGACGCCATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCgg                                                                                                                                             <  7:276858/149‑1 (MQ=255)
        aGGCTACCTGGCCGTTGTGTGTCAGGGGCATTTCGACGCCATGGCGTGCCTGAGCGCGTACCTGGCCCTGCCAAAAACGGAGGTTCCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCCGCAGgccgcccg                                                                                                                                             >  7:162365/1‑147 (MQ=255)
              ccTGGCCGTTGTGTGTCAGGGGCATTTCGACGCCATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCGGCCTGCCGGCGaa                                                                                                                                         <  4:83482/147‑1 (MQ=255)
              ccTGGCCGTTGTGTGTCAGGGGCATTTCGACGCCATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCGGCCTGCCGGCGaa                                                                                                                                         <  2:52452/147‑1 (MQ=255)
                tGGCCGTTGTGTGTCAGGGGCATTTCGACGCCATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGc                                                                                                                                     >  5:332102/1‑149 (MQ=255)
                  gCCGTTGTGTGTCAGGGGCATTTCGACGCCATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCg                                                                                                                                    <  7:198387/148‑1 (MQ=255)
                  gCCGTTGTGTGTCAGGGGCATTTCGACGCCATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGt                                                                                                                                   <  5:191414/149‑1 (MQ=255)
                    cGTTGTGTGTCAGGGGCATTTCGACGCCATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGa                                                                                                                                 >  4:144728/1‑149 (MQ=255)
                    cGTTGTGTGTCAGGGGCATTTCGACGCCATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGa                                                                                                                                 <  3:300229/149‑1 (MQ=255)
                             tCAGGGGCATTTCGACGCCATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCgcgc                                                                                                                        >  6:342012/1‑149 (MQ=255)
                               aGGGGCATTTCGACGCCATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGccc                                                                                                                      >  5:145039/1‑149 (MQ=255)
                                ggggCATTTCGACGCCATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAg                                                                                                                                        >  8:8784/1‑130 (MQ=255)
                                ggggCATTTCGACGCCATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAg                                                                                                                                        <  7:8784/130‑1 (MQ=255)
                                ggggCATTTCGACGCCATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCt                                                                                                                     >  6:205780/1‑149 (MQ=255)
                                 gggCATTTCGACGCCATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTg                                                                                                                    <  4:25599/149‑1 (MQ=255)
                                  ggCATTTCGACGCCATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTgg                                                                                                                   <  6:332102/149‑1 (MQ=255)
                                    cATTTCGACGCCATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTggg                                                                                                                  >  5:368194/1‑148 (MQ=255)
                                      tttCGACGCCATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCaa                                                                                                               >  5:172908/1‑149 (MQ=255)
                                             gCCATGGCGTGCGTGAGGGCGTACCTGGCCCTGCCAAAATCGGATGTTGCAGCATGACCTACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGc                                                                                                        <  8:162365/149‑1 (MQ=255)
                                             gCCATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGCCCTTCGCGCCCTGGGCAACCACGGc                                                                                                        >  2:279436/1‑149 (MQ=255)
                                             gCCATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGc                                                                                                        >  7:329050/1‑149 (MQ=255)
                                             gCCATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGc                                                                                                        >  8:317736/1‑149 (MQ=255)
                                               cATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCt                                                                                                      >  8:157489/1‑149 (MQ=255)
                                                 tGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGt                                                                                                    <  7:157489/149‑1 (MQ=255)
                                                       gCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTcc                                                                                              >  5:311291/1‑149 (MQ=255)
                                                          tGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCt                                                                                           >  3:210694/1‑149 (MQ=255)
                                                           gAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGCCCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCtt                                                                                          >  4:18484/1‑149 (MQ=255)
                                                                    ccTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTc                                                                                 >  1:283721/1‑149 (MQ=255)
                                                                     cTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCTACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTcc                                                                                >  1:243760/1‑149 (MQ=255)
                                                                               cAAAAACGGAGGTTGCAGAATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCGTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTGCATCGCTCCAGcgcgcgcg                                                                      >  8:189057/1‑149 (MQ=255)
                                                                                aaaaaCGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCgg                                                                     >  4:46552/1‑149 (MQ=255)
                                                                                 aaaaCGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCggg                                                                    >  6:197354/1‑149 (MQ=255)
                                                                                 aaaaCGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCggg                                                                    <  5:197354/149‑1 (MQ=255)
                                                                                  aaaCGGAGGCTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGt                                                                   <  8:326946/149‑1 (MQ=255)
                                                                                     cGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGAGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGt                                                                >  7:149911/1‑149 (MQ=255)
                                                                                      ggAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTg                                                               <  6:311291/149‑1 (MQ=255)
                                                                                      ggAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTg                                                               <  3:144728/149‑1 (MQ=255)
                                                                                       gAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTgg                                                              <  5:205780/149‑1 (MQ=255)
                                                                                        aGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGt                                                             <  2:243760/149‑1 (MQ=255)
                                                                                         ggTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTa                                                            >  5:70801/1‑149 (MQ=255)
                                                                                             gcaggaTGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCt                                                         <  7:317736/148‑1 (MQ=255)
                                                                                                 caTGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATcc                                                    <  8:329050/149‑1 (MQ=255)
                                                                                                       cAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAgg                                              <  5:342012/149‑1 (MQ=255)
                                                                                                       cAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAgg                                              >  3:357328/1‑149 (MQ=255)
                                                                                                       cAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAgg                                              >  8:213993/1‑149 (MQ=255)
                                                                                                            ccgccgCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGt                                                                   <  5:334206/123‑1 (MQ=255)
                                                                                                            ccgccgCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGt                                                                   >  6:334206/1‑123 (MQ=255)
                                                                                                                       aaaaTCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGa                               <  6:145039/148‑1 (MQ=255)
                                                                                                                       aaaaTCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGa                               <  6:70801/148‑1 (MQ=255)
                                                                                                                       aaaaTCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGa                               <  8:149911/148‑1 (MQ=255)
                                                                                                                         aaTCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGt                                        >  4:273514/1‑137 (MQ=255)
                                                                                                                         aaTCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGt                                        <  3:273514/137‑1 (MQ=255)
                                                                                                                          aTCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCgg                           <  3:46552/149‑1 (MQ=255)
                                                                                                                              gCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCGGCTc                        >  2:166954/1‑148 (MQ=255)
                                                                                                                              gCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCGGCTc                        <  1:166954/148‑1 (MQ=255)
                                                                                                                                gATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCGGCTCATa                     <  3:18484/149‑1 (MQ=255)
                                                                                                                                 aTGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATTGGCTCATAc                    <  4:210694/149‑1 (MQ=255)
                                                                                                                                     agagATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCGGCTCATACCctg                >  3:216515/1‑149 (MQ=255)
                                                                                                                                         aTTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCGGCTCATACCctgagct            >  4:70494/1‑149 (MQ=255)
                                                                                                                                         aTTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCATGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCGGCTCATACCctgagct            >  2:305642/1‑149 (MQ=255)
                                                                                                                                          ttggCCGGCGGCCTGCCGGCGAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCGGCTCATACCctgagct            >  3:87748/4‑148 (MQ=255)
                                                                                                                                          ttggCCGGCGGCCTGCCGGCGAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCGGCTCATACCctgagct            >  1:3031/4‑148 (MQ=255)
                                                                                                                                           tCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCGGCTCATACCctgagct            >  4:299352/1‑147 (MQ=255)
                                                                                                                                            ggCCGGCGGCCTGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCGGCTCATACCctg                <  2:199244/141‑1 (MQ=255)
                                                                                                                                            ggCCGGCGGCCTGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCGGCTCATACCctg                >  1:199244/2‑142 (MQ=255)
                                                                                                                                             gCCGGCATGCCGCCGGCCAAGTTCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCGGCTCATACCCtgagctgagc        <  5:331273/149‑1 (MQ=255)
                                                                                                                                             gCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCGGCTCATACCCtgagctgagc        >  8:106608/1‑149 (MQ=255)
                                                                                                                                               cGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCGGCTCATAcc                   <  2:293038/136‑1 (MQ=255)
                                                                                                                                               cGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCGGCTCATAcc                   >  1:293038/1‑136 (MQ=255)
                                                                                                                                                ggCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCGGCTCATACCctgagctga          >  7:319623/1‑144 (MQ=255)
                                                                                                                                                ggCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCGGCTCATACCctgagctga          <  8:319623/144‑1 (MQ=255)
                                                                                                                                                  cAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCGGCTCATACCCTGAGCTGAGCGCTGc   >  1:339756/1‑149 (MQ=255)
                                                                                                                                                   aGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCGGCTCATACCCTGAGCTGAGCGCTGcc  <  3:299352/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AACGACCAAGGCTACCTGGCCGTTGTGTGTCAGGGGCATTTCGACGCCATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCGGCTCATACCCTGAGCTGAGCGCTGCC  >  NC_002947/1752915‑1753210

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