Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 189,893 A→C 5.6% intergenic (+55/+102) PP_0165 → / ← PP_0166 GGDEF domain‑containing protein/HlyD family secretion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947189,8930AC 5.6% 879.1 / 25.3 274intergenic (+55/+102)PP_0165/PP_0166GGDEF domain‑containing protein/HlyD family secretion protein
Reads supporting (aligned to +/- strand):  ref base A (144/113);  new base C (5/12);  total (149/125)
Fisher's exact test for biased strand distribution p-value = 4.34e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 277 total aligned reads displayed.

GAGTCGAGACGGAAGGGGAGCGGTTGGTGTTGCTGGAGATGGGAGTAGGCGGGATTCAGGGGCAACTGGTGGGTGAGCCTGCGCCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAATGGCGAAACCCGCAGGTTGGCTATCACCGCTCGCGCAGCGCTTCCGAACGCGCCTTCATGATCGGCTTGAG  >  NC_002947/189744‑190042
                                                                                                                                                     |                                                                                                                                                     
   tCGAGACGGAAGGGGAGCGGTTGGTGTTGCTGGAGATGGGAGTAGGCGGGATTCAGGGGCAACTGGTGGGTGAGCCTGCGCCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCg                                                                                                                                                    >  2:5266268/1‑150 (MQ=255)
      agaCGGAAGGGGAGCGGTTGGTGTTGCTGGAGATGGGAGTAGGCGGGATTCAGGGGCAACTGGTGGGTGAGCCTGCGCCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGcc                                                                                                                                                 <  1:5055759/150‑1 (MQ=255)
        aCGGAAGGGGAGCGGTTGGTGTTGCTGGAGATGGGAGTAGGCGGGATTCAGGGGCAACTGGTGGGTGAGCCTGCGCCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCTACGCGGccct                                                                                                                                               <  1:4597509/150‑2 (MQ=255)
        aCGGAAGGGGAGCGGTTGGTGTTGCTGGAGATGGGAGTAGGCGGGATTCAGGGGCAACTGGTGGGTGAGCCTGCGCCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCTACGCGGccct                                                                                                                                               <  1:964402/150‑2 (MQ=255)
         cGGAAGGGGAGCGGTTGGTGTTGCTGGAGATGGGAGTAGGCGGGATTCAGGGGCAACTGGTGGGTGAGCCTGCGCCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCTACGCGGccctt                                                                                                                                              <  1:1378616/150‑3 (MQ=255)
          ggAAGGGGAGCGGTTGGTGTTGCTGGAGATGGGAGTAGGCGGGATTCAGGGGCAACTGGTGGGTGAGCCTGCGCCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCTACGCGGcccttt                                                                                                                                             <  2:5337012/150‑4 (MQ=255)
          ggAAGGGGAGCGGTTGGTGTTGCTGGAGATGGGAGTAGGCGGGATTCAGGGGCAACTGGTGGGTGAGCCTGCGCCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCa                                                                                                                                             >  2:33342/1‑150 (MQ=255)
          ggAAGGGGAGCGGTTGGTGTTGCTGGAGATGGGAGTAGGCGGGATTCAGGGGCAACTGGTGGGTGAGCCTGCGCCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCa                                                                                                                                             >  1:2695810/1‑150 (MQ=255)
            aaGGGGAGCGGTTGGTGTTGCTGGAGATGGGAGTAGGCGGGATTCAGGGGCAACTGGTGGGTGAGCCTGCGCCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAgg                                                                                                                                           <  2:1621290/150‑1 (MQ=255)
                 gAGCGGTTGGTGTTGCTGGAGATGGGAGTAGGCGGGATTCAGGGGCAACTGGTGGGTGAGCCTGCGCCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCAc                                                                                                                                      >  1:1001847/1‑150 (MQ=255)
                     ggTTGGTGTTGCTGGAGATGGGTGTAGGGGGGATTCAGGGGCAACTGGTGGGTGAGCCTGCGCGTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAg                                                                                                                                  <  1:6355075/150‑1 (MQ=255)
                          gtgtTGCTGGAGATGGGAGTAGGCGGGATTCAGGGGCAACTGGTGGGTGAGCCTGCGCCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGAtt                                                                                                                             >  2:2088617/1‑150 (MQ=255)
                            gtTGCTGGAGATGGGAGTAGGCGGGATTCAGGGGCAACTGGTGGGTGAGCCTGCGCCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTg                                                                                                                           >  2:6771271/1‑150 (MQ=255)
                                cTGGAGATGGGAGTAGGCGGGATTCAGGGGCAACTGGTGGGTGAGCCTGCGCCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAgg                                                                                                                       <  1:2226283/150‑1 (MQ=255)
                                cTGGAGATGGGAGTAGGCGGGATTCAGGGGCAACTGGTGGGTGAGCCTGCGCCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTACCGGCCCCCAGGTTCACGTAGCGATTTGCAgg                                                                                                                       >  2:1400977/1‑150 (MQ=255)
                                 tGGAGATGGGAGTAGGCGGGATTCAGGGGCAACTGGTGGGTGAGCCTGCGCCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAggg                                                                                                                      >  2:4261496/1‑150 (MQ=255)
                                 tGGAGATGGGAGTAGGCGGGATTCAGGGGCAACTGGTGGGTGAGCCTGCGCCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAggg                                                                                                                      >  1:4196738/1‑150 (MQ=255)
                                     gttGGGAGTCGGCGGGTTTCCGGGGCAAGGGGTGGGTGAGCCTGCGCCTTGGCGCTGAAGGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAACCGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGgctgc                                                                                                                  <  1:1532284/148‑1 (MQ=255)
                                     gaTGGGAGTAGGCGGGATTCAGGGGCAACTGGTGGGTGAGCCTGCGCCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGgctgc                                                                                                                  <  2:818168/150‑1 (MQ=255)
                                     gaTGGGAGTAGGCGGGATTCAGGGGCAACTGGTGGGTGAGCCTGCGCCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGgctgc                                                                                                                  >  1:5386852/1‑150 (MQ=255)
                                        gggAGTAGGCGGGATTCAGGGGCAACTGGTGGGTGAGCCTGCGCCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTg                                                                                                               >  2:3773448/1‑150 (MQ=255)
                                          gAGTAGGCGGGATTCAGGGGCAACTGGTGGGTGAGCCTGCGCCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTa                                                                                                             <  2:4860035/150‑1 (MQ=255)
                                            gTAGGCGGGATTCAGGGGCAACTGGTGGGTGAGCCTGCGCCTTGGCGCTGAATGGCCTCTCGCGTGCCTCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACa                                                                                                           >  1:4676157/1‑150 (MQ=255)
                                               ggggggTTTCTGGGGCAAGTGGTGGGTGAGCCTGCGCCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCATGGGCTAGAGCCGCAAACGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGcc                                                                                                        <  1:5874178/147‑1 (MQ=255)
                                               ggCGGGTTTCAGGGGCAACTGGTGGGTGAGCCGGCGCCTTGGCGCTGAAGGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGCGCCGCAAACGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGcc                                                                                                        <  1:4302457/150‑1 (MQ=255)
                                               ggCGGGATTCAGGGGCAACTGGTGGGTGAGCCTGCGCCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGcc                                                                                                        >  2:198429/1‑150 (MQ=255)
                                               ggCGGGATTCAGGGGCAACTGGTGGGTGAGCCTGCGCCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGcc                                                                                                        <  1:1938800/150‑1 (MQ=255)
                                                  gggATTCAGGGGCAACTGGTGGGTGAGCCTGCGCCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTg                                                                                                     >  1:6241404/1‑150 (MQ=255)
                                                      ttCTGGGGCAACTGGTGGGTGAGCCTGCGCCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTcgc                                                                                                 <  2:361676/150‑1 (MQ=255)
                                                      ttCAGGGGCAAGGGGTGGGTGAGCCTGCGCCTTGGCGGTGAAGGGCCTCGCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTCGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTcgc                                                                                                 <  1:66774/150‑1 (MQ=255)
                                                      ttCAGGGGCAACTGGTGGGTGAGCCTGCGCCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTcgc                                                                                                 <  1:4493448/150‑1 (MQ=255)
                                                       tCAGGGGCAACTGGTGGGTGAGCCTGCGCCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGAAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTcgcg                                                                                                >  2:5417593/1‑150 (MQ=255)
                                                        cAGGGGCAACTGGTGGGTGAGCCTGCGCCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGa                                                                                               >  1:1197354/1‑150 (MQ=255)
                                                          ggggCAACTGGTGGGTGAGCCTGCGCCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACa                                                                                             >  1:4686141/1‑150 (MQ=255)
                                                           gggCAACTGGTGGGTGAGCCTGCGCCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACAt                                                                                            >  1:3477033/1‑150 (MQ=255)
                                                             gCAACTGGTGGGTGAGCCTGCGCCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTCGGCGTGGGCTAGAGCCGCAAACGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGa                                                                                          <  1:5914789/150‑1 (MQ=255)
                                                                aCTGGTGGGTGAGCCTGCGCCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGc                                                                                       <  2:1826139/150‑1 (MQ=255)
                                                                  tggtggGTGAGCCTGCGCCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCg                                                                                     >  2:2399441/1‑150 (MQ=255)
                                                                  tggtggGTGAGCCTGCGCCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCg                                                                                     >  1:3242631/1‑150 (MQ=255)
                                                                    gtgggtgAGCCTGCGCCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTAc                                                                                                            <  1:5758854/125‑1 (MQ=255)
                                                                    gtgggtgAGCCTGCGCCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCt                                                                                   >  1:6159501/1‑150 (MQ=255)
                                                                      gggtgAGCCTGCGCCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAggg                                                                                                                      >  2:3347790/1‑113 (MQ=255)
                                                                        gtgAGCCTGCGCCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTa                                                                               >  2:4039356/1‑150 (MQ=255)
                                                                        gtgAGCCTGCGCCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTa                                                                               >  2:1717699/1‑150 (MQ=255)
                                                                         tgAGCCTGCGCCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGACCTGTCGCGACATGAGGCCGCTCCTaa                                                                              <  1:6130457/150‑1 (MQ=255)
                                                                          gAGCCTGCGCCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTaaa                                                                             <  2:6609858/150‑1 (MQ=255)
                                                                          gAGCCTGCGCCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGATAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTaaa                                                                             <  2:2883565/150‑1 (MQ=255)
                                                                            gCCTGCGCCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGc                                                                           >  1:1772912/1‑150 (MQ=255)
                                                                             ccTGCGCCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCa                                                                          >  1:2552772/1‑150 (MQ=255)
                                                                               tGCGCCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTATAGCCGAAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAAt                                                                        <  1:2290139/150‑1 (MQ=255)
                                                                               tGCGCCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAAt                                                                        >  1:2159731/1‑150 (MQ=255)
                                                                               tGCGCCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAAt                                                                        <  1:2011557/150‑1 (MQ=255)
                                                                               tGCGCCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAAt                                                                        >  2:4594904/1‑150 (MQ=255)
                                                                                  gcCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAATGGc                                                                     <  2:3977958/150‑1 (MQ=255)
                                                                                  gcCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAATGGc                                                                     <  1:5274409/150‑1 (MQ=255)
                                                                                  gcCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAATGGc                                                                     <  1:4440326/150‑1 (MQ=255)
                                                                                   ccTTGGCGCTGAATGGCCTCTCGCGTGCCGCGGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAATGGCg                                                                    >  2:5511365/1‑150 (MQ=255)
                                                                                    cTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAATGGCGa                                                                   >  2:4725215/1‑150 (MQ=255)
                                                                                         cgcTGACTGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAATGGCGAAAccc                                                              <  2:4425553/150‑1 (MQ=255)
                                                                                             gAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAATGGCGAAACCCGCAg                                                          >  1:3555775/1‑150 (MQ=255)
                                                                                                   ccTCTCGCGTGCCGCTGTAGGCGTGGGCTAGCGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAATGGCGAAACCCGCAGGTTGGc                                                    <  2:4286874/150‑1 (MQ=255)
                                                                                                   ccTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAATGGCGAAACCCGCAGGTTGGc                                                    >  1:813334/1‑150 (MQ=255)
                                                                                                   ccTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAATGGCGAAACCCGCAGGTTGGc                                                    <  2:4320103/150‑1 (MQ=255)
                                                                                                   ccTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAATGGCGAAACCCGCAGGTTGGc                                                    <  2:476844/150‑1 (MQ=255)
                                                                                                   ccTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAATGGCGAAACCCGCAGGTTGGc                                                    >  1:5608504/1‑150 (MQ=255)
                                                                                                    ctctCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAATGGCGAAACCCGCAGGTTGGCt                                                   <  1:4887893/150‑1 (MQ=255)
                                                                                                      ctcGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAATGGCGAAACCCGCAGGTTGGCTAt                                                 <  2:6211196/150‑1 (MQ=255)
                                                                                                       tcGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAATGGCGAAACCCGCAGGTTGGCTATc                                                >  1:4336321/1‑150 (MQ=255)
                                                                                                         gcgTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAATGGCGAAACCCGCAGGTTGGCTATCAc                                              >  2:6180632/1‑150 (MQ=255)
                                                                                                         gcgTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAATGGCGAAACCCGCAGGTTGGCTATCAc                                              <  2:5136068/150‑1 (MQ=255)
                                                                                                             gCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAATGGCGAAACCCGCAGGTTGGCTATCACCGCt                                          >  2:613818/1‑150 (MQ=255)
                                                                                                               cGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAATGGCGAAACCCGCAGGTTGGCTATCACCGCTcg                                        >  2:2216725/1‑150 (MQ=255)
                                                                                                                       gCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAATGGCGAAACCCGCAGGTTGGCTATCACCGCTCGCGCAgcgc                                >  1:3046858/1‑150 (MQ=255)
                                                                                                                          tGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAATGGCGAAACCCGCAGGTTGGCTATCACCGCTCGCGCAGCGCTTc                             <  1:1866206/150‑1 (MQ=255)
                                                                                                                           gggCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAATGGCGAAACCCGCAGGTTGGCTATCACCGCTCGCGCAGCGCTTcc                            <  2:35221/150‑1 (MQ=255)
                                                                                                                             gCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAATGGCGAAACCCGCAGGTTGGCTATCACCGCTCGCGCAGCGCTTCCGa                          >  1:1135440/1‑150 (MQ=255)
                                                                                                                             gCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAATGGCGAAACCCGCAGGTTGGCTATCACCGCTCGCGCAGCGCTTCCGa                          >  1:358129/1‑150 (MQ=255)
                                                                                                                              cTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAATGGCGAAACCCGCAGGTTGGCTATCACCGCTCGCGCAGCGCTTCCGaa                         <  2:6069777/150‑1 (MQ=255)
                                                                                                                               tAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAATGGCGAAACCCGCAGGTTGGCTATCACCGCTCGCGCAGCGCTTCCGAAc                        >  2:3178791/1‑150 (MQ=255)
                                                                                                                                 gagCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAATGGCGAAACCCGCAGGTTGGCTATCACCGCTCGCGCAGCGCTTCCGAAcgc                      <  2:5744163/150‑1 (MQ=255)
                                                                                                                                  agCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAATGGCGAAACCCGCAGGTTGGCTATCACCGCTCGCGCAGCGCTTCCGAAcgcg                     >  2:4503436/1‑150 (MQ=255)
                                                                                                                                    ccGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAATGGCGAAACCCGCAGGTTGGCTATCACCGCTCGCGCAGCGCTTCCGAACGCGcc                   >  2:3818540/1‑150 (MQ=255)
                                                                                                                                       cAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAATGGCGAAACCCGCAGGTTGGCTATCACCGCTCGCGCAGCGCTTCCGAACGCGCCTTc                <  1:3113553/150‑1 (MQ=255)
                                                                                                                                       cAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAATGGCGAAACCCGCAGGTTGGCTATCACCGCTCGCGCAGCGCTTCCGAACGCGCCTTc                <  1:4050984/150‑1 (MQ=255)
                                                                                                                                        ctgggggCCGCTACGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAATGGCGAAACCCGCAGGTTGGCTATCACCGCTCGCGCAGCGCTTCCGAACGCGCCTTCa               >  2:1660096/5‑150 (MQ=255)
                                                                                                                                        ctgggggCCGCTACGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAATGGCGAAACCCGCAGGTTGGCTATCACCGCTCGCGCAGCGCTTCCGAACGCGCCTTCa               >  2:2139646/5‑150 (MQ=255)
                                                                                                                                        ctgggggCCGCTACGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAAAGGCGAAACCCGCAGGTTGGCTATCACCGCTCGCGCAGCGCTTCCGAACGCGCCTTCa               >  1:721732/5‑150 (MQ=255)
                                                                                                                                        aaaaGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAATGGCGAAACCCGCAGGTTGGCTATCACCGCTCGCGCAGCGCTTCCGAACGCGCCTTCa               >  1:5795794/1‑150 (MQ=255)
                                                                                                                                        aaaaGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAATGGCGAAACCCGCAGGTTGGCTATCACCGCTCGCGCAGCGCTTCCGAACGCGCCTTCa               <  1:2989385/150‑1 (MQ=255)
                                                                                                                                        aaaaGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAATGGCGAAACCCGCAGGTTGGCTATCACCGCTCGCGCAGCGCTTCCGAACGCGCCTTCa               <  1:2088617/150‑1 (MQ=255)
                                                                                                                                        aaaaGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAATGGCGAAACCCGCAGGTTGGCTATCACCGCTCGCGCAGCGCTTCCGAACGCGCCTTCa               >  1:1424108/1‑150 (MQ=255)
                                                                                                                                         tgggggCCGCTACGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAATGGCGAAACCCGCAGGTTGGCTATCACCGCTCGCGCAGCGCTTCCGAACGCGCCTTCAt              >  2:1166579/4‑150 (MQ=255)
                                                                                                                                         aaaGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAATGGCGAAACCCGCAGGTTGGCTATCACCGCTCGCGCAGCGCTTCCGAACGCGCCTTCAt              <  1:3880670/150‑1 (MQ=255)
                                                                                                                                         aaaGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAATGGCGAAACCCGCAGGTTGGCTATCACCGCTCGCGCAGCGCTTCCGAACGCGCCTTCAt              <  2:2054604/150‑1 (MQ=255)
                                                                                                                                                cgcgTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAATGGCGAAACCCGCAGGTTGGCTATCACCGCTCGCGCAGCGCTTCCGAACGCGCCTTCATGATCGGc       >  1:3667189/1‑150 (MQ=255)
                                                                                                                                                 gcgTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAATGGCGAAACCCGCAGGTTGGCTATCACCGCTCGCGCAGCGCTTCCGAACGCGCCTTCATGATCGGCt      >  1:6537848/1‑150 (MQ=255)
                                                                                                                                                 gcgTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAATGGCGAAACCCGCAGGTTGGCTATCACCGCTCGCGCAGCGCTTCCGAACGCGCCTTCATGATCGGCt      <  2:1197354/150‑1 (MQ=255)
                                                                                                                                                  cgTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAATGGCGAAACCCGCAGGTTGTCTATCACCGCTCGCGCAGCGCTTCCGAACGCGCCTTCATGATCGGCtt     <  2:311434/150‑1 (MQ=255)
                                                                                                                                                   gTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAATGGCGAAACCCGCAGGTTGGCTATCACCGCTCGCGCAGCGCTTCCGAACGCGCCTTCATGATCGGCTTg    <  2:5670007/150‑1 (MQ=255)
                                                                                                                                                    tAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAATGGCGAAACCCGCAGGTTGGCTATCACCGCTCGCGCAGCGCTTCCGAACGCGCCTTCATGATCGGCTTGa   <  2:2695810/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAGTCGAGACGGAAGGGGAGCGGTTGGTGTTGCTGGAGATGGGAGTAGGCGGGATTCAGGGGCAACTGGTGGGTGAGCCTGCGCCTTGGCGCTGAATGGCCTCTCGCGTGCCGCTGTAGGCGTGGGCTAGAGCCGCAAAAGGGCCGCGTAGCGGCCCCCAGGTTCACGTAGCGATTTGCAGGGCTGCTTGTACAGCCCTGTCGCGACATGAGGCCGCTCCTAAAGCAATGGCGAAACCCGCAGGTTGGCTATCACCGCTCGCGCAGCGCTTCCGAACGCGCCTTCATGATCGGCTTGAG  >  NC_002947/189744‑190042

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: