Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,728,844 A→T 5.8% intergenic (‑207/+117) opgG ← / ← dtd glucans biosynthesis protein G/D‑tyrosyl‑tRNA(Tyr) deacylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,728,8440AT 5.8% 208.4 / 2.8 69intergenic (‑207/+117)opgG/dtdglucans biosynthesis protein G/D‑tyrosyl‑tRNA(Tyr) deacylase
Reads supporting (aligned to +/- strand):  ref base A (34/31);  new base T (2/2);  total (36/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 = 8.65e-01

GAGCAACCCTCGCCGAAAACAGATCATGAAACCAACGGCCAGCGACAAATGCCGGGTTGGCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCTGTGAACAACGGGTTTTTTGACCTCATATTTGTAACATAAATGTTACCGGGC  >  NC_002947/5728698‑5728989
                                                                                                                                                  |                                                                                                                                                 
gAGCAACCCTCGCCGAAAACAGATCATGAAACCAACGGCCAGCGACAAATGCCGGGTTGGCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAAct                                                                                                                                                 >  6:282081/1‑148 (MQ=255)
  gCAACCCTCGCCGAAAACAGATCATGAAACCAACGGCCAGCGACAAATGCCGGGTTGGCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACg                                                                                                                                                 >  6:262268/1‑147 (MQ=255)
         tCGCCGAAAACAGATCATGAAACCAACGGCCAGCGACAAATGCCGGGTTGGCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGc                                                                                                                                                >  4:143253/1‑141 (MQ=255)
         tCGCCGAAAACAGATCATGAAACCAACGGCCAGCGACAAATGCCGGGTTGGCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGc                                                                                                                                                <  3:143253/141‑1 (MQ=255)
         tCGCCGAAAACAGATCATGAAACCAACGGCCAGCGACAAATGCCGGGTTGGCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCAt                                                                                                                                        >  3:312236/1‑149 (MQ=255)
           gCCGAAAACAGATCATGAAACCAACGGCCAGCGACAAATGCCGGGTTGGCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTGCGTTCGTAGGTccgttt                                                                                                                                      <  8:285900/149‑5 (MQ=255)
                    aGATCATGAAACCAACGGCCAGCGACAAATGCCGGGTTGGCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGa                                                                                                                             <  2:107647/149‑1 (MQ=255)
                       tCATGAAACCAACGGCCAGCGACAAATGCCGGGTTGGCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAgg                                                                                                                                             <  5:200632/130‑1 (MQ=255)
                       tCATGAAACCAACGGCCAGCGACAAATGCCGGGTTGGCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAgg                                                                                                                                             >  6:200632/1‑130 (MQ=255)
                         aTGAAACCAACGGCCAGCGACAAATGCCGGGTTGGCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATc                                                                                                                        <  5:320038/149‑1 (MQ=255)
                                  aCGGCCAGCGACAAATGCCGGGTTGGCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGtt                                                                                                               <  3:66215/149‑1 (MQ=255)
                                              aaaTGCCGGGTTGGCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACACTGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACg                                                                                                    >  6:192758/1‑148 (MQ=255)
                                              aaaTGCCGGGTTGGCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTGCGTTCGTAGGt                                                                                                                                            <  6:162775/108‑1 (MQ=255)
                                              aaaTGCCGGGTTGGCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTGCGTTCGTAGGt                                                                                                                                            >  5:162775/1‑108 (MQ=255)
                                               aaTGCCGGGTTGGCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCa                                                                                                  <  8:144066/149‑1 (MQ=255)
                                               aaTGCCGGGTTGGCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCa                                                                                                  <  6:145709/149‑1 (MQ=255)
                                                aTGCCGGGTTGGCGCGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAATCGGATCCAGACAATCCTCGAGTCAGAGTGAAACGGACCTACGAACGCCGGTCCATGCAGCGCGCGATTATACAGAAGGTTGCGTGACAACCCa                                                                                                  >  5:325598/1‑148 (MQ=255)
                                                aTGCCGGGTTGGCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGAtt                                                                                                                           >  6:73822/1‑123 (MQ=255)
                                                aTGCCGGGTTGGCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGAtt                                                                                                                           <  5:73822/123‑1 (MQ=255)
                                                   ccGGGTTGGCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCa                                                                                               <  5:282081/148‑1 (MQ=255)
                                                   ccGGGTTGGCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCa                                                                                               <  6:98236/148‑1 (MQ=255)
                                                   ccGGGTTGGCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAg                                                                                              >  1:220770/1‑149 (MQ=255)
                                                         tGGCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCTTGCAGCGCGCGATTATCCAGCAGGTTGCg                                                                                                            >  1:11187/1‑129 (MQ=255)
                                                         tGGCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCTTGCAGCGCGCGATTATCCAGCAGGTTGCg                                                                                                            <  2:11187/129‑1 (MQ=255)
                                                          ggCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGc                                                                                                   <  8:289135/137‑1 (MQ=255)
                                                          ggCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGc                                                                                                   >  7:289135/1‑137 (MQ=255)
                                                            cGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGc                                                                                      >  5:90059/1‑148 (MQ=255)
                                                             gAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCGAGACAATGCTCGAGTCAGAGTGAAACGGAGCTACGCACGCAGGTGGATGGAGCGCGGGGTTTTCGAGCAGGTTGCGTGCGAGCGGAGCAGGCTgggggag                                                                                    >  8:209706/1‑144 (MQ=255)
                                                                    gTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTa                                                                             >  6:218879/1‑149 (MQ=255)
                                                                     ttCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTa                                                                             <  5:101830/148‑1 (MQ=255)
                                                                     ttCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTa                                                                             >  6:101830/1‑148 (MQ=255)
                                                                     ttCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAAGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATc                                                                           >  6:133643/1‑149 (MQ=255)
                                                                        ggACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGGAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGc                                                                                                             >  1:101932/1‑113 (MQ=255)
                                                                        ggACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGc                                                                                                             <  2:101932/113‑1 (MQ=255)
                                                                        ggACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCg                                                                          <  5:262268/148‑1 (MQ=255)
                                                                          aCCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTa                                                                        >  4:76204/1‑148 (MQ=255)
                                                                                    ggCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCtt                                                             >  7:199400/1‑149 (MQ=255)
                                                                                       gTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAgc                                                                                                                                   <  8:276614/76‑1 (MQ=255)
                                                                                       gTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAgc                                                                                                                                   >  7:276614/1‑76 (MQ=255)
                                                                                          aTGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGtttt                                                        <  3:76204/148‑1 (MQ=255)
                                                                                          aTGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTc                                                       <  6:90059/149‑1 (MQ=255)
                                                                                             aTTCGCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCt                                                     <  5:7648/148‑1 (MQ=255)
                                                                                              ttCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCtgt                                                   >  8:186669/1‑149 (MQ=255)
                                                                                              ttCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCtgt                                                   <  2:178237/149‑1 (MQ=255)
                                                                                                 cccAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAAcc                                                                 >  3:73818/1‑132 (MQ=255)
                                                                                                 cccAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAAcc                                                                 <  4:73818/132‑1 (MQ=255)
                                                                                                 cccAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCTGTGa                                                 >  5:286938/1‑148 (MQ=255)
                                                                                                  ccAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCTGTGaa                                                >  7:5885/1‑148 (MQ=255)
                                                                                                  ccAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCTGTGaa                                                >  1:81953/1‑148 (MQ=255)
                                                                                                  ccAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCTGTGaa                                                >  6:2446/1‑148 (MQ=255)
                                                                                                  ccAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCTGTGaa                                                >  3:102276/1‑148 (MQ=255)
                                                                                                    aaaCGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCTGTGaacaa                                             >  3:149809/1‑149 (MQ=255)
                                                                                                       cGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCTGTGAACAACgg                                          <  2:145776/149‑1 (MQ=255)
                                                                                                          aTCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCg                                                                          <  6:153774/114‑1 (MQ=255)
                                                                                                          aTCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCg                                                                          >  5:153774/1‑114 (MQ=255)
                                                                                                              aGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCTGTGAACAACGGGtttttt                                   <  6:286938/149‑1 (MQ=255)
                                                                                                                                agTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCtgt                                                   <  7:218087/115‑1 (MQ=255)
                                                                                                                                agTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCtgt                                                   >  8:218087/1‑115 (MQ=255)
                                                                                                                                    aaaCGGACCTACTAACGAAGGTCCATGCAGCGCGCGATTATCCAGAAGGTTGAGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACACCTTGTTTTCCTGTGAAAAACCGTTTTTTTGACCTCATATTTGTAACATaaa             >  2:263372/1‑149 (MQ=255)
                                                                                                                                    aaaCGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCTGTGAACAACGGGTTTTTTGACCTCATATTTGTAACATaaa             >  4:66212/1‑149 (MQ=255)
                                                                                                                                     aaCGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAAccc                                                                <  6:174063/97‑1 (MQ=255)
                                                                                                                                     aaCGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAAccc                                                                >  5:174063/1‑97 (MQ=255)
                                                                                                                                     aaCGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCTGTGAACAACGGGTTTTTTGACCTCATATTTGTAACATAAAt            <  8:5885/149‑1 (MQ=255)
                                                                                                                                      atggACCTGCGTTCGTAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCTGTGAACAACGGGTTTTTTGACCTCATATTTGTAACATAAATg           >  6:287353/3‑149 (MQ=255)
                                                                                                                                          aCCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCTGTGAACAACGGGTTTTTTGACCTCATATTTGTAACATAAATGTTAc       <  5:218879/149‑1 (MQ=255)
                                                                                                                                              aCGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCTGTGAACAACGGGtttttt                                   <  7:135855/117‑1 (MQ=255)
                                                                                                                                              aCGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCTGTGAACAACGGGtttttt                                   >  8:135855/1‑117 (MQ=255)
                                                                                                                                               cGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCTGTGAACAACGGGTTTTTTGACCTCATATTTGTAACATAAATGTTACCggg   <  5:2446/148‑1 (MQ=255)
                                                                                                                                               cGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCTGTGAACAACGGGTTTTTTGACCTCATATTTGTAACATAAATGTTACCGGGc  <  4:149809/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GAGCAACCCTCGCCGAAAACAGATCATGAAACCAACGGCCAGCGACAAATGCCGGGTTGGCGAGCACTGTTCGGACCCCGAGAGGGCGTAATGATTCCCCAAACGGATCCAGACAATGCTCGAGTCAGAGTGAAACGGACCTACGAACGCAGGTCCATGCAGCGCGCGATTATCCAGCAGGTTGCGTGACAACGCATCAGGCTGGAGCAACTATTTATCGTACAAAACCCCTTGTTTTCCTGTGAACAACGGGTTTTTTGACCTCATATTTGTAACATAAATGTTACCGGGC  >  NC_002947/5728698‑5728989

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: