Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,400,104 C→A 6.9% A156S (GCC→TCC)  PP_3006 ← ECF family RNA polymerase sigma‑70 factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,400,1040CA 6.9% 183.3 / 3.9 58A156S (GCC→TCC) PP_3006ECF family RNA polymerase sigma‑70 factor
Reads supporting (aligned to +/- strand):  ref base C (28/26);  new base A (2/2);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.60e-01

CGACAAGGTTCCCAGCACGTATTCCCCGGCGAGGCTGTCCAGTTCTTCGGCAGAATCGGGCTTCATTCCAGACACTCCCGTAGCGCTTTGAGGCTGCGTTTGATCCAGGCCTTGACGGTGCCCAGTGGCGCACCCAATAAACCGGCTATCTCGGCGTGGCTGCAGCCGTCCACGTAGGCGTGCAGGATGCAGCGTTGCGGTGCGGGCTGCAGTTGTTCCAGGCAACGGTAAATCTTTGCTGAGCCGGACCAGAGGTCAAGGTCATCAACCCTCCCGTGCAATGTTGATGCA  >  NC_002947/3399959‑3400249
                                                                                                                                                 |                                                                                                                                                 
cGACAAGGTTCCCAGCACGTATTCCCCGGCGAGGCTGTCCAGTTCTTCGGCAGAATCGGGCTTCATTCCAGACACTCCCGTAGCGCTTTGAGGCTGCGTTTGATCCAGGCCTTGACGGTGCCCAGTGGCGCACCCAATAAACCGGCt                                                                                                                                                  <  8:288107/147‑1 (MQ=255)
  aCAAGGTTCCCAGCACGTATTCCCCGGCGAGGCTGTCCAGTTCTTCGGCAGAATCGGGCTTCATTCCAGACACTCCCGTAGCGCTTTGAGGCTGCGTTTGATCCAGGCCTTGACGGTGCCCAGTGGCGCACCCAATAAACCGGCTAtct                                                                                                                                              >  2:9151/1‑149 (MQ=255)
        ttCCCAGCACGTATTCCCCGGCGAGGCTGTCCAGTTCTTCGGCAGAATCGGGCTTCATTCCAGACACTCCCGTAGCGCTTTGAGGCTGCGTTTGATCCAGGCCTTGACGGTGCCCAGTGGCGCACCCAATAAACCGGCTATCTCGGCGt                                                                                                                                        >  7:61252/1‑149 (MQ=255)
           ccAGCACGTATTCCCCGGCGAGGCTGTCCAGTTCTTCGGCAGAATCGGGCTTCATTCCAGACACTCCCGTAGCGCTTTGAGGCTGCGTTTGATCCAGGCCTTGACGGTGCCCAGTGGCGCACCCAATAAACCGGCTATCTCGGCGTgg                                                                                                                                      <  5:123678/148‑1 (MQ=255)
            cAGCACGTATTCCCCGGCGAGGCTGTCCAGTTCTTCGGCAGAATCGGGCTTCATTCCAGACACTCCCGTAGCGCTTTGAGGCTGCGTTTGATCCAGGCCTTGACGGTGCCCAGTGGCGCACCCAATAAACCGGCTATCTCGGCGTGGCt                                                                                                                                    >  7:102511/1‑149 (MQ=255)
             aGCACGTATTCCCCGGCGAGGCTGTCCAGTTCTTCGGCAGAATCGAGCTTCATTCCAGACACTCCCGTAGCGCTTTGAGGCTGCGTTTGATCCAGGCCTTGACGGTGCCCAGTGGCGCACCCAATAAACCGGCTATCTCGGCGTGGCTg                                                                                                                                   >  1:242260/1‑149 (MQ=255)
                aCGTATTCCCCGGCGAGGCCGTCCAGTTCTTCGGCAGAATCGGGCTTCATTCCAGACACTCCCGTAGCGCTTTGAGGCTGCGTTTGATCCAGGCCTTGACGGTGCCCAGTGGCGCACCCAATAAACCGGCTATCTCGGCGTGGCTGCAg                                                                                                                                <  8:292682/149‑1 (MQ=255)
                   tATTCCCCGGCGAGGCTGTCCAGTTCTTCGGCAGAATCGGGCTTCATTCCAGACACTCCCGTAGCGCTTTGAGGCTGCGTTTGATCCAGGCCTTGACGGTGCCCAGTGGCGCACCCAATAAACCGGCTATCTCGGCGTGGCTGCAGCCg                                                                                                                             >  5:334128/1‑149 (MQ=255)
                        cccGGCGAGGCTGTCCAGTTCTTCGGCAGAATCGGGCTTCATTCCAGACACTCCCGTAGCGCTTTGAGGCTGCGTTTGATCCAGGCCTTGACGGTGCCCAGTGGCGCACCCAATAAACCGGCTATCTCGGCGTGGCTGCAGCCGTCCAc                                                                                                                        <  6:334128/149‑1 (MQ=255)
                         ccGGCGAGGCTGTCCAGTTCTTCGGCAGAATCGGGCTTCATTCCAGACACTCCCGTAGCGCTTTGAGGCTGCGTTTGATCCAGGCCTTGACGGTGCCCAGTGGCGCACCCAATAAACCGGCTATCTCGGCGTGGCTGCAGCCGTCAACg                                                                                                                       >  1:285158/1‑149 (MQ=255)
                           ggCGAGGCTGTCCATGTCTTCGGCAGAATCGGGCTTCATTCCAGACACTCCCGTAGCGCTTTGAGGCTGCGTTTGATCCAGGCCTTGACGGTGCCCAGTGGCGCACCCAATAAACCGGCTATCTCGGCGTGGCTGCAGCCGTCCACGTa                                                                                                                     <  6:232746/149‑1 (MQ=255)
                                  cTGTCCAGTTCTTCGGCAGAATCGGGCTTCATTCCAGACACTCCCGTAGCGCTTTGAGGCTGCGTTTGATCCAGGCCTTGACGGTGCCCAGTGGCGCACCCAATAAACCGGCTATCTCGGCGTGGCTGCAGCCGTCCACGTAGGCGTGc                                                                                                              <  1:149392/149‑1 (MQ=255)
                                   tGTCCAGTTCTTCGGCAGAATCGGGCTTCATTCCAGACACTCCCGTAGCGCTTTGAGGCTGCGTTTGATCCAGGCCTTGACGGTGCCCAGTGGCGCACCCAATAAACCGGCt                                                                                                                                                  <  2:119814/112‑1 (MQ=255)
                                     tCCAGTTCTTCGGCAGAATCGGGCTTCATTCCAGACACTCCCGTAGCGCTTTGAGGCTGCGGTTGATCCAGGCCTTGACGGTGCCCAGTGGCGCACCCACTAAACCGGCTATCTCGGCGTGGCTGCAGCCGTCCACGTAGGCGTGCa                                                                                                             <  7:57177/147‑1 (MQ=255)
                                         gTTCTTCGGCAGAATCGGGCTTCATTCCAGACACTCCCGTAGCGCTTTGAGGCTGCGTTTGATCCAGGCCTTGACGGTGCCCAGTGGCGCACCCAATAAACCGGCTATCTCGGCGTGGCTGCAGCCGTCCACGTAGGCGTGCAGGATg                                                                                                        >  4:205026/1‑148 (MQ=255)
                                              tcGGCAGAATCGGGCTTCATTCCAGACACTCCCGTAGCGCTTTGAGGCTGCGTTTGATCCAGGCCTTGACGGTGCCCAGTGGCGCACCCAATAAACCGGCTATCTCGGCGTGGCTGCAGCCGTCCACGTAGGCGTGCAGGATGCAGCGt                                                                                                  >  1:149341/1‑149 (MQ=255)
                                                                aTTCCAGACACTCCCGTAGCGCTTTGAGGCTGCGTTTGATCCAGGGCTTGACGGTGCCCAGTGGCGCACCCAATAAACCGGCTATCTCGGCGTGGCTGCAGCCGTCCACGTAGGCGTGCAGGATGCAGCGTTGCGGTGCGGGCTGCAg                                                                                 <  3:247904/148‑1 (MQ=255)
                                                                aTTCCAGACACTCCCGTAGCGCTTTGAGGCTGCGTTTGATCCAGGCCTTGACGGTGCCCAGTGGCGCACCCAATAAACCGGCTATCTAGGCGTGGCTGCAGCCGTCCACCTAGGCGTGCa                                                                                                             >  7:207890/1‑120 (MQ=255)
                                                                 ttCCAGACACTCCCGTAGCGCTTTGAGGCTGCGTTTGATCCAGGCCTTGACGGTGCCCAGTGGCGCACCCAATAAACCGGCTATCTCGGCGTGGCTGCAGCCGTCCACGTAGGCGTGCAGGATGCAGCGTTGCGGTGCGGGCTGCAgtt                                                                               >  4:283066/1‑149 (MQ=255)
                                                                     aGACACTCCCGTAGCGCTTTGAGGCTGCGTTTGATCCAGGCCTTGACGGTGCCCAGTGGCGCACCCAATAAACCGGCTATCTCGGCGTGGCTGCAGCCGTCCACGTAGGCGTGCAGGATGCAGCGTTGCGGTGCGGGCTGCAGTTGTTc                                                                           >  2:198651/1‑149 (MQ=255)
                                                                         acTCCCGTAGCGCTTTGAGGCTGCGTTTGATCCAGGCCTTGACGGTGCCCAGTGGCGCACCCAATAAACCGGCTATCTCGGCGTGGCTGCAGCCGTCCACGTAGGCGTGCAGGATGCAGCGTTGCGGTGCGGGCTGCAGTTGTTCCAgg                                                                       <  8:102511/149‑1 (MQ=255)
                                                                             ccGTAGCGCTTTGAGGCTGCGTTTGATCCAGGCCTTGACGGTGCCCAGTGGCGCACCCAATAAACCGGCTATCTCGGCGTGGCTGCAGCCGTCCACGTAGGCGTGCAGGATGCAGCGTTGCGGTGCGGGCTGCAGTTGTTCCAGGCaa                                                                    <  3:205026/148‑1 (MQ=255)
                                                                             ccGTAGCGCTTTGAGGCTGCGTTTGATCCAGGCCTTGACGGTGCCCAGTGGCGCACCCAATAAACCGGCTATCTCGGCGTGGCTGCAGCCGTCCACGTAGGCGTGCAGGATGCAGCGTTGCGGTGCGGGCTGCAGTTGTTCCAGGCAAc                                                                   >  6:341893/1‑149 (MQ=255)
                                                                             ccGTAGCGCTTTGAGGCTGCGTTTGATCCAGGCCTTGACGGTGCCCAGTGGCGCACCCAATAAACCGGCTATCTCGGCGTGGCTGCAGCCGTCCACGTAGGCGTGCAGGATGCAGCGTTGCGGTGCGGGCTGCAGTTGTTCCAGGCAAc                                                                   >  6:337293/1‑149 (MQ=255)
                                                                              cGTAGCGCTTTGAGGCTGCGTTTGATCCAGGCCTTGACGGTGCCCAGTGGCGCACCCAATAAACCGGCTATCTCGGCGTGGCTGCAGCCGTCCACGTAGGCGTGCAGGATGCAGCGTTGCGGTGCGGGCTGCAGTTGTTCCAGGCAACg                                                                  >  3:156607/1‑149 (MQ=255)
                                                                               gTAGCGCTTTGAGGCTGCGTTTGATCCAGGCCTTGACGGTGCCCAGTGGCGCACCCAATAAACCGGCTATCTCGGCGTGGCTGCAGCCGTCCACGTAGGCGTGCAGGATGCAGCGTTGCGGTGCGGGCTGCAGTTGTTCCAGGCAACgg                                                                 <  1:198651/149‑1 (MQ=255)
                                                                                         gAGGCTGCGTTTGATCCAGGCCTTGACGGTGCCCAGTGGCGCACCCAATAAACCGGCTATCTCGGCGTGGCTGCAGCCGTCCACGTAGGCGTGCAGGATGCAGCGTTGCGGTGCGGGCTGCAGTTGTTCCAGGCAACGGTAAATCttt                                                        <  2:149341/148‑1 (MQ=255)
                                                                                         gAGGCTGCGTTTGATCCAGGCCTTGACGGTGCCCAGTGGCGCACCCAATAAACCGGCTATCTCGGCGTGGCTGCAGCCGTCCACGTAGGCGTGCAGGATGCAGCGTTGCGGTGCGGGCTGCAGTTGTTCCAGGCAACGGTAAATCttt                                                        <  5:341893/148‑1 (MQ=255)
                                                                                         gAGGCTGCGTTTGATCCAGGCCTTGACGGTGCCCAGTGGCGCACCCAATAAACCGGCTATCTCGGCGTGGCTGCAGCCGTCCACGTAGGCGTGCAGGATGCAGCGTTGCGGTGCGGGCTGCAGTTGTTCCAGGCAACGGTAAATCttt                                                        <  5:337293/148‑1 (MQ=255)
                                                                                          aGGCTGCGTTTGATCCAGGCCTTGACGGTGCCCAGTGGCGCACCCAATAAACCGGCTATCTCGGCGTGGCTGCAGCCGTCCACGTAGGCGTGCAGGATGCAGCGTTGCGGTGCGGGCTGCAGTTGTTCCAGGCAACGGTAAATCTTTGc                                                      <  8:61252/149‑1 (MQ=255)
                                                                                          aGGCTGCGTTTGATCCAGGCCTTGACGGTGCCCAGTGACGCACCCAATAAACCGGCTATCTCGGCGTGGCTGCAGCCGTCCACGTAGGCGTGCAGGATGCAGCGTTGCGGTGCGGGCTGCAGTTGTTCCAGGCAACGGTAAATCTTTGc                                                      >  2:341291/1‑149 (MQ=255)
                                                                                              tGCGTTTGATCCAGGCCTTGACGGTGCCCAGTGGCGCACCCAATAAACCGGCTATCTCGGCGTGGCTGCAGCCGTCCACGTAGGCGTGCAGGATGCAGCGTTGCGGTGCGGGCTGCAGTTGTTCCAGGCAACGGTAAATCTTTGCTGAg                                                  >  8:182416/1‑149 (MQ=255)
                                                                                                cGTTTGATCCAGGCCTTGACGGTGCCCAGTGGCGCACCCAATAAACCGGCTATCTCGGCGGGGCTGCAGCCGTCCACGTAGGCGTGCAGGATGCAGCGTtgcggtgcgg                                                                                        <  2:245203/109‑1 (MQ=255)
                                                                                                  tttGATCCAGGCCTTGACGGTGCCCAGTGGCGCACCCAATAAACCGGCTATCTCGGCGTGGCTGCAGCCGTCCACGTAGGCGTGCAGGATGCAGCGTTGCGGTGCGGGCTGCAGTTGTTCCAGGCAACGGTAAATCTTTGCTGAGCCg                                               >  4:74391/1‑148 (MQ=255)
                                                                                                      aTCCAGGCCTTGACGGTGCCCAGTGGCGCACCCAATAAACCGGATATCTCGGCGTGGCTGCAGCCGTCCACGTAGGCGTGCAGGATGCAGCGTTGCGGTGCGGGCTGCAGTTGTTCCAGGCAACGGTAAATCTTTGCTGAGCCGGAc                                            >  1:180880/1‑147 (MQ=255)
                                                                                                      aTCCAGGCCTTGACGGTGCCCAGTGGCGCACCCAATAAACCGGATATCTCGGCGTGGCTGCAGCCGTCCACGTAGGCGTGCAGGATGCAGCGTTGCGGTGCGGGCTGCAGTTGTTCCAGGCAACGGTAAATCTTTGCTGAGCCGGAc                                            >  1:186356/1‑147 (MQ=255)
                                                                                                       tCCAGGCCTTGACGGTGCCCAGTGGCGCACCCAATAAACCGGCTATCTCGGCGTGGCTGCAGCCGTCCACGTAGGCGTGCAGGATGCAGCGTTGCGGTGCGGGCTGCAGTTGTTCCAGGCAACGGTAAATCTTTGCTGAGCCGGAc                                            >  1:113315/1‑146 (MQ=255)
                                                                                                         cAGGCCTTGACGGTGCCCAGTGGCGCACCCAATAAACCGGCTATCTCGGCGTGGCTGCAGCCGTCCACGTAGGCGTGCAGGATGCAGCGTTGCGGTGCGGGCTGCAGTTGTTCCAGGCAACGGTAAATCTTTGCTGAGCCGGACCaga                                        <  7:182416/148‑1 (MQ=255)
                                                                                                           ggCCTTGACGGTGCCCAGTGGCGCACCCAATAAACCGGCTATCTCGGCGTGGCTGCAGCCGTCCACGTAGGCGTGCAGGATGCAGCGTTGCGGTGCGGGCTGCAGTTGTTCCAGGCAACGGTAAATCTTTg                                                       >  8:216586/1‑131 (MQ=255)
                                                                                                           ggCCTTGACGGTGCCCAGTGGCGCACCCAATAAACCGGCTATCTCGGCGTGGCTGCAGCCGTCCACGTAGGCGTGCAGGATGCAGCGTTGCGGTGCGGGCTGCAGTTGTTCCAGGCAACGGTAAATCTTTg                                                       <  7:216586/131‑1 (MQ=255)
                                                                                                               ttGACGGTGCCCAGTGGCGCACCCAATAAACCGGCTATCTCGGCGTGGCTGCAGCCGTCCACGTAGGCGTGCAGGATGCAGCGTTGCGGTGCGGGCTGCAGTTGTTCCAGGCAACGGTAAATCTTTGCTGAGCCGGACCAGaggtcaag                                 >  4:253582/1‑149 (MQ=255)
                                                                                                                 gACGGTGCCCAGTGGCGCACCCAATAAACCGGCTATCTCGGCGTGGCTGCAGCCGTCCACGTAGGCGTGCAGGATGCAGCGTTGCGGTGCGGGCTGCAGTTGTTCCAGGCAACGGTAAATCTTTGCTGAGCCGGACCAGaggtcaag                                 <  2:113315/147‑1 (MQ=255)
                                                                                                                 gACGGTGCCCAGTGGCGCACCCAATAAACCGGATATCTCGGCGTGGCTGCAGCCGTCCACGTAGGCGTGCAGGATGCAGCGTTGCGGTGCGGGCTGCAGTTGTTCCAGGCAACGGTAAATCTTTGCTGAGCCGGACCAGaggtcaag                                 <  2:180880/147‑1 (MQ=255)
                                                                                                                 gACGGTGCCCAGTGGCGCACCCAATAAACCGGATATCTCGGCGTGGCTGCAGCCGTCCACGTAGGCGTGCAGGATGCAGCGTTGCGGTGCGGGCTGCAGTTGTTCCAGGCAACGGTAAATCTTTGCTGAGCCGGACCAGaggtcaag                                 <  2:186356/147‑1 (MQ=255)
                                                                                                                          cAGTGGCGCACCCAATAAACCGGCTATCTCGGCGTGGCTGCAGCCGTCCACGTAGGCGTGCAGGATGCAGCGTTGCGGTGCGGGCTGCAGTTGTTCCAGGCAACGGTAAATCTTTGCTGAGCCGGACCAGAGGTCAAGGTCATCAAccc                      >  1:157841/1‑149 (MQ=255)
                                                                                                                               gcgcACCCAATAAACCGGCTATCTCGGCGTGGCTGCAGCCGTCCACGTAGGCGTGCAGGATGCAGCGTTGCGGTGCGGGCTGCAGTTGTTCCAGGCAACGGTAAATCTTTGCTGAGCCGGACCagag                                       <  6:200289/127‑1 (MQ=255)
                                                                                                                               gcgcACCCAATAAACCGGCTATCTCGGCGTGGCTGCAGCCGTCCACGTAGGCGTGCAGGATGCAGCGTTGCGGTGCGGGCTGCAGTTGTTCCAGGCAACGGTAAATCTTTGCTGAGCCGGACCagag                                       >  5:200289/1‑127 (MQ=255)
                                                                                                                               gcgcACCCAATAAACCGGCTATCTCGGCGTGGCTGCAGCCGTCCACGTAGGCGTGCAGGATGCAGCGTTGCGGTGCGGGCTGCAGTTGTTCCAGGCAACGGTAAATCTTTGCTGAGCCGGACCAGAGGTCAAGGTCATCAACCCTccc                  <  3:253582/148‑1 (MQ=255)
                                                                                                                                   aCCCAATAAACCGGCTATCTCGGCGTGGCTGCAGCCGTCCACGTAGGCGTGCAGGATGCAGCGTTGCGGTGCGGGCTGCAGTTGTTCCAGGCAACGGTAAATCTTTGCTGAGACGGACCAGAGGTCAAGGTCATCAACCCTCCCGTGc              >  5:7594/1‑148 (MQ=255)
                                                                                                                                     ccAATAAACCGGCTATCTCGGCGTGGCTGCAGCCGTCCACGTAGGCGTGCAGGATGCAGCGTTGCGGTGCGGGCTGCAGTTGTTCCAGGCAACGGTAAATCTTTGCTGAGCCGGACCAGAGGTCAAGGTCATCAACCCTCCCGTGCAAt           <  3:74391/149‑1 (MQ=255)
                                                                                                                                      cAATAAACCGGCTATCTCGGCGTGGCTGCAGCCGTCCACGTAGGCGTGCAGGATGCAGCGTTGCGGTGCGGGCTGCAGTTGTTCCAGGCAACGGTAAATCTTTGCTGAGCCGGACCAGAGGTCAAGGTCATCAACCCTCCCGTGCAATg          <  2:285158/149‑1 (MQ=255)
                                                                                                                                       aaaaaaCAGGCTATCTCGGCGTGGCTGCAGCCGTCCACGTAGGCGTGCAGGATGCAGCGTTGCGGTGCGGGCTGCAGTTGTTCCAGGCAACGGTAAATCTTTGCTGAGCCGGACCAGAGGTCAAGGTCATCAACCCTCCCGTGCAATGt         >  7:155846/4‑149 (MQ=255)
                                                                                                                                          aaaCCGGCTATCTCGGCGTGGCTGCAGCCGTCCACGTAGGCGTGCAGGATGCAGCGTTGCGGTGCGGGCTGCAGTTGTTCCAGGCAACGGTAAATCTTTGCTGAGCCGGACCAGAGGTCAAGGTCATCAACCCTCCCGTGCAATGTTGa      >  1:70820/1‑149 (MQ=255)
                                                                                                                                          aaaCCGGCTATCTCGGCGTGGCTGCAGCCGTCCACGTAGGCGTGCAGGATGCAGCGTTGCGGTGCGGGCTGCAGTTGTTCCAGGCAACGGTAAATCTTTGCTGAGCCGGACCAGAGGTCAAGGTCATCAACCCTCCCGTGCAATGTTGa      >  1:75574/1‑149 (MQ=255)
                                                                                                                                           caCCGGCTATCTCGGCGGGGCTGCAGCCGTCCCCGTAGGCGTGCAGGATGCATCGTTGCGGTGCGGGCTGCAGTTGTTCCAGGCAACGGTAAATCTTTGCTGAGCCGGACCAGAGGTCAAGGTCATCAACCCTCCCGTGCAATGTTGAt     <  2:242260/148‑1 (MQ=255)
                                                                                                                                           aaCCGGCTATCTCGGCGTGGCTGCAGCCGTCCACGTAGGCGTGCAGGATGCAGCGTTGCGGTGCGGGCTGCAGTTGTTCCAGGCAACGGTAAATCTTTGCTGAGCCGGACCAGAGGTCAAGGTCATCAACCCTCCCGTGCAATGTTGAt     >  5:34199/1‑149 (MQ=255)
                                                                                                                                            aCCGGCTATCTCGGCGTGGCTGCAGCCGTCCACGTAGGCGTGCAGGATGCAGCGTTGCGGTGCGGGCTGCAGTTGTTCCAGGCAACGGTAAATCTTTGCTGAGCCGGACCAGAGGTCAAGGTCATCAACCCTCCCGTGCAATGTTGATg    >  1:42995/1‑149 (MQ=255)
                                                                                                                                              cGGCTATCTCGGCGTGGCTGCAGCCGTCCACGTAGGCGTGCAGGATGCAGCGTTGCGGTGCGGGCTGCAGTTGTTCCAGGCAACGGTAAATCTTTGCTGAGCCGGACCAGAGGTCAAGGTCATCAACCCTCCCGTGCAATGTTGATGCa  <  1:9151/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
CGACAAGGTTCCCAGCACGTATTCCCCGGCGAGGCTGTCCAGTTCTTCGGCAGAATCGGGCTTCATTCCAGACACTCCCGTAGCGCTTTGAGGCTGCGTTTGATCCAGGCCTTGACGGTGCCCAGTGGCGCACCCAATAAACCGGCTATCTCGGCGTGGCTGCAGCCGTCCACGTAGGCGTGCAGGATGCAGCGTTGCGGTGCGGGCTGCAGTTGTTCCAGGCAACGGTAAATCTTTGCTGAGCCGGACCAGAGGTCAAGGTCATCAACCCTCCCGTGCAATGTTGATGCA  >  NC_002947/3399959‑3400249

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