Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,904,885 G→T 5.7% intergenic (‑52/+183) PP_4312 ← / ← PP_4313 hypothetical protein/peptidylprolyl isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,904,8850GT 5.7% 256.1 / 3.1 79intergenic (‑52/+183)PP_4312/PP_4313hypothetical protein/peptidylprolyl isomerase
Reads supporting (aligned to +/- strand):  ref base G (38/36);  new base T (2/3);  total (40/39)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTCGACAGCCTCCAGCGCCAACGGCGCGGCGATGGCAGTGAGCACCAGGCTCAGCAGCAAGTGACGACGACGGCGAAAGAATGTTGGCAACTTCATCAAGGCAGGCTCGTGTGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATT  >  NC_002947/4904737‑4905029
                                                                                                                                                      |                                                                                                                                                     
tttCGACAGCCTCCAGCGCCAACGGCGCGGCGATGGCAGTGAGCACCAGGCTCAGCAGCAAGTGACGACGACGGCGAAAGAATGTTGGCAACTTCATCAAGGCAGGCTCGTGTGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCg                                                                                                                                                       >  4:218979/1‑149 (MQ=255)
tttCGACAGCCTCCAGCGCCAACGGCGCGGCGATGGCAGTGAGCACCAGGCTCAGCAGCAAGTGACGACGACGGCGAAAGAATGTTGGCAACTTCATCAAGGCAGGCTCGTGTGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCg                                                                                                                                                       >  1:346614/1‑149 (MQ=255)
  tCGACAGCCTCCAGCGCCAACGGCGCGGCGATGGCAGTGAGCACCAGGCTCAGCAGCAAGTGACGACGACGGCGAAAGAATGTTGGCAACTTCATCAAGGCAGGCTCGTGTGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCGGc                                                                                                                                                     >  3:179055/1‑149 (MQ=255)
   cGACAGCCTCCAGCGCCAACGGCGCGGCGATGGCAGTGAGCACCAGGCTCAGCAGCAAGTGACGACGACGGCGAAAGAATGTTGGCAACTTCATCAAGGCAGGCTCGTGTGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCGGc                                                                                                                                                     <  5:169255/148‑1 (MQ=255)
    gACAGCCTCCAGCGCCAACGGCGCGGCGATGGCAGTGAGCACCAGGCTCAGCAGCAAGTGACGACGACGGCGAAAGAATGTTGGCAACTTCATCAAGGCAGGCTCGTGTGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCGGc                                                                                                                                                     >  4:48609/1‑147 (MQ=255)
         ccTCCAGCGCCAACGGCGCGGCGATGGCAGTGAGCACCAGGCTCAGCAGCAAGTGACGACGACGGCGAAAGAATGTTGGCAACTTCATCAAGGCAGGCTCGTGTGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGc                                                                                                                                              >  5:186951/1‑149 (MQ=255)
         ccTCCAGCGCCAACGGCGCGGCGATGGCAGTGAGCACCAGGCTCAGCAGCAAGTGACGACGACGGCGAAAGAATGTTGGCAACTTCATCAAGGCAGGCTCGTGTGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGc                                                                                                                                              <  2:120301/149‑1 (MQ=255)
             cAGCGCCAACGGCGCGGCGATGGCAGTGAGCACCAGGCTCAGCAGCAAGTGACGACGACGGCGAAAGAATGTTGGCAACTTCATCAAGGCAGGCTCGTGTGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAgg                                                                                                                                           <  3:276023/148‑1 (MQ=255)
               ctgcCATCGCCGCGCCGTTGGCAGTGAGCACCAGGCTCAGCAGCAAGTGACGACGACGGCGAAAGAATGTTGGCAACTTCATCAAGGCAGGCTCGTGTGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGgcg                                                                                                                                        >  5:184887/3‑149 (MQ=255)
                 gcCAACGGCGCGGCGATGGCAGTGAGCACCAGGCTCAGCAGCAAGTGACGACGACGGCGAAAGAATGTTGGCAACTTCATCAAGGCAGGCTCGTGTGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCTGCAACTAGCCGGgcgat                                                                                                                                      >  4:207077/1‑147 (MQ=255)
                 gcCAACGGCGCGGCGATGGCAGTGAGCACCAGGCTCAGCAGCAAGTGACGACGACGGCGAAAGAATGTTGGCAACTTCATCAAGGCAGGCTCGTGTGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCTGCAACTAGCCGGgc                                                                                                                                         >  1:318733/1‑146 (MQ=255)
                 gcCAACGGCGCGGCGATGGCAGTGAGCACCAGGCTCAGCAGCAAGTGACGACGACGGCGAAAGAATGTTGGCAACTTCATCAAGGCAGGCTCGTGTGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCTGCAACTAGCCGGgc                                                                                                                                         <  2:318733/146‑1 (MQ=255)
                 gcCAACGGCGCGGCGATGGCAGTGAGCACCAGGCCCAGCAGCAAGTGACGACGACGGCGAAAGAATGTTGGCAACTTCATCAAGGCAGGCTCGTGTGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGcc                                                                                                                                      >  2:68873/1‑149 (MQ=255)
                  ccAACGGCGCGGCGATGGCAGTGAGCACCAGGCTCAGCAGCAAGTGACGACGACGGCGAAAGAAAGTTGGCAACTTCATCAAGGCAGGCTCGTGTGGCAGAGAATATCGGGTGTTACCCGCCACATCGCCCTGCAACTAGCCGGgcgat                                                                                                                                      <  1:61352/149‑3 (MQ=255)
                   cAACGGCGCGGCGATGGCAGTGAGCACCAGGCTCAGCAGCAAGTGACGACGACGGCGAAAGAATGTTGGCAACTTCATCAAGGCAGGCTCGTGTGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCTGCAACTAGCCGGgcgatgt                                                                                                                                    <  3:207077/149‑5 (MQ=255)
                   cAACGGCGCGGCGATGGCAGTGAGCACCAGGCTCAGCAGCAAGTGACGACGACGGCGAAAGAATGTTGGCAACTTCATCAAGGCAGGCTCGTGTGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCg                                                                                                                                    >  4:113171/1‑149 (MQ=255)
                     aCGGCGCGGCGATGGCAGTGAGCACCAGGCCCAGCAGCAAGTGACGACGACGGCGAAAGAATGTTGGCAACTTCATCAAGGCAGGCTCGTGTGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTg                                                                                                                                  <  1:68873/149‑1 (MQ=255)
                       ggcgcggcgATGGCAGTGAGCACCAGGCTCAGCAGCAAGTGACGACGACGGCGAAAGAATGTTGGCAACTTCATCAAGGCAGGCTCGTGTGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCa                                                                                                                                <  3:48609/149‑1 (MQ=255)
                              cgATGGCAGTGAGCACCAGGCTCAGCAGCAAGTGACGACGACGGCGAAAGAATGTTGGCAACTTCATCAAGGCAGGCTCGTGTGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATg                                                                                                                         <  4:179055/149‑1 (MQ=255)
                                       tGAGCACCAGGCTCAGCAGCAAGTGACGACGACGGCGAAAGAATGTTGGCAACTTCATCAAGGCAGGCTCGTGTGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTg                                                                                                                <  6:184887/149‑1 (MQ=255)
                                         aGCACCAGGCTCAGCAGCAAGTGACGACGACGGCGAAAGAATGTTGGCAACTTCATCAAGGCAGGCTCGTGTGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTgg                                                                                                               <  6:203862/148‑1 (MQ=255)
                                          gCACCAGGCTCAGCAGCAAGTGACGACGACGGCGAAAGAATGTTGGCAACTTCATCAAGGCAGGCTCGTGTGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGc                                                                                                              >  5:294184/1‑148 (MQ=255)
                                                   tCAGCAGCAAGTGACGACGACGGCGAAAGAATGTTGGCAACTTCATCAAGGCAGGCTCGTGTGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGtt                                                                                                                                                >  3:21365/1‑105 (MQ=255)
                                                   tCAGCAGCAAGTGACGACGACGGCGAAAGAATGTTGGCAACTTCATCAAGGCAGGCTCGTGTGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGtt                                                                                                                                                <  4:21365/105‑1 (MQ=255)
                                                    cagcagCAAGTGACGACGACGGCGAAAGAATGTTGGCAACTTCATCAAGGCAGGCTCGTGTGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑Agg                                                                                                 <  2:124345/149‑1 (MQ=255)
                                                          caAGTGACGACGACGGCGAAAGAATGTTGGCAACTTCATCAAGGCAGGCTCGTGTGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCaga                                                                                           <  6:186951/149‑1 (MQ=255)
                                                                 cgacgacgGCGAAAGAATGTTTGCAACTTCATCAAGGCAGGCTCGTGTGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGcccc                                                                                    >  4:6808/1‑149 (MQ=255)
                                                                  gacgacgGCGAAAGAATGTTGGCAACTTCATCAAGGCAGGCTCGTGTGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGcccc                                                                                    >  2:354553/1‑148 (MQ=255)
                                                                  gacgacgGCGAAAGAATGTTGGCAACTTCATCAAGGCAGGCTCGTGTGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGcccc                                                                                    >  3:94324/1‑148 (MQ=255)
                                                                     gacgGCGAAAGAATGTTGGCAACTTCATCAAGGCAGGCTCGTGTGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAAc                                                                             <  3:218979/149‑1 (MQ=255)
                                                                     gacgGCGAAAGAATGTTGGCAACTTCATCAAGGCAGGCTCGTGTGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAAc                                                                             <  7:269635/149‑1 (MQ=255)
                                                                        ggCGAAAGAATGTTGGCAACTTCATCAAGGCAGGCTCGTGTGGCAGAGAATATCGGGTTTACCAGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTGAGAGTGCCCC‑‑‑GAACACg                                                                          >  8:329365/1‑149 (MQ=255)
                                                                          cGAAAGAATGTTTGCAACTTCATCAAGGCAGGCTCGTGTGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGtt                                                                        <  3:6808/149‑1 (MQ=255)
                                                                          cGAAAGAATGTTGGCAACTTCATCAAGGCAGGCTCGTGTGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGtt                                                                        <  6:200314/149‑1 (MQ=255)
                                                                           gAAAGAATGTTGGCAACTTCATCAAGGCAGGGCTCGTGTGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGAAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGTTCAGAGTGCCCC‑‑‑GAACCCGtt                                                                        >  1:215388/1‑149 (MQ=255)
                                                                           gAAAGAATGTTGGCAACTTCATCAAGGCAGGCTCGTGTGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGttt                                                                       >  2:299767/1‑149 (MQ=255)
                                                                             aaGAATGTTGGCAACTTCATCAAGGCAGGCTCGTGTGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTgg                                                                     <  1:354553/149‑1 (MQ=255)
                                                                              aGAATGTTGGCAACTTCATCAAGGCAGGCTCGTGTGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGt                                                                    <  6:294184/149‑1 (MQ=255)
                                                                                      ggCAACTTCATCAAGGCAGGCTCGTGTGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCt                                                            >  3:98668/1‑149 (MQ=255)
                                                                                      ggCAACTTCATCAAGGCAGGCTCGTGTGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGCGCCCC‑‑‑GAACCCGTTTGGTtcct                                                                <  2:279876/145‑1 (MQ=255)
                                                                                      ggCAACTTCATCAAGGCAGGCTCGTGTGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGCGCCCC‑‑‑GAACCCGTTTGGTtcct                                                                >  1:279876/1‑145 (MQ=255)
                                                                                         aaCTTCATCAAGGCAGGCTCGTGTGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTc                                                                                              <  3:274953/115‑1 (MQ=255)
                                                                                         aaCTTCATCAAGGCAGGCTCGTGTGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTc                                                                                              >  4:274953/1‑115 (MQ=255)
                                                                                          aCTTCATCAAGGCAGGCTCGTGTGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAg                                                        >  2:285321/1‑149 (MQ=255)
                                                                                               atcaAGGCAGGCTCGTGTGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGaa                                                   >  1:156949/1‑149 (MQ=255)
                                                                                                tcaAGGCAGGCTCGTGTGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAAt                                                  <  2:215388/149‑1 (MQ=255)
                                                                                                  aaGGCAGGCTCGTGTGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGc                                                <  3:113171/149‑1 (MQ=255)
                                                                                                   aggcaggcTCGTGTGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGcc                                               <  4:98668/149‑1 (MQ=255)
                                                                                                   aggcaggcTCGTGTGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGcc                                               >  7:246954/1‑149 (MQ=255)
                                                                                                    ggcaggcTCGTGTGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCg                                              >  5:156876/1‑149 (MQ=255)
                                                                                                    ggcaggcTCGTGTGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCg                                              >  2:81166/1‑149 (MQ=255)
                                                                                                       aggcTCGTGTGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGtta                                           >  5:302030/1‑148 (MQ=255)
                                                                                                             cggggggCAGAGAATATCGGTGTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAg                                      <  4:203576/144‑1 (MQ=255)
                                                                                                              gtgtGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAgg                                     >  1:218316/1‑149 (MQ=255)
                                                                                                              gtgtGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAg                                      <  6:156876/148‑1 (MQ=255)
                                                                                                                 tgGCAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGc                                    >  7:97753/1‑147 (MQ=255)
                                                                                                                   gCAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAg                                 <  6:302030/148‑1 (MQ=255)
                                                                                                                    cAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGc                               >  8:221186/1‑149 (MQ=255)
                                                                                                                     agagAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGt                                                                                                                              <  5:40376/58‑1 (MQ=255)
                                                                                                                     agagAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGt                                                                                                                              >  6:40376/1‑58 (MQ=255)
                                                                                                                      gagaATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGcc                              >  2:231455/1‑148 (MQ=255)
                                                                                                                      gagaATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGcc                              >  6:93793/1‑148 (MQ=255)
                                                                                                                        gaATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAg                           <  8:97753/149‑1 (MQ=255)
                                                                                                                         aaTATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAgg                          >  5:338336/1‑149 (MQ=255)
                                                                                                                          atatCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCCTTGCCTAgg                          <  1:118651/148‑1 (MQ=255)
                                                                                                                           tatCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGtt                        >  2:277399/1‑149 (MQ=255)
                                                                                                                           tatCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGt                         >  6:346006/1‑148 (MQ=255)
                                                                                                                           tatCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGt                         >  1:309361/1‑148 (MQ=255)
                                                                                                                             tCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTTCAGTTAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTa                       <  8:246954/148‑1 (MQ=255)
                                                                                                                             tCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGt                                                                    >  3:221840/1‑103 (MQ=255)
                                                                                                                             tCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGt                                                                    <  4:221840/103‑1 (MQ=255)
                                                                                                                               gggTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATTCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCCTTGCCTAgggaacg                     <  7:329365/148‑6 (MQ=255)
                                                                                                                                ggTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAgggg                   <  2:346614/149‑1 (MQ=255)
                                                                                                                                 gTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGc                  >  7:113355/1‑149 (MQ=255)
                                                                                                                                      aCCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTgg               <  8:113355/148‑1 (MQ=255)
                                                                                                                                         cGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGgtgt           >  6:369211/1‑149 (MQ=255)
                                                                                                                                             acaTCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTgcg       >  6:321394/1‑149 (MQ=255)
                                                                                                                                               aTCGCCCGGCTAGTTTCAGGGCGCCCGTGCAGGGCATGCCATCCTTGGCAAAGCGTGTCTGGGTCAGAGTGCCCCCCAGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGtt          <  1:299767/149‑1 (MQ=255)
                                                                                                                                                  gCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTg         >  7:101157/1‑142 (MQ=255)
                                                                                                                                                  gCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTg         <  8:101157/142‑1 (MQ=255)
                                                                                                                                                  gCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAAtt  <  4:94324/149‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
TTTCGACAGCCTCCAGCGCCAACGGCGCGGCGATGGCAGTGAGCACCAGGCTCAGCAGCAAGTGACGACGACGGCGAAAGAATGTTGGCAACTTCATCAAGGCAGGCTCGTGTGGCAGAGAATATCGGGTTTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATT  >  NC_002947/4904737‑4905029

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