Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 607,608 T→G 5.3% L171R (CTG→CGG)  ribA → GTP cyclohydrolase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947607,6080TG 5.3% 220.3 / 2.0 70L171R (CTG→CGG) ribAGTP cyclohydrolase II
Reads supporting (aligned to +/- strand):  ref base T (36/30);  new base G (2/2);  total (38/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGATGCTCGAGCACCTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGTGCCCCAGGCGTGACCCGCGCCGAGGTCAAGCGGCGCCTGGCGCTGGCCTGGTGG  >  NC_002947/607461‑607755
                                                                                                                                                   |                                                                                                                                                   
gCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGATGCTCGAGCACCTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTg                                                                                                                                                    <  4:122624/149‑1 (MQ=255)
      cTGCCGACCAGCGTGACTACGCCATATGCCTGCCGATGCTCGAGCACCTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACAcc                                                                                                                                              >  3:85287/1‑149 (MQ=255)
         ccGACCAGCGTGACTACGCCATATGCCTGCCGATGCTCGAGCACCTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCggg                                                                                                                                           <  4:85287/149‑1 (MQ=255)
          cGACCAGCGTGACTACGCCATATGCCTGCCGATGCTCGAGCACCTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCAGCGGCacacgctcg                                                                                                                                          <  6:283910/149‑6 (MQ=255)
               aGCGTGACTACGCCATATGCCTGCCGATGCTCGAGCACCTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCa                                                                                                                                         <  4:247237/145‑1 (MQ=255)
               aGCGTGACTACGCCATATGCCTGCCGATGCTCGAGCACCTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCa                                                                                                                                         >  3:247237/1‑145 (MQ=255)
               aGCGTGACTACGCCATATGCCTGCCGATGCTCGAGCACCTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAAc                                                                                                                                     <  3:141109/149‑1 (MQ=255)
                     aCTACGCCATATGCCTGCCGATGCTCGAGCACCTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGca                                                                                                                                <  4:177277/148‑1 (MQ=255)
                          gCCATATGCCTGCCGATGCTCGAGCACCTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATa                                                                                                                                       >  5:280085/1‑136 (MQ=255)
                          gCCATATGCCTGCCGATGCTCGAGCACCTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATa                                                                                                                                       <  6:280085/136‑1 (MQ=255)
                           ccATATGCCTGCCGATGCTCGAGCACCTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCAGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACCCCGGGCATAACCCGCACAACCGc                                                                                                                         >  7:232107/1‑149 (MQ=255)
                                 gcctgccGATGCTCGAGCACCTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCtactac                                                                                                                   >  6:330159/1‑149 (MQ=255)
                                  cctgccGATGCTCGAGCACCTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCtactac                                                                                                                   <  5:293961/148‑1 (MQ=255)
                                  cctgccGATGCTCGAGCACCTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTAcc                                                                                                                  >  7:171716/1‑149 (MQ=255)
                                        gATGCTCGAGCACCTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGa                                                                                                            >  6:215791/1‑149 (MQ=255)
                                           gCTCGAGCACCTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCg                                                                                                                          <  7:317708/132‑1 (MQ=255)
                                           gCTCGAGCACCTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCg                                                                                                                          >  8:317708/1‑132 (MQ=255)
                                               gAGCACCTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAg                                                                                                       <  8:171716/147‑1 (MQ=255)
                                                aGCACCTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGcc                                                                                                    >  6:60677/1‑149 (MQ=255)
                                                   aCCTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCgg                                                                                                  >  7:104845/1‑148 (MQ=255)
                                                     cTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGc                                                                                                 >  5:149976/1‑147 (MQ=255)
                                                      tGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAg                                                                                              >  3:120419/1‑149 (MQ=255)
                                                      tGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCAGCGGCacacgctc                                                                                                                                           <  5:125495/104‑5 (MQ=255)
                                                      tGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCAGCGGCacacgctc                                                                                                                                           >  6:125495/1‑100 (MQ=255)
                                                          cGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTgg                                                                                          <  6:149976/149‑1 (MQ=255)
                                                           gTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTggg                                                                                         <  1:2383/149‑1 (MQ=255)
                                                           gTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACA‑CGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTCGGc                                                                                        >  7:167567/1‑149 (MQ=255)
                                                                gTCGCTGCGCCTTATGACCAACAATCCGCGCAAGGGCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCcaca                                                                                    <  8:232107/149‑1 (MQ=255)
                                                                gTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCcaca                                                                                    >  2:26572/1‑149 (MQ=255)
                                                                 tCGCTGCGCCTGATGACCATCAATCCGCGCAAGGTCAATGGCCTGACCGACATGATCAGCGTGGTTGCCGAGCGTGTGCCGCTGCACA‑CGGGCATAACCCGCACAACCGGTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACAt                                                                                   <  8:167567/148‑1 (MQ=255)
                                                                         ccTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGc                                                                            >  3:79484/1‑148 (MQ=255)
                                                                          cTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGAc                                                                                                           >  1:354254/1‑116 (MQ=255)
                                                                          cTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGAc                                                                                                           <  2:354254/116‑1 (MQ=255)
                                                                           tgatgaCCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAAc                                                                         >  2:239136/1‑149 (MQ=255)
                                                                               gaCCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGc                                                                                                                                          <  8:181190/80‑1 (MQ=255)
                                                                               gaCCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGc                                                                                                                                          >  7:181190/1‑80 (MQ=255)
                                                                               gaCCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCg                                                                                                                          >  7:317784/1‑96 (MQ=255)
                                                                               gaCCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCg                                                                                                                          <  8:317784/96‑1 (MQ=255)
                                                                                   aacaaTCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTg                                                                                           <  3:2246/123‑1 (MQ=255)
                                                                                   aacaaTCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTg                                                                                           >  4:2246/1‑123 (MQ=255)
                                                                                   aacaaTCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAacac                                                                   >  1:144423/1‑147 (MQ=255)
                                                                                    acaaTCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAg                                                                >  5:219147/1‑149 (MQ=255)
                                                                                    acaaTCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGGGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAg                                                                <  3:45537/149‑1 (MQ=255)
                                                                                        tCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTgg                                                                                          <  5:159304/119‑1 (MQ=255)
                                                                                        tCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTgg                                                                                          >  6:159304/1‑119 (MQ=255)
                                                                                          cgcgCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGa                                                           <  6:219147/148‑1 (MQ=255)
                                                                                          cgcgCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGa                                                           <  5:330159/148‑1 (MQ=255)
                                                                                            cgcAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGt                                                        <  8:104845/149‑1 (MQ=255)
                                                                                                  gTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGTGccc                                                    >  4:223782/1‑147 (MQ=255)
                                                                                                       aGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACAAACACCAGGGCGAGGTGCCCCAGGCGt                                             >  2:307924/1‑149 (MQ=255)
                                                                                                                    aTGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCGACGAACACCAGGGCGAGGTGCCCCAGGCGTGACCCGCGCCGAg                                >  3:349093/1‑149 (MQ=255)
                                                                                                                      gAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGTGCCCCAGGCGTGACCCGCGCCGAGGt                              >  2:120146/1‑149 (MQ=255)
                                                                                                                      gAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGTGCCCCAGGCGTGACCCGCGCCGAGGt                              >  4:307308/1‑149 (MQ=255)
                                                                                                                         cATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGTGCCCCAGGCGTGACCCGCGCCGAGGTc                             >  3:200813/1‑147 (MQ=255)
                                                                                                                            cGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGTGCCCCAGGCGTGACCCGCGCCGAGGTCAAgcg                        >  5:176010/1‑149 (MQ=255)
                                                                                                                             gTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGTGCCCCAGGCGTGACCCGCGCCGAGGTCAAgcgg                       >  3:89697/1‑149 (MQ=255)
                                                                                                                                    ccGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGTGCCCCAGGCGTGACCCGCGCCGAGGTCAAGCGGCGCCTg                 <  5:60677/148‑1 (MQ=255)
                                                                                                                                      gAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGTGCCCCAGGCGTGACCCGCGCCGAGGTCAAGCGGCGCCTGgcg              >  2:323387/1‑149 (MQ=255)
                                                                                                                                      gAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGTGCCCCAGGCGTGACCCGCGCCGAGGTCAAGCGGCGCCTGgcg              >  4:336529/1‑149 (MQ=255)
                                                                                                                                       ccggtgtGCAGCGGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACCCCAGGGCGAGGTGCCCCAGGCGTGACCCGCGCCGAGGTCAAGCGGCGCCTGgc               >  4:286220/4‑147 (MQ=255)
                                                                                                                                         cGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGTGCCCCAGGCGTGACCCGCGCCGAGGTCAAGCGGCGCCTGGCGCTg           <  4:79484/149‑1 (MQ=255)
                                                                                                                                           tgtgCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGTGCCCCAGGCGTGACCCGCGCCGAGGTCAAGCGGCGCCTGGCGCTGGc         <  1:323387/149‑1 (MQ=255)
                                                                                                                                            gtgCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGt                                                                                             <  5:233151/64‑1 (MQ=255)
                                                                                                                                            gtgCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGt                                                                                             >  6:233151/1‑64 (MQ=255)
                                                                                                                                            gtgCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGTGCCCCAGGCGTGACCCGCGCCGAGGTCAAGCGGCGCCTGGCGCTGGcc        >  6:147456/1‑149 (MQ=255)
                                                                                                                                              gCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGTGCCCCAGGCGTGACCCGCGCCGAGGTCAAGCGGCGCCTGGCGCTGGCCtg      <  4:200813/149‑1 (MQ=255)
                                                                                                                                               ccGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGTGCCCCAGGCGTGACCCGCGCCGAGGTCAAGCGGCGCCTGGCGCTGGCCtg      <  4:120419/148‑1 (MQ=255)
                                                                                                                                                  cTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGTGCCCCAGGCGTGACCCGCGCCGAGGTCAAGCGGCGCCTGGCGCTGGCCtggtgg  <  1:239136/149‑1 (MQ=255)
                                                                                                                                                  cTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGTGCCCCAGGCGTGACCCGCGCCGAGGTCAAGCGGCGCCTGGCGCTGGCCt       <  6:51499/144‑1 (MQ=255)
                                                                                                                                                  cTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGTGCCCCAGGCGTGACCCGCGCCGAGGTCAAGCGGCGCCTGGCGCTGGCCt       >  5:51499/1‑144 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCTTCGCTGCCGACCAGCGTGACTACGCCATATGCCTGCCGATGCTCGAGCACCTGGGCGTCAAGTCGCTGCGCCTGATGACCAACAATCCGCGCAAGGTCAAAGCCCTGACCGACATGAACATCGTGGTTGCCGAGCGTGTGCCGCTGCACACCGGGCATAACCCGCACAACCGCTACTACCTGGCGACCAAGGCCGGCAAGTTGGGCCACATGCTGGGCAACGAACACCAGGGCGAGGTGCCCCAGGCGTGACCCGCGCCGAGGTCAAGCGGCGCCTGGCGCTGGCCTGGTGG  >  NC_002947/607461‑607755

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