Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,686,417 C→G 5.2% intergenic (‑134/+74) PP_4986 ← / ← PP_4987 hemolysin III family channel protein/putative Chemotaxis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,686,4170CG 5.2% 408.4 / 5.9 124intergenic (‑134/+74)PP_4986/PP_4987hemolysin III family channel protein/putative Chemotaxis protein
Reads supporting (aligned to +/- strand):  ref base C (61/56);  new base G (2/5);  total (63/61)
Fisher's exact test for biased strand distribution p-value = 2.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 129 total aligned reads displayed.

TCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAATGCGACCAGGTCCGGTACCCGCGCCAC  >  NC_002947/5686270‑5686565
                                                                                                                                                   |                                                                                                                                                    
tCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCgg                                                                                                                                                    <  2:1807658/150‑1 (MQ=255)
tCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCgg                                                                                                                                                    <  1:643738/150‑1 (MQ=255)
   ccGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGgc                                                                                                                                                 <  2:2294050/150‑1 (MQ=255)
    cGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGgct                                                                                                                                                <  1:80417/150‑1 (MQ=255)
    cGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGgct                                                                                                                                                <  1:80437/150‑1 (MQ=255)
      tAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGgctgc                                                                                                                                              <  1:454444/150‑1 (MQ=255)
      tAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCGTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGgctgc                                                                                                                                              >  2:377697/1‑150 (MQ=255)
         ggtgcgtgacgcggatctgcaacctgaccttggttaggatgtgaggcgcagggcctggtggtctggaGTGGGAGCGGCCTTGTGTCGCGATAGGGCTGCGCAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCtgctgc                                                                                                                                           >  1:5903667‑M1/68‑150 (MQ=255)
           cATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCGTGTGGCGCGCATGGGGCGCGAAGCGGCCCCAGGCTTTCAATGGCACGGCCAGCACCCCCGGGGCGGCGgcgc                                                                                                                                         >  1:3394502/1‑150 (MQ=255)
            aTGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCa                                                                                                                                        <  1:3899763/150‑1 (MQ=255)
            aTGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCa                                                                                                                                        <  1:4031416/150‑1 (MQ=255)
             tGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTGGCGAAAGGGGCGCGAAGCGGCCCCACGCGTTAATTGCCCCTGCATGAAACGCCGGGGGTGCTtcggcc                                                                                                                                       >  2:3113108/1‑147 (MQ=38)
             tGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCCCTGCAAGACCCGCCGGGGCTGCTgcgccg                                                                                                                                       >  2:4502727/1‑148 (MQ=255)
             tGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCCAAAGGGGCGCGAAGCGGCCCCACGCATTCATTTCCACAGCATGGACCGCCCGGGCTGGTGCGCAg                                                                                                                                       >  2:2690942/1‑150 (MQ=38)
              gCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGc                                                                                                                                      <  2:6682378/150‑1 (MQ=255)
              gCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGc                                                                                                                                      <  2:3584555/150‑1 (MQ=255)
              gCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGc                                                                                                                                      <  2:3934530/150‑1 (MQ=255)
              gCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCATGTGTCGCGAAACGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCAATGCATGAAACGCCGGGGCTGCTGcgccgc                                                                                                                                      >  2:5131201/1‑150 (MQ=255)
                  acacAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGcccc                                                                                                                                   <  1:2352879/149‑2 (MQ=255)
                   ccccgccctcctGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTAt                                                                                                                                 <  2:1146075/144‑1 (MQ=255)
                   cacaGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTAt                                                                                                                                 <  2:1146260/150‑1 (MQ=255)
                    cggatctgcaacctgaccttggttaggatgtgaggcgcagggcctggtggtctggaGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATc                                                                                                                                >  1:1136482‑M1/57‑150 (MQ=255)
                     caGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATcg                                                                                                                               <  1:3891785/150‑1 (MQ=255)
                      gatctgcaacctgaccttggttaggatgtgaggcgcagggcctggtggtctggaGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATcgc                                                                                                                              >  1:5217361‑M1/55‑150 (MQ=255)
                      gatctgcaacctgaccttggttaggatgtgaggcgcagggcctggtggtctggaGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATcgc                                                                                                                              >  2:3895182‑M1/55‑150 (MQ=255)
                      aGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATcgc                                                                                                                              <  2:716935/150‑1 (MQ=255)
                      aGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCGATcgc                                                                                                                              >  2:3950339/1‑150 (MQ=255)
                        aCCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTTCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTAACGCGa                                                                                                                            >  2:458252/1‑150 (MQ=255)
                        aCCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGa                                                                                                                            <  1:4458014/150‑1 (MQ=255)
                          ctcctGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGaca                                                                                                                          <  2:319394/150‑1 (MQ=255)
                           tcctGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCGTTCATTGCCACAGCATGAACCGCCGGGGCGGCAGCGCAGCCCAAACGCGacac                                                                                                                         >  2:3095539/1‑150 (MQ=38)
                           gcaacctgaccttggttaggatgtgaggcgcagggcctggtggtctggaGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGacac                                                                                                                         >  2:1855923‑M1/50‑150 (MQ=255)
                             ctGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACaa                                                                                                                       <  2:1418658/150‑1 (MQ=255)
                             ctGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACaa                                                                                                                       <  2:5083577/150‑1 (MQ=255)
                              tGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAg                                                                                                                      <  1:1802018/150‑1 (MQ=255)
                              tGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAg                                                                                                                      <  1:5487424/150‑1 (MQ=255)
                               gAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAgg                                                                                                                     <  2:2253995/150‑1 (MQ=255)
                                aGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGc                                                                                                                    <  1:6462152/150‑1 (MQ=255)
                                         gttaggatgtgaggcgcagggcctggtggtctggaGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATc                                                                                                                                <  1:3785095‑M1/94‑1 (MQ=255)
                                         gttaggatgtgaggcgcagggcctggtggtctggaGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATc                                                                                                                                >  2:3785095‑M1/36‑129 (MQ=255)
                                          ggTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATcgc                                                                                                                              <  2:4747537/130‑1 (MQ=255)
                                          ggTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATcgc                                                                                                                              >  1:4747537/1‑130 (MQ=255)
                                                tgttgtGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTTTGTCGCGAAAGGGGCGCGGAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCGGCTGCGCCGCCCGATCGCGGACCAAGGCCGCGCCCCCACCAGac                                                                                                    >  2:5265632/1‑150 (MQ=255)
                                                        tGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCATGTTTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACAGCATGAAACGCCCGGGCTGCTGCGCCGCCCCAAAGCCACACAAGGCCCCACCCACACCACAACACCACGc                                                                                            >  2:6069380/1‑150 (MQ=21)
                                                                  ttCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCCCGAAGCGGCCCCAGGCTTTCCTTGGCACTGCATTAAACGCCGGGGCTGCTGCGCCGCCCTATAGCCACACAAGGCCGCTCCCACTACAGAACACCAGGCCCAGCGccac                                                                                  >  2:456667/1‑148 (MQ=25)
                                                                  ttCGTTCTTTGTGGGAGCGGCCGTGTTTGGCGGAATGGGCGCGAAGCCGCCCCAGGCGTTCATTGCCCATGCAAGTACCGCCCGGGGGGGTTGGCAGCCCCATAGCCGCCCAAGGGCGCACCCACGCCAGAACCCCCGCCCCAGCGccgc                                                                                  >  2:1973644/1‑148 (MQ=255)
                                                                   tCGTTCTTTGTGGGAGCGGCCGTGTGTGGCGCAAAGGGCGCGAAGCCGCCCCACCCGTTCAATTCCACAGCAAGAACACCCCGGGCGGCGGCGCCGCCCCAAACCCCCACAACGCCCCGCCCACAACACACCAACACGCCCAGCGccgcc                                                                                 >  2:656979/1‑147 (MQ=255)
                                                                      ttCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACAGCATGAACCGCCCGGGCTGCTGCGCAGCCCCATCGCCACACAACGCCCCGACCCCAACACAACACCACGCCCAGCCCCGcacaa                                                                              >  2:4958373/1‑149 (MQ=18)
                                                                      ttCTTTGTGGGAGCGGCCATGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCAATGCCACAGCATGGACCCCCGGGGCGGCTGCGCAGCCCAAACCCCACACAACGCCCCACCCCAACCAGAACACCACGCCCAGCCCCCcacaa                                                                              >  2:4958345/1‑149 (MQ=16)
                                                                           tgtgGGAGCGGCCTTGTGTCGCGAAAGGGGCGCCAAGCGGCCCCAGGCTTTCATTGGCACAGCATTAACCGCCCGGGCTGCTGCGCCGCCCGATCGCGACACAAGGCCGCGCCCCCGCCAGACCCCCAGGCCCAGCGCCGCACAACCaaa                                                                         >  2:6695627/1‑150 (MQ=21)
                                                                                     GCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAg                                                               <  1:1657262‑M1/150‑1 (MQ=255)
                                                                                      CCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCAATAGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAgg                                                              <  1:2789308‑M1/150‑1 (MQ=255)
                                                                                        TTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGtt                                                            <  2:5021863‑M1/150‑1 (MQ=255)
                                                                                        TTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACAACCTAACCAAGGTCAGGtt                                                            <  1:1411602‑M1/150‑1 (MQ=255)
                                                                                          GTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGc                                                          >  2:3796319‑M1/1‑150 (MQ=255)
                                                                                            GTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCGATCGCGCCACAAGGCCGCTCCCACACCAGACCACCCGGCCCTGCGCCACACATCCGAACCAAGGTCACGTTGCAg                                                        <  1:5893398‑M1/150‑1 (MQ=255)
                                                                                             TCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGGCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGa                                                       <  1:6232777‑M1/150‑1 (MQ=255)
                                                                                             TCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGGCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGa                                                       <  1:911410‑M1/150‑1 (MQ=255)
                                                                                             TCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGGCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGa                                                       <  1:911467‑M1/150‑1 (MQ=255)
                                                                                             TCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGa                                                       <  2:1482814‑M1/150‑1 (MQ=255)
                                                                                             TCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCGCACATCCTAACCAAGGTCAGGTTGCAGa                                                       <  1:5633820‑M1/150‑1 (MQ=255)
                                                                                               GCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCGATCGCGACACAACGCCGCTCCCCCACCACACCACCCGGCCCGGCGCCGCACCACCGAACCAAGGTCACGGTGCAGATc                                                     <  1:2360793‑M1/150‑1 (MQ=255)
                                                                                                CGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATcc                                                    <  1:2208055‑M1/150‑1 (MQ=255)
                                                                                                CGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATcc                                                    <  1:5674004‑M1/150‑1 (MQ=255)
                                                                                                 GAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCcg                                                   <  2:3512736‑M1/150‑1 (MQ=255)
                                                                                                 GAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCGATCGCGACACAAGGCCGCTCCCACACCAGACCACCAGGCCCAGCGCCTCACATCCTAAACAAGGTCACGTTGCAGATCcg                                                   <  1:4127569‑M1/150‑1 (MQ=255)
                                                                                                       ggCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCAc                                             >  2:5210962/1‑150 (MQ=255)
                                                                                                       ggCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCAc                                             >  2:5211011/1‑150 (MQ=255)
                                                                                                       ggCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCAc                                             >  2:6676937/1‑150 (MQ=255)
                                                                                                        gcgcGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAg                                                                    <  2:2190449/126‑1 (MQ=255)
                                                                                                        gcgcGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCGATCGCGACACAAGGCCGCGCCCCCTCCAGACCACCAGGCCCTGCGCCTCACATCCAAACCAAGGTCAGGGTGCAGATCCGCGTCACg                                            >  1:4649622/1‑150 (MQ=255)
                                                                                                         cgcgAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGc                                           >  2:5396059/1‑150 (MQ=255)
                                                                                                          gcgAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCAATCGCGACACAAGGCCGCGCCCACACCAGAACACCACGCCCAGCGCCGCACATCCAAAACAACGTCAGGTTTCAGATCCGCGTCACGCa                                          >  1:22167/1‑150 (MQ=37)
                                                                                                            gAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCAcc                                        >  2:4696810/1‑150 (MQ=255)
                                                                                                            gAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACACCAGACCACCACGCCCTGCGCCTCACATCCAAACCAAGGTCAGGTTGCAGATCCGCGTCACGCAcc                                        >  1:5002857/1‑150 (MQ=255)
                                                                                                                  gCCCCCGGCTTTCTTGGCCACTGCATGACCCGCGGGGGGTGCTGCGCAGCCCTAGCGCGACACACGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCtt                                  <  2:4936531/150‑1 (MQ=255)
                                                                                                                  gCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCtt                                  >  2:5562277/1‑150 (MQ=255)
                                                                                                                  gCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCtt                                  <  1:3512736/150‑1 (MQ=255)
                                                                                                                  gCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCtt                                  >  2:6446844/1‑150 (MQ=255)
                                                                                                                    cccTGGCTTTCTTGGCCACTGAATGAACCGCGGGGGCTGCTGCGCAGCCATATCGCTACACATGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTc                                <  2:513558/150‑1 (MQ=255)
                                                                                                                     ccAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTcc                               >  2:2826662/1‑150 (MQ=255)
                                                                                                                     ccAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTcc                               >  2:562022/1‑150 (MQ=255)
                                                                                                                     ccAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCAATCGCGACACCAGGCCGCTCCCACACCAGACCACCAGGCCCTGCGCCTCACATCCAAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCATTcc                               >  1:5483523/1‑150 (MQ=39)
                                                                                                                      cAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCa                              >  2:4489235/1‑150 (MQ=255)
                                                                                                                        ggCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAAt                            <  1:2826662/150‑1 (MQ=255)
                                                                                                                           tttCATTGCCACTGCATGAACCGCCGGGGCTTCTTCGCAGACCTATCGCGAAAAAAAGCCCCTCCCACACCACACCACCAAGCCCTGCGCCTCACAAACAAACAAAGGTAAGGTTTAAGATCCCCGTCAAGCAACACCCATTACAATgag                         >  1:3817262/1‑148 (MQ=255)
                                                                                                                           tttCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAATGCg                         >  2:2574037/1‑150 (MQ=255)
                                                                                                                           tttCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCGATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCAAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTATCCAATGCg                         >  1:3817356/1‑150 (MQ=255)
                                                                                                                           tttCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCGATCGCCACACAAGGCCGCTCCCACTCCCGACCACCAGGCCCAGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAATGCg                         >  2:7377/1‑150 (MQ=255)
                                                                                                                               aTTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAATGCGACCa                     >  2:6758345/1‑150 (MQ=255)
                                                                                                                                   ccACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAATGCGACCAGGTc                 >  2:3443354/1‑150 (MQ=255)
                                                                                                                                     aCTGCTTTAACCGCGGGGGCTGCTGCGCAGCCCTATCGCGACACACGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAATGCGACCAGGTCCg               <  2:3635919/150‑1 (MQ=255)
                                                                                                                                     aCTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAAAGCGACCAGGTCCg               >  2:2596548/1‑150 (MQ=255)
                                                                                                                                        gCATGAACCGCCGGGGCTGCTGCGCAGCCCAATCGCGACACAAGGCCCCGCCCACTCCAGACCACCAGGCCCAGCGCCGCACAACCAAACCCAGGTCACGTTGCAGATCCGCGTCACGCACCAGCCATTCCAAAGCGACCAGGTCCGGTa            >  1:6146741/1‑150 (MQ=255)
                                                                                                                                          aTGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAACGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAATGCGACCAGGTCCGGTAcc          >  2:6469224/1‑150 (MQ=255)
                                                                                                                                           tGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCAGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCATTCCAATGCGACCAGGTCCGGTAccc         >  2:518090/1‑150 (MQ=255)
                                                                                                                                               ccgccgGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAATGCGACCAGGTCCGGTACCCgcgc     <  2:5674004/150‑1 (MQ=255)
                                                                                                                                                  ccgGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAg                                                                    <  2:2222254/84‑1 (MQ=255)
                                                                                                                                                  ccgGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAg                                                                    >  1:2222254/1‑84 (MQ=255)
                                                                                                                                                  ccgGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAATGCGACCAGGTCCGGTACCCGCGCCAc  <  1:6602698/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TCACCGTAATACATGCAAACACAGACCTCCTGAGGTGGGCTGGGTTCCTGTTGTGATGCTGATGGCTTCGTTCTTTGTGGGAGCGGCCTTGTGTCGCGAAAGGGGCGCGAAGCGGCCCCAGGCTTTCATTGCCACTGCATGAACCGCCGGGGCTGCTGCGCAGCCCTATCGCGACACAAGGCCGCTCCCACTCCAGACCACCAGGCCCTGCGCCTCACATCCTAACCAAGGTCAGGTTGCAGATCCGCGTCACGCACCAGCCTTTCCAATGCGACCAGGTCCGGTACCCGCGCCAC  >  NC_002947/5686270‑5686565

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