Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,753,073 C→G 6.9% Q15E (CAA→GAA)  PP_1555 → putative Phage integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,753,0730CG 6.9% 230.8 / 5.6 72Q15E (CAA→GAA) PP_1555putative Phage integrase
Reads supporting (aligned to +/- strand):  ref base C (36/31);  new base G (2/3);  total (38/34)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.78e-01

CCTGGCCGTTGTGTGTCAGGGGCATTTCGACGCCATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCGGCTCATACCCTGAGCTGAGCGCTGCCAAGGTGGCCG  >  NC_002947/1752929‑1753220
                                                                                                                                                |                                                                                                                                                   
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                                                                         <  3:144728/149‑1 (MQ=255)
                                                                        ggAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTg                                                                         <  6:311291/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                                         <  8:149911/148‑1 (MQ=255)
                                                                                                         aaaaTCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGa                                         <  6:70801/148‑1 (MQ=255)
                                                                                                         aaaaTCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGa                                         <  6:145039/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                                  <  1:166954/148‑1 (MQ=255)
                                                                                                                gCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCGGCTc                                  >  2:166954/1‑148 (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                      >  1:3031/4‑148 (MQ=255)
                                                                                                                            ttggCCGGCGGCCTGCCGGCGAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCGGCTCATACCctgagct                      >  3:87748/4‑148 (MQ=255)
                                                                                                                             tCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCGGCTCATACCctgagct                      >  4:299352/1‑147 (MQ=255)
                                                                                                                              ggCCGGCGGCCTGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCGGCTCATACCctg                          >  1:199244/2‑142 (MQ=255)
                                                                                                                              ggCCGGCGGCCTGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCGGCTCATACCctg                          <  2:199244/141‑1 (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)
                                                                                                                                           ccgGCGAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCGGCTCATACCCTGAGCTGAGCGCTGCCAAGGTg      <  4:87748/149‑1 (MQ=255)
                                                                                                                                               ccAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCGGCTCATACCCTGAGCTGAGCGCTGCCAAGGTGGCcg  <  6:172908/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
CCTGGCCGTTGTGTGTCAGGGGCATTTCGACGCCATGGCGTGCCTGAGGGCGTACCTGGCCCTGCCAAAAACGGAGGTTGCAGCATGACCAACACCGCCGCTGTGAAAATCAGCGATGCAGAGATTCGCCGGCAGGCCGCCGGCCAAGTGCGTGACCTTCGCGCCCTGGGCAACCACGGCCTGTATTTCCGCTTTCATCGCTCCCGCGAGCGCGGGTCGTGGTACCTGATCCACAAGGGCAAGTGGAACCTGATCGGCTCATACCCTGAGCTGAGCGCTGCCAAGGTGGCCG  >  NC_002947/1752929‑1753220

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: ‑