Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 604,469 C→T 7.1% P101S (CCA→TCA) ‡ nusB → transcription antitermination protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947604,4690CT 7.1% 203.8 / 4.3 71P101S (CCA→TCA) ‡nusBtranscription antitermination protein
Reads supporting (aligned to +/- strand):  ref base C (34/31);  new base T (3/2);  total (38/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.30e-01

AAGCGCAATTCCGGGTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTA  >  NC_002947/604326‑604612
                                                                                                                                               |                                                                                                                                                
aaGCGCAATTCCGGGTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTc                                                                                                                                             <  1:5827/149‑1 (MQ=255)
     cAATTCCGGGTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAgct                                                                                                                                        <  2:35884/149‑1 (MQ=255)
        ttCCGGGTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAgctggc                                                                                                                                     <  5:343568/149‑1 (MQ=255)
        ttCCGGGTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAgctggc                                                                                                                                     >  5:103192/1‑149 (MQ=255)
           cGGGTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACTGGGTCGAgctcttcca                                                                                                                                  <  1:309721/149‑7 (MQ=255)
              gTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGtgct                                                                                                                               >  3:112293/1‑149 (MQ=255)
               cgttctaCGATTTCACCGATGTAGACTGTTCCTATTTCCGCTAGATCCTGCATGGTGTCCCGTCAATCATTGGGGTAATCGACAAATCGCTGTTGCCTTGCATGTCCATCGCTCTGGAAGAGCTCGACCCAGTCGTGCTGTCTGtgctg                                                                                                                              <  4:163917/143‑1 (MQ=255)
                  aTAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTgcgc                                                                                                                           <  7:22832/149‑1 (MQ=255)
                      cGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTAGACCCAGTCGAGCTGGCTGtgctgcagatcg                                                                                                                       >  3:163917/1‑143 (MQ=255)
                             aCCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTg                                                                                                                >  1:153513/1‑149 (MQ=255)
                                      gACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTggt                                                                                                       >  2:235637/1‑147 (MQ=255)
                                       aCGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGAGCTGCGCCTGTCCACCTGGGAGCTGAt                                                                                                       >  6:83550/1‑148 (MQ=255)
                                            gCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAAcg                                                                                                 >  3:352071/1‑149 (MQ=255)
                                              cTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAAcgcg                                                                                               >  3:228378/1‑149 (MQ=255)
                                                   tCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGAGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGa                                                                                           <  5:83550/148‑1 (MQ=255)
                                                        gagaTCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCaaa                                                                                                   >  8:107141/1‑135 (MQ=255)
                                                        gagaTCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCaaa                                                                                                   <  7:107141/135‑1 (MQ=255)
                                                         agaTCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCt                                                                                    >  6:42860/1‑149 (MQ=255)
                                                          gaTCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACTGGGTCGAgctct                                                                                                                                      >  5:107653/1‑96 (MQ=255)
                                                          gaTCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACTGGGTCGAgctct                                                                                                                                      <  6:107653/98‑3 (MQ=255)
                                                          gaTCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCt                                                                                    >  6:379700/1‑148 (MQ=255)
                                                              cTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTAcc                                                                                <  4:228378/148‑1 (MQ=255)
                                                              cTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACcg                                                                               >  2:43793/1‑149 (MQ=255)
                                                                gCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACcgcg                                                                             >  6:122887/1‑149 (MQ=255)
                                                                      ggTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGcgtcgt                                                                         <  1:167955/147‑1 (MQ=255)
                                                                      ggTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGcgtcgt                                                                         >  2:167955/1‑147 (MQ=255)
                                                                        tCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATc                                                                     <  2:153513/149‑1 (MQ=255)
                                                                                  cAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGgtg                                                           >  2:35427/1‑149 (MQ=255)
                                                                                   aaGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGTAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGgtgt                                                          <  4:112293/149‑1 (MQ=255)
                                                                                   aaGGGCGAAATCGACAAAGCGCTGGAGCCTTGCATGACCATCGCACTGGAAGAGCTCGAAACAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCTTCGACGTGCCTTACCGCGTCGTCATCAACGAAGgtgt                                                          >  5:154512/1‑149 (MQ=255)
                                                                                    aGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGtgtg                                                         <  1:330023/149‑1 (MQ=255)
                                                                                    aGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGtgtg                                                         <  1:35427/149‑1 (MQ=255)
                                                                                           aaTCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAg                                                              >  4:160730/1‑137 (MQ=255)
                                                                                           aaTCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTgg                                                 <  6:103192/149‑1 (MQ=255)
                                                                                                aCAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAg                                            >  3:296595/1‑149 (MQ=255)
                                                                                                aCAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAg                                            <  4:177111/149‑1 (MQ=255)
                                                                                                 cAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAgg                                           <  4:352071/149‑1 (MQ=255)
                                                                                                  aaaGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGtgtg                                                         <  3:135214/135‑1 (MQ=255)
                                                                                                  aaaGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGtgtg                                                         >  4:135214/1‑135 (MQ=255)
                                                                                                  aaaGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGt                                          >  1:185574/1‑149 (MQ=255)
                                                                                                  aaaGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGt                                          >  5:255763/1‑149 (MQ=255)
                                                                                                    aGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCt                                        >  6:148210/1‑149 (MQ=255)
                                                                                                      cgcTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTc                                      <  5:42860/149‑1 (MQ=255)
                                                                                                           gTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGgcgc                                 >  6:111872/1‑149 (MQ=255)
                                                                                                           gTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGgcgc                                 >  7:343654/1‑149 (MQ=255)
                                                                                                             gCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCa                               >  2:149152/1‑149 (MQ=255)
                                                                                                                  gCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGa                           <  5:379700/148‑1 (MQ=255)
                                                                                                                  gCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGAc                          >  7:113928/1‑149 (MQ=255)
                                                                                                                     gcaCCAGCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGc                       <  3:197936/147‑1 (MQ=255)
                                                                                                                      gACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTgg                                                       >  6:258386/1‑116 (MQ=255)
                                                                                                                      gACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTgg                                                       <  5:258386/116‑1 (MQ=255)
                                                                                                                      gACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATAAACGAAGGTGTGGAACTGGACAAGGTCTACGGCGCCACTTACGGc                       >  7:231619/1‑148 (MQ=255)
                                                                                                                      gACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCATCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGc                       >  2:217105/1‑148 (MQ=255)
                                                                                                                       aCCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGCGAACTGGCCAAGGTCTTCGGCGCCACTGACGGc                       <  1:43793/148‑1 (MQ=255)
                                                                                                                        ccATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCcac                    <  4:296595/149‑1 (MQ=255)
                                                                                                                          aTCGCACTGGAAGAGCTCGACCCAGTCGACCAGGCTGTGCTGCCCCTGTCCACCACGGAGCTGAACAAACCCGACGACGAGCCTTACCGCGTCGGGATCAACGAAGGTGTGGAACTGGCCAAGGTCGTCGGCGCCACTGCCGGCcacac                  >  8:283707/1‑148 (MQ=255)
                                                                                                                              cACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTc              <  8:113928/149‑1 (MQ=255)
                                                                                                                              cACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTc              <  5:122887/149‑1 (MQ=255)
                                                                                                                               aCTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCg             >  3:240750/1‑149 (MQ=255)
                                                                                                                                   gAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACaa                  <  6:189238/140‑1 (MQ=255)
                                                                                                                                   gAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACaa                  >  5:189238/1‑140 (MQ=255)
                                                                                                                                   gAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTc           <  2:185574/147‑1 (MQ=255)
                                                                                                                                      gagCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACaa                  >  5:215252/1‑137 (MQ=255)
                                                                                                                                      gagCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACaa                  <  6:215252/137‑1 (MQ=255)
                                                                                                                                       agCTCGACTGGGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGc     >  3:322627/1‑149 (MQ=255)
                                                                                                                                       agCTCGACTGGGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGc     >  5:292600/1‑149 (MQ=255)
                                                                                                                                        gCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGGTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCg    <  6:154512/149‑1 (MQ=255)
                                                                                                                                        gCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGaa                                                     <  4:19857/100‑1 (MQ=255)
                                                                                                                                        gCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGaa                                                     >  3:19857/1‑100 (MQ=255)
                                                                                                                                         cTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCg    >  2:278165/1‑148 (MQ=255)
                                                                                                                                            gACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGGTCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTa  <  1:235637/147‑1 (MQ=255)
                                                                                                                                              cccAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCa                               >  5:163573/1‑116 (MQ=255)
                                                                                                                                              cccAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCa                               <  6:163573/116‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                
AAGCGCAATTCCGGGTGGATAACGATTTCACCGATGTCGACGGTGCCTATTTCCGCGAGATCCTGCATGGGGTCCCGGCAATCAAGGGCGAAATCGACAAAGCGCTGGTGCCTTGCATGACCATCGCACTGGAAGAGCTCGACCCAGTCGAGCTGGCTGTGCTGCGCCTGTCCACCTGGGAGCTGATCAAACGCGTCGACGTGCCTTACCGCGTCGTGATCAACGAAGGTGTGGAACTGGCCAAGGTCTTCGGCGCCACTGACGGCCACAAGTTCGTCAACGGCGTA  >  NC_002947/604326‑604612

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: