Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 604,470 C→G 5.4% P101R (CCA→CGA)  nusB → transcription antitermination protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947604,4700CG 5.4% 279.2 / 3.7 86P101R (CCA→CGA) nusBtranscription antitermination protein
Reads supporting (aligned to +/- strand):  ref base C (45/36);  new base G (3/2);  total (48/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.30e-01

CGAAGCGCAATTCCGGGTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTGG  >  NC_002947/604324‑604616
                                                                                                                                                  |                                                                                                                                                  
cGAAGCGCAATTCCGGGTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAg                                                                                                                                                  >  7:7309/1‑149 (MQ=255)
  aaGCTCAATTCAGGGTTGATAACGATTTCACCGATGTCGACGGTGGCTATTTCCGCGAGATCCTGCATGGGGTCCGGGCAATCAAGGTCGAAATCGACAAAGCGCTGGTGCCTTGCATGACGATCGCACTGGAAGAGCTCGACCCAGTc                                                                                                                                                <  7:93416/149‑1 (MQ=255)
  aaGCGCAATTCCGGGTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTc                                                                                                                                                >  6:226217/1‑149 (MQ=255)
      gcAATTCCGGGTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAgc                                                                                                                                            <  6:34265/149‑1 (MQ=255)
       cAATTCCGGGTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAgct                                                                                                                                           <  8:244315/149‑1 (MQ=255)
         aTTCCGGGTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACTGGGTCGAgctct                                                                                                                                         <  6:162677/149‑3 (MQ=255)
              gggTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCtgt                                                                                                                                     >  7:44735/1‑148 (MQ=255)
              gggTGGATAACGATTTCACCGATGCCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCtgtg                                                                                                                                    >  2:34899/1‑149 (MQ=255)
                      aaCGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTgcgc                                                                                                                              >  1:310957/1‑147 (MQ=255)
                      aaCGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTgcgc                                                                                                                              <  2:310957/147‑1 (MQ=255)
                               aCCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTc                                                                                                                         <  4:251580/143‑1 (MQ=255)
                               aCCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTc                                                                                                                         >  3:251580/1‑143 (MQ=255)
                               aCCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCt                                                                                                                    <  1:60651/148‑1 (MQ=255)
                               aCCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTg                                                                                                                   >  2:366750/1‑149 (MQ=255)
                               aCCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTg                                                                                                                   >  6:169171/1‑149 (MQ=255)
                                   aTGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAg                                                                                                               <  5:226217/149‑1 (MQ=255)
                                        gACGGTGCCTATTTGCGCGAGATCCTGGATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGGGCTGGTGGCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCGACCTGGGAGCTGa                                                                                                           <  1:10883/148‑1 (MQ=255)
                                        gACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGAt                                                                                                          >  5:348652/1‑149 (MQ=255)
                                        gACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGGGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGa                                                                                                           <  8:156344/148‑1 (MQ=255)
                                              gCCTATTTCCGCGAGCTCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTCGTGCCTTGCATGACCATCGCACTGCAAGAGCTCGACCCAGTCGAGCTGCCTGTGCTGCGACTGTCCACATGGGAGCTGATCAAAcg                                                                                                    >  2:112195/1‑149 (MQ=255)
                                               ccTATTTCCGCGAGATCCTGCATGGGTTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGTTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTTCTGCGCCTGTGCACCTGGGAGCTGATCAAAcgc                                                                                                   <  7:61774/149‑1 (MQ=255)
                                                cTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCa                                                                                                                       <  2:103601/128‑1 (MQ=255)
                                                cTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCa                                                                                                                       >  1:103601/1‑128 (MQ=255)
                                                  aTTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTc                                                                                                >  5:372646/1‑149 (MQ=255)
                                                  aTTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGACTTGCATGACCATCGCACTCGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTc                                                                                                >  7:182599/1‑149 (MQ=255)
                                                     tCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGtgctgc                                                                                                                                >  7:198014/1‑114 (MQ=255)
                                                     tCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGtgctgc                                                                                                                                <  8:198014/114‑1 (MQ=255)
                                                        gcgAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTg                                                                                          <  7:145052/149‑1 (MQ=255)
                                                                cTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTAcc                                                                                   <  6:135766/148‑1 (MQ=255)
                                                                cTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTAcc                                                                                   <  8:44735/148‑1 (MQ=255)
                                                                cTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACcg                                                                                  >  3:265670/1‑149 (MQ=255)
                                                                  gCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTg                                                                                          <  1:391709/139‑1 (MQ=255)
                                                                  gCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTg                                                                                          >  2:391709/1‑139 (MQ=255)
                                                                  gCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACcgcg                                                                                >  8:313682/1‑149 (MQ=255)
                                                                      ggtgTCCCGGCTATCAAGGGCGAAATCGACAAAGCGCTTGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGcgtcg                                                                             <  8:182599/145‑1 (MQ=255)
                                                                        ggTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGa                                                                          <  4:383798/149‑1 (MQ=255)
                                                                          tCCCGGCAATCAAGGGCGAAATTGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGAt                                                                         >  6:7596/1‑148 (MQ=255)
                                                                            ccGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCaa                                                                      >  3:327767/1‑149 (MQ=255)
                                                                            ccGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCaa                                                                      >  2:356325/1‑149 (MQ=255)
                                                                              ggCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACg                                                                    <  1:34899/149‑1 (MQ=255)
                                                                                 aaTCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAg                                                                 >  7:234262/1‑149 (MQ=255)
                                                                                 aaTCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAg                                                                 >  2:88192/1‑149 (MQ=255)
                                                                                     aaGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGgtgt                                                             >  2:319033/1‑149 (MQ=255)
                                                                                      aGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGtgtg                                                            >  3:259639/1‑149 (MQ=255)
                                                                                       gggCGAAATCGACAAAGCGCTGGTGCGTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTTCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTTATCAACGAAGGTGTgg                                                           <  1:112195/149‑1 (MQ=255)
                                                                                          cGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTgg                                                                                                                  <  5:170188/91‑1 (MQ=255)
                                                                                          cGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTgg                                                                                                                  >  6:170188/1‑91 (MQ=255)
                                                                                             aaTCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTgg                                                     <  1:153083/149‑1 (MQ=255)
                                                                                              aTCGACAAAGCGCTGGTGCCTTGCATGTCCTTCGCACTGGAAGAGCTCGACCGAGTCGTGCTGGCTGTGCTGCGCCTGTACACCTGGGAGCTGATCAAACGCGTCGAGGTGCGTTAACGCGTCGTGATCAACGATGGTGTGGAACTGGc                                                    <  4:268847/149‑1 (MQ=255)
                                                                                              aTCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGc                                                    <  1:366750/149‑1 (MQ=255)
                                                                                              aTCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGc                                                    >  5:250662/1‑149 (MQ=255)
                                                                                                  aCAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAg                                                >  1:49791/1‑149 (MQ=255)
                                                                                                  aCAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAg                                                <  2:363836/149‑1 (MQ=255)
                                                                                                            ggTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGgcg                                      >  5:236795/1‑149 (MQ=255)
                                                                                                                  ttGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACt                                 >  8:146923/1‑148 (MQ=255)
                                                                                                                  ttGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACt                                 >  4:199754/1‑148 (MQ=255)
                                                                                                                  ttGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACt                                 >  6:289475/1‑148 (MQ=255)
                                                                                                                   tGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCg                                                                                               >  4:349616/1‑85 (MQ=255)
                                                                                                                   tGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCg                                                                                               <  3:349616/85‑1 (MQ=255)
                                                                                                                    gCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGAc                              >  1:198240/1‑149 (MQ=255)
                                                                                                                     cATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATc                                                                        <  7:39740/106‑1 (MQ=255)
                                                                                                                     cATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATc                                                                        >  8:39740/1‑106 (MQ=255)
                                                                                                                     cATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACg                             <  1:356325/149‑1 (MQ=255)
                                                                                                                        gACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGc                           >  2:214801/1‑148 (MQ=255)
                                                                                                                         aCCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTgg                                                     <  5:323879/121‑1 (MQ=255)
                                                                                                                         aCCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTgg                                                     >  6:323879/1‑121 (MQ=255)
                                                                                                                          ccATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCca                         <  4:259639/148‑1 (MQ=255)
                                                                                                                           cATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAAc                                                                                                     <  1:349238/71‑1 (MQ=255)
                                                                                                                           cATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAAc                                                                                                     >  2:349238/1‑71 (MQ=255)
                                                                                                                            aTCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTTCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACaa                      <  6:176445/149‑1 (MQ=255)
                                                                                                                            aTCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCcac                        >  6:343180/1‑147 (MQ=255)
                                                                                                                                 aCTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGgcgc                                     <  3:191741/129‑1 (MQ=255)
                                                                                                                                 aCTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGgcgc                                     >  4:191741/1‑129 (MQ=255)
                                                                                                                                 aCTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCg                 >  3:119568/1‑149 (MQ=255)
                                                                                                                                   tGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGt                >  8:110463/1‑148 (MQ=255)
                                                                                                                                    cagccagCTCGACTGGGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTc               >  3:194115/6‑148 (MQ=255)
                                                                                                                                     gAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCaa             >  6:308455/1‑149 (MQ=255)
                                                                                                                                     gAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCaa             >  5:52214/1‑149 (MQ=255)
                                                                                                                                     gAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCaa             <  5:308455/149‑1 (MQ=255)
                                                                                                                                     agccagCTCGACTGGGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCaa             >  2:237200/5‑149 (MQ=255)
                                                                                                                                       ccagCTCGACTGGGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACg           >  5:233516/3‑149 (MQ=255)
                                                                                                                                          gCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCg        <  5:7596/149‑1 (MQ=255)
                                                                                                                                          gCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCg        >  1:190968/1‑149 (MQ=255)
                                                                                                                                               aCCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTg   <  2:49791/149‑1 (MQ=255)
                                                                                                                                               aCCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTg   <  8:234262/149‑1 (MQ=255)
                                                                                                                                               aCCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTg   >  3:386225/1‑149 (MQ=255)
                                                                                                                                                cccAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAATGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTg   >  5:354406/1‑148 (MQ=255)
                                                                                                                                                cccAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCGCAAGTTCGTCAACGGCGTACTgg  <  6:250662/149‑1 (MQ=255)
                                                                                                                                                cccAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTgg  <  2:198240/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CGAAGCGCAATTCCGGGTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTACTGG  >  NC_002947/604324‑604616

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: