Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,424,890 C→G 12.8% intergenic (+478/‑34) PP_3039 → / → PP_3040 hypothetical protein/holin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,424,8900CG12.8% 194.1 / 19.3 77intergenic (+478/‑34)PP_3039/PP_3040hypothetical protein/holin
Reads supporting (aligned to +/- strand):  ref base C (32/35);  new base G (3/7);  total (35/42)
Fisher's exact test for biased strand distribution p-value = 3.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.94e-01

AACAGACATTCCCCGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGC‑‑CACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGCCGGCGT  >  NC_002947/3424746‑3425033
                                                                                                                                                  |                                                                                                                                               
aaCAGACATTCCCCGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGC‑‑CACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTgg                                                                                                                                             >  6:285515/1‑149 (MQ=255)
aaCAGACATTCCCCGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGC‑‑CACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTgg                                                                                                                                             >  4:60059/1‑149 (MQ=255)
 aCAGACATTCCCCGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGCAACACCGGAATCCAGCTTTAAAATGTACACATATTCGAGAAGAGCGCGCGCACAAAGCCGCCCGCAAAACACTCAAAACCCCGGCCGTg                                                                                                                                              >  1:108214/1‑149 (MQ=255)
 aCAGACATTCCCCGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGC‑‑CACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCCCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGgc                                                                                                                                            >  7:161658/1‑149 (MQ=255)
       aTTCCCCGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGC‑‑CACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCGCCAGGGCCggg                                                                                                                                      >  4:153964/1‑149 (MQ=255)
       aTTCCCCGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGC‑‑CACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGAACAAAGCGGCCCGCCAAACACTGAAAACCCCCGCCCTGGCGCCggg                                                                                                                                      >  7:325533/1‑149 (MQ=255)
        ttCCCCGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGT‑‑CACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCgggg                                                                                                                                     >  6:363826/1‑149 (MQ=255)
        ttCCCCGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGC‑‑CACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCGCCAGGGCCgggg                                                                                                                                     <  7:46253/149‑1 (MQ=255)
          ccccGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGC‑‑CACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCGCCAGGGCCGGGGt                                                                                                                                    <  7:325710/148‑1 (MQ=255)
           cccGGACGTTTAACACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGC‑‑CACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCGCCAGGGCCGGGGttt                                                                                                                                  <  2:347413/149‑1 (MQ=255)
           cagggACGTTTAACACTTCTTGCCGGGTGGCATTAAAACTCGCTTGCTGC‑‑CACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGttt                                                                                                                                  <  4:31631/146‑1 (MQ=255)
            ccGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGC‑‑CACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGtttt                                                                                                                                 >  5:188431/1‑149 (MQ=255)
             cGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGC‑‑CACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGttttt                                                                                                                                <  1:199080/149‑1 (MQ=255)
              ggACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGC‑‑CACCGGAATCGACCTGTAAAAAGTACCCATCTTGGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCGCCAGGGCCGGGGttttca                                                                                                                               <  1:82402/149‑3 (MQ=255)
              ggACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGC‑‑CACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTg                                                                                                                               >  7:282624/1‑149 (MQ=255)
               gACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGC‑‑CACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCGCCAGGGCCGGGGttttcag                                                                                                                              <  3:153964/149‑4 (MQ=255)
               gACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGC‑‑CACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGc                                                                                                                              <  5:322469/149‑1 (MQ=255)
               caCGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGC‑‑CACCGGAATCGACCTTTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCCGCGCCAGGGCCGGGGttttcag                                                                                                                              <  7:56008/148‑4 (MQ=255)
                aCGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGC‑‑CACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCGCCAGGGCCGGGGttttcagt                                                                                                                             <  6:300244/149‑5 (MQ=255)
                     tACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGC‑‑CACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTg                                                                                                                         <  2:352104/148‑1 (MQ=255)
                       ccACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGC‑‑CACCGGAATCGACCTGTAAAATGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTggg                                                                                                                       <  2:73217/148‑1 (MQ=255)
                       ccACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGC‑‑CACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTggg                                                                                                                       <  3:60059/148‑1 (MQ=255)
                        cACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGC‑‑CACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTggggg                                                                                                                     <  7:190540/149‑1 (MQ=255)
                          cTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGC‑‑CACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTgg                                                                                                                        <  6:157387/144‑1 (MQ=255)
                          cTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGC‑‑CACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTgg                                                                                                                        >  5:157387/1‑144 (MQ=255)
                           ttATTGCCGGGTGGCATTAAAACTCGCTTGCTGC‑‑CACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGAt                                                                                                                   >  8:97371/1‑148 (MQ=255)
                              ttGCCGGGTGGCATTAAAACTCGCTTGCTGC‑‑CACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGa                                                                                                               >  6:30448/1‑149 (MQ=255)
                                gCCGGGTGGCATTAAAACTCGCTTGCTGC‑‑CACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGAt                                                                                                              <  6:124432/148‑1 (MQ=255)
                                  cGGGTGGCATTAAAACTCGCTTGCTGC‑‑CACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGaa                                                                                                           <  8:161658/149‑1 (MQ=255)
                                         cgttAAAACTCGCTTGCTGC‑‑CACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGag                                                                                                    <  3:176466/147‑1 (MQ=255)
                                          aTTAAAACTCGCTTGCTGC‑‑CACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCgcg                                                                                                    >  2:159958/1‑146 (MQ=255)
                                          aTTAAAACTCGCTTGCTGC‑‑CACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGag                                                                                                    >  8:180340/1‑148 (MQ=255)
                                           ttAAAACTCGCTTGCTGC‑‑CACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGagca                                                                                                  <  5:30448/149‑1 (MQ=255)
                                                aCTCGCTTACTGC‑‑CACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAg                                                                                             >  5:260397/1‑149 (MQ=255)
                                                  tCGCTTGCTGC‑‑CACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCg                                                                                           >  2:326834/1‑149 (MQ=255)
                                                    gCTTGCTGC‑‑CACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGt                                                                                          <  3:349385/148‑1 (MQ=255)
                                                    gCTTGCTG‑‑CACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGc                                                                                                       >  1:290924/1‑134 (MQ=255)
                                                    gCTTGCTG‑‑CACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGc                                                                                                       <  2:290924/134‑1 (MQ=255)
                                                     cTTACTGC‑‑CACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTa                                                                                        <  6:260397/149‑1 (MQ=255)
                                                        gctgc‑‑CACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGCGCAGCAAGCGTTAGtt                                                                                     >  2:24149/1‑149 (MQ=255)
                                                         ctgc‑‑CACCGGAATCTAACAGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGtt                                                                                     <  8:251576/148‑1 (MQ=255)
                                                            c‑‑cACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCg                                                                                                      >  8:197848/1‑128 (MQ=255)
                                                            c‑‑cACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCg                                                                                                      <  7:197848/128‑1 (MQ=255)
                                                            c‑‑cACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTg                                                                                    <  6:188431/146‑1 (MQ=255)
                                                                   ggAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGcc                                                                             >  1:140657/1‑148 (MQ=255)
                                                                    gAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCaa                                                                           >  1:223306/1‑149 (MQ=255)
                                                                          aCCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCtggt                                                                     >  4:132270/1‑149 (MQ=255)
                                                                           ccTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCtggtt                                                                    <  7:97371/149‑1 (MQ=255)
                                                                            cTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCtggttg                                                                   <  5:285515/149‑1 (MQ=255)
                                                                               tAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAGGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGAt                                                                                >  3:285836/1‑133 (MQ=255)
                                                                               tAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAGGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGAt                                                                                <  4:285836/133‑1 (MQ=255)
                                                                               tAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCtggttggt                                                                 <  1:24149/148‑1 (MQ=255)
                                                                                  aaaGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATc                                                             <  2:223306/149‑1 (MQ=255)
                                                                                  aaaGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCATGGCGCCGGGGTTTTTGCCTTTGGTGATTCGATTAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCt                                                            <  8:177769/149‑1 (MQ=255)
                                                                                   aaGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCt                                                            <  1:326834/149‑1 (MQ=255)
                                                                                       aCCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTgaga                                                                                 >  3:61696/1‑124 (MQ=255)
                                                                                       aCCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTgaga                                                                                 <  4:61696/124‑1 (MQ=255)
                                                                                          cATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTc                                                          <  5:60003/144‑1 (MQ=255)
                                                                                           aTCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATc                                                             >  5:38141/1‑140 (MQ=255)
                                                                                           aTCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATc                                                             <  6:38141/140‑1 (MQ=255)
                                                                                              ttCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAAACCCGGCCCAGGCGCCGGGGGTTTTGGCGTTGGGGGTTTGGATGAACAGCGAGCCGCAAGGGTTATTTGAGATGCCAATCTGGTTGGTGTTTTTCCTGTCCCTg                                                 >  2:353597/1‑149 (MQ=255)
                                                                                               tCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACCCTGAAAACCCCGGCCCTGGGGCCGGGGGTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGGTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCGGTCCCTgg                                                >  7:299361/1‑149 (MQ=255)
                                                                                                cGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTg                                                 <  7:59379/147‑1 (MQ=255)
                                                                                                      gTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTcggcgg                                         >  7:269475/1‑149 (MQ=255)
                                                                                                       tGCGGGCTCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCggg                                        <  3:204659/149‑1 (MQ=255)
                                                                                                       tGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCggg                                        <  7:180340/149‑1 (MQ=255)
                                                                                                            gcgcACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTc                                   <  1:181936/149‑1 (MQ=255)
                                                                                                                        gCCCGCCAAAGACTGAAACCCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGAtgtg                       <  8:325533/149‑1 (MQ=255)
                                                                                                                         cccGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTgg                      >  6:123089/1‑149 (MQ=255)
                                                                                                                              cAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCTGGGTGTCAGGCGAGATGTGGCGGGc                 >  4:253087/1‑149 (MQ=255)
                                                                                                                                     tGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATg          >  1:241835/1‑149 (MQ=255)
                                                                                                                                     caaaaaCCCCGGCGCCAGGGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATg          >  8:225456/3‑149 (MQ=255)
                                                                                                                                      gAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATgg         >  5:342956/1‑149 (MQ=255)
                                                                                                                                       aaaaCCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATg          <  2:241835/147‑1 (MQ=255)
                                                                                                                                           ccccGGCGCCAGGGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTACCTGGTCGGAGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGCCGGc    >  7:336133/1‑149 (MQ=255)
                                                                                                                                           ccccGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGCCGGc    <  8:282624/149‑1 (MQ=255)
                                                                                                                                            cccGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTCAGATGCCATTCTGGTTGGTTAACTTCCTGTCCCTGGTCGGCGGGTTGTCAGGCGAGATGTGGCGGGCCGACATGGACGGCg   >  1:59051/1‑149 (MQ=255)
                                                                                                                                             ccctcccTGTGGCCGTGGTTTTTGCCTTTGGGGGAATCGATGAACAGCGAGCATCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTTGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGGCGACATGGCCGGCGt  <  6:257988/145‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
AACAGACATTCCCCGGACGTTTACCACTTATTGCCGGGTGGCATTAAAACTCGCTTGCTGC‑‑CACCGGAATCGACCTGTAAAAAGTACCCATCTTCGAGACGTGCGGGCGCACAAAGCGGCCCGCCAAACACTGAAAACCCCGGCCCTGGCGCCGGGGTTTTTGCCTTTGGGGGATTCGATGAACAGCGAGCAGCAAGCGTTAGTTGAGATGCCAATCTGGTTGGTGATCTTCCTGTCCCTGGTCGGCGGGGTGTCAGGCGAGATGTGGCGGGCCGACATGGCCGGCGT  >  NC_002947/3424746‑3425033

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