Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 475,314 T→C 5.5% Q64R (CAG→CGG)  rpsU ← 30S ribosomal protein S21

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947475,3140TC 5.5% 270.5 / 3.6 91Q64R (CAG→CGG) rpsU30S ribosomal protein S21
Reads supporting (aligned to +/- strand):  ref base T (38/47);  new base C (3/2);  total (42/49)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

ATGATGCAGTTCGCATGAAGCTTGTGCAGGGCCTTGAGCTAGCTGCTCAAGTTTGAAACGACCACTGCCATCGCCCGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCTTTGCGCTTACGCTCTGCGGTCGGCTTCTCGTAAAAC‑TCGCGGCTACGAACTTCAGCCAGTACACCGGCTTTTTCGCAGGAGCGCTTGAAACGACGCAGA  >  NC_002947/475168‑475460
                                                                                                                                                  |                                                                                                                                                   
atgatgCAGTTCGCATGAAGCTTGTGCAGGGCCTTGAGCTAGCTGCTCAAGTTTGAAACGACCACTGCCATCGCCCGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGt                                                                                                                                                   >  4:146057/1‑149 (MQ=255)
        gTTCGCATGAAGCTTGTGCAGGGCCTTGAGCTAGCTGCTCAAGTTTGAAACGACCACTGCCATCGCCCGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTcgc                                                                                                                                               <  8:302261/145‑1 (MQ=255)
        gTTCGCATGAAGCTTGTGCAGGGCCTTGAGCTAGCTGCTCAAGTTTGAAACGACCACTGCCATCGCCCGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTcgc                                                                                                                                               >  7:302261/1‑145 (MQ=255)
        gTTCGCATGAAGCTTGTGCAGGGCCTTGAGCTAGCTGCTCAAGTTTGAAACGACCACTGCCATCGCCCGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACAGCGCCGCGCCGTTcgtctg                                                                                                                                            >  3:49132/1‑144 (MQ=255)
         ttagCATGAAGCTTGTGCAGGGCCTTGAGCTAGCTGCTCAAGTTTGAAAAGACCACTGCCATCGCCCGGGGGTTTGACGGGCAGGGCAGAAGGTTTGCGTTGAACGCCTGGATTATTAAAGACGAACGGCGCGGCGCGGTTCGCGGTg                                                                                                                                           >  6:52495/4‑148 (MQ=255)
         ttCGCATGAAGCTTGTGCAGGGCCTTGAGCTAGCTGCTCAAGTTTGAAACGACCACTGCCATCGCCCGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACAGCGCCGCGCCGTTcgt                                                                                                                                               >  5:60518/1‑143 (MQ=255)
         ttCGCATGAAGCTTGTGCAGGGCCTTGAGCTAGCTGCTCAAGTTTGAAACGACCACTGCCATCGCCCGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACAGCGCCGCGCCGTTcgt                                                                                                                                               <  6:60518/144‑2 (MQ=255)
            gCATGAAGCTTGTGCAGGGCCTTGAGCTAGCTGCTCAAGTTTGAAACGACCACTGCCATCGCCCGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACt                                                                                                                                       <  8:242779/149‑1 (MQ=255)
              aTGAAGCTTGTGCAGGGCCTTGAGCTAGCTGCTCAAGTTTGAAACGACCACTGCCATCGCCCGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACtt                                                                                                                                      >  8:55933/1‑148 (MQ=255)
               tGAAGCTTGTGCAGGGCCTTGAGCTAGCTGCTCAAGTTTGAAACGACCACTGCCATCGCCCGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTc                                                                                                                                    <  6:329579/149‑1 (MQ=255)
               tGAAGCTTGTGCAGGGCCTTGAGCTAGCTGCTCAAGTTTGAAACGACCACTGCCATCGCCCGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTc                                                                                                                                    <  2:226834/149‑1 (MQ=255)
                gAAGCTTGTGCATGGCCTTGAGCTAGCTGCTCAAGTTTGAAACTACCACTGCCATCTCCCGGGGGTTAGACGGGCAGGGCAGAAGCTTTGCGTTGAACGCATGTATTAGTACAGACGAATGTCGCGGCGCTGTTAGCTCTGAACTTTCt                                                                                                                                   >  8:111363/1‑149 (MQ=255)
                gAAGCTTGTGCAGGGCCTTGAGCTAGCTGCTCAAGTTTGAAACGACCACTGCCATCGCCCGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCt                                                                                                                                   >  6:20694/1‑149 (MQ=255)
                gAAGCTTGTGCAGGGCCTTGAGCTAGCTGCTCAAGTTTGAAACGACCACTGCCATCGCCCGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCt                                                                                                                                   >  5:191203/1‑149 (MQ=255)
                gAAGCTTGTGCAGGGCCTTGAGCTAGCTGCTCAAGTTTGAAACGACCACTGCCATCGCCCGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCt                                                                                                                                   >  4:308053/1‑149 (MQ=255)
                 aaGCTTGTGCAGGGCCTTGAGCTAGCTGCTCAAGTTTGAAACGACCACTGCCATCGCCCGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCt                                                                                                                                   <  3:115884/148‑1 (MQ=255)
                 aaGCTTGTGCAGGGCCTTGAGCTAGCTGCTCAAGTTTGAAACGACCACTGCCATCGCCAGGGGGTTAGCGGGGCAGGGCAGAAGCTTTGCGTTGAACTCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCt                                                                                                                                   <  5:152005/148‑1 (MQ=255)
                    cTTGTGCAGGGCCTTGAGCTAGCTGCTCAAGTTTGAAACGACCACTGCCATCGCCCGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGc                                                                                                                               <  3:146057/149‑1 (MQ=255)
                     ttGTGCAGGGCCTTGAGCTAGCTGCTCAAGTTTGAAACGACCACTGCCATCGCCCGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCg                                                                                                                              >  8:143486/1‑149 (MQ=255)
                        tgCAGGGCCTTGAGCTAGCTGCTCAAGTTTGAAACGACCACTGCCATCGCCCGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGa                                                                                                                           <  1:350458/149‑1 (MQ=255)
                          cAGGGCCTTGAGCTAGCTGCTCAAGTTTGAAACGACCACTGCCATCGCCCGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGAc                                                                                                                          <  6:342152/148‑1 (MQ=255)
                          cAGGGCCTTGAGCTAGCTGCTCAAGTTTGAAACGACCACTGCCATCGCCCGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACg                                                                                                                         >  2:179606/1‑149 (MQ=255)
                          cAGGGCCTTGAGCTAGCTGCTCAAGTTTGAAACGACCACTGCCATCGCCCGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACAGCGCCGCGCCGTTcgt                                                                                                                                               >  8:6027/1‑126 (MQ=255)
                          cAGGGCCTTGAGCTAGCTGCTCAAGTTTGAAACGACCACTGCCATCGCCCGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACAGCGCCGCGCCGTTcgt                                                                                                                                               <  7:6027/127‑2 (MQ=255)
                            gggCCTTGAGCTAGCTGCTCAAGTTTGAAACGACCACTGCCATCGCCCGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACttt                                                                                                                                     <  7:51036/135‑1 (MQ=255)
                            gggCCTTGAGCTAGCTGCTCAAGTTTGAAACGACCACTGCCATCGCCCGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACttt                                                                                                                                     >  8:51036/1‑135 (MQ=255)
                               ccTTGAGCTAGCTGCTCAAGTTTGAAACGACCACTGCCATCGCCCGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTaa                                                                                                                    <  5:323711/149‑1 (MQ=255)
                                 ttGAGCTAGCTGCTCAAGTTTGAAACGACCACTGCCATCGCCCGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTaa                                                                                                                    >  3:333996/1‑147 (MQ=255)
                                       tagctGCTCAAGTTTGAAACGACCACTGCCATCGCCCGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAAcagcag                                                                                                              >  5:38722/1‑147 (MQ=255)
                                         gctgctCAAGTTTGAAACGACCACTGCCATCGCCCGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACagcagcagc                                                                                                          >  3:52737/1‑149 (MQ=255)
                                          ctgctCAAGTTTGAAACGACCACTGCCATCGCCCGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAgct                                                                                                         >  3:132523/1‑149 (MQ=255)
                                           tgctCAAGTTTGAAACGACCACTGCCATCGCCCGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAgctg                                                                                                        >  8:69125/1‑149 (MQ=255)
                                            gctCAAGTTTGAAACGACCACTGCCATCGCCCGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCGTGTATTAGTACAGACGACCGGCGCGGCGCTGTTCGCGCGGAACTTTCTTGGCGTGACGCTTAACAGCAGCCgctga                                                                                                       >  6:92779/1‑148 (MQ=255)
                                            gctCAAGTTTGAAACGACCACTGCCATCGCCCGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGATGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAgctgc                                                                                                       >  6:87509/1‑149 (MQ=255)
                                              tCAAGTTTGAAACGACCACTGCCATCGCCCGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCtt                                                                                                     >  8:29862/1‑149 (MQ=255)
                                              tCAAGTTTGAAACGACCACTGCCATCGCCCGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCtt                                                                                                     >  3:340716/1‑149 (MQ=255)
                                                     tGAAACGACCACTGCCATCGCCCGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCTTTGCGCt                                                                                               >  3:34943/1‑148 (MQ=255)
                                                       aaaCGACCACTGCCATCGCCCGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCTTTGCGCTTa                                                                                             <  5:52495/148‑1 (MQ=255)
                                                        aaCGACCACTGCCATCGCCCGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCTTTGCGCTTACg                                                                                           <  7:22343/149‑1 (MQ=255)
                                                        aaCGACCACTGCCATCGCCCGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCTTTGCGCTTACg                                                                                           >  8:200190/1‑149 (MQ=255)
                                                            aCCACTGCCATCGCCCGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCTTTGCGCTTACGctct                                                                                       <  7:143486/149‑1 (MQ=255)
                                                             cctcTGCCCTCTCCCGGGGGTTAGCCGGGCAGGGCAGAAGCGTTGCGGTGAACGCCTGTCTTAGTACAGACGATCGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCTTTGCGCTTACGCTCTg                                                                                      <  8:234027/146‑1 (MQ=255)
                                                              cACTGCCATCGCCCGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCTTTGCGCTTACGCTCTGc                                                                                     <  4:52737/149‑1 (MQ=255)
                                                                   ccATCGCCCGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCTTTGCGCTTACGCTCTGCGGTCg                                                                                >  8:76879/1‑149 (MQ=255)
                                                                    cATCGCCCGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTCTAATAGTACAGACGAACGCCGCGGCGCTGTTCGCGCTGAACGTTCTTGGCTTGACGCTTAACAGCAGCAGCTGCTTTGCGCTTACCCTCTCCGGTAgg                                                                               >  6:274182/1‑149 (MQ=255)
                                                                     aTCGCCCGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCTTTGCGCTTACGCTCTGCGGTCGGc                                                                              <  6:191203/149‑1 (MQ=255)
                                                                     aTCGCCCGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCTTTGCGCTTACGCTCTGCGGTCGGc                                                                              <  1:179606/149‑1 (MQ=255)
                                                                     aTCGCCCGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCTTTGCGCTTACGCTCTGCGGTCGGc                                                                              >  8:46242/1‑149 (MQ=255)
                                                                        gCCCGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGAc                                                                                                                          >  6:34647/1‑102 (MQ=255)
                                                                        gCCCGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGAc                                                                                                                          <  5:34647/102‑1 (MQ=255)
                                                                        gCCCGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCTTTGCGCTTACGCTCTGCGGTCGGCTtc                                                                           <  7:76879/149‑1 (MQ=255)
                                                                         cccGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAAcagcagc                                                                                                             <  8:28011/114‑1 (MQ=255)
                                                                         cccGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAAcagcagc                                                                                                             >  7:28011/1‑114 (MQ=255)
                                                                          ccGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCTTTGCGCTTACTCTCTGCGGTCGGCTtctc                                                                         <  5:32999/149‑1 (MQ=255)
                                                                          ccGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCTTTGCGCTTACGCTCTGCg                                                                                    >  5:3741/1‑138 (MQ=255)
                                                                          ccGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCTTTGCGCTTACGCTCTGCg                                                                                    <  6:3741/138‑1 (MQ=255)
                                                                              gggTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCTTTGCGCTTACGCTCTGAGGTCGGCTTCTCGTa                                                                      >  8:157535/1‑148 (MQ=255)
                                                                               gggtatccGTGCATGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAATGGCGCGGCGCTGTTCGCTCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCTTTGCGCTTACGCTCTGCGGTCGGCTTCTCGTaaa                                                                    <  7:111363/143‑1 (MQ=255)
                                                                               ggTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCTTTGCGCTTACGCTCTGCGGTCGGCTTCTCGTaaa                                                                    <  5:20694/149‑1 (MQ=255)
                                                                               ggTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCTTTGCGCTTACGCTCTGCGGTCGGCTTCTCGTaa                                                                     <  2:322623/148‑1 (MQ=255)
                                                                               ggTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCTTTGCGCTTACGCTCTGCGGTCGGCTTCTCGTaa                                                                     >  1:322623/1‑148 (MQ=255)
                                                                                     ccGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCTTTGCGCTTACGCTCTGCGGTCGGCTTCTCGTAAAAC‑Tcgc                                                             <  4:34943/149‑1 (MQ=255)
                                                                                         gcagggcagAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCTTTGCGCTTACGCTCTGCGGTCGGCTTCTCGTAAAAC‑TCGCGGc                                                          >  7:240129/1‑148 (MQ=255)
                                                                                              gcagAAGCTTTGCGTTGCACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCTTTGCGCTTACGCTCTGCGGGCGGCTTCGAGTAAAAC‑TCGCGGCTACGa                                                     >  7:328214/1‑148 (MQ=255)
                                                                                              gcagAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCTTTGCGCTTACGCTCTGCGGTCGGCTTCTCGTAAAAC‑TCGCGGCTACGaa                                                    >  2:93988/1‑149 (MQ=255)
                                                                                               cagAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCTTTGCGCTTACGCTCTGCGGTCGGCTTCTCGTAAAAC‑TCGCGGCTACGAAc                                                   <  7:55933/149‑1 (MQ=255)
                                                                                                  aaGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCTTTGCGCTTACGCTCTGCGGTCGGCTTCTCGTAAAAC‑TCGCGGCTACGAACtt                                                 >  2:1610/1‑148 (MQ=255)
                                                                                                   aGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCTTTGCGCTTACGCTCTGCGGTCGGCTTCTCGTAAAAC‑TCGCGGCTACGAACTTCa                                               <  4:132523/149‑1 (MQ=255)
                                                                                                        tGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCTTTGCGCTTACGCTCTGCGGTCGGCTTCTCGTAAAAC‑TCGCGGCTACGAACTTCAGCCa                                           <  6:38722/148‑1 (MQ=255)
                                                                                                        tGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCTTTGCGCTTACGCTCTGCGGTCGGCTTCTCGTAAAAC‑TCGCGGCTACGAACTTCAGCCa                                           <  7:69125/148‑1 (MQ=255)
                                                                                                         gCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCTTTGCGCTTACGCTCTGCGGTCGGCTTCTCGTAAAAC‑TCGCGGCTACGAACTTCAGCCAg                                          <  7:29862/148‑1 (MQ=255)
                                                                                                           gTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCTTTGCGCTTACGCTCTGCGGTCGGCTTCTCGTAAAAC‑TCGCGGCTACGAACTTCAGCCAGTa                                        >  6:334724/1‑148 (MQ=255)
                                                                                                              gAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCTTTGCGCTTACGCTCTGCGGTCGGCTTCTCGTAAACC‑TCGCGGCTACGAACTTCAGCCAGTACAcc                                    >  6:280739/1‑149 (MQ=255)
                                                                                                              gAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCTTTGCGCTTACGCTCTGCGGTCGGCTTCTCGTAAAAC‑TCGCGGCTACGAACTTCAGCCAGTacac                                     <  4:340716/148‑1 (MQ=255)
                                                                                                               aaCGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCTTTGCGCTTACGCTCTGCGGTCGGCTTCTCGTAAAAC‑TCGCGGCTACGAACTTCAGc                                             <  1:313731/139‑1 (MQ=255)
                                                                                                               aaCGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCTTTGCGCTTACGCTCTGCGGTCGGCTTCTCGTAAAAC‑TCGCGGCTACGAACTTCAGc                                             >  2:313731/1‑139 (MQ=255)
                                                                                                                       tATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCTTTGCGCTTACGCTCTGCGGTCGGCTTCTCGTAAAAC‑TCGCGGCTACGAACTTCAGCCAGTACACCGGCttttt                            <  8:121537/148‑1 (MQ=255)
                                                                                                                       tATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCTTTGCGCTTACGCTCTGCGGTCGGCTTCTCGTAAAAC‑TCGCGGCTACGAACTTCAGCCAGTACACCGGCTTTTTc                           <  3:176977/149‑1 (MQ=255)
                                                                                                                         ttAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCTTTGCGCTTACGCTCTGCGGTCGGCTTCTCGTAAAAC‑TCGCGGCTACGAACTTCAGCCAGTACACCGGCTTTTTCg                          <  5:87509/148‑1 (MQ=255)
                                                                                                                         ttAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCTTTGCGCTTACGCTCTGCGGTCGGCTTCTCGTAAAAC‑TCGCGGCTACGAACTTCAGCCAGTACACCGGCTTTTTCg                          <  5:92779/148‑1 (MQ=255)
                                                                                                                           aGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCTTTGCGCTTACGCTCTGCGGTCGGCTTCTCGTAAAAC‑TCGCGGCTACGAACTTCAGCCAGTACACCGGCTTTTTCGCa                        <  8:240129/148‑1 (MQ=255)
                                                                                                                            gTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCTTTGCGCTTACGCTCTGCGGTCGGCTTCTCGTAAAAC‑TCGCGGCTACGAACTTCAGCCAGTACACCGGCTTTTTCGCAg                       >  7:39460/1‑148 (MQ=255)
                                                                                                                            gTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCTTTGCGCTTACGCTCTGCGGTCGGCTTCTCGTAAAAC‑TCGCGGCTACGAACTTCAGCCAGTACACCGGCTTTTTCGCAg                       <  5:274182/148‑1 (MQ=255)
                                                                                                                              aCAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCTTTGCGCTTACGCTCTGCGGTCGGCTTCTCGTAAAAC‑TCGCGGCTACGAACTTCAGCCa                                           >  5:287380/1‑126 (MQ=255)
                                                                                                                                   cGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCTTTGCGCTTACGCTCTGCGGTCGGCTTCTCGTAAAAC‑TCGCGGCTACGAACTTCAGCCAGTACACCGGCTTTTTCGCAGGAGCGCtt               <  3:118592/149‑1 (MQ=255)
                                                                                                                                      aCGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCTTTGCGCTTACGCTCTGCGGTCGGCTTCTCGTAAAAC‑TCGCGGCTACGAACTTCAGCCAGTACACCGGCTTTTTCGCAGGAGCGCTTGaa            <  7:46242/149‑1 (MQ=255)
                                                                                                                                       cggcgcggcgCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCTTTGCGCTTACGCTCTGCGGTCGGCTTCTCGTAAAACTTCGCGGCTACGAACTTCAGCCAGTACACCGGCTTTTTCGCAGGAGCGCTTGaa            <  1:93988/149‑1 (MQ=255)
                                                                                                                                       cggcgcggcgCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCTTTGCGCTTACGCTCTGCGGTCGGCTTCTCGTAAAAC‑TCGCGGCTACGAACTTCAGCCAGTACACCGGCTTTTTCGCAGGAGCGCTTGaaa           <  4:351911/149‑1 (MQ=255)
                                                                                                                                             ggcgcTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCTTTGCGCTTACGCTCTGCGGTCGGCTTCTCGTAAAAC‑TCGCGGCTACGAACTTCAGCCAGTACACCGGCTTTTTCGCAGGAGCGCTTGAAACGACGc     <  4:333996/149‑1 (MQ=255)
                                                                                                                                                gcTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCTTTGCGCTTACGCTCTGCGGTCGGCTTCTCGTAAAAC‑TCGCGGCTACGAACTTCAGCCAGTACACCGGCTTTTTCGCAGGAGCGCTTGAAACGACGCaga  <  7:200190/149‑1 (MQ=255)
                                                                                                                                                gcTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCTTTGCGCTTACGCTCTGCGGTCGGCTTCTCGTAAAAC‑TCGCGGCTACGAACTTCAGCCAGTACACCGGCTTTTTCGCAGGAGCGCTTGAAACGACGCaga  <  1:1610/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
ATGATGCAGTTCGCATGAAGCTTGTGCAGGGCCTTGAGCTAGCTGCTCAAGTTTGAAACGACCACTGCCATCGCCCGGGGGTTAGCCGGGCAGGGCAGAAGCTTTGCGTTGAACGCCTGTATTAGTACAGACGAACGGCGCGGCGCTGTTCGCGCTGAACTTTCTTGGCGTGACGCTTAACAGCAGCAGCTGCTTTGCGCTTACGCTCTGCGGTCGGCTTCTCGTAAAAC‑TCGCGGCTACGAACTTCAGCCAGTACACCGGCTTTTTCGCAGGAGCGCTTGAAACGACGCAGA  >  NC_002947/475168‑475460

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