Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 655,063 G→A 8.3% A428V (GCC→GTC)  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,0630GA 8.3% 200.6 / 7.5 73A428V (GCC→GTC) PP_0563two‑component system response regulator
Reads supporting (aligned to +/- strand):  ref base G (33/33);  new base A (3/3);  total (37/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.15e-01

TCCAGACTTGAGGTATGTGAAATGCCTGGGGAAGCGGGCGAGCCTGTGAAGAAGCCGGCACCGTGAGTGCGAATTTCAGGGCTCTACACGAACGACGCAATTGCGCCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGCAAATCCCGCGGCTGCGCCGGGTAGAGAATTTCGGCAAAGCGCCGGCGAATTTCGTCCAGCACCTTGTGTGCGGCATCCAGCGC  >  NC_002947/654915‑655208
                                                                                                                                                    |                                                                                                                                                 
tCCAGACTTGAGGTATGTGAAATGCCTGGGGAAGCGGGCGAGCCTGTGAAGAAGCCGGCACCGTGAGTGCGAATTTCAGGGCTCTACACGAACGACGCAATTGCGCCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGc                                                                                                                                                  >  2:1255702/1‑150 (MQ=255)
tCCAGACTTGAGGTATGTGAAATGCCTGGGGAAGCGGGCGAGCCTGTGAAGAAGCCGGCACCGTGAGTGCGAATTTCAGGGCTCTACACGAACGACGCAATTGCGCCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGc                                                                                                                                                  >  2:698488/1‑150 (MQ=255)
tCCAGACTTGAGGTATGTGAAATGCCTGGGGAAGCGGGCGAGCCTGTGAAGAAGCCGGCACCGTGAGTGCGAATTTCAGGGCTCTACACGAACGACGCAATTGCGCCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGc                                                                                                                                                  <  1:189643/150‑1 (MQ=255)
       ttGAGGTATGTGAAATGCCTGGGGAAGCGGGCGAGCCTGTGAAGAAGCCGGCACCGTGAGTGCGAATTTCAGGGCTCTACACGAACGACGCAATTGCGCCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGACCATggcc                                                                                                                                           >  2:1296004/1‑148 (MQ=255)
        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                                                                                                                               >  1:768460/1‑148 (MQ=255)
                     aTGCCTGTGGAAGCGGGCGAGCCTGTGAAGAAGCCGGCACCGTGAGTGCGAATTTCAGGGCTCTACACGAACGACGCAATTGCGCCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATCCAGc                                                                                                                               <  2:768460/148‑1 (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                                                                                                                                   <  2:524415/108‑3 (MQ=255)
                                                         gCACCGTGAGTGCGAATTTCAGGGCTCTACACGAACGACGCAATTGCGCCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGACCATGGCCAGTGCagc                                                                                                                                   >  1:524415/1‑106 (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                                                                                                                                  >  1:1137795/1‑76 (MQ=255)
                                                                                          aacgaCGCAATTGCGCCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATcc                                                                                                                                  <  2:1137795/76‑1 (MQ=255)
                                                                                             gacgCAATTGCGCCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGCAAATCCCGCGGCTGCGCCGGGTAGAGAATTTc                                                     <  1:686525/150‑1 (MQ=255)
                                                                                             gacgCAATTGCGCCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGCAAATCCCGCGGCTGCGCCGGGTAGAGAATTTc                                                     >  2:686525/1‑150 (MQ=255)
                                                                                             gacgCAATTGCGCCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGCAAATCCCGCGGCTGCGCCGGGTAGAGAATTTc                                                     >  1:622372/1‑150 (MQ=255)
                                                                                                  aaTTGCGCCCTGCATGCTTGGCCCGGTACAACGCCTCGTCCGCTGCACTGTCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGCAAATCACGCGGCTGCGCCGGGTAGAGAATTTCGGCaa                                                >  1:629354/1‑150 (MQ=255)
                                                                                                   aTTGCGCCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCgg                                                                                                      <  2:1416453/95‑1 (MQ=255)
                                                                                                   aTTGCGCCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCgg                                                                                                      >  1:1416453/1‑95 (MQ=255)
                                                                                                        gcCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGCAAATCCCGCGGCTGCGCCgg                                                                 <  2:1406426/127‑1 (MQ=255)
                                                                                                        gcCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGCAAATCCCGCGGCTGCGCCgg                                                                 >  1:1406426/1‑127 (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:1184157/1‑45 (MQ=255)
                                                                                                                   ttGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGt                                                                                                                                        <  1:1184083/45‑1 (MQ=255)
                                                                                                                   ttGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGt                                                                                                                                        >  2:1184083/1‑45 (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                             >  1:1025972/1‑125 (MQ=255)
                                                                                                                                              gCACTGACCATGGCCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGCAAATCCCGCGGCTGCGCCGGGTAGAGAATTTCGGCAAAGCGCCGGCGAATTTCGTc                             <  2:1025972/125‑1 (MQ=255)
                                                                                                                                                aCTGGCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGCAAATCCCGCGGCTGCGCCGGGTAGAGAATTTCGGCAAAGCGCCGGCGAATTTCGTCCAGCACCTTGTGTGCGGCATCCAgcgc  <  2:1351660/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TCCAGACTTGAGGTATGTGAAATGCCTGGGGAAGCGGGCGAGCCTGTGAAGAAGCCGGCACCGTGAGTGCGAATTTCAGGGCTCTACACGAACGACGCAATTGCGCCCTGCATGCTTGGCCCGGTACAGCGCCTCGTCCGCTGCACTGGCCATGGTCAGTGCATCCAGCCCCTCGTCCAGTTGCACCACCCCGGCACTGAAGGTGCACTGCAAATCCCGCGGCTGCGCCGGGTAGAGAATTTCGGCAAAGCGCCGGCGAATTTCGTCCAGCACCTTGTGTGCGGCATCCAGCGC  >  NC_002947/654915‑655208

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: