Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,233,549 G→A 11.1% G138D (GGC→GAC)  PP_4611 → ECF family RNA polymerase sigma‑70 factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,233,5490GA11.1% 218.2 / 13.1 82G138D (GGC→GAC) PP_4611ECF family RNA polymerase sigma‑70 factor
Reads supporting (aligned to +/- strand):  ref base G (37/35);  new base A (4/5);  total (41/41)
Fisher's exact test for biased strand distribution p-value = 7.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.21e-01

CTGCAAGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGTGCATTGGT  >  NC_002947/5233401‑5233694
                                                                                                                                                     |                                                                                                                                                 
cTGCAAGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACgg                                                                                                                                                   >  6:140641/1‑149 (MQ=255)
cTGCAAGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACgg                                                                                                                                                   >  2:173848/1‑149 (MQ=255)
 tGCAAGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGc                                                                                                                                                  >  8:17810/1‑149 (MQ=255)
  gCAAGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGcc                                                                                                                                                 >  5:61202/1‑149 (MQ=255)
     aGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCGTTCCTGCTGGCCCAGGTCGACGGCCTgg                                                                                                                                              <  1:182711/149‑1 (MQ=255)
     aGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTgg                                                                                                                                              <  1:187727/149‑1 (MQ=255)
     aGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTgg                                                                                                                                              <  5:2929/149‑1 (MQ=255)
     aGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGCCGTCGACCTgg                                                                                                                                              <  1:290913/149‑1 (MQ=255)
       cgcTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGc                                                                                                                                            <  6:135983/149‑1 (MQ=255)
       cgcTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGCCGTCGACCTGGGc                                                                                                                                            <  4:165256/149‑1 (MQ=255)
          tCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGCCGTCGACCTGGGcc                                                                                                                                           <  6:27862/147‑1 (MQ=255)
          tCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGCCGTCGACCTGGGcc                                                                                                                                           >  5:27862/1‑147 (MQ=255)
            gCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTTGACGGCCTGGGCCAAgg                                                                                                                                       >  1:307017/1‑149 (MQ=255)
               caGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGa                                                                                                                                    <  6:69538/149‑1 (MQ=255)
                        gAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGaa                                                                                                                                   <  5:275204/141‑1 (MQ=255)
                        gAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGaa                                                                                                                                   >  6:275204/1‑141 (MQ=255)
                             ggTAGCGCCGTGCGATGAACAGAATGCTATAATCTTCGAGACATTGCTAGCGCTCGACGACTTGCTGGACGGCTTTCCAGCGTTGGTAAAAGGCGCGTTCCTTCTGGCCCAGGTCGACGTCATTGGCGTAGTCGAAATCGGTGGCGaa                                                                                                                       <  8:69322/148‑1 (MQ=255)
                                 gcgcCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTggg                                                                                                                  <  5:140641/149‑1 (MQ=255)
                                 gcgcCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGCCGTCGACCTGGGc                                                                                                                                            <  7:130540/123‑1 (MQ=255)
                                 gcgcCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGCCGTCGACCTGGGc                                                                                                                                            >  8:130540/1‑123 (MQ=255)
                                       tGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTc                                                                                                            >  8:32578/1‑149 (MQ=255)
                                        gCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAg                                                                                                                                        >  4:131187/1‑120 (MQ=255)
                                        gCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAg                                                                                                                                        <  3:131187/120‑1 (MQ=255)
                                          gAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGc                                                                                                          >  1:85394/1‑149 (MQ=255)
                                          gAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCt                                                                                                         >  2:180427/1‑149 (MQ=255)
                                             gAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGc                                                                                                      >  6:265106/1‑149 (MQ=255)
                                                      gCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGACAAGGCGAAATCGCTCGCGAACTCGGCATTTCGCTGGCCACGGTc                                                                                               >  3:142510/1‑147 (MQ=255)
                                                        gATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAgc                                                                                           <  6:188817/149‑1 (MQ=255)
                                                           catcTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAgcgc                                                                                         >  6:61333/1‑148 (MQ=255)
                                                            atcTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTTGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTa                                                                                       <  2:307017/149‑1 (MQ=255)
                                                                 cGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTAcc                                                                                     <  5:265106/146‑1 (MQ=255)
                                                                 cGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTg                                                                                   >  2:182604/1‑148 (MQ=255)
                                                                   agaCATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAAc                                                                                <  4:296582/149‑1 (MQ=255)
                                                                    gaCATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAAc                                                                                >  3:265066/1‑148 (MQ=255)
                                                                     aCATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAAc                                                                                >  4:67440/1‑147 (MQ=255)
                                                                     aCATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACaa                                                                              <  1:153258/149‑1 (MQ=255)
                                                                       aTTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAgg                                                                            >  8:158340/1‑149 (MQ=255)
                                                                       aTTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAgg                                                                            >  5:118467/1‑149 (MQ=255)
                                                                             aGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCa                                                                      >  7:125945/1‑149 (MQ=255)
                                                                             aGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCa                                                                      >  3:233468/1‑149 (MQ=255)
                                                                              gcgcTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCAt                                                                     >  3:66364/1‑149 (MQ=255)
                                                                                gcTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCTAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATgc                                                                   >  6:249722/1‑149 (MQ=255)
                                                                                gcTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATgc                                                                   >  2:197591/1‑149 (MQ=255)
                                                                                gcTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATgc                                                                   >  2:191754/1‑149 (MQ=255)
                                                                                            gCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTc                                                        <  5:87080/148‑1 (MQ=255)
                                                                                                 gACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCt                                                   >  8:41085/1‑149 (MQ=255)
                                                                                                   cGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGc                                                                                                                 <  6:295806/85‑1 (MQ=255)
                                                                                                       cTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAAc                                             <  6:118467/149‑1 (MQ=255)
                                                                                                          ccACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACg                                            <  2:186943/147‑1 (MQ=255)
                                                                                                          ccACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGc                                           >  4:215661/1‑148 (MQ=255)
                                                                                                           cACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCt                                         >  7:300509/1‑149 (MQ=255)
                                                                                                           cACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCt                                         <  7:179514/149‑1 (MQ=255)
                                                                                                                    tCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTccc                                >  8:106265/1‑149 (MQ=255)
                                                                                                                    tCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTccc                                >  4:244124/1‑149 (MQ=255)
                                                                                                                        aCGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCt                            <  8:125945/149‑1 (MQ=255)
                                                                                                                         cgcgCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCtg                           <  1:191754/149‑1 (MQ=255)
                                                                                                                          gcgcCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTc                                  <  2:106279/141‑1 (MQ=255)
                                                                                                                          gcgcCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTc                                  >  1:106279/1‑141 (MQ=255)
                                                                                                                            gcCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTcgcg                                                     >  6:124882/1‑120 (MQ=255)
                                                                                                                            gcCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTcgcg                                                     <  5:124882/120‑1 (MQ=255)
                                                                                                                            gcCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTcc                                 >  2:23869/1‑139 (MQ=255)
                                                                                                                            gcCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTcc                                 <  1:23869/139‑1 (MQ=255)
                                                                                                                              cTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGc                      >  3:102630/1‑149 (MQ=255)
                                                                                                                              cTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGc                      <  4:233468/149‑1 (MQ=255)
                                                                                                                              cTTCCTGAAGGCCCAGGTCGCCCGCCAGGGCAAAGGCCAAATCTCTCGCGAAATCGGCATATCGCTGGCAACGGTCAATCGCTACCTGAACAAGGCAGCCATCCGCTGCTACTTCGCCATGAGAACCAGTCGTTTTCCCCGCTAGTGGc                      >  4:214071/1‑149 (MQ=255)
                                                                                                                                 ccTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCg                     <  5:61333/147‑1 (MQ=255)
                                                                                                                                 ccTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCg                     <  4:265066/147‑1 (MQ=255)
                                                                                                                                    cttgGCCCAGGCCGTCGACCTGGGCCAGTGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCgctgc                                                              <  1:266625/101‑1 (MQ=38)
                                                                                                                                    cttgGCCCAGGCCGTCGACCTGGGCCAGAGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCgctgc                                                              >  2:266625/3‑103 (MQ=38)
                                                                                                                                    cttgGCCCAGGCCGTCGACCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGaa                  >  3:77407/3‑147 (MQ=255)
                                                                                                                                     ctgGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAg               >  6:138888/1‑149 (MQ=255)
                                                                                                                                        gCCCTGGTCTACGGCCTGGGCCAAGGCGAAATCGCTCGCTAACTGGGCAGTTCGCTGGCCACGGCCAAGCGCTCCCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGcc            <  3:214071/149‑1 (MQ=255)
                                                                                                                                        gCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGcc            <  6:61202/149‑1 (MQ=255)
                                                                                                                                        gCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGcc            >  4:315423/1‑149 (MQ=255)
                                                                                                                                        gCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGc             <  7:41085/148‑1 (MQ=255)
                                                                                                                                         cccAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCg           <  2:85394/149‑1 (MQ=255)
                                                                                                                                          ccAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGt          >  2:176172/1‑149 (MQ=255)
                                                                                                                                          ccAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGt          <  1:197591/149‑1 (MQ=255)
                                                                                                                                           cAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTATCCCCGCTGCTGGCGGAACAGCCCGTg         >  8:275951/1‑149 (MQ=255)
                                                                                                                                              gTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCg                              <  3:275482/125‑1 (MQ=255)
                                                                                                                                              gTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCg                              >  4:275482/1‑125 (MQ=255)
                                                                                                                                                cGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGTGCATTg    <  5:249722/149‑1 (MQ=255)
                                                                                                                                                 gACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGTGCATTg    <  3:315423/148‑1 (MQ=255)
                                                                                                                                                  aCGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGTGCATTGGt  <  1:180427/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
CTGCAAGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGTGCATTGGT  >  NC_002947/5233401‑5233694

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: