Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,674,660 T→G 9.9% intergenic (+22/+118) rhlE‑II → / ← yceI ATP‑dependent RNA helicase/polyisoprene‑binding acidic stress response factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,674,6600TG 9.9% 161.5 / 11.0 61intergenic (+22/+118)rhlE‑II/yceIATP‑dependent RNA helicase/polyisoprene‑binding acidic stress response factor
Reads supporting (aligned to +/- strand):  ref base T (29/26);  new base G (2/4);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 4.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.56e-01

GAAAAACAGCCAGTGATCATTCGCAAGGAGTCCAAGCTGGACCGTTTCCCGACGCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGC‑GCGCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTG  >  NC_002947/5674513‑5674805
                                                                                                                                                   |                                                                                                                                                  
gAAAAACAGCCAGTGATCATTCGCAAGGAGTCCAAGCTGGACCGTTTCCCGACGCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCtt                                                                                                                                                   >  6:369413/1‑149 (MQ=255)
     aCAGCCAGTGATCATTCGCAAGGAGTCCAAGCTGGACCGTTTCCCGACGCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAg                                                                                                                                              >  6:299252/1‑149 (MQ=255)
       aGCCAGTGATCATTCGCAAGGAGTCCAAGCTGGACCGTTTCCCGACGCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCGGCGAAGCAGcc                                                                                                                                            <  8:55895/149‑1 (MQ=255)
        gCCAGTGATCATTCGCAAGGAGTCCAAGCTGGACCGTTTCCCGACGCCTGAGCCACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCGAAGCAGccc                                                                                                                                           <  5:368531/149‑1 (MQ=255)
        gCCAGTGATCATTCGCAAGGAGTCCAAGCTGGACCGTTTCCCGACGCCTGAGCAACTGGCTGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCGAAGCAGccc                                                                                                                                           <  1:258449/149‑1 (MQ=255)
         ccAGTGATCATTCGCAAGGAGTCCAAGCTGGACCGTTTCCCGACGCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCGAAGCAGcccc                                                                                                                                          <  6:169696/149‑2 (MQ=255)
               aTCATTCGCAAGGAGTCCAAGCTGGACCGTTTCCCGACGCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCAc                                                                                                                                     <  2:340690/148‑1 (MQ=255)
                tCATTCGCAAGGAGTCCAAGCTGGACCGTTTCCCGACGCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGa                                                                                                                                   <  1:304137/149‑1 (MQ=255)
                tCATTCGCAAGGAGTCCAAGCTGGACCGTTTCCCGACGCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGa                                                                                                                                   <  5:258790/149‑1 (MQ=255)
                 cATTCGCAAGGAGTCCAAGCTGGACCGTTTCCCGACGCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGa                                                                                                                                   >  5:329664/1‑148 (MQ=255)
                       cAATGAGTCCTAGCTGGACCGTTTCCCGACGCCTGAGCAACTGGATGACCTGCCAAGCCTCCCACGCGGTGAGCGACCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTc                                                                                                                                                  <  6:368824/127‑1 (MQ=255)
                       cAAGGAGTCCAAGCTGGACCGTTTCCCGACGCCTGACCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCAAGAGCTGCTGACCCGCAAGGCCTGATAGAAGAACCGAGGGGGCTGCTTc                                                                                                                                                  >  5:368824/1‑127 (MQ=255)
                             gTCCAAGCTGGACCGTTTCCCGACGCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTc                                                                                                                      >  5:200340/1‑149 (MQ=255)
                                 aaGCTGGACCGTTTCGCGACGACTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGa                                                                                                                                   <  6:15147/132‑1 (MQ=255)
                                       gACCGTTTCCCGACGCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCCGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGaa                                                                                                            >  8:231060/1‑149 (MQ=255)
                                        aCCGTTTCCCGACGCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAAt                                                                                                           >  8:202137/1‑149 (MQ=255)
                                         ccGTTTCCCGACGCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCAc                                                                                                                                     <  5:353990/122‑1 (MQ=255)
                                         ccGTTTCCCGACGCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCAc                                                                                                                                     >  6:353990/1‑122 (MQ=255)
                                            tttCCCGACGCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCt                                                                                                       <  5:369413/149‑1 (MQ=255)
                                                 cGACGCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGAtc                                                                                                  >  4:139492/1‑149 (MQ=255)
                                                  gACGCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCTCGACACAAGGCCGCTCCTGCAGAGAATTTCTTGAtct                                                                                                 >  8:165647/1‑149 (MQ=255)
                                                    cGCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTAGCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATctct                                                                                               <  5:80887/149‑1 (MQ=255)
                                                     gCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACaa                                                                                                                              >  6:11607/1‑117 (MQ=255)
                                                     gCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACaa                                                                                                                              <  5:11607/117‑1 (MQ=255)
                                                       cTGAGCAACTGGATGACCTGCCAAGCCGCTCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCt                                                                                                       <  8:75395/138‑1 (MQ=255)
                                                       cTGAGCAACTGGATGACCTGCCAAGCCGCTCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCt                                                                                                       >  7:75395/1‑138 (MQ=255)
                                                       cTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGacac                                                                                                                                >  4:220115/1‑113 (MQ=255)
                                                       cTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGacac                                                                                                                                <  3:220115/113‑1 (MQ=255)
                                                       cTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACa                                                                                            <  6:329664/149‑1 (MQ=255)
                                                           gCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAg                                                                                        >  1:304712/1‑149 (MQ=255)
                                                            cAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCAATAAACACACAAAGCCCCACCAGCACAGAAATTATTGATGAGTAAAGGAGc                                                                                       >  2:151613/1‑149 (MQ=255)
                                                                 ggATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCt                                                                                  >  6:202256/1‑149 (MQ=255)
                                                                     gACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTtgtg                                                                              >  3:303331/1‑149 (MQ=255)
                                                                       ccTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTc                                                                            <  7:202137/149‑1 (MQ=255)
                                                                        cTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCg                                                                           <  1:151613/149‑1 (MQ=255)
                                                                        cTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCg                                                                           >  1:82439/1‑149 (MQ=255)
                                                                        cTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCg                                                                           >  1:334354/1‑149 (MQ=255)
                                                                        cTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAACGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCg                                                                           >  2:96486/1‑149 (MQ=255)
                                                                         tGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGt                                                                          >  5:326119/1‑149 (MQ=255)
                                                                           ccAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGt                                                                          <  2:334354/147‑1 (MQ=255)
                                                                                 cgccCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGc                                                                  >  7:333272/1‑149 (MQ=255)
                                                                                    ccACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCACTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCg                                                                >  3:301645/1‑148 (MQ=255)
                                                                                    ccACGCGGTGAGCGCCCAGCGCTGCCGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCg                                                                >  1:307318/1‑148 (MQ=255)
                                                                                          ggTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCTCGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAgc‑g                                                        <  7:165647/149‑1 (MQ=255)
                                                                                              aGCGCCCAGCGCTGCTGACCCACAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCCCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAgc‑ggccc                                                    >  6:304227/1‑145 (MQ=255)
                                                                                                   ccAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCGCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCg                                                <  5:299252/149‑1 (MQ=255)
                                                                                                    cAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCgg                                               <  4:303331/149‑1 (MQ=255)
                                                                                                         ctgctgACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTc                                          <  6:200340/149‑1 (MQ=255)
                                                                                                               aCCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCt                                      <  2:304712/147‑1 (MQ=255)
                                                                                                               aCCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTgg                                    >  1:70210/1‑149 (MQ=255)
                                                                                                               aCCCACAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCCCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCt                                      <  5:304227/147‑1 (MQ=37)
                                                                                                                 ccGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTgg                                    >  6:372468/1‑147 (MQ=255)
                                                                                                                        ggCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTa                           <  1:96486/149‑1 (MQ=255)
                                                                                                                         gCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTAc                          <  8:333272/149‑1 (MQ=255)
                                                                                                                             aTAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACtt                        <  1:20545/147‑1 (MQ=255)
                                                                                                                              tAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCt                      <  2:70210/148‑1 (MQ=255)
                                                                                                                                 aaGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCtt                  >  7:217893/1‑149 (MQ=255)
                                                                                                                                   gTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAg                                        <  3:353786/125‑1 (MQ=255)
                                                                                                                                   gTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAg                                        >  4:353786/1‑125 (MQ=255)
                                                                                                                                   gTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTca                >  7:238232/1‑149 (MQ=255)
                                                                                                                                       tgatagggCTGCGAAGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCCCAcc            >  7:30777/6‑149 (MQ=25)
                                                                                                                                          tagggCTGCGAAGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTc         >  4:294647/3‑149 (MQ=37)
                                                                                                                                               ctgctTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGTATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGacac    <  5:202256/149‑1 (MQ=255)
                                                                                                                                                  ctTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTg  <  2:307318/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GAAAAACAGCCAGTGATCATTCGCAAGGAGTCCAAGCTGGACCGTTTCCCGACGCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGC‑GCGCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTG  >  NC_002947/5674513‑5674805

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