Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 655,070 T→C 6.9% T426A (ACC→GCC)  PP_0563 ← two‑component system response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947655,0700TC 6.9% 216.6 / 4.6 72T426A (ACC→GCC) PP_0563two‑component system response regulator
Reads supporting (aligned to +/- strand):  ref base T (34/33);  new base C (2/3);  total (36/36)
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.45e-01

TGAGGTATGTGAAATGCCTGGGGAAGCGGGCGAGCCTGTGAAGAAGCCGGCACCGTGAGTGCGAATTTCAGGGCTCTACACGAACGACGCAATTGCGCCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGCAAATCCCGCGGCTGCGCCGGGTAGAGAATTTCGGCAAAGCGCCGGCGAATTTCGTCCAGCACCTTGTGTGCGGCATCCAGCGCGGTATTGGGCA  >  NC_002947/654923‑655219
                                                                                                                                                   |                                                                                                                                                     
tGAGGTATGTGAAATGCCTGGGGAAGCGGGCGAGCCTGTGAAGAAGCCGGCACCGTGAGTGCGAATTTCAGGGCTCTACACGAACGACGCAATTGCGCCCTTCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCa                                                                                                                                                     <  1:537149/150‑1 (MQ=255)
   ggTATGTGAAATGCCTGGGGAAGCGGGCGAGCCTGTGAAGAAGCCGGCACCGTGAGTGCGAATTTCAGGGCTCTACACGAACGACGCAATTGCGCCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGACCATGGCCAGTg                                                                                                                                                  <  1:1296004/150‑1 (MQ=255)
     tATGTGAAATGCCTGGGGAAGCGGGCGAGCCTGTGAAGAAGCCGGCACCGTGAGTGCGAATTTCAGGGCTCTACACGAACGACGCAATTGCGCCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCa                                                                                                                                                >  1:1126016/1‑150 (MQ=255)
           aaaTGCCTGGGGAAGCGGGCGAGCCTGTGAAGAAGCCGGCACCGTGAGTGCGAATTTCAGGGCTCTACACGAACGACGCAATTGCGCCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATCCAGc                                                                                                                                          <  2:1150812/150‑1 (MQ=255)
             aTGCCTGTGGAAGCGGGCGAGCCTGTGAAGAAGCCGGCACCGTGAGTGCGAATTTCAGGGCTCTACACGAACGACGCAATTGCGCCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATCCAGc                                                                                                                                          <  2:768460/148‑1 (MQ=255)
             aTGCCTGTGGAAGCGGGCGAGCCTGTGAAGAAGCCGGCACCGTGAGTGCGAATTTCAGGGCTCTACACGAACGACGCAATTGCGCCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATCCAGc                                                                                                                                          >  1:768460/1‑148 (MQ=255)
              tGCCTGGGGAAGCGGGCGAGCCTGTGAAGAAGCCGGCACCGTGAGTGCGAATTTCAGGGCTCTACACGAACGACGCAATTGCGCCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATCCAGcccc                                                                                                                                       <  2:419038/150‑1 (MQ=255)
                ccTGGGGAAGCGGGCGAGCCTGTGAAGAAGCCGGCACCGTGAGTGCGAATTTCAGGGCTCTACACGAACGACGCAATTGCGCCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATCCAGCCCCTc                                                                                                                                     >  1:997801/1‑150 (MQ=255)
                        aGCGGGCGAGCCTGTGAAGAAGACGGCACCGTGAGTGCGAATTTCAGGGCTCTACACGAACGACGCAATTGCGCCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATCCAGCCCCTCGTCCAGtt                                                                                                                             <  2:611947/150‑1 (MQ=255)
                               gAGCCTGTGAAGAAGCCGGCACCGTGAGTGCGAATTTCAGGGCTCTACACGAACGACGCAATTGCGCCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGAACTGGCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGcaccac                                                                                                                      <  2:1126016/150‑1 (MQ=255)
                                   cTGTGAAGAAGCCGGCACCGTGAGTGCGAATTTCAGGGCTCTACACGAACGACGCAATTGCGCCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCGTGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCg                                                                                                                  >  2:343092/1‑150 (MQ=255)
                                    tgtgAAGAAGCCGGCACCGTGAGTGCGAATTTCAGGGCTCTACACGAACGACGCAATTGCGCCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCgg                                                                                                                 >  2:1186660/1‑150 (MQ=255)
                                        aagaagCCGGCACCGTGAGTGCGAATTTCAGGGCTCTACACGAACGACGCAATTGCGCCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCagcgga                                                                                                                                           <  2:538960/120‑6 (MQ=255)
                                        aagaagCCGGCACCGTGAGTGCGAATTTCAGGGCTCTACACGAACGACGCAATTGCGCCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCagcgga                                                                                                                                           >  1:538960/1‑115 (MQ=255)
                                        aagaagCCGGCACCGTGAGTGCGAATTTCAGGGCTCTACACGAACGACGCAATTGCGCCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATc                                                                                                                                              >  2:1104154/1‑117 (MQ=255)
                                        aagaagCCGGCACCGTGAGTGCGAATTTCAGGGCTCTACACGAACGACGCAATTGCGCCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATc                                                                                                                                              <  1:1104154/117‑1 (MQ=255)
                                          gaagCCGGCACCGTGAGTGCGAATTTCAGGGCTCTACACGAACGACGCAATTGCGACCCGAATGCTTGGCCTGGAACAGCTCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGa                                                                                                           >  1:531763/1‑150 (MQ=255)
                                             gCCGGCACCGTGAGTGCGAATTTCAGGGCTCTACACGAACGACGCAATTGCGCCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAgg                                                                                                        >  1:539014/1‑150 (MQ=255)
                                              ccGGCACCGTGAGTGCGAATTTCAGGGCTCTACACGAACGACGCAATTGCGCCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCGTGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGt                                                                                                       <  1:343092/150‑1 (MQ=255)
                                                ggCACCGTGAGTGCGAATTTCAGGGCTCTACACGAACGACGCAATTGCGCCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATc                                                                                                                                              <  2:1349227/109‑1 (MQ=255)
                                                ggCACCGTGAGTGCGAATTTCAGGGCTCTACACGAACGACGCAATTGCGCCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATc                                                                                                                                              >  1:1349227/1‑109 (MQ=255)
                                                 gCACCGTGAGTGCGAATTTCAGGGCTCTACACGAACGACGCAATTGCGCCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGACCATGGCCAGTGCagc                                                                                                                                              >  1:524415/1‑106 (MQ=255)
                                                 gCACCGTGAGTGCGAATTTCAGGGCTCTACACGAACGACGCAATTGCGCCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGACCATGGCCAGTGCagc                                                                                                                                              <  2:524415/108‑3 (MQ=255)
                                                             cGAATTTCAGGGCTCTACACGAACGACGCAATTGCGCCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGCAAATCCcg                                                                                        >  1:1324822/1‑150 (MQ=255)
                                                                aTTTCAGGGCTCTACACGAACGACGCAATTGCGCCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCCGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTg                                                                                                      <  2:557832/133‑1 (MQ=255)
                                                                aTTTCAGGGCTCTACACGAACGACGCAATTGCGCCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTg                                                                                                      >  1:557832/1‑133 (MQ=255)
                                                                   tCAGGGCTCTACACGAACGACGCAATTGCGCCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGCAAATCCCGCGGCTg                                                                                  <  2:1384601/150‑1 (MQ=255)
                                                                            tACACGAACGACGCAATTGCGCCCTGCATGCTTGGCCTACCGCAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGCAAATCCCGCGGCTGCGCCGGGTa                                                                         <  1:727833/150‑1 (MQ=255)
                                                                            tACACGAACGACGCAATTGCGCCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGCAAATCCCGCGGCTGCGCCGGGTa                                                                         >  1:1351660/1‑150 (MQ=255)
                                                                            tACACGAACGACGCAATTGCGCCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGCAAATCCCGCGGCTGCGCCGGGTa                                                                         <  2:196115/150‑1 (MQ=255)
                                                                              cacGAACGACGCAATTGCGCCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGCAAATCCCGCGGCTGCGCCGGGTaga                                                                       >  2:510701/1‑150 (MQ=255)
                                                                                  aacgaCGCAATTGCGCCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATcc                                                                                                                                             <  2:1137795/76‑1 (MQ=255)
                                                                                  aacgaCGCAATTGCGCCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATcc                                                                                                                                             >  1:1137795/1‑76 (MQ=255)
                                                                                     gacgCAATTGCGCCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGCAAATCCCGCGGCTGCGCCGGGTAGAGAATTTc                                                                <  1:686525/150‑1 (MQ=255)
                                                                                     gacgCAATTGCGCCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGCAAATCCCGCGGCTGCGCCGGGTAGAGAATTTc                                                                >  1:622372/1‑150 (MQ=255)
                                                                                     gacgCAATTGCGCCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGCAAATCCCGCGGCTGCGCCGGGTAGAGAATTTc                                                                >  2:686525/1‑150 (MQ=255)
                                                                                          aaTTGCGCCCTGCATGCTTGGCCCGGTACAACGCCTCGTCCGCTGCACTGTCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGCAAATCACGCGGCTGCGCCGGGTAGAGAATTTCGGCaa                                                           >  1:629354/1‑150 (MQ=255)
                                                                                           aTTGCGCCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCgg                                                                                                                 >  1:1416453/1‑95 (MQ=255)
                                                                                           aTTGCGCCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCgg                                                                                                                 <  2:1416453/95‑1 (MQ=255)
                                                                                                gcCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGCAAATCCCGCGGCTGCGCCgg                                                                            >  1:1406426/1‑127 (MQ=255)
                                                                                                gcCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGCAAATCCCGCGGCTGCGCCgg                                                                            <  2:1406426/127‑1 (MQ=255)
                                                                                                gcCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGCAAATCCCGCGGCTGCGCCGGGTAGAGAATTTCGGCAAAGCGcc                                                     <  1:1186660/150‑1 (MQ=255)
                                                                                                 cccTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGCAAATCCCGCGGCTGCGCCGGGTAGAGAATTTCGGCAAAGCGCCg                                                    <  1:1255702/150‑1 (MQ=255)
                                                                                                     gCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGCAAATCCCGCGGCTGCGCCGGGTAGAGAATTTCGGCAAAGCGCCGGCGa                                                >  2:412846/1‑150 (MQ=255)
                                                                                                      cATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGCAAATCCCGCGGCTGCGCCGGGTAGAGAATTTCGGCAAAGCGCCGGCGaa                                               >  1:344329/1‑150 (MQ=255)
                                                                                                         gCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGCAAATCCCGCGGCTGCGCCGGGTAGAGAATTTCGGCAAAGCGCCGGCGAAttt                                            >  1:1300524/1‑150 (MQ=255)
                                                                                                           ttGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGt                                                                                                                                                   >  2:1184083/1‑45 (MQ=255)
                                                                                                           ttGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGt                                                                                                                                                   >  2:1184157/1‑45 (MQ=255)
                                                                                                           ttGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGt                                                                                                                                                   <  1:1184083/45‑1 (MQ=255)
                                                                                                           ttGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGt                                                                                                                                                   <  1:1184157/45‑1 (MQ=255)
                                                                                                            tGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGc                                                                                                <  2:631511/95‑1 (MQ=255)
                                                                                                            tGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGc                                                                                                >  1:631511/1‑95 (MQ=255)
                                                                                                            tGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGCAAATCCCGCGGCTGCGCCGGGTAGAGAATTTCGGCAAAGCGCCGGCGAATTTCGt                                         <  2:622372/150‑1 (MQ=255)
                                                                                                             ggCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGCACATCCCGCGGCTGCGCCGGGTAGAGAATTTCGGCAAAGCGCCGGCGAATTTCGTc                                        <  2:1122440/150‑1 (MQ=255)
                                                                                                               cccGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGCAAATCCCGCGGCTGCGCCGGATAGAGAATTTCGGCAAAGCGCCGGCGAATTTCGTCCa                                      >  2:844331/1‑150 (MQ=255)
                                                                                                                ccGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATCCCGAACCTGCTGCAGGTGCACCACCCCGGCACTGAAGGTGCACTGCAAATCCCGCGGCTGCGCCGGGTAGAGAATTTCGGCAAAGCGCCGGCGAATTTCGTCCAg                                     >  2:1058956/1‑150 (MQ=255)
                                                                                                                        gggcCTCGTCCGCTGCACTGGCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGCAAATCCCGCGGCTGCGCCGGGTAGAGAATTTCGGCAAAGCGCCGGCGAATTTCGTCCAGCACCTtgt                             >  2:239250/3‑150 (MQ=255)
                                                                                                                        gcgcCTCGTCCGCTGCACTGGCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGCAAATCCCGCGGCTGCGCCGGGTAGAGAATTTCGGCAAAgcgc                                                      <  1:570265/125‑1 (MQ=255)
                                                                                                                        gcgcCTCGTCCGCTGCACTGGCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGCAAATCCCGCGGCTGCGCCGGGTAGAGAATTTCGGCAAAgcgc                                                      >  2:570265/1‑125 (MQ=255)
                                                                                                                        gcgcCTCGTCCGCTGCACTGGCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGCAAATCCCGCGGCTGCGCCGGGTAGAGAATTTCGGCAAAGCGCCGGCGAATTTCGTCCAGCACCTtgt                             <  2:1300524/150‑1 (MQ=255)
                                                                                                                        gcgcCTCGTCCGCTGCACTGGCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGCAAATCCCGCGGCTGCGCCGGGTAGAGAATTTCGGCAAAGCGCCGGCGAATTTCGTCCAGCACCTtgt                             >  1:701640/1‑150 (MQ=255)
                                                                                                                            cTCGTCCGCTGCACTGGCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGCAAATCCCGCGGCTGCGCCGGGTAGAGAATTTCGGCAAAGCGCCGGCGAATTTCGTCCAGCACCTTGTGTGc                         <  2:1324822/150‑1 (MQ=255)
                                                                                                                             tCGTCCGCTGCACTGGCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGCAAATCCCGCGGCTGCGCCGGGTAGAGAATTTCGGCAAAGCGCCGGCGAATTTCGTCCAGCACCTTGTGTGCg                        <  2:344329/150‑1 (MQ=255)
                                                                                                                               gTCCGCTGCACTGGCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGCAAATCCcgc                                                                                       >  2:1027407/1‑85 (MQ=255)
                                                                                                                               gTCCGCTGCACTGGCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGCAAATCCcgc                                                                                       <  1:1027407/85‑1 (MQ=255)
                                                                                                                                  cGCTGCACTGGCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGCAAATCCCGCGGCTGCGCCGGGTAGAGAATTTCGGCAAAGCGCCGGCGAATTTCGTCCAGCACCTTGTGTGCGGCATc                   <  1:412846/150‑1 (MQ=255)
                                                                                                                                      gCACTGACCATGGCCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGCAAATCCCGCGGCTGCGCCGGGTAGAGAATTTCGGCAAAGCGCCGGCGAATTTCGTc                                        <  2:1025972/125‑1 (MQ=255)
                                                                                                                                      gCACTGACCATGGCCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGCAAATCCCGCGGCTGCGCCGGGTAGAGAATTTCGGCAAAGCGCCGGCGAATTTCGTc                                        >  1:1025972/1‑125 (MQ=255)
                                                                                                                                        aCTGGCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGCAAATCCCGCGGCTGCGCCGGGTAGAGAATTTCGGCAAAGCGCCGGCGAATTTCGTCCAGCACCTTGTGTGCGGCATCCAgcgc             <  2:1351660/150‑1 (MQ=255)
                                                                                                                                              cATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGCAAATCCCGCGGCTGCGCCGGGTAGAGAATTTCGGCAAAGCGCCGGCGAATTTCGTCCAGCACCTTGTGTGCGGCATCCAGCGCGGTAtt       >  1:906660/1‑150 (MQ=255)
                                                                                                                                               aTGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGCAAATCCCGCGGCTGCGCCGGGTAGAGAATTTCGGCAAAGCGCCGGCGAATTTCGTCCAGCACCTTGTGTGc                         >  2:1240031/1‑131 (MQ=255)
                                                                                                                                               aTGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGCAAATCCCGCGGCTGCGCCGGGTAGAGAATTTCGGCAAAGCGCCGGCGAATTTCGTCCAGCACCTTGTGTGc                         <  1:1240031/131‑1 (MQ=255)
                                                                                                                                                   tCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGCAAATCCCGCGGCTGCGCCGGGTAGAGAATTTCGGCAAAGCGCCGGCGAATTTCGTCCAGCACCTTGTGTGc                         >  1:64909/1‑127 (MQ=255)
                                                                                                                                                   tCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGCAAATCCCGCGGCTGCGCCGGGTAGAGAATTTCGGCAAAGCGCCGGCGAATTTCGTCCAGCACCTTGTGTGc                         <  2:64909/127‑1 (MQ=255)
                                                                                                                                                   tCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGCAAATCCCGCGGCTGCGCCGGGTAGAGAATTTCGGCAAAGCGCCGGCGAATTTCGTCCAGCACCTTGTGTGCGGCATCCAGCGCGGTATTGGGCa  >  1:1295191/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TGAGGTATGTGAAATGCCTGGGGAAGCGGGCGAGCCTGTGAAGAAGCCGGCACCGTGAGTGCGAATTTCAGGGCTCTACACGAACGACGCAATTGCGCCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGCAAATCCCGCGGCTGCGCCGGGTAGAGAATTTCGGCAAAGCGCCGGCGAATTTCGTCCAGCACCTTGTGTGCGGCATCCAGCGCGGTATTGGGCA  >  NC_002947/654923‑655219

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