Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,127,163 A→T 7.3% V241E (GTG→GAG) ‡ PP_5375 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,127,1630AT 7.3% 183.4 / 5.3 69V241E (GTG→GAG) ‡PP_5375LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (29/34);  new base T (2/3);  total (32/37)
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.82e-01

CGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCA‑CGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCGGGCCAACCTACCCAGTCGCGCTGGGAGTCTTCAAGCCAGAGCC  >  NC_002947/6127022‑6127310
                                                                                                                                             |                                                                                                                                                    
cGCACTCTGTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCTGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCA‑CGCCCTg                                                                                                                                              <  7:206067/149‑1 (MQ=255)
  cACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCGT‑GGCCCTGGc                                                                                                                                            <  2:263204/149‑1 (MQ=255)
  cACGCTTTCACGAGGACCACGCTCCGTGGACGACTCCAGCAGGCAGAACTTGGAGTCGCTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAAAGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCA‑CGCCCTGGc                                                                                                                                            >  2:105005/1‑149 (MQ=255)
   aCGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCGT‑GGCCCTGGcc                                                                                                                                           <  3:73700/149‑1 (MQ=255)
     gCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCA‑CGCCCTGGCCAc                                                                                                                                         <  6:172251/149‑1 (MQ=255)
        ttCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCA‑CGCCCTGGCCACTg                                                                                                                                       >  2:154923/1‑148 (MQ=255)
          cACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCA‑CGCCCTGGCCACTGAt                                                                                                                                     <  4:280451/148‑1 (MQ=255)
             gAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCA‑CGCCCTGGCCACTGATCGcc                                                                                                                                 <  7:73493/149‑1 (MQ=255)
              aGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCA‑CGCCCTGGCCACTGATCGcc                                                                                                                                 >  3:161465/1‑148 (MQ=255)
               ggACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGGTGCCCCAGGCCAGGGCCC‑CGCCCTGGCCACTGAGCGCCGa                                                                                                                               >  8:103342/1‑149 (MQ=255)
                   cACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCA‑CGCCCTGGCCACTGATCGCCGACTg                                                                                                                            >  2:313611/1‑148 (MQ=255)
                      gCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGACCCAGGCCAGGGCCA‑CGCCCTGGCCACTGATCGCCGACTGGAt                                                                                                                         >  5:148451/1‑148 (MQ=255)
                       gtCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGTGCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCA‑CGCCCTGGCCACTGATCGCCGACTGGATCa                                                                                                                       <  4:161465/148‑1 (MQ=255)
                          ggTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCA‑CGCCCTGGCCACTGATCGCCGACTGGATCAGc                                                                                                                     <  7:103342/148‑1 (MQ=255)
                           gTGGACGACTCCAGCAGGCAGAAGTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGATCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCA‑CGCCCTGGCCACTGATCGCCGACTGGATCAGCAt                                                                                                                   >  7:201667/1‑149 (MQ=255)
                            tGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCGT‑GGCCCTGGcctg                                                                                                                                         >  1:7124/1‑124 (MQ=255)
                            tGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCGT‑GGCCCTGGcctg                                                                                                                                         <  2:7124/126‑3 (MQ=255)
                              gacgacTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCA‑CGCCCTGGCCACTGATCGCCGCCTGGATCAGCATAcc                                                                                                                >  5:317177/1‑149 (MQ=255)
                                  acTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCA‑CGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAg                                                                                                            <  7:208115/149‑1 (MQ=255)
                                  acTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCA‑CGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAg                                                                                                            <  6:247857/149‑1 (MQ=255)
                                   cTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCA‑CGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGc                                                                                                           <  5:122826/149‑1 (MQ=255)
                                      cagcagGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCA‑CGCCCt                                                                                                                                               <  4:176805/110‑1 (MQ=255)
                                      cagcagGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCA‑CGCCCt                                                                                                                                               >  3:176805/1‑110 (MQ=255)
                                      cagcagGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCA‑CGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGt                                                                                                        <  2:260211/149‑1 (MQ=255)
                                      cagcagGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCA‑CGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGt                                                                                                        >  7:176708/1‑149 (MQ=255)
                                            gcagAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCA‑CGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGatat                                                                                                  >  8:207841/1‑149 (MQ=255)
                                            gcagAACTTGGCGTCGGTCAGGACTACAACAACGTTGGGCATGAACAGCGCACAGGTGATAAACGAATCGACCACCAGGTTCCGCCAGGACAGCGCGA‑CGCCCATGCCACGGATCCCCGACTGTATCAGCATACAGTAGCTGTTGatat                                                                                                  >  4:217971/1‑149 (MQ=255)
                                                 aCTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCAAGGCCAGGGCCATCGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGAc                                                                                              <  1:88982/149‑1 (MQ=255)
                                                           ggTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCA‑CGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTgg                                                                                   <  1:154923/149‑1 (MQ=255)
                                                                ggAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCA‑CGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACACGATGCTTGGGGGt                                                                               >  5:70899/1‑148 (MQ=255)
                                                                 gAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCCGGGCCA‑CGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCa                                                                             >  5:204696/1‑149 (MQ=255)
                                                                   gTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGTGCCA‑CGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATg                                                                           <  7:207841/149‑1 (MQ=255)
                                                                        acaTCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCA‑CGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTc                                                                      >  2:247223/1‑149 (MQ=255)
                                                                        acaTCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCA‑CGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACACGATGCTTGGGGGTCATGGGCTc                                                                      <  6:70899/149‑1 (MQ=255)
                                                                                 ggCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCA‑CGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAg                                                             <  1:313611/149‑1 (MQ=255)
                                                                                    ttGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCA‑CGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGgcg                                                          >  2:79072/1‑149 (MQ=255)
                                                                                         cAGCGCACCGGTGCTCAAGGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCA‑CGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATTGGCTCGTATCCCAGGCGCTTGa                                                     <  1:105005/149‑1 (MQ=255)
                                                                                         cAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCA‑CGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCt                                                                       <  8:28983/131‑1 (MQ=255)
                                                                                         cAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCA‑CGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCt                                                                       >  7:28983/1‑131 (MQ=255)
                                                                                          aGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCA‑CGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGaa                                                    <  6:204696/149‑1 (MQ=255)
                                                                                           gcgcACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCA‑CGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAAc                                                   <  8:176708/149‑1 (MQ=255)
                                                                                              cACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCA‑CGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCAt                                                <  1:294180/149‑1 (MQ=255)
                                                                                                 cGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCA‑CGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCg                                             >  3:305836/1‑149 (MQ=255)
                                                                                                  ggTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCA‑CGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCgg                                            <  8:201667/149‑1 (MQ=255)
                                                                                                   gTGCTCAACGCATCGTCCACCAGGTTGCCCGAGGCCAGGGCCA‑CGCCCTGGGCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCGTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCggg                                           <  7:6258/149‑1 (MQ=255)
                                                                                                   gTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCA‑CGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCggg                                           >  1:213882/1‑149 (MQ=255)
                                                                                                       tCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCA‑CGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCg                                             <  4:32686/143‑1 (MQ=255)
                                                                                                       tCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCA‑CGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCg                                             >  3:32686/1‑143 (MQ=255)
                                                                                                         aaCGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCA‑CGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCGGGCCAAc                                      <  6:317177/148‑1 (MQ=255)
                                                                                                               tCGTCCACCAGGTTGCCCCAGGCCAGGGCCA‑CGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCGGGCCAACCTACCCa                               >  2:34852/1‑149 (MQ=255)
                                                                                                                 gTCCACCAGGTTGCCCCAGGCCAGGGCCA‑CGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCGGGCCAACCTACCCAGt                             <  1:247223/149‑1 (MQ=255)
                                                                                                                    caccaGGTTGCCCCAGGCCAGGGCCA‑CGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCGGGCCAACCTACCCAGTcgc                          >  5:282103/1‑149 (MQ=255)
                                                                                                                     accaGGTTGCCCCAGGCCAGGGCCA‑CGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCGGGCCAACCTACCCAGTcgcg                         >  4:285645/1‑149 (MQ=255)
                                                                                                                       caGGTTGCCCCAGGCCAGGGCCA‑CGCCCTGGCAACTGATCGCCGACTGGATCAGCATACCGTAGATGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTGGGGCCAACCTACCCAGTcgcggt                       >  7:202450/1‑147 (MQ=255)
                                                                                                                           ttGCCCCAGGCCAGGGCCA‑CGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCGGGCCAACCTACCCAGTCGCGCTGGGa                   <  3:285645/149‑1 (MQ=255)
                                                                                                                            tGCCCCAGGCCAGGGCCA‑CGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCGGGCCAACCTACCCAGTCGCGCTGGGAg                  >  6:311374/1‑149 (MQ=255)
                                                                                                                             gCCCCAGGCCAGGGCCA‑CGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCGGGCCAACCTa                                   >  7:16008/1‑131 (MQ=255)
                                                                                                                             gCCCCAGGCCAGGGCCA‑CGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCGGGCCAACCTa                                   <  8:16008/131‑1 (MQ=255)
                                                                                                                               cccAGGCCAGGGCCA‑CGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCGGGCCAACCTACCCAGTCGCGCTGGGAGTCt               >  5:155348/1‑149 (MQ=255)
                                                                                                                               cccAGGCCAGGGCCA‑CGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCGGGCCAACCTACCCAGTCGCGCTGGGAGTCt               >  4:78920/1‑149 (MQ=255)
                                                                                                                               cagtggccaggGCGT‑GGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCGGGCCAACCTACCCAGTCGCGCTGGGAGTCt               >  5:280282/5‑149 (MQ=255)
                                                                                                                                 caggccaggGCCA‑CGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCGGGCCAACCTACCCAGTCGCGCTGGGAGTCtt              >  6:266951/1‑148 (MQ=255)
                                                                                                                                  aggccaggGCCA‑CGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACTCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACGATTCGGGCCAACCTACCCAGTCGCGCTGGGAGTCTTCa            <  6:148451/149‑1 (MQ=255)
                                                                                                                                   ggccaggGCCA‑CGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCGGGCCAACCTACCCAGTCGCGCTGGGAGTCTTCaa           <  6:282103/149‑1 (MQ=255)
                                                                                                                                      atttgCCA‑CGCCCTTGCCACTGATCGCAGACTGGATCATGATTCCGTAGCTGTTGATATTTACTCGATTCTTGGGGGTCATGGGCTCGTATCCCAGGCTCTTGAACTATTGGGGGCAACCTACCCATTCGCGCTGGGAGTCTTCaa           <  2:300099/142‑1 (MQ=255)
                                                                                                                                         ggCCA‑CGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCGGGCCAACCTACCCAGTCGCGCTGGGAGTCTTCAagccaga     <  1:79072/149‑1 (MQ=255)
                                                                                                                                           ccA‑CGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCGGGCCAACCTACCAAGTCGCGCTGGGAGTCTTCaa           >  1:300099/1‑141 (MQ=255)
                                                                                                                                            cA‑CGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCggg                                           <  8:78129/108‑1 (MQ=255)
                                                                                                                                            cA‑CGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCGGGCCAACCTACCCAGTCGCGCTGGGAGTCTTCAagccagagcc  <  4:305836/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
CGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCA‑CGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCGGGCCAACCTACCCAGTCGCGCTGGGAGTCTTCAAGCCAGAGCC  >  NC_002947/6127022‑6127310

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