Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 655,691 A→T 7.3% C218S (TGC→AGC)  PP_0563 ← two‑component system response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947655,6910AT 7.3% 236.1 / 7.2 82C218S (TGC→AGC) PP_0563two‑component system response regulator
Reads supporting (aligned to +/- strand):  ref base A (37/39);  new base T (4/2);  total (41/41)
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.98e-01

GCTTGGTCAAAAAGTCATCACCGCCTTCGCTCATGGCATCCAGTTGCTTGTCCAGGTCGTCCTCGGCAGACAGGTAGATGATCGGCACACTCACATAACGGTCACTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGCTGGAAGTCGGCAAGCTCGGCCATGGTGCGGATCGGGTCGGTCAGGCTGCGGGTGAGCATGCCGGCGCCGGCCAGTACACGCTCGGTGTACATGGCCTGTGTGCGCGAATCG  >  NC_002947/655550‑655839
                                                                                                                                             |                                                                                                                                                    
gCTTGGTCAAAAAGTCATCACCGCCTTCGCTCATGGCATCCAGTTGCTTGTCCAGGTCGTCCTCGGCAGACAGGTAGATGATCGGCACACTCACATAACGGTCACTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGc                                                                                                                                               <  1:117018/149‑1 (MQ=255)
gCTTGGTCAAAAAGTCATCACCGCCTTCGCTCATGGCATCCAGTTGCTTGTCCAGGTCGTCCTCGGCAGACAGGTAGATGATCGGCACACTCACATAACGGTCACTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGc                                                                                                                                               <  6:248466/149‑1 (MQ=255)
     gTCAAAAAGTCATCACCGCCTTCGCTCATGGCATCCAGTTGCTTGTCCAGGTCGTCCTCGGCAGACAGGTAGATGATCGGCACACTCACATAACGGTCACTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCAt                                                                                                                                             >  5:57676/1‑146 (MQ=255)
     gTCAAAAAGTCATCACCGCCTTCGCTCATGGCATCCAGTTGCTTGTCCAGGTCGTCCTCGGCAGACAGGTAGATGATCGGCACACTCACATAACGGTCACTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCAt                                                                                                                                             <  6:57676/146‑1 (MQ=255)
     gTCAAAAAGTCATCACCGCCTTCGCTCATGGCATCCAGTTGCTTGTCCAGGTCGTCCTCGGCAGACAGGTAGATGATCGGCACACTCACATAACGGTCACTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTa                                                                                                                                          >  5:90681/1‑149 (MQ=255)
       cAAAAAGTCATCACCGCCTTCGCTCATGGCATCCAGTTGCTTGTCCAGGTCGTCCTCGGCAGACAGGTAGATGATCGGCACACTCACATAACGGTCACTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTa                                                                                                                                          >  7:341805/1‑147 (MQ=255)
        aaaaaGTCATCACCGCCTTCGCTCATGGCATCCAGTTGCTTGTCCAGGTCGTCCTCGGCAGACAGGTAGATGATCGGCACACTCACATAACGGTCACTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAg                                                                                                                                       >  7:293060/1‑149 (MQ=255)
           aaGTCATCACCGCCTTCGCTCATGGCATCCAGTTGCTTGTCCAGGTCGTCCTCGGCAGACAGGTAGATGATCGGCACACTCACATAACGGTCACTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTc                                                                                                                                    <  2:341819/149‑1 (MQ=255)
            aGTCATCACCGCCTTCGCTCATGGCATCCAGTTGCTTGTCCAGGTCGTCCTCGGCAGACAGGTAGATGATCGGCACACTCACATAACGGTCACTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCg                                                                                                                                   <  8:341805/149‑1 (MQ=255)
             gTCATCACCGCCTTCGCTCATGGCATCCAGTTGCTTGTCCAGGTCGTCCTCGGCAGACAGGTAGATGATCGGCACACTCACATAACGGTCACTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGa                                                                                                                                  <  6:335233/149‑1 (MQ=255)
                atcaCCGCCTTCGCTCATGGCATCCAGTTGCTTGTCCAGGTCGTCCTCGGCAGACAGGTAGATGATCGGCACACTCACATAACGGTCACTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAgg                                                                                                                                <  7:27007/148‑1 (MQ=255)
                 tcaCCGCCTTCGCTCATGGCATCCAGTTGCTTGTCCAGGTCGTCCTCGGCAGACAGGTAGATGATCGGCACACTCACATAACGGTCACTGTGGCGGATTACCTTGGCCTGCACCGGCCCGGAGCTGGccaag                                                                                                                                               <  6:265499/132‑4 (MQ=255)
                 tcaCCGCCTTCGCTCATGGCATCCAGTTGCTTGTCCAGGTCGTCCTCGGCAGACAGGTAGATGATCGGCACACTCACATAACGGTCACTGTGGCGGATTACCTTGGCCTGCACCGGCCCGGAGCTGGccaag                                                                                                                                               >  5:265499/1‑129 (MQ=255)
                  caCCGCCTTCGCTCATGGCATCCAGTTGCTTGTCCAGGTCGTCCTCGGCAGACAGGTAGATGATCGGCACACTCACATAACGGTCACTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGgat                                                                                                                              >  6:74848/1‑148 (MQ=255)
                  caCCGCCTTCGCTCATGGCATCCAGTTGCTTGTCCAGGTCGTCCTCGGCAGACAGGTAGATGATCGGCACACTCACATAACGGTCACTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGga                                                                                                                               <  6:134237/147‑1 (MQ=255)
                  caCCGCCTTCGCTCATGGCATCCAGTTGCTTGTCCAGGTCGTCCTCGGCAGACAGGTAGATGATCGGCACACTCACATAACGGTCACTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGga                                                                                                                               >  5:134237/1‑147 (MQ=255)
                   aCCGCCTTCGCTCATGGCATCCAGTTGCTTGTCCAGGTCGTCCTCGGCAGACAGGTAGATGATCGGCACACTCACATAACGGTCACTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGgatga                                                                                                                            >  6:38547/1‑149 (MQ=255)
                   aCCGCCTTCGCTCATGGCATCCAGTTGCTTGTCCAGGTCGTCCTCGGCAGACAGGTAGATGATCGGCACACTCACATAACGGTCACTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGgatg                                                                                                                             <  8:119269/148‑1 (MQ=255)
                   aCCGCCTTCGCTCATGGCATCCAGTTGCTTGTCCAGGTCGTCCTCGGCAGACAGGTAGATGATCGGCACACTCACATAACGGTCACTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGgatg                                                                                                                             <  8:124388/148‑1 (MQ=255)
                   aCCGCCTTCGCTCATGGCATCCAGTTGCTTGTCCAGGTCGTCCTCGGCAGACAGGTAGATGATCGGCACACTCACATAACGGTCACTGTGGCGGATTACCTAGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAgcaag                                                                                                                             >  1:145864/1‑144 (MQ=255)
                    cccccTTCGCTCATGGCATCCAGTTGCTTGTCCAGGTCGTCCTCGGCACACAGGAAGATGATCGGCACACTCAAATAACGGTCACTGAGGCGGATTACCTTGGACAGCTACGGGCCGGTGCAGGCCGGCATGTACAGGTAGAGgaaga                                                                                                                            >  2:5144/4‑148 (MQ=255)
                     cGCCTTCGCTCATGGCATCCAGTTGCTTGTCCAGGTCGTCCTCGGCAGACAGGTAGATGATCGGCACACTCACATAACGGTCACTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGgatgat                                                                                                                           >  3:321691/1‑148 (MQ=255)
                      gCCTTCGCTCATGGCATCCAGTTGCTTGTCCAGGTCGTCCTCGGCAGACAGGTAGATGATCGGCACACTCACATAACGGTCACTGTGGCGGATTACCTTGGCCAGCTCCTGGCCGGTGCAGGCCGGCATGTACAGGTCGAGgatgat                                                                                                                           >  5:102617/1‑147 (MQ=255)
                       ccTTCGCTCATGGCATCCAGTTGCTTGTCCAGGTCGTCCTCGGCAGACAGGTAGATGATCGGCACACTCACATAACGGTCACTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGGATGATCAg                                                                                                                        <  5:292050/149‑1 (MQ=255)
                             cTCATGGCATCCAGTTGCTTGTCCAGGTCGTCCTCGGCAGACAGGTAGATGATCGGCACACTCACATAACGGTCACTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGGATGATCAGGTCggg                                                                                                                  >  4:28482/1‑149 (MQ=255)
                              tCATGGCATCCAGTTGCTTGTCCAGGTCGTCCTCGGCAGACAGGTAGATGATCGGCACACTCACATAACGGTCACTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGc                                                                                                                 >  3:253448/1‑149 (MQ=255)
                                     aTCCAGTTGCTTGTCCAGGTCGTCCTCGGCAGACAGGTAGATGATCGGCACACTCACATAACGGTCACTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCgg                                                                                                                                                <  8:290544/111‑1 (MQ=255)
                                      tCCAGTTGCTTGTCCAGGTCGTCCTCGGCAGACAGGTAGATGATCGGCACACTCACATAACGGTCACTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGCTGGAAGTc                                                                                                         <  4:270417/149‑1 (MQ=255)
                                       ccAGTTGCTTGTCCAGGTCGTCCTCGGCAGACAGGTAGATGATCGGCACACTCACATAACGGTCACTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGCTGGAAGTCg                                                                                                        >  7:46150/1‑149 (MQ=255)
                                        cAGTTGGTTGTGCAGGTCGTCCTCGGCAGACAGGTTGATGATCGGCACACTCACATAAGGGTCAGTGTGGCGGATTAGCTTGGCCAGCTCCGTGCCGGTGCAGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGCTGGAAGtag                                                                                                        <  1:5144/148‑3 (MQ=255)
                                        cAGTTGCTTGTCCAGGTCGTCCTCGGCAGACAGGTAGATGATCGGCACACTCACATAACGGTCACTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGCTGGAAGTCg                                                                                                        <  4:253448/148‑1 (MQ=255)
                                           ttGCTTGTCCAGGTCGTCCTCGGCAGACAGGTAGATGATCGGCACACTCACATAACGGTCACTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGCTGGAAGTCGGCaa                                                                                                    <  6:90681/149‑1 (MQ=255)
                                                          gtcCTCGGCATACATGTAGATGATCGGCACACTCTCATAACGGTGACTGTGGCGGATTTCCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGCTGGAAGTCGGCAAGCTCGGCCATGGTg                                                                                      <  2:270752/148‑1 (MQ=255)
                                                          gtcCTCCGCAGACAGCTAGATGATCGGCACACTCACATAACGGTCACTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGCTGGAAGTCGGCAAGCTCGGCCATGGTGc                                                                                     <  2:132404/149‑1 (MQ=255)
                                                             cTCGGCAGACAGGTAGATGATCGGCACACTCACATAACGGTCACTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGACTGGAAGTCGGCAAGCTCGGCCATGGTGCGGa                                                                                  >  6:257728/1‑149 (MQ=255)
                                                                ggCAGACAGGTAGATGATCGGCACACTCACATAACGGTCACTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGCTGGAAGTCGGCAAGCTCGGCCATGGTGCGGATCg                                                                               <  8:293060/149‑1 (MQ=255)
                                                                 gCAGACAGGTTGATGATCGGCACTCTCAGATAACTGTCAGTGTGGCGGATTACCTTGGCCAGCTCGGGGCCGGTGCAGGGCTGCATGTAGAGGTCGAGGATGATCAGGTCGGGCTGGAAGTCGGCAAGCTCGGCCATGGTGCGGATCgg                                                                              <  2:145864/149‑1 (MQ=255)
                                                                    gACAGGTAGATGATCGGCACACTCACATAACGGTCACTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGCTGGAAGTCGGCAAGCTCGGCCATGGTGCGGATCGggtc                                                                           >  7:256950/1‑149 (MQ=255)
                                                                     aCAGGTAGATGATCGGCACACTCACATAACGGTCACTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGCTGGAAGTCGGCAAGCTCGGCCATGGTGCGGATCGggtcc                                                                          >  3:310966/1‑148 (MQ=255)
                                                                        ggTAGATGATCGGCACACTCACATAACGGTCACTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGCTGGAAGTCGGCAAGCTCGGCCATGGTGCGGATCGggtcggt                                                                        >  5:271416/1‑148 (MQ=255)
                                                                         gTAGATGATCGGCACACTCACATAACGGTCACTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAgg                                                                                                                                <  8:240994/91‑1 (MQ=255)
                                                                         gTAGATGATCGGCACACTCACATAACGGTCACTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAgg                                                                                                                                >  7:240994/1‑91 (MQ=255)
                                                                         gTAGATGATCGGCACACTCACATAACGGTCACTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGCTGGAAGTCGGCAAGCTCGGCCATGGTGCGGATCGggtcggt                                                                        >  2:89763/1‑147 (MQ=255)
                                                                         gTAGATGATCGGCACACTCACATAACGGTCACTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACACGTCGAGGATGATCAGGTCGGGCTGGAAGTCGGCAAGCTCGGCCATGGTGCGGATCGggtcggt                                                                        >  1:280274/1‑147 (MQ=255)
                                                                          tAGATGATCGGCACACTCACATAACGGTCACTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGCTGGAAGTCGGCAAGCTCGGCCATGGTGCGGATCGGGTCGGTCAg                                                                     >  6:27378/1‑149 (MQ=255)
                                                                           aGATGATCGGCACACTCACATAACGGTCACTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGCTGGAAGTCGGCAAGCTCGGCCATGGTGCGGATCGGGTCGGTCAgg                                                                    >  7:245392/1‑149 (MQ=255)
                                                                           aGATGATCGGCACAATCACATAACGGTCACTGTGGCGGATTACCTTGGCCAGATCCGCGCCGGTGCAGGCCGGCATGTACAGCTCCAGGATGATCAGGTCGGCCTGGAAGTCGGCAAGCTCGGCCATGGTGCGGATCGGGTCGGTCAgg                                                                    >  5:169811/1‑149 (MQ=255)
                                                                            gatgatCGGCACACTCACATAACGGTCACTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGCTGGAAGTCGGCAAGCTCGGCCATGGTGCGGATCGGGTCGGTCAGGc                                                                   >  8:57279/1‑149 (MQ=255)
                                                                                    gCACACTCACATAACGGTCACTGTGGCGGATTACCTTGGCCAGCTCTGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGCTGGAAGTCGGCAAGCTCGGCCATGGTGCGGATCGGGTCGGTCAGGCTGCGGGt                                                            >  6:277946/1‑148 (MQ=255)
                                                                                     cacaCTCACATAACGGTCACTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGCTGGAAGTCGGCAAGCTCGGCCATGGTGCGGATCGGGTCGGTCAGGCTGCGGGTGa                                                          <  8:46150/149‑1 (MQ=255)
                                                                                     cacaCTCACATAACGGTCACTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGCTGGAAGTCGGCAAGCTCGGCCATGGTGCGGATCGGGTCGGTCAGGCTGCGGGTGa                                                          <  5:27378/149‑1 (MQ=255)
                                                                                              aTAACGGTCACTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGCTGGAAGTCGGCAAGCTCGGCCATGGTGCGGAt                                                                                 <  6:326564/117‑1 (MQ=255)
                                                                                              aTAACGGTCACTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGCTGGAAGTCGGCAAGCTCGGCCATGGTGCGGAt                                                                                 >  5:326564/1‑117 (MQ=255)
                                                                                                aaCGGTCACTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGCTGGAAGTCGGCAAGCTCGGCCATGGTGCGGATCGGGTCGGTCAGGCTGCGGGTGAGCATgccggcg                                               >  4:15077/1‑149 (MQ=255)
                                                                                                   ggTCACTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGCTGGAAGTCGGCAAGCTCGGCCATGGTGCGGATCggg                                                                             >  6:301393/1‑116 (MQ=255)
                                                                                                   ggTCACTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGCTGGAAGTCGGCAAGCTCGGCCATGGTGCGGATCggg                                                                             <  5:301393/116‑1 (MQ=255)
                                                                                                   ggTCACTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGCTGGAAGTCGGCAAGCTCGGCCATGGTGCGGATCGGGTCGGTCAGGCTGCGGGTGAGCATgccggcgccg                                            >  7:115189/1‑149 (MQ=255)
                                                                                                       aCTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGCTGGAAGTCGGCAAGCTCGGCCATGGTGCGGATCGTGTCGGTCAGGCTGCGGGTGAGCATGCCGGCGCCGGCCa                                        <  2:280274/149‑1 (MQ=255)
                                                                                                       aCTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGCTGGAAGTCGGCAAGCTCGGCCATGGTGCGGATCGGGTCGGTCAGGCTGCGGGTGAg                                                         <  6:136077/132‑1 (MQ=255)
                                                                                                       aCTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGCTGGAAGTCGGCAAGCTCGGCCATGGTGCGGATCGGGTCGGTCAGGCTGCGGGTGAg                                                         >  5:136077/1‑132 (MQ=255)
                                                                                                       aCTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGCTGGAAGTCGGCAAGCTCGGCCATGGTGCGGATCGGGTCGGTCAGGCTGCGGGTGAGCATgccggcgccggc                                          >  1:49095/1‑147 (MQ=255)
                                                                                                            ggCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGCTGGAAGTCGGCAAGCTCGGCCATGGTGCGGATCGGGTCGGTCAGGCTGCGGGTGAGCATGCCGGCGCCGGCCAGTaca                                   <  1:89763/149‑1 (MQ=255)
                                                                                                             gCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGCTGGAAGTCGGCAAGCTCGGCCATGGTGCGGATCGGGTCGGTCAGGCTGCGGGTGAGCATGCCGGCGCCGGCCAGTacac                                  <  2:49095/149‑1 (MQ=255)
                                                                                                              cGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGCTGGAAGTCGGCAAGCTCGGCCATGGTGCGGATCGGGTCGGTCAGGCTGCGGGTGAGCATGCCGGCGCCGGCCAGTacac                                  <  5:38547/148‑1 (MQ=255)
                                                                                                               ggATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGCTGGAAGTCGGCAAGCTCGGCCATGGTGCGGATCGGGTCGGTCAGGCTGCGGGTGAGCATGCCGGCGCCGGCCAGTACACGc                                >  8:184634/1‑149 (MQ=255)
                                                                                                                 aTTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGCTGGAAGTCGGCAAGCTCGGCCATGGTGCGGATCGGGTCGGTCAGGCTGCGGGTGAGCATGCCGGCGCCGGCCAGTACACGCTc                              <  5:74848/149‑1 (MQ=255)
                                                                                                                  ttACCTTGGCCAGCTCTGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGCTGGAAGTCGGCAAGCTCGGCCATGGTGCGGATCGGGTCGGTCAGGCTGCGGGTGAGCATGCCGGCGCCGGCCAGTACACGCTCg                             <  5:277946/149‑1 (MQ=255)
                                                                                                                    atgccggCCTGCACCGGCCCGGAGCTGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGCTGGAAGTCGGCAAGCTCGGCCATGGTGCGGATCGGGTCGGTCAGGCTGCGGGTGAGCATGCCGGCGCCGGCCAGTACACGCTCGgt                           <  5:195609/144‑1 (MQ=255)
                                                                                                                    atgccggCCTGCACCGGCCCGGAGCTGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGCTGGAAGTCGGCAAGCTCGGCCATGGTGCGGATCGGGTCGGTCAGGCTGCGGGTGAGCATGCCGGCGCCGGCCAGTACACGCTCGgt                           >  6:195609/6‑149 (MQ=255)
                                                                                                                     ccTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGCTGGAAGTCGGCAAGCTCGGCCATGGTGCGGATCGGGTCGGTCAGGCTGCGGGTGAGCATGCCGGCGCCGGCCAGTACACGCTCGgt                           <  8:245392/148‑1 (MQ=255)
                                                                                                                      gccggCCTGCACCGGCCCGGAGCTGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGCTGGAAGTCGGCAAGCTCGGCCATGGTGCGGATCGGGTCGGTCAGGCTGCGGGTGAGCATGCCGGCGCCGGCCAGTACACGCTCGgtg                          >  6:73358/4‑148 (MQ=255)
                                                                                                                        cggCCTGCACCGGCCCGGAGCTGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGCTGGAAGTCGGCAAGCTCGGCCATGGTGCGGATCGGGTCGGTCAGGCTGCGGGTGAGCATGCCGGCGCCGGCCAGTACACGCTCGGTGTa                        >  1:71572/2‑148 (MQ=255)
                                                                                                                             aGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGCTGGAAGTCGGCAAGCTCGGCCATGGTGCGGATCGGGTCGGTCAGGCTGCGGGTGAGCATGCCGGCGCCGGCCAGTACACGCTCGGTGTACATGGc                  <  8:256950/149‑1 (MQ=255)
                                                                                                                               cTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGCTGGAAGTCGGCAAGCTCGGCCATGGTGCGGATCGGGTCGGTCAGGCTGCGGGTGAGCATGCCGGCGCCGGCCAGTACACTCTCGGTGTACATGGCCt                <  4:310966/149‑1 (MQ=255)
                                                                                                                               cTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGCTGGAAGTCGGCAAGCTCGGCCATGGTGCGGATCGGGTCGGTCAGGCTGCGGGTGAGCATGCCGGCGCCGGCCAGTACACGCTCGGTGTACATGGCCt                <  7:319559/149‑1 (MQ=255)
                                                                                                                               cTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGCTGGAAGTCGGCAAGCTCGGCCATGGTGCGGATCGGGTCGGTCAGGCTGCGGGTGAGCATGCCGGCGCCGGCCAGTACACGCTCGGTGTACATGGCCt                >  8:319559/1‑149 (MQ=255)
                                                                                                                                 ccGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGCTGGAAGTCGGCAAGCTCGGCCATGGTGCGGATCGGGTCGGTCAGGCTGCGGGTGAGCATGCCGGCGCCGGCAAGTACACGCTCgg                            <  2:34276/135‑1 (MQ=255)
                                                                                                                                 ccGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGCTGGAAGTCGGCAAGCTCGGCCATGGTGCGGATCGGGTCGGTCAGGCTGCGGGTGAGCATGCCGGCGCCGGCAAGTACACGCTCgg                            >  1:34276/1‑135 (MQ=255)
                                                                                                                                     gCCGGTGCAGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGCTGGAAGTCGGCAAGCTCGGCCATGGTGCGGATCGGGTCGGTCAGGCTGCGGGTGAGCATGCCGGCGCCGGCCAGTACACGCTCGGTGTACATGGCCTGTGTgc          <  7:57279/149‑1 (MQ=255)
                                                                                                                                       cGGTGCAGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGCTGGAAGTCGGCAAGCTCGGCCATGGTGCGGATCGGGTCGGTCAGGCTGCGGGTGAGCATGCCGGCGCCGGCCAGTACACGCTCGGTGTACATGGCCTGTGTgcgc        <  8:115189/149‑1 (MQ=255)
                                                                                                                                       cGGTGCAGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGCTGGAAGTCGGCAAGCTCGGCCATGGTGCGGATCGGGTCGGTCAGGCTGCGGGTGAGCATGCCGGCGCCGGCCAGTACACGCTCGGTGTACATGGCCTGTGTgcgc        >  8:94251/1‑149 (MQ=255)
                                                                                                                                       cGGTGCAGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGCCGGAAGTCGGCAAGCTAGGCCATGGTGCGGATCGGGTCGGTCAGGCTGCGGGTGAGCATGCCGGCGCCGGCCAGTACACGCTCGGTGTACATGGCCTGTGTgcgc        >  1:295098/1‑149 (MQ=255)
                                                                                                                                             aGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGCTGGAAGTCGGCAAGCTTGGCCATGGTGCGGATCGGGTCGGTCAGGCTGCGGGTGAGCATGCCGGCGCCGGCCAGTACACGCTCGGTGTACATGGCCTGTGTGCGCGAAtcg  >  3:104631/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
GCTTGGTCAAAAAGTCATCACCGCCTTCGCTCATGGCATCCAGTTGCTTGTCCAGGTCGTCCTCGGCAGACAGGTAGATGATCGGCACACTCACATAACGGTCACTGTGGCGGATTACCTTGGCCAGCTCCGGGCCGGTGCAGGCCGGCATGTACAGGTCGAGGATGATCAGGTCGGGCTGGAAGTCGGCAAGCTCGGCCATGGTGCGGATCGGGTCGGTCAGGCTGCGGGTGAGCATGCCGGCGCCGGCCAGTACACGCTCGGTGTACATGGCCTGTGTGCGCGAATCG  >  NC_002947/655550‑655839

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: