Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,992,862 C→G 7.2% D71H (GAT→CAT)  PP_2616 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,992,8620CG 7.2% 213.8 / 6.8 69D71H (GAT→CAT) PP_2616hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (32/32);  new base G (2/3);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.33e-01

TGAGCCGGCTCGTAACCGTGCTGGCGCATTTCTTCCTGAATGTTGTCGATCTGGCTGCGTACTTGCTGGTGCACGTCCTGGATGTTCGGCAGCTCGTCCAGCGTA‑CGCAGGCGTATGACCAAGCCGAACAACGGTGATGCGGCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGAAATTCCGGGTCAGCGGGGTAGCCTTCGGAAACGGGCTCTGGCTCAGGTTTGTCCGTGGATGCTTGCCCATTGGTTGCTGCAGCCCGGGGGGCGGGTACTGGCTGTTGAAG  >  NC_002947/2992717‑2993006
                                                                                                                                                  |                                                                                                                                                   
tGAGCCGGCTCGTAACCGTGCTGGCGCATTTCTTCCTGAATGTTGTCGATCTGGCTGCGTACTTGCTGGTGCACGTCCTGGATGTTCGGCAGCTCGTCCAGCGTA‑CGCAGGCGTATGACCAAGCCGAACAACGGTGATGCGGCATCCAg                                                                                                                                                  >  8:287934/1‑149 (MQ=255)
        cTCGTAACCGTGCTGGCGCATTTCTTCCTGAATGTTGTCGATCTGGCTGCGTACTTGCTGGTGCACGTCCTGGATGTTCGGCAGCTCGTCCAGCGTA‑CGCAGGCGTATGACCAAGCCGAACAACGGTGATGCTGGATGCCGCATCaccg                                                                                                                                          <  6:182596/149‑4 (MQ=255)
         tCGTAACCGTTCTGGCGGATTTCTTGCTGAATGTTGTCGATCTGGCTGCGTACTTGCTGGTGCACGTCCTGGATGTTCGGCAGCTCGTCGAGCGTA‑CGCAGGCGTATGACCAAGCCGAACAAGGGTGATGGGGCATCCAGCATCAGGtt                                                                                                                                         <  2:45248/149‑1 (MQ=255)
         tCGTAACCGTGCTGGCGCATTTCTTCCTGAATGTTGTCGATCTGGCTGCGTACTTGCTGGTGCACGTCCTGGATGTTCGGCAGCTCGTCCAGCGTA‑CGCAGGCGTATGACCAAGCCGAACAACGGTGATGCTGGATGCCGCATCa                                                                                                                                             <  6:368629/145‑1 (MQ=255)
         tCGTAACCGTGCTGGCGCATTTCTTCCTGAATGTTGTCGATCTGGCTGCGTACTTGCTGGTGCACGTCCTGGATGTTCGGCAGCTCGTCCAGCGTA‑CGCAGGCGTATGACCAAGCCGAACAACGGTGATGCTGGATGCCGCATCa                                                                                                                                             >  5:368629/1‑145 (MQ=255)
         tCGTAACCGTGCTGGCGCATTTCTTCCTGAATGTTGTCGATCTGGCTGCGTACTTGCTGGTGCACGTCCTGGATGTTCGGCAGCTCGTCCAGCGTA‑CGCAGGCGTATGACCAAGCCGAACAACGGTGATGCGGCATCCAGCATCAGGtt                                                                                                                                         >  6:228020/1‑149 (MQ=255)
          cGTAACCGTGCTGGCGCATTTCTTCCTGAATGTTGTCGATCTGGCTGCGTACTTGCTGGTGCACGTCCTGGATGTTCGGCAGCTCGTCCAGCGTA‑CGCAGGCGTATGACCAAGCCGAACAACGGTGATGCGGCATCCAGCATCAGGTTg                                                                                                                                        >  1:179305/1‑149 (MQ=255)
              aCCGTGCTGGCGCATTTCTTCCTGAATGTTGTCGATCTGGCTGCGTACTTGCTGGTGCACGTCCTGGATGTTCGGCAGCTCGTCCAGCGTA‑CGCAGGCGTATGACCAAGCCGAACAACGGTGATGCGGCATCCAGCATCAGGTTGGCg                                                                                                                                     <  5:228020/148‑1 (MQ=255)
                         gcATTTATTCCTGAATGTTGTCGATCTGGCTGCGTACTTGCTGGTGCACGTCCTGGATGTACGGCAGCTCGTCCAGCGTA‑CGCAGGCGTATGACCAAGCCGAACAACGGTGATGCGGCATCCAGAAACAGGTTGGCGAACCCGC‑‑‑CCCg                                                                                                                       >  8:165300/1‑148 (MQ=255)
                              tcttcCTGAATGTTGTCGATCTGGCTGCGTACTTGCTGGTGCACGTCCTGGATGTTCGGCAGCTCGTCCAGCGTA‑CGCAGGCGTATGACCAAGCCGAACAACGGTGATGCGGCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCt                                                                                                                  <  2:179305/148‑1 (MQ=255)
                                 tcCTGAATGTTGTCGATCTGGCTGCGTACTTGCTGGTGCACGTCCTGGATGTTCGGCAGCTCGTCCAGCGTA‑CGCAGGCGTATGACCAAGCCGAACAACGGTGATGCGGCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGaaa                                                                                                              >  5:24965/1‑149 (MQ=255)
                                    tGAATGTTGTCGATCTGGCTGCGTACTTGCTGGTGCACGTCCTGGATGTTCGGCAGCTCGTCCAGCGTA‑CGCAGGCGTATGACCAAGCCGAACAACGGTGATGCTGCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCGGAAATTc                                                                                                           >  3:301440/1‑149 (MQ=255)
                                        tgttgtCGATCTGGCTGCGTACTTGCTGGTGCACGTCCTGGATGTTCGGCAGCTCGTCCAGCGTA‑CGCAGGCGTATGACCAAGCCGAACAACGGTGATGCGGCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGAAATTCCggg                                                                                                       <  8:115098/149‑1 (MQ=255)
                                          ttgtCGATCTGGCTGCGTACTTGCTGGTGCACGTCCTGGATGTTCGGCAGCTCGTCCAGCGTATGGCAGGCGTATGACCAAGCCGAACAACGGTGATGCGGCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGAAATTCCGGGt                                                                                                      >  1:296567/1‑149 (MQ=255)
                                            gtCGATCTGGCTGCGTACTTGCTGGTGCACGTCCTGGATGTTCGGCAGCTCGTCCAGCGTA‑CGCAGGCGTATGACCAAGCCGAACAACGGTGATGCGGCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGAAATTCCGGGTCAg                                                                                                   >  8:209617/1‑149 (MQ=255)
                                                  cTGGCTGCGTACTTGCTGGTGCATGTCCTGGATGTTCGGCAGCTCGTCCAGCGTA‑CGCAGGCGTATGACCAAGCCGAACAACGGTGATGCGGCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGAAATTCCGGGTCAGCgggg                                                                                              >  2:124440/1‑148 (MQ=255)
                                                     gCTGCGTACTTGCTGGTGCACGTCCTGGATGTTCGGCAGCTCGTCCAGCGTA‑CGCAGGCGTATGACCAAGCCGAACAACGGTGATGCGGCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGAAATTCCGGGTCAGCGGGGTa                                                                                            <  5:261500/147‑1 (MQ=255)
                                                     gCTGCGTACTTGCTGGTGCACGTCCTGGATGTTCGGCAGCTCGTCCAGCGTA‑CGCAGGCGTATGACCAAGCCGAACAACGGTGATGCGGCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGAAATTCCGGGTCAGCGGGGTAGc                                                                                          >  8:158552/1‑149 (MQ=255)
                                                       tGCGTACTTGCTGGTGCACGTCCTGGATGTTCGGCAGCTCGTCCAGCGTA‑CGCAGGCGTATGACCAAGCCGAACAACGGTGATGCGGCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGAAATTCCGGGTCAGCGGGGTAGCCt                                                                                        >  3:6684/1‑149 (MQ=255)
                                                          gTACTTGCTGGTGCACGTCCTGGATGTTCGGCAGCTCGTCCAGCGTA‑CGCAGGCGTATGACCAAGCCGATCAACGGTTATGCGGCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCATCTGAAATTCCGGGTCAGCGGGGTAGCCTTc                                                                                      <  7:165300/148‑1 (MQ=255)
                                                          gTACTTGCGGGTGCACGTCCTGGATGTTCGGCAGCTCGTCCAGCGTA‑CGCAGGCGTATGACCAAGCCGAACAACGGTGATGCGGCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGAAATTCCGGGTCAGCGGGGTAGCCTTCg                                                                                     <  2:296567/149‑1 (MQ=255)
                                                              ttGCTGGTGCACGTCCTGGATGTTCGGCAGCTCGTCCAGCGTA‑CGCAGGCGTATGACCAAGCCGAACAACGGTGATGCGGCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGAAATTCCGGGTCAGCGGGGTAGCCTTCGGaaa                                                                                 >  1:69917/1‑149 (MQ=255)
                                                                gCTGGTGCACGTCCTGGATGTTCGGCAGCTCGTCCAGCGTA‑CGCAGGCGTATGACCAAGCCGAACAACGGTGATGCGGCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGAAATTCCGGTTCAGCGGGGTAGCCTTCGGaaa                                                                                 <  7:287934/147‑1 (MQ=255)
                                                                gCTGGTGCACGTCCTGGATGTTCGGCAGCTCGTCAAGCGTA‑CGCAGGCGTATGACCAAGCCGAACAACGGTGATGCGGCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGAAATTCCGGGTCAGCGGGGTAGCCTTCGGAAACg                                                                               >  1:290400/1‑149 (MQ=255)
                                                                 cTGGTGCACGTCCTGGATGTTCGGCAGCTCGTCCAGCGTA‑CGCAGGCGTATGACCAAGCCGAACAACGGTGATGCGGCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGAAATTCCGGGTCAGCGGGGTAGCCTTCGGAAACgg                                                                              <  6:24965/149‑1 (MQ=255)
                                                                 cTGGTGCACGTCCTGGATGTTCGGCAGCTCGTCCAGCGTA‑CGCAGGCGTATGACCAAGCCGAACAACGGTGATGCGGCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGAAATTCCGGGTCAGCGGGGTAGCCTTCGGAAACgg                                                                              <  7:158552/149‑1 (MQ=255)
                                                                  tGGTGCACGTCCTGGATGTTCGGCAGCTCGTCCAGCGTA‑CGCAGGCGTATGACCAAGCCGAACAACGGTGATGCGGCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGAAATTCCGGGTCAGCGGGGTAGCCTTCGGAAACggg                                                                             >  8:106620/1‑149 (MQ=255)
                                                                  tGGTGCACGTCCTGGATGTTCGGCAGCTCGTCCAGCGTA‑CGCAGGCGTATGACCAAGCCGAACAACGGTGATGCGGCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGAAATTCCGGGTCAGCGGGGTAGCCTTCGGAAACggg                                                                             >  1:328890/1‑149 (MQ=255)
                                                                    gTGCACGTCCTGGATGTTCGGCAGCTCGTCCAGCGTA‑CGCAGGCGTATGACCAAGCCGAACAACGGTGATGCGGCATCCAGCATCAGGTTGGCGaa                                                                                                                                   <  1:234612/96‑1 (MQ=255)
                                                                    gTGCACGTCCTGGATGTTCGGCAGCTCGTCCAGCGTA‑CGCAGGCGTATGACCAAGCCGAACAACGGTGATGCGGCATCCAGCATCAGGTTGGCGaa                                                                                                                                   >  2:234612/1‑96 (MQ=255)
                                                                       cACGTCCTGGATGTTCGGCAGCTCGTCCAGCGTA‑CGCAGGCGTATGACCAAGCCGAACAACGGTGATGCGGCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGAAATTCCGGGTCAGCGGGGTAGCCTTCGGAAACGGGctct                                                                         <  2:202115/148‑1 (MQ=255)
                                                                        aCGTCCTGGATGTTCGGCAGCTCGTCCAGCGTA‑CGCAGGCGTATGACCAAGCCGAACAACGGTGATGCGGCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGAAATTCCGGGTCAGCGGGGTAGCCTTCGGAAACGGGCTCTgg                                                                       >  3:259909/1‑149 (MQ=255)
                                                                          gTCCTGGATGTTCGGCAGCTCGTCCAGCGTA‑CGCAGGCGTATGACCAAGCCGAACAACGGTGATGCGGCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGAAATTCCGGGTCAGCGGGGTAGCCTTCGGAAACGGGCTCTGGCt                                                                     >  1:140302/1‑149 (MQ=255)
                                                                              tGGATGTTCGGCAGCTCGTCCAGCGTA‑CGCAGGCGTATGACCAAGCCGAACAACGGTGATGCGGCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGAAATTCCGGGTCAGCGGGGTAGCCTTCGGAAACGGGCTCTGGCTCAgg                                                                 >  1:78491/1‑149 (MQ=255)
                                                                               ggATGTTCGGCAGCTCGTCCAGCGTA‑CGCAGGCGTATGACCAAGCCGAACAACGGTGATGCGGCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGAAATTCCGGGTCAGCGGGGTAGCCTTCGGAAACGGGCTCTGGCTCAGGt                                                                >  2:281391/1‑149 (MQ=255)
                                                                               ggATGTTCGGCAGCTCGTCCAGCGTA‑CGCAGGCGTATGACAAAGCCGAACAACGGTGATGCGGCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGAAATTCCGGGTCAGCGGGGTAGCCTTCGGAAACGGGCTCTGGCTCAGGt                                                                >  3:118802/1‑149 (MQ=255)
                                                                                gATGTTCGGCAGCTCGTCCAGCGTA‑CGCAGGCGTATGACCAAGCCGAACAACGGTGATGCGGCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGAAATTCCGGGTCAGCGGGGTAGCCTTCGGAAACGGGCTCTGGCTCAGGt                                                                <  2:69917/148‑1 (MQ=255)
                                                                                       ggCAGCTCGTCCAGCGTA‑CGCAGGCGTATGACCAAGCCGAACAACGGTGATGCGGCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGAAATTCCGGGTCAGCGGGGTAGCCTTCGGAAACGGGCTCTGGCTCAGGTTTGTCCGt                                                        <  2:140302/149‑1 (MQ=255)
                                                                                            cTCGTCCAGCGTA‑CGCAGGCGTATGACCAAGCCGAACAACGGTGATGCGGCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGAAATTCCGGGTCAGCGGGGTAGCCTTCGGAAACGGGCTCTGGCTCAGGTTTGTCCGTGGATg                                                   >  2:87356/1‑149 (MQ=255)
                                                                                             tCGTCCAGCGTA‑CGCAGGCGTATGACCAAGCCGAACAACGGTGATGCGGCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGAAATTCCGGGTCAGCGGGGTAGCCTTCg                                                                                     >  5:210742/1‑114 (MQ=255)
                                                                                             tCGTCCAGCGTA‑CGCAGGCGTATGACCAAGCCGAACAACGGTGATGCGGCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGAAATTCCGGGTCAGCGGGGTAGCCTTCg                                                                                     <  6:210742/114‑1 (MQ=255)
                                                                                                tCCAGCGTA‑CGCAGGCGTATGACCAAGCCGAACAACGGTGATGCGGCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGAAATTCCGGGTCAGCGGGGTAGCCTTCGGAAACGGGCTCTGGCTCAGGt                                                                <  1:280657/132‑1 (MQ=255)
                                                                                                tCCAGCGTA‑CGCAGGCGTATGACCAAGCCGAACAACGGTGATGCGGCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGATGAAATTCCGGGTCAGCGGGGTAGCCTTCGGAAACGGGCTCTGGCTCAGGt                                                                >  2:280657/1‑132 (MQ=255)
                                                                                                  cAGCGTA‑CGCAGGCGTATGACCAAGCCGAACAACGGTGATGCGGCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGAAATTCCGGGTCAGCGGGGTAGCCTTCGGaaa                                                                                 >  4:315331/1‑113 (MQ=255)
                                                                                                  cAGCGTA‑CGCAGGCGTATGACCAAGCCGAACAACGGTGATGCGGCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGAAATTCCGGGTCAGCGGGGTAGCCTTCGGaaa                                                                                 <  3:315331/113‑1 (MQ=255)
                                                                                                        a‑CGCAGGCGTATGACCAAGCCGAACAACGGTGATGCGGCATCCAGGATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGAAATTCCGGGTCAGCGGGGTAGCCTTCGGAAACGGGCTCTGGCTCAGGTTTGTCCGTGGATGCTTGCCCATTgg                                       >  7:39914/1‑149 (MQ=255)
                                                                                                                       ccAAGCCGAACAACGGTGATGCGGCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGAAATTCCGGGTCAGCGGGGTAGCCTTCGGAAACGGGCTCTGGCTCAGGTTTGTCCGTGGATGCTTGCCCATTGGTTGCTGCAGCCCgg                         >  6:84897/1‑149 (MQ=255)
                                                                                                                       ccAAGCCGAACAACGGTGATGCGGCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGAAATTCCGGGTCAGCGGGGTAGCCTTCGGAAACGGGCTCTGGCTCAGGTTTGTCCGTGGATGCTTGCCCATTGGTTGCTGCAGCCCgg                         >  7:218781/1‑149 (MQ=255)
                                                                                                                        cAAGCCGAACAACGGTGATGCGGCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGAAATTCCGGGTCAGCGGGGTAGCCTTCGGAAACGGGCTCTGGCTCAGGTTTGTCCGTGGATGCTTGCCCATTGGTTGCTGCAGCCCggg                        <  7:106620/149‑1 (MQ=255)
                                                                                                                         aaGCCGAACAACGGTGATGCGGCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGAAATTCCGGGTCAGCGGGGTAGCCTTCGGAAACGGGCTCTGGCTCAGGTTTGTCCGTGGATGCTTGCCCATTGGTTGCTGCAGCCCgggg                       <  4:6684/149‑1 (MQ=255)
                                                                                                                          aGCCGAACAACGGTGATGTGGCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGAAATTCCGGGTCAGCGGGGTAGCCTTCGGAAACGGGCTCTGGCTCAGGTTTGTCCGTGGATGCTTGCCCATTGGTTGCTGCAGCCCggggg                      >  7:115305/1‑149 (MQ=255)
                                                                                                                            ccGAACAACGGTGATGCGGCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGAAATTCCGGGTCAGCGGGGTAGCCTTCGGAAACGGGCTCTGGCTCATGTTTGTCCGTGGATGCTTGCCCATTGGTTGCTGCAGCCCGGGGGGc                    <  2:290400/149‑1 (MQ=255)
                                                                                                                             cGAACAACGGTGATGCGGCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGAAATTCCGGGTCAGCGGGGTAGCCTTCGGAAACGGGCTCTGGCTCAGGTTTGTCCGTGGATGCTTGCCCATTGGTTGCTGCAGCCCGGGGGGCg                   <  5:84897/149‑1 (MQ=255)
                                                                                                                             cGAACAACGGTGATGCGGCATCCAGCATCAGGGTGGCGAAGCCCCGTACTTGCAGCTGAAATTCCGGGTCAGCGGGGTAGCCTTCGGAAACGGGCTCTGGCTCAGGTTTGTCCGTGGATGCTTGCCCATTGGTTGCTGCAGCCCggggg                      <  1:124440/149‑1 (MQ=255)
                                                                                                                              gAACAACGGTGATGCGGCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGAAATTCCGGGTCAGCGGGGTAGCCTTCGGAAACGGGCTCTGGCTCAGGTTTGTCCGTGGATGCTTGCCCATTGGTTGCTGCAGCCCGGGGGGCgg                  <  6:278768/149‑1 (MQ=255)
                                                                                                                              gAACAACGGTGATGCGGCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGAAATTCCGGGTCAGCGGGGTAGCCTTCGGAAACGGGCTCTGGCTCAGGTTTGTCCGTGGATGCTTGCCCATTGGTTGCTGCAGCCCGGGGGGCgg                  <  4:259909/149‑1 (MQ=255)
                                                                                                                               aacaacGGTGATGCGGCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGAAATTCCGGGTCAGCGGGGTAGCCTTCGGAAACGGGCTCTGGCTCAGGTTTGTCCGTGGATGGTTGCCCATTGGTTGCTGCAGCCCGGGGGGCggg                 <  1:281391/149‑1 (MQ=255)
                                                                                                                               aacaacGGTGATGCGGCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGAAATTCCGGGTCAGCGGGGTAGCCTTCGGAAACGGGCTCTGGCTCAGGTTTGTCCGTGGATGCTTGCCCATTGGTTGCTGCAGCCCGGGGGGCggg                 <  7:209617/149‑1 (MQ=255)
                                                                                                                                acaacGGTGATGCGGCATCCAGCATCAGGTTGGCGAAGCCGT‑‑‑CCCGCAGCTGAAATTCCGGGTCAGCGGGGTAGCCTTCGGAAACGGGCTCTGGCTCAGGTTTGTCCGTGGATGCTTGCCCATTgg                                       <  4:236631/126‑1 (MQ=255)
                                                                                                                                acaacGGTGATGCGGCATCCAGCATCAGGTTGGCGAAGCCGT‑‑‑CCCGCAGCTGAAATTCCGGGTCAGCGGGGTAGCCTTCGGAAACGGGCTCTGGCTCAGGTTTGTCCGTGGATGCTTGCCCATTgg                                       >  3:236631/1‑126 (MQ=255)
                                                                                                                                  caacctGATGCTGGATGCCGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGAAATTCCGGGTCAGCGGGGTAGCCTTCGGAAACGGGCTCTGGCTCAGGTTTGTCCGTGGATGCTTGCCCATTGGTtgctg                                 >  3:366795/6‑130 (MQ=255)
                                                                                                                                  caacctGATGCTGGATGCCGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGAAATTCCGGGTCAGCGGGGTAGCCTTCGGAAACGGGCTCTGGCTCAGGTTTGTCCGTGGATGCTTGCCCATTGGTtgctg                                 <  4:366795/125‑1 (MQ=255)
                                                                                                                                     ggTGATGCGGCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGAAATTCCGGGTCAGCGGGGTAGCCTTCGGAAACGGGCTCTGGCTCAGGTTTGTCCGTGGATGCTTGCCCATTGGTTGCTGCAGCCCGGGGGGCGGGTACTgg           <  2:78491/149‑1 (MQ=255)
                                                                                                                                     ggTGATGCGGCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGAAATTCCGGGTCAGCGGGGTAGCCTTCGGAAACGGGCTCTGGCTCAGGTTTGTCCGTGGATGCTTGCCCATTGGTTGCTGCAGCCCGGGGGGCGGGTACTgg           <  2:328890/149‑1 (MQ=255)
                                                                                                                                      gTGATGCGGCATCCAGGATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGAAATTCCGGGTCAGCGGGGTAGCCTTCGGAAACGGGCTCTGGCTCAGGTTTGTCCGTGGATGCTTGCCCATTGGTTGCTGCAGCCCGGGGGGCGGGTACTGGc          <  8:39914/149‑1 (MQ=255)
                                                                                                                                             ggCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGAAATTCCGGGTCAGCGGGGTAGc                                                                                          >  7:203229/1‑62 (MQ=255)
                                                                                                                                             ggCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGAAATTCCGGGTCAGCGGGGTAGc                                                                                          <  8:203229/62‑1 (MQ=255)
                                                                                                                                             ggCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGAAATTCCGGGTCAGCGGGGTAGCCTTCGGAAACGGGCTCTGGCTCAGGTTTGTCCGTGGATGCTTGCCCATTGGTTGCTGCAGCCCGGGGGGCGGGTACTGGCTGTTGaa   <  8:218781/149‑1 (MQ=255)
                                                                                                                                              gCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGAAATTCCGGGTCAGCGCGGTAGCCTTCGGAAACGGGCTCTGGCTCAGGTTTGTCCGTGGATGCTTGCCCATTGGTTGCTGCAGCCGGGGGGGCGGGTACTGGCTGTTGAAg  >  6:37793/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TGAGCCGGCTCGTAACCGTGCTGGCGCATTTCTTCCTGAATGTTGTCGATCTGGCTGCGTACTTGCTGGTGCACGTCCTGGATGTTCGGCAGCTCGTCCAGCGTA‑CGCAGGCGTATGACCAAGCCGAACAACGGTGATGCGGCATCCAGCATCAGGTTGGCGAAGCCGC‑‑‑CCCGCAGCTGAAATTCCGGGTCAGCGGGGTAGCCTTCGGAAACGGGCTCTGGCTCAGGTTTGTCCGTGGATGCTTGCCCATTGGTTGCTGCAGCCCGGGGGGCGGGTACTGGCTGTTGAAG  >  NC_002947/2992717‑2993006

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