Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 451,272 G→A 8.6% R201C (CGT→TGT) ‡ PP_0371 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947451,2720GA 8.6% 162.3 / 4.9 58R201C (CGT→TGT) ‡PP_0371LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (23/30);  new base A (2/3);  total (25/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

TCTTCGGCGGCCTCCTTGGCGCGCCCCGGCAGAATGGTCGCCAGCGCTGTGTGCGCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGCGGGTGGGTATCGGGTATGACCAGTTGCAGCTCATCTTCGTACAACAGCAGGCCATGCAGCCCTGGCTGG  >  NC_002947/451124‑451420
                                                                                                                                                    |                                                                                                                                                    
tcTTCGGCGGCCTCCTTGGCGCGCCCCGGCAGAATGGTCGCCAGCGCTGTGTGCGCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACg                                                                                                                                                      >  1:1219/1‑149 (MQ=255)
   tCGGCGGCCTCCTTGGCGCGCCCCGGCAGAATGGTCGCCAGCGCTGTGTGCGCGAGGCTATCGAGAATCCCTGCCATATTATTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGAt                                                                                                                                                   >  4:246106/1‑149 (MQ=255)
    cggcggCCTCCTTGGCGCGCCCCGGCAGAATGGTCGCCAGCGCTGTGTGCGCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGAt                                                                                                                                                   <  5:275277/148‑1 (MQ=255)
      gcggcCTCCTTGGCGCGCCCCGGCAGAATGGTCGCCAGCGCTGTGTGCGCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTg                                                                                                                                                <  5:50377/149‑1 (MQ=255)
       cggcCTCCTTGGCGCGCCCCGGCAGAATGGTCGCCAGCGCTGTGTGCGCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTg                                                                                                                                                >  2:192595/1‑148 (MQ=255)
         gcCTCCTTGGCGCGCCCCGGCAGAATGGTCGCCAGCGCTGTGTGCGCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCCGTCAGATCTggca                                                                                                                                             <  4:113796/149‑3 (MQ=255)
           ctcctTGGCGCGCCCCGGCAGAATGGTCGCCAGCGCTGTGTGCGCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAAct                                                                                                                                           >  1:61896/1‑149 (MQ=255)
            tcctTGGCGCGCCCCGGCAGAATGGTCGCCAGCGCTGTGTGCGCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCCGTCAGATCTggca                                                                                                                                             >  5:110901/1‑144 (MQ=255)
            tcctTGGCGCGCCCCGGCAGAATGGTCGCCAGCGCTGTGTGCGCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCCGTCAGATCTggca                                                                                                                                             <  6:110901/146‑3 (MQ=255)
             cctTGGCGCGCCCCGGCAGAATGGTCGCCAGCGCTGTGTGCGCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAActcc                                                                                                                                         >  3:84359/1‑149 (MQ=255)
                  gcgcgcCCCGGCAGAATGGTCGCCAGCGCTGTGTGCGCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGCTCTGACGGATCTGGCACTCCTCGcc                                                                                                                                    >  4:225012/1‑149 (MQ=255)
                  gcgcgcCCCGGCAGAATGGTCGCCAGCGCTGTGTGCGCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGc                                                                                                                                     <  5:96275/148‑1 (MQ=255)
                    gcgcCCCGGCAGAATGGTCGCCAGCGCTGTGTGCGCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGcc                                                                                                                                    >  5:102872/1‑147 (MQ=255)
                       ccccGGCAGAATGGTCGCCAGCGCTGTGTGCGCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCcagca                                                                                                                               >  2:47051/1‑149 (MQ=255)
                        cccGGCAGAATGGTCGCCAGCGCTGTGTGCGCGAGGCTATCGAGAATCCCTGCCAGATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCcatgct                                                                                                                              >  1:276041/1‑145 (MQ=255)
                            gCAGAATGGTCGCCAGCGCTGTGTGCGCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCcagcatcagc                                                                                                                          <  2:6067/149‑1 (MQ=255)
                               gAATGGTCGCCAGCGCTGTGTGCGCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATg                                                                                                                       <  1:192595/149‑1 (MQ=255)
                                  tGGTCGCCAGCGCTGTGTGCGCGTGGCTATCGAGAATCCCTGCCAGATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGAGCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATggg                                                                                                                     >  5:77357/1‑148 (MQ=255)
                                  tGGTCGCCAGCGCTGTGTGCCCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGc                                                                                                                    <  6:42862/149‑1 (MQ=255)
                                    gTCGCCAGCGCTGTGTGCGCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCaa                                                                                                                  <  5:316797/149‑1 (MQ=255)
                                            cgcTGTGTGCGCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTcggcgg                                                                                                          <  2:57764/149‑1 (MQ=255)
                                                   tgCGCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACg                                                                                                   <  7:254174/149‑1 (MQ=255)
                                                     cgcgAGGCTATCGAGAATCCCTGCCATATTATTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACgg                                                                                                  <  3:246106/148‑1 (MQ=255)
                                                      gcgAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCTGCGGCCTGACGGAt                                                                                                >  7:258656/1‑149 (MQ=255)
                                                       cgAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATg                                                                                               <  2:296408/149‑1 (MQ=255)
                                                          ggCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATg                                                                                               <  3:298206/146‑1 (MQ=255)
                                                          ggCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATg                                                                                               >  4:298206/1‑146 (MQ=255)
                                                               tCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGGCCTGACTGGCCAGTTGCGCCTGACAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACAcc                                                                                         <  4:257321/147‑1 (MQ=255)
                                                                     aTCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGaa                                                                                 <  6:102872/149‑1 (MQ=255)
                                                                     aTCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGaa                                                                                 <  1:47051/149‑1 (MQ=255)
                                                                     aTCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGaa                                                                                 <  2:1219/149‑1 (MQ=255)
                                                                         cTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGa                                                                                                     >  4:57306/1‑125 (MQ=255)
                                                                         cTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGa                                                                                                     <  3:57306/125‑1 (MQ=255)
                                                                         cTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGAc                                                                                                    >  2:257380/1‑126 (MQ=255)
                                                                         cTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGAc                                                                                                    <  1:257380/126‑1 (MQ=255)
                                                                         cTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTcc                                                                             >  6:9506/1‑149 (MQ=255)
                                                                          tGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGa                                                                                                 <  6:304571/128‑1 (MQ=255)
                                                                          tGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGa                                                                                                 >  5:304571/1‑128 (MQ=255)
                                                                          tGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCt                                                                            >  6:49062/1‑149 (MQ=255)
                                                                                         ctctGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGCGGGTGGGTAt                                                             >  6:31201/1‑149 (MQ=255)
                                                                                              cctgcACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGCGGGTGGGTATCgg                                                          <  5:9506/147‑1 (MQ=255)
                                                                                              cctgcACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGCGGGTGGGTATCgg                                                          <  2:61896/147‑1 (MQ=255)
                                                                                                 gcACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGCGGGTGGGTATCGGGTAt                                                      >  5:246028/1‑148 (MQ=255)
                                                                                                  cACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGCGGGTGGGTATCGGGTAt                                                      <  5:31201/147‑1 (MQ=255)
                                                                                                   aCCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGCGGGTGGGTATCGGGTATGAc                                                   >  7:116941/1‑149 (MQ=255)
                                                                                                     cTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGc                                                                       >  1:43443/1‑127 (MQ=255)
                                                                                                     cTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGc                                                                       <  2:43443/127‑1 (MQ=255)
                                                                                                     cTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGCGGGTGGGTATCGGGTATGAc                                                   >  5:203930/1‑147 (MQ=255)
                                                                                                     cTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGCGGGTGGGTATCGGGTATGACCa                                                 <  2:314280/149‑1 (MQ=255)
                                                                                                        cGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGCGGGTGGGTATCGGGTATGACCAg                                                <  6:203930/147‑1 (MQ=255)
                                                                                                            cGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGCGGGTGGGTATCGGGTATGACCAGTTGCAg                                          <  5:49062/149‑1 (MQ=255)
                                                                                                                ggCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGCGGGTGGGTATCGGGTATGACCAGTTGCAGCTCa                                      >  3:267192/1‑149 (MQ=255)
                                                                                                                            ccAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGCgggtgg                                                                 <  4:111295/110‑1 (MQ=255)
                                                                                                                               gTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGCGGGTGGGTATCGGGTATGACCAGTTGCAGCTCATCTTCGTACAAcagc                       <  6:246028/149‑1 (MQ=255)
                                                                                                                                         ccAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGCGGGTGGGTATCGGGTATGACCAGTTGCAGCTCATCTTCGTACAACAGCAGGCCATGc              >  6:118667/1‑148 (MQ=255)
                                                                                                                                         ccAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGCGGGTGGGTATCGGGTATGACCAGTTGCAGCTCATCTTCGTACAACAGCAGGCCATGc              <  5:118667/148‑1 (MQ=255)
                                                                                                                                          cAGATCCGTCAGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGCGGGTGGGTATCGGGTATGACCAGTTGCAGCTCATCTTCGTACAACAGCAgg                    <  6:255921/141‑1 (MQ=255)
                                                                                                                                          cAGATCCGTCAGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGCGGGTGGGTATCGGGTATGACCAGTTGCAGCTCATCTTCGTACAACAGCAgg                    >  5:255921/1‑141 (MQ=255)
                                                                                                                                                  aCGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGCGGGTGGGTATCGGGTATGACCAGTTGCAGCTCATCTTCGTACAACAGCAGGCCATGCAGCCctggc     <  8:116941/148‑1 (MQ=255)
                                                                                                                                                    ggATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGCGGGTGGGTATCGGGTATGACCAGTTGCAGCTCATCTTCGTACAACAGCAGGCCATGCAGCCctggctgg  >  4:17375/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCTTCGGCGGCCTCCTTGGCGCGCCCCGGCAGAATGGTCGCCAGCGCTGTGTGCGCGAGGCTATCGAGAATCCCTGCCATATTGTTCATCTCTGCCTGCACCTGCGGCCGTCGGCCCTGACTGGCCAGTTGCGCCTGCCAGATCTGACGGATCTGGAACTCCTCGCCCAGCATCAGCATGGGCAACTCGGCGGCCTGACGGATGGACACCTTCTTGAAGTCCTTCAGCGGGTGGGTATCGGGTATGACCAGTTGCAGCTCATCTTCGTACAACAGCAGGCCATGCAGCCCTGGCTGG  >  NC_002947/451124‑451420

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