Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 853,763 T→C 9.1% R62R (CGT→CGC) ‡ prmC → release factor‑(glutamine‑N5) methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947853,7630TC 9.1% 176.9 / 7.5 66R62R (CGT→CGC) ‡prmCrelease factor‑(glutamine‑N5) methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (30/30);  new base C (2/4);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.52e-01

GCCTGAATCGCCCACTGAGCGCCTGGATGCCGAACTGCTGCTGGCTGCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGCGCCGCATACCCTGATCCCGCGGCCGGATACCGAGCTGCTGGTCGAGACGGCGCTTGAACTGCAACCCGCTTCGC  >  NC_002947/853616‑853902
                                                                                                                                                   |                                                                                                                                           
gCCTGAATCGCCCACTGAGCGCCTGGATGCCGAACTGCTGCTGGCTGCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTg                                                                                                                                            >  3:335446/1‑149 (MQ=255)
  cTGAATCGCCCACTGAGCGCCTGGATGCCGAACTGCTGCTGGCTGCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTggc                                                                                                                                          >  6:26792/1‑149 (MQ=255)
     aaTCGCCCACTGAGCGCCTGGATGCCGAACTGCTGCTGGCTGCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTggcggc                                                                                                                                       >  4:286843/1‑149 (MQ=255)
       tCGCCCACTGAGCGCCTGGATGCCGAACTGCTCCTGGCTGCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGa                                                                                                                                     <  6:264315/149‑1 (MQ=255)
         gCCCACTGAGCGCCTGGATGCCGAACTGCTGCTGGCTGCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGc                                                                                                                                   >  3:343703/1‑149 (MQ=255)
          cccACTGAGCGCCTGGATGCCGAACTGCTGCTGGCTGCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGc                                                                                                                                   >  3:37483/1‑148 (MQ=255)
          cccACTGAGCGCCTGGATGCCGAACTGCTGCTGGCTGCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCACGgcggcgcttc                                                                                                                                  <  8:104326/149‑5 (MQ=255)
            cACTGAGCGCCTGGATGCCGAACTGCTGCTGGCTGCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCACGgcggcgcttcag                                                                                                                                <  4:93423/149‑7 (MQ=255)
                     ccTGGATGCCGAACTGCTGCTGGCTGCGGCCATCGGTAAATCGCGCAGCTACATGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCAGTGGCGGCGAGCCAGTGGCCtata                                                                                                                       >  3:367865/1‑149 (MQ=255)
                      cTGGATGCCGAACTGCTGCTGGCTGCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTatat                                                                                                                      >  4:250063/1‑149 (MQ=255)
                          aTGCCGAACTGCTGCTGGCTGCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAATATGCCGAGCCATATGCCGGCTACCTGAATCGCCTCCGTGGCGGCGAGCCAGTGGCCTCTATCAtt                                                                                                                  >  8:164313/1‑149 (MQ=255)
                             ccGAAGTGGTGCTGGCTGCGGCCATGGTTATATCTCGCTGTTACTTCCATACATGGCCCGAGCGTATTGTCAGCTGCGAAGATGCCGAGACCTATGCGGGCTACCTGAAGCGCGGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTgg                                                                                                                <  4:367865/148‑1 (MQ=255)
                              cGAACTGCTGCTGGCTGCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGc                                                                                                              >  5:141970/1‑149 (MQ=255)
                                aaCTGCTGCTGGCTGCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTg                                                                                                            <  5:26792/149‑1 (MQ=255)
                                   tgctgctgGCTGCGGCCATCGGTAAATCCCTCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAg                                                                                                         <  7:164313/149‑1 (MQ=255)
                                      tgctgGCTGCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGc                                                                                                      >  2:196769/1‑149 (MQ=255)
                                        ctggctgCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCtt                                                                                                    >  8:289760/1‑149 (MQ=255)
                                           gctgCGGCCATCGGTAAATCGGGCAGCTACCTGCATACCTGGCCCTATCGAATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTAGCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGTCCTATATCCTTGGGCTGCAGCGCTTc                                                                                                   <  4:263918/147‑1 (MQ=255)
                                               cGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGa                                                                                               <  4:335446/147‑1 (MQ=255)
                                               cGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAg                                                                                             <  3:250063/149‑1 (MQ=255)
                                                    aTCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGa                                                                                        <  6:141970/149‑1 (MQ=255)
                                                     tCGGTAACTCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGAc                                                                                       <  8:348711/149‑1 (MQ=255)
                                                       ggTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCACGgcggcgcttca                                                                                                                                 >  1:29479/1‑100 (MQ=255)
                                                       ggTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCACGgcggcgcttca                                                                                                                                 <  2:29479/105‑6 (MQ=255)
                                                         tAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTgg                                                                                                <  7:275921/136‑1 (MQ=255)
                                                         tAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTgg                                                                                                >  8:275921/1‑136 (MQ=255)
                                                          aaaTCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCACGgcggcgc                                                                                                                                     <  8:168144/98‑2 (MQ=255)
                                                          aaaTCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCACGgcggcgc                                                                                                                                     >  7:168144/1‑97 (MQ=255)
                                                           aaTCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGa                                                                                  <  4:343703/148‑1 (MQ=255)
                                                           aaTCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAg                                                                                 >  7:276059/1‑149 (MQ=255)
                                                                     cTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGcgccg                                                                        >  2:365735/1‑148 (MQ=255)
                                                                      tACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGCGCCGCa                                                                      >  3:95161/1‑149 (MQ=255)
                                                                        ccTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGCGCCGCAt                                                                     <  3:286843/148‑1 (MQ=255)
                                                                                 cTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGCGCCGCATACCCTGATcc                                                           >  4:27055/1‑149 (MQ=255)
                                                                                        gAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCT‑GAAGATCGACCTGGAGGTCGCGCCGCATACCCTGATCCCGCGGCCg                                                   >  4:264509/1‑149 (MQ=255)
                                                                                         aGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGGGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGCGCCGCATACCCTGATCCCGCGGCCg                                                   <  4:74119/149‑1 (MQ=255)
                                                                                              aTTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGCGCCGCATACCCTGATCCCGCGGCCGGATAc                                              >  3:14479/1‑149 (MQ=255)
                                                                                                  tCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTg                                                                                                 >  2:348041/1‑94 (MQ=255)
                                                                                                  tCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTg                                                                                                 <  1:348041/94‑1 (MQ=255)
                                                                                                  tCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGaa                                                                                              <  7:88291/97‑1 (MQ=255)
                                                                                                  tCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGaa                                                                                              >  8:88291/1‑97 (MQ=255)
                                                                                                  tCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGCGCCGCATACCCTGATCCCGCGGCCGGATACCGAg                                          >  1:187596/1‑149 (MQ=255)
                                                                                                  tCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGCGCCGCATACCCTGATCCCGCGGCCGGATACCGAg                                          >  4:118895/1‑149 (MQ=255)
                                                                                                  gcagcagCGAAGATGCCGAGACCTATGCCGGCTACCGGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGCGCCGCATACCCTGATCCCGCGGCCGGATACCGAg                                          <  8:155289/148‑1 (MQ=255)
                                                                                                     gcagcGAAGATGCTGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGa                                                                                  <  5:37976/106‑1 (MQ=255)
                                                                                                     gcagcGAAGATGCAGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGa                                                                                  >  6:37976/1‑106 (MQ=255)
                                                                                                      cagcGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTggg                                                                                                               <  8:178864/76‑1 (MQ=255)
                                                                                                      cagcGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTggg                                                                                                               >  7:178864/1‑76 (MQ=255)
                                                                                                      cagcGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGCGCCGCATACCCTGATCCCGCGGCCGGATACCGAgctgc                                      <  8:276059/149‑1 (MQ=255)
                                                                                                        gcGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGCGCCGCATACCCTGATCCCGCGGCCGGATACCGAGctgctg                                    <  2:187596/149‑1 (MQ=255)
                                                                                                        gcGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGCGCCGCATACCCTGATCCCGCGGCCGGATACCGAGctgctg                                    >  2:193537/1‑149 (MQ=255)
                                                                                                             gATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGCGCCGCATACCCTGATCCCGCGGCCGGATACCGAGCTGCTGGTCga                               <  7:289760/149‑1 (MQ=255)
                                                                                                             gATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGCGCCGCATACCCTGATCCCGCGGCCGGATACCGAGCTGCTGGTCga                               <  4:95161/149‑1 (MQ=255)
                                                                                                                  cGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGCGCCGCATACCCTGATCCCGCGGCTGGATACCGAGCTGCTGGTCGAGACgg                          <  1:196769/149‑1 (MQ=255)
                                                                                                                       ccTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGCGCCGCATACCCTGATCCCGCGGCCGGATACCGAGCTGCTGGTCGAGACGGCGCtt                     <  4:37483/149‑1 (MQ=255)
                                                                                                                       ccTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCT‑GAAGATCGACCTGGAGGTCGCGCCGCATACCCTGATCCCGCGGCCGGATACCGAGCTGCTGGTCGAGACGGCGCtt                     <  3:264509/148‑1 (MQ=255)
                                                                                                                           tGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGCGCCGCATACCCTGATCCCGCGGCCGGATACCGAGCTGCTGGTCGAGACGGCGCTTGaa                  <  1:193537/148‑1 (MQ=255)
                                                                                                                           tGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGCGCCGCATACCCTGATCCCGCGGCCGGATACCGAGCTGCTGGTCGAGACGGCGCTTGAAc                 >  7:42496/1‑149 (MQ=255)
                                                                                                                              cGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGCGCCGCATACCCTGATCCCGCGGCCGGATACCGAGCTGCTGGTCGAGACGGCGCTTGAACTGc              <  1:10100/149‑1 (MQ=255)
                                                                                                                                gCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGCGCCGCATACCCTGATCCCGCGGCCGGATACCGAGCTGCTGGTCGAGACGGCGCTTGAACTGCaa            <  1:365735/149‑1 (MQ=255)
                                                                                                                                  ctgcTGTAGCGCCGCCGTGGCGGCTCGCCTGTGTCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGCGCCTCATACCCTGATCGCGCGGCCGGATACCGAGCTGCTGGTCGAGACTGCGCTTGAACTGCAAc           <  6:235433/145‑1 (MQ=255)
                                                                                                                                     cTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCATTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGCGCCGCATACCCTGATCCCGCGGCCGGATACCGAGCTGCTGGTCGAGACGGCGCTTGAACTGCAAc           >  5:235433/1‑145 (MQ=255)
                                                                                                                                       gAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGCCCTGGAGGTCGCGCCGCATACCCTGATCCCGCGGCCGGATACCGAGCTGCTGGTCGAGACGGCGCTTGAACTGCAACCCGCtt     >  7:261323/1‑149 (MQ=255)
                                                                                                                                          gcgcCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGCGCCGCATACCCTGATCCCGCGGCCGGATACCGAGCTGCTGGTCGAGACGGCGCTTGAACTGCAACCCGCTTcgc  >  4:158554/1‑149 (MQ=255)
                                                                                                                                                ccgTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTcgc                                                                            >  6:12958/1‑69 (MQ=255)
                                                                                                                                                ccgTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTcgc                                                                            <  5:12958/69‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                           
GCCTGAATCGCCCACTGAGCGCCTGGATGCCGAACTGCTGCTGGCTGCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGCGCCGCATACCCTGATCCCGCGGCCGGATACCGAGCTGCTGGTCGAGACGGCGCTTGAACTGCAACCCGCTTCGC  >  NC_002947/853616‑853902

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑