Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,437,398 A→T 6.3% V54E (GTG→GAG)  PP_4775 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,437,3980AT 6.3% 186.5 / 3.0 64V54E (GTG→GAG) PP_4775hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (27/33);  new base T (2/2);  total (29/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 = 4.35e-01

CCAGCCATTGCCATGGCCATCCGGGCGCCCCTGCCAGTCGCCACCGGGCCGTTGCCCCGGGCCACCGGTTTGCCAGTGGTGATCATCACCGTGGCGGCGGGGGATGTCCTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAACCAACACCTGCCAACTGC  >  NC_002947/5437253‑5437541
                                                                                                                                                 |                                                                                                                                               
ccAGCCATTGCCATGCCCATCCGGGCGCCCCTGCCAGTCGCCACCGGGCCGTTGCCCCGGGCCACCGGTTTGCCAGTGGTGATCATCACCGTGGCGGCGGGGGATGTCCTGGTAATACCCCTTGCGCGGCGGTTGGGTCTGGCGCACtt                                                                                                                                              <  6:231188/149‑1 (MQ=255)
 cAGCCATTGCCATGGCCATCCGGGCGCCCCTGCCAGTCGCCACCGGGCCGTTGCCCCGGGCCACCGGTTTGCCAGTGGTGATCATCACCGTGGCGGCGGGGGATGTCCTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTtc                                                                                                                                             <  4:4145/149‑1 (MQ=255)
  aGCCATTGCCATGGCCATCCGGGCGCCCCTGCGAGTCGCCACCGGGCCGTTGCCCCGGGCCACCGGTTTGCCAGTGGTTATCATCACCGTGGCGGCGGGGGATGTCCTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGCCGAACTtcg                                                                                                                                            >  2:317469/1‑149 (MQ=255)
          ccATGGCCATCCGGGCGCCCCTGCCAGTCGCCACCGGGCCGTTGCCCCGGGCCACCGGTTTGCCAGTGGTGATCATCACCGTGGCGGCGGGGGATGTCCTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTcgcgc                                                                                                                                      >  3:362261/1‑147 (MQ=255)
            aTGGCCATCCGGGCGCCCCTGCCAGTCGCCACCGGGCCGTTGCCCCGGGCCACCGGTTTGCCAGTGGTGATCATCACCGTGGCGGCGGGGGATGTCCTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGt                                                                                                                                  >  3:156488/1‑149 (MQ=255)
                tcATCCGGGCGCCCCTGCCAGTCGCCACCGGGCCGTTGCCCCGGGCCACCGGTTTGCCAGTGGTGATCATCACCGTGGCGGCGGGGGATGTCCTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTgag                                                                                                                              >  1:142085/2‑149 (MQ=255)
                ccATCCGGGCGCCCCTGCCAGTCGCCACCGGGCCGTTGCCCCGGGCCACCGGTTTGCCAGTGGTGATCATCACCGTGGCGGCGGGGGATGTCCTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTcgcgc                                                                                                                                      >  6:90034/1‑141 (MQ=255)
                ccATCCGGGCGCCCCTGCCAGTCGCCACCGGGCCGTTGCCCCGGGCCACCGGTTTGCCAGTGGTGATCATCACCGTGGCGGCGGGGGATGTCCTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTcgcgc                                                                                                                                      <  5:90034/141‑1 (MQ=255)
                 cATCCGGGCGCCCCTGCCAGTCGCCACCGGGCCGTTGCCCCGGGCCACCGGTTTGCCAGTGGTGATCATCACCGTGGCGGCGGGGGATGTCCTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTgagc                                                                                                                             >  4:198202/1‑149 (MQ=255)
                            cccTGCCAGTCGCCACCGGGCCGTTGCCCCGGGCCACCGGTTTGCCAGTGGTGATCATCACCGTGGCGGCGGGGGATGTCCTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCgg                                                                                                                  >  8:155814/1‑149 (MQ=255)
                            cccTGCCAGTCGCCACCGGGCCGTTGCCCCGGGCCACCGGTTTGCCAGTGGTGATCATCACCGTGGCGGCGGGGGATGTCCTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCgg                                                                                                                  >  6:254318/1‑149 (MQ=255)
                            cccTGCCAGTCGCCACCGGGCCGTTGCCCCGGGCCACAGGTTTGCCAGTGGTGATCATCACCGTGGCGGCGGGGGATGTCCTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCgg                                                                                                                  >  7:171333/1‑149 (MQ=255)
                             ccTGCCAGTCGCCACCGGGCCGTTGCCCCGGGCCACCGGTTTGCCAGTGGTGATCATCACCGTGGCGGCGGGGGATGTCCTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGt                                                                                                                 <  4:362261/149‑1 (MQ=255)
                                gCCAGTCGCCACCGGGCCGTTGCCCCGGGCCACCGGTTTGCCAGTGGTGATCATCACCGTGGCGGCGGGGGATGTCCTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGTGAGTTGAGCGAGCGCCCCGGTCCa                                                                                                              >  3:71853/1‑149 (MQ=255)
                                 ccAGTCGCCACCGGGCCGTTGCCCCGGGCCACCGGTTTGCCAGTGGTGATCATCACCGTGGCGGCGGGGGATGTCCTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTcgc                                                                                                                                        <  7:331130/122‑1 (MQ=255)
                                 ccAGTCGCCACCGGGCCGTTGCCCCGGGCCACCGGTTTGCCAGTGGTGATCATCACCGTGGCGGCGGGGGATGTCCTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTcgc                                                                                                                                        >  8:331130/1‑122 (MQ=255)
                                 ccAGTCGCCACCGGGCCGTTGCCCCGGGCCACCGGTTTGCCAGTGGTGATCATCACCGTGGCGGCGGGGGATGTCCTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCa                                                                                                              <  5:254318/148‑1 (MQ=255)
                                  cAGTCGCCACCGGGCCGTTGCCCCGGGCCACCGGTTTGCCAGTGGTGATCATCACCGTGGCGGCGGGGGATGTCCTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCAcc                                                                                                            <  2:51704/149‑1 (MQ=255)
                                   aGTCGCCACCGGGCCGTTGCCCCGGGCCACCGGTTTGCCAGTGGTGATCATCACCGTGGCGGCGGGGGATGTCCTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCAccc                                                                                                           <  2:142085/149‑1 (MQ=255)
                                        ccACCGGGCCGTTGCCCCGGGCCACCGGTTTGCCAGTGGTGATCATCACCGTGGCGGCGGGGGATGTCCTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCAcc                                                                                                      >  1:71981/1‑149 (MQ=255)
                                        ccACCGGGCCGTTGCCCCGGGCCACCGGTTTGCCAGTGGTGATCATCACCGTGGCGGCGGGGGATGTCCTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCAcc                                                                                                      <  1:356325/149‑1 (MQ=255)
                                          aCCGGGCCGTTGCCCCGGGCCACCGGTTTGCCAGTGGTGATCATCACCGTGGCGGCGGGGGATGTCCTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGTGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCtt                                                                                                    <  4:71853/149‑1 (MQ=255)
                                           ccGGGCCGTTGCCCCGGGCCACCGGTTTGCCAGTGGTGATCATCACCGTGGCGGCGGGGGATGTCCTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTg                                                                                                   <  7:310341/149‑1 (MQ=255)
                                             gggCCGTTGCCCCGGGCCACCGGTTTGCCAGTGGTGATCATCACCGTGGCGGCGGGGGATGTCCTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCa                                                                                                 <  2:158513/149‑1 (MQ=255)
                                              ggCCGTTGCCCCGGGCCACCGGTTTGCCAGTGGTGATCATCACCGTGGCGGCGGGGGATGTCCTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCAc                                                                                                <  6:27961/149‑1 (MQ=255)
                                              ggCCGTTGCCCCGGGCCACCGGTTTGCCAGTGGTGATCATCACCGTGGCGGCGGGGGATGTCCTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCAc                                                                                                <  1:151535/149‑1 (MQ=255)
                                                 cGTTGCCCCGGGCCACCGGTTTGCCAGTGGTGATCATCACCGTGGCGGCGGGGGATGTCCTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCTCTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACctg                                                                                             >  3:271249/1‑149 (MQ=255)
                                                 cGTTGCCCCGGGCCACCGGTTTGCCAGTGGTGATCATCACCGTGGCGGCGGGGGATGTCCTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCTCTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACctg                                                                                             >  3:266259/1‑149 (MQ=255)
                                                    cggcccGGGCCACCGGTTTGCCAGTGGTGATCATCACCGTGGCGGCGGGGGATGTCCTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACctgctg                                                                                          <  5:147196/146‑1 (MQ=255)
                                                       cccGGGCCACCGGTTTGCCAGTGGTGATCATCACCGTGGCGGCGGGGGATGTCCTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGcc                                                                                       <  6:14254/149‑1 (MQ=255)
                                                       cccGGGCCACCGGTTTGCCAGTGGTGATCATCACCGTGGCGGCGGGGGATGTCCTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGc                                                                                        >  6:43381/1‑148 (MQ=255)
                                                                   gTTTGCCAGTGGTGATCATCACCGTGGCGGCGGGGGATGTCCTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCTCTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTg                                                                           <  4:266259/149‑1 (MQ=255)
                                                                   gTTTGCCAGTGGTGATCATCACCGTGGCGGCGGGGGATGTCCTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCTCTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTg                                                                           <  4:271249/149‑1 (MQ=255)
                                                                    tttGCCAGTGGTGATCATCACCGTGGCGGCGGGGGATGTCCTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGc                                                                          <  5:43381/149‑1 (MQ=255)
                                                                       gCCAGTGGTGATCATCACCGTGGCGGCGGGGGATGTCCTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGccc                                                                        <  3:198202/148‑1 (MQ=255)
                                                                              gtgATCATCACCGTGGCGGCGGGGGATGTCCTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACctgc                                                                                            <  8:5485/121‑1 (MQ=255)
                                                                              gtgATCATCACCGTGGCGGCGGGGGATGTCCTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACctgc                                                                                            >  7:5485/1‑121 (MQ=255)
                                                                                     tcaCCGTGGCGGCGGGGGATGTCCTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCtt                                                          >  2:286310/1‑148 (MQ=255)
                                                                                             gcggcgGGGGATGTCCTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTAGCCGGGCCCTTCAGCCTAgg                                                 >  7:217668/1‑149 (MQ=255)
                                                                                              cggcggGGGATGTCCTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGc                                                <  2:71981/149‑1 (MQ=255)
                                                                                                   ggggATGTCCTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCaa                                           <  8:171333/149‑1 (MQ=255)
                                                                                                     ggATGTCCTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGTGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGc                                                                          >  2:157088/1‑116 (MQ=255)
                                                                                                     ggATGTCCTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGTGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGc                                                                          <  1:157088/116‑1 (MQ=255)
                                                                                                     ggATGTCCTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGACGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAgg                                         >  5:228438/1‑149 (MQ=255)
                                                                                                        tGTCCTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTg                                      <  1:286310/149‑1 (MQ=255)
                                                                                                            cTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGTGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCa                                                                                <  1:77584/103‑1 (MQ=255)
                                                                                                            cTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCa                                                                                >  2:77584/1‑103 (MQ=255)
                                                                                                            cTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAgg                                                 >  6:276895/1‑134 (MQ=255)
                                                                                                            cTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAgg                                                 <  5:276895/134‑1 (MQ=255)
                                                                                                            cTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATg                                   <  3:208678/148‑1 (MQ=255)
                                                                                                            cTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGc                                  >  3:221349/1‑149 (MQ=255)
                                                                                                             tGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGCGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGCCCCTTCAGCCAAGGCCCCAAGGCTGATggt                                 >  7:173637/1‑147 (MQ=255)
                                                                                                               gTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCt                                 <  7:210563/147‑1 (MQ=255)
                                                                                                                tAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCaa                              <  7:155814/149‑1 (MQ=255)
                                                                                                                  aTACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGCCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCaaa                             <  6:228438/148‑1 (MQ=255)
                                                                                                                     ccccTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCa                         <  4:156488/149‑1 (MQ=255)
                                                                                                                           gggcggcgGTTTGGTCTGGAGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCACTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAAc                   <  1:317469/147‑1 (MQ=255)
                                                                                                                             gcggcgGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAAcc                  >  7:165057/1‑148 (MQ=255)
                                                                                                                                    ttGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAACCAACACCTg          >  5:232117/1‑149 (MQ=255)
                                                                                                                                      gggTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAACCAACACCTGc         >  6:137231/1‑148 (MQ=255)
                                                                                                                                        gTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAACCAACACCTGCCaa      >  2:265174/1‑149 (MQ=255)
                                                                                                                                          cTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAACCAacac             >  3:1895/1‑140 (MQ=255)
                                                                                                                                          cTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAACCAacac             <  4:1895/140‑1 (MQ=255)
                                                                                                                                             gcgcACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAACCAACACCTGCCAACTGc  <  5:137231/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
CCAGCCATTGCCATGGCCATCCGGGCGCCCCTGCCAGTCGCCACCGGGCCGTTGCCCCGGGCCACCGGTTTGCCAGTGGTGATCATCACCGTGGCGGCGGGGGATGTCCTGGTAATACCCCTGGCGCGGCGGTTGGGTCTGGCGCACTTCGTCGCGCGAGTTGAGCGAGCGCCCCGGTCCACCCTCACCTTGCACCTGCTGGCCACGTTCACCGTGCCCCTCGCCGGGCCCTTCAGCCAAGGCCCCAAGGCTGATGCTCAAACCCAGTAAACCAACACCTGCCAACTGC  >  NC_002947/5437253‑5437541

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: ‑