Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 535,819 A→C 6.1% E85A (GAA→GCA) ‡ rplA → 50S ribosomal protein L1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947535,8190AC 6.1% 208.5 / 4.6 66E85A (GAA→GCA) ‡rplA50S ribosomal protein L1
Reads supporting (aligned to +/- strand):  ref base A (30/32);  new base C (2/2);  total (32/34)
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.19e-01

CCAAGTTCGTCGAGTCCTACGACATCGCCGTTAACCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTGGGTCAGGTTCTGGGTCCTCG  >  NC_002947/535672‑535966
                                                                                                                                                   |                                                                                                                                                   
ccAAGTTCGTCGAGTCCTACGACATCGCCGTTAACCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGa                                                                                                                                                     <  2:179129/148‑1 (MQ=255)
 cAAGTTCGTCGAGTCCTACGACATCGCCGTTAACCTCGGTGTTGACCCGCGTAAATCCGCCCAGGTCGTACGTAGCGCTACCGTGCAGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAg                                                                                                                                                   >  7:237554/1‑149 (MQ=255)
  aaGTTCGTCGAGTCCTACGACATCGCCGTTAACCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAgc                                                                                                                                                  >  2:59434/1‑149 (MQ=255)
       cgtcgCGTCCTACGACATCGCCGTTAACCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTTCAgctgctg                                                                                                                                             <  4:47104/149‑2 (MQ=255)
         tcgAGTCCTACGACATCGCCGTTAACCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTg                                                                                                                                           >  6:62934/1‑149 (MQ=255)
            aGTCCTACGACATCGCCGTTAACCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCt                                                                                                                                        >  4:7370/1‑149 (MQ=255)
              tCCTACGACATCGCCGTTAACCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGc                                                                                                                                      <  8:201941/149‑1 (MQ=255)
                    gACATCGCCGTTAACCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTAGCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGgcgc                                                                                                                                <  5:76328/149‑1 (MQ=255)
                       aTCGCCGTTAACCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGa                                                                                                                             >  5:24605/1‑149 (MQ=255)
                       aTCGCCGTTAACCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGa                                                                                                                             <  7:286044/149‑1 (MQ=255)
                            cGTTAACCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATcgcg                                                                                                                        <  3:7370/149‑1 (MQ=255)
                                aaCCTCGGTGTTGACCCGCGTAAATCCGAGCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAg                                                                                                                     <  6:95387/148‑1 (MQ=255)
                                  ccTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTa                                                                                                                  >  2:37691/1‑149 (MQ=255)
                                            gACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCt                                                                                                        >  7:305337/1‑149 (MQ=255)
                                             aCCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTg                                                                                                       >  3:241300/1‑149 (MQ=255)
                                               ccGCGTAAATCCCACCAGGTCGTACGTAGCGCTACCGTGCTGCCCCACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTa                                                                                                                      <  4:35384/132‑1 (MQ=255)
                                                  cgTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGc                                                                                                  <  1:37691/149‑1 (MQ=255)
                                                   gTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGcc                                                                                                 >  6:7572/1‑149 (MQ=255)
                                                      aaTCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGa                                                                                               >  6:13318/1‑148 (MQ=255)
                                                       aTCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGaaa                                                                                             <  4:274377/149‑1 (MQ=255)
                                                        tCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAAt                                                                                            <  4:241300/149‑1 (MQ=255)
                                                         ccGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAAt                                                                                            >  2:34126/1‑148 (MQ=255)
                                                            aCCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGaaa                                                                                        <  1:59434/149‑1 (MQ=255)
                                                             ccAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAg                                                                                       >  1:125735/1‑149 (MQ=255)
                                                              cAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTTCAgctgctgga                                                                                                                                           <  5:59476/96‑4 (MQ=255)
                                                              cAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTTCAgctgctgga                                                                                                                                           >  6:59476/1‑93 (MQ=255)
                                                                       cgtaGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCt                                                                             <  7:68068/149‑1 (MQ=255)
                                                                        gtaGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCt                                                                             >  3:223420/1‑148 (MQ=255)
                                                                        gtaGCGCTACCGTGCAGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGGCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTg                                                                            <  8:237554/149‑1 (MQ=255)
                                                                          aGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATcgc                                                                                                                         <  7:308448/102‑1 (MQ=255)
                                                                          aGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATcgc                                                                                                                         >  8:308448/1‑102 (MQ=255)
                                                                           gcgcTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGaa                                                                          <  1:34126/148‑1 (MQ=255)
                                                                            cgcTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAAc                                                                         >  4:264332/1‑148 (MQ=255)
                                                                            cgcTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAAc                                                                         >  5:273176/1‑148 (MQ=255)
                                                                             gcTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTa                                                                       <  3:264332/149‑1 (MQ=255)
                                                                                 ccGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTTAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGa                                                                    <  8:305337/148‑1 (MQ=255)
                                                                                 ccGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATg                                                                     <  4:223420/147‑1 (MQ=255)
                                                                                 ccGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGa                                                                    <  6:273176/148‑1 (MQ=255)
                                                                                  cGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTg                                                                                                       >  1:145184/1‑112 (MQ=255)
                                                                                  cGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTg                                                                                                       <  2:145184/112‑1 (MQ=255)
                                                                                    tgctgcCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGAcgtc                                                                <  3:171803/149‑1 (MQ=255)
                                                                                      ctgaCACACGGCACTGGCAAGACCCTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGAcgtcgt                                                              <  4:42762/149‑1 (MQ=255)
                                                                                       tgcCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGAcgtcgt                                                              >  5:238505/1‑148 (MQ=255)
                                                                                         ccACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTa                                                            <  5:7572/148‑1 (MQ=255)
                                                                                           acacGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACt                                                                        <  8:178624/134‑1 (MQ=255)
                                                                                           acacGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACt                                                                        >  7:178624/1‑134 (MQ=255)
                                                                                             acGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTatcgca                                                       >  2:195438/1‑149 (MQ=255)
                                                                                                      gCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTatcgcatc                                                     <  1:208755/142‑1 (MQ=255)
                                                                                                      gCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTatcgcatc                                                     >  2:208755/1‑142 (MQ=255)
                                                                                                       cAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATg                                             <  2:125735/149‑1 (MQ=255)
                                                                                                                    gTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCACCGCCTGATGCCATGCGCGTTGt                                <  7:206192/149‑1 (MQ=255)
                                                                                                                        cTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAgg                             <  6:24605/148‑1 (MQ=255)
                                                                                                                             ttCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTTATCGCGTAGGTATGGACGATCTGGCTGCCGAAATg                                                                                           <  7:165162/81‑1 (MQ=255)
                                                                                                                             ttCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATg                                                                                           >  8:165162/1‑81 (MQ=255)
                                                                                                                             ttCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATCGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCTCGTTGTAGGTCAGCt                       >  8:92015/1‑149 (MQ=255)
                                                                                                                                  ccAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTatcgcatcgc                                                   >  2:52183/1‑116 (MQ=255)
                                                                                                                                  ccAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTatcgcatcgc                                                   <  1:52183/116‑1 (MQ=255)
                                                                                                                                  ccAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTGGGt                   >  4:213431/1‑148 (MQ=255)
                                                                                                                                  ccAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTGGGt                   >  5:112978/1‑148 (MQ=255)
                                                                                                                                   cAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTGGGt                   >  6:267224/1‑147 (MQ=255)
                                                                                                                                       ccagagcagcTTCAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTGGGTCAGGt              >  4:172617/5‑148 (MQ=255)
                                                                                                                                        tCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTGGGTCAGGTTc            >  5:29568/1‑149 (MQ=255)
                                                                                                                                           agcagcTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTGGGTCAGGTTCTgg         <  5:62934/149‑1 (MQ=255)
                                                                                                                                           agcagcTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTGGGTCAGGTTCTgg         <  1:195438/149‑1 (MQ=255)
                                                                                                                                            gcagcTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTa                               <  3:213631/126‑1 (MQ=255)
                                                                                                                                            gcagcTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTa                               >  4:213631/1‑126 (MQ=255)
                                                                                                                                                  gAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTGGGTCAGGTTGTGGGTCCTcg  >  7:57235/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CCAAGTTCGTCGAGTCCTACGACATCGCCGTTAACCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTACCGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTGGGTCAGGTTCTGGGTCCTCG  >  NC_002947/535672‑535966

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