Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 939,935 C→G 7.0% G1897A (GGC→GCC)  PP_0806 ← surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947939,9350CG 7.0% 683.0 / 28.0 247G1897A (GGC→GCC) PP_0806surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (96/133);  new base G (13/5);  total (109/138)
Fisher's exact test for biased strand distribution p-value = 1.46e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 261 total aligned reads displayed.

GTTGCCTGCCGCATCGGTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAATGGCCACCACGTCCAGCAGTTCG  >  NC_002947/939787‑940084
                                                                                                                                                     |                                                                                                                                                     
     ctgccGCATCGGTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGtt                                                                                                                                                 <  1:1394479/150‑1 (MQ=255)
      tgccGCATCGGTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGTTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGGCAGCACCGTACGCGTGCTGGCCgccga                                                                                                                                                <  1:6013561/150‑5 (MQ=14)
        ccGCATCGGTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCCGCCCCCGGACCGTGCTGCCCCGTTCgg                                                                                                                                              >  1:4221850/1‑150 (MQ=14)
            aTCGGTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGcgcg                                                                                                                                          <  2:4922660/150‑1 (MQ=255)
                 tGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCgg                                                                                                                                     <  2:2491642/150‑1 (MQ=255)
                 tGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCgg                                                                                                                                     <  1:5915552/150‑1 (MQ=255)
                 tGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCgg                                                                                                                                     <  2:4740504/150‑1 (MQ=255)
                 tGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCACCACCCCGAACGTGCTGCCCCGTTCGGCGCGGcccg                                                                                                                                     >  1:5428299/1‑148 (MQ=14)
                    ccGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCACCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGGGa                                                                                                                                  >  1:1044409/1‑150 (MQ=255)
                     cGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCCGCACGCCTACGGTGCTGCCCCGTTCGGCGCGGCCGGTgtc                                                                                                                                 >  1:5509237/1‑148 (MQ=255)
                       tGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCCGCACCCGTACGGTGCGGCCCCGGGCGGCGCGGCCCGGGAccc                                                                                                                               >  1:4842789/1‑149 (MQ=14)
                        gACTTGGAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCAc                                                                                                                              <  2:2077551/150‑1 (MQ=255)
                           ttGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCt                                                                                                                           <  2:1958602/150‑1 (MQ=255)
                           ttGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGGACGGTGCTGCCCGGTTCGGCGCGGCCCGTGACCAcccg                                                                                                                           >  2:2577825/1‑147 (MQ=255)
                              aaGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGcc                                                                                                                        <  1:2587535/150‑1 (MQ=40)
                              aaGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCGACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGcc                                                                                                                        <  2:5735648/150‑1 (MQ=255)
                                  aCCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCg                                                                                                                    <  1:4413539/150‑1 (MQ=39)
                                     tCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTTCGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGGCGGCa                                                                                                                 >  2:1005426/1‑150 (MQ=255)
                                     tCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACCCGTACGGTGCGGCCCGGTTCGGCGCGGCCGGTGACCACCCGGCCGGCGGCa                                                                                                                 >  1:2152895/1‑150 (MQ=11)
                                      cGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCAt                                                                                                                <  1:3153515/150‑1 (MQ=35)
                                      cGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCAt                                                                                                                <  2:113627/150‑1 (MQ=35)
                                           tCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATc                                                                                                           <  1:1743733/150‑1 (MQ=34)
                                             tATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCAc                                                                                                         <  1:4917538/150‑1 (MQ=34)
                                             tATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCAc                                                                                                         <  1:1859654/150‑1 (MQ=34)
                                              aTCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGGGGGCAGCACGCCGACGGGGCTGCCCGGTGGGGCCCGGCCCGTGGCCACCCGGCCGTCGGCAGTGGTAAcc                                                                                                        >  1:1497989/1‑150 (MQ=11)
                                               tCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCa                                                                                                       <  2:3473782/150‑1 (MQ=34)
                                                     gCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCCGCACCCCTACCGTGCTGCCCCGGTCGGCGCCGCCCGTGACCACCCTGCCCGCGGCAGTGATCCCCACTtcgg                                                                                                 >  2:2407304/1‑150 (MQ=11)
                                                      cGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTtcggt                                                                                                <  1:276796/150‑1 (MQ=34)
                                                      cGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTtcggt                                                                                                <  2:5044707/150‑1 (MQ=34)
                                                        ggTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTtcggtcg                                                                                              <  1:5950013/150‑1 (MQ=34)
                                                        ggTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTtcggtcg                                                                                              <  1:5859239/150‑1 (MQ=34)
                                                          tCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTcggtcggt                                                                                            <  2:4378257/150‑1 (MQ=34)
                                                          tCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTcggtcggt                                                                                            <  2:3441615/150‑1 (MQ=34)
                                                            gAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGc                                                                                          <  2:1813034/150‑1 (MQ=34)
                                                            gAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGc                                                                                          <  2:4516836/150‑1 (MQ=34)
                                                               ggTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCCGCACGCGTACGGTGCTGCCCCGTTTGGCGCGGCCGGGGACCACCCAGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGc                                                                                       >  2:717378/1‑150 (MQ=11)
                                                               ggTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGGACGGTGCTGCCCGGGTCGGCGCGGCCGGGGGCCACCCAGCCGGCGGCATTTATGACCACTTCGGTCGGTGGGGc                                                                                       >  1:1401929/1‑150 (MQ=11)
                                                                    tGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCCCTTCGGTCGGTGCGGCAgggg                                                                                  >  1:6819064/1‑150 (MQ=21)
                                                                      cTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCCGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCg                                                                                >  2:1609610/1‑150 (MQ=21)
                                                                       tGAAGCTGCCGGTAGCGCCATCTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCgg                                                                               <  1:3391142/150‑1 (MQ=21)
                                                                       tGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGGGGCCAGCCCGCGGACGGTGCTGCCCCGGTGGGCGCGGCCCGGGACCACCCGGCCCTCCGCCGTGATAAACCCGTAGGTCGGGGGGGCAGGGGCgg                                                                               >  1:3592308/1‑150 (MQ=2)
                                                                       tGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCgg                                                                               <  1:5686603/150‑1 (MQ=34)
                                                                       tGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCgg                                                                               <  2:6022956/150‑1 (MQ=34)
                                                                       tGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGGTCGGCGCGGGCGGTGACCCCGCGGCCGGCGGCAGTGTTCACCAGTTCGGTCGGGGCGGCGGGGGCgg                                                                               >  1:1024111/1‑150 (MQ=11)
                                                                        gAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCCGCCCCCGGACGGGGCGGCCCGGGTCGGCGCGGCCCGGGCCCACCCCGCCGGCGGCCGTGATCCCCCGGTCGGGCGGGGCGGCCGGGGCggg                                                                              >  1:460233/1‑149 (MQ=2)
                                                                        gAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGt                                                                              <  1:2547845/150‑1 (MQ=34)
                                                                           gCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGat                                                                           <  1:2936864/150‑1 (MQ=34)
                                                                              gCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATc                                                                        >  2:1544647/1‑150 (MQ=34)
                                                                              gCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGGGCCCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATc                                                                        >  2:2198013/1‑150 (MQ=21)
                                                                                  gTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGGACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGGCGGCATTGATCACCAGGTCGGTCGGGGCGGCAGGGGCGGTGATATCcggg                                                                    >  2:3928991/1‑149 (MQ=255)
                                                                                    aGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCcggtgc                                                                  <  2:6843099/150‑1 (MQ=34)
                                                                                     gcgcCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCcggtgcg                                                                 <  1:4138838/150‑1 (MQ=34)
                                                                                     gcgcCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCCTACGGTGCTGCCCGGGTCGGCGCGGCCGGTGACCACGCAGCCGGCGGCCTTGATCACCAGTGCGGTCGGGGCGGCAGGGGCGGTGATATCcggtgcg                                                                 >  1:2165874/1‑150 (MQ=11)
                                                                                      cgcCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCcggtgcgg                                                                <  1:5154359/150‑1 (MQ=34)
                                                                                       gcCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGcc                                                                                                                                                      >  2:882478/1‑63 (MQ=35)
                                                                                       gcCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGcc                                                                                                                                                      <  1:882478/63‑1 (MQ=35)
                                                                                        ccAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACCCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACCCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATAcggtgcggtg                                                              >  1:5642368/1‑150 (MQ=255)
                                                                                         cAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGa                                                             <  2:6456386/150‑1 (MQ=34)
                                                                                         cAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGa                                                             <  2:5326082/150‑1 (MQ=34)
                                                                                          aaCTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGGCAGCGCGCGTACGGTGGTGGCCGGTTGGGGGCGGCCGGTGACCACGCTGCCGTGGGCAATGATAACCAGTTAGGTGGGTGCGGCAGGGGCGGTGATATACGGTTCGGTGAt                                                            >  2:1377774/1‑150 (MQ=11)
                                                                                            cTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGAttt                                                          <  2:5499260/150‑1 (MQ=34)
                                                                                             tATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCCGCCCGCGGGCGGTGCTGGCCGGTTCGGCGCGGCCGGTGGCCACGCTGCCGGCGGCCTTGATCACCAGTTCGGTCGGGGCGGCCGGGGCGGGGGGATCcggggcggggggttg                                                         >  2:2285876/1‑143 (MQ=11)
                                                                                             tATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCg                                                                                                                                      <  1:3826373/73‑1 (MQ=39)
                                                                                             tATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGGGATATCcggtgcggtgggttg                                                         >  2:3812768/1‑145 (MQ=19)
                                                                                              aTCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTtgg                                                        <  2:6494367/150‑1 (MQ=33)
                                                                                              aTCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTtgg                                                        <  2:3614733/150‑1 (MQ=33)
                                                                                              aTCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTtgg                                                        <  2:6594127/150‑1 (MQ=33)
                                                                                                 aCGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCAGTGATCACCAGTTCGGTCGGGGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTtgggtg                                                     >  2:6208184/1‑150 (MQ=17)
                                                                                                 aCGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGGTCGGCGCGGCCGGGGACCACGCGGCCGTCGGCATTGATCACCAGTTCGGGCGGGGCGGCAGGGGCGGTGGTATCCGGTGCGGTGATTtgggtg                                                     >  1:3920488/1‑150 (MQ=11)
                                                                                                  cGCTACCCAGCACAGTGGTGCCGTCGGCGGCCCGCCCGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCACTTCGGTCGGTGCGGCAGGGGCGGGGATATACGGTGCGGTGATTtgggtgg                                                    >  1:2817910/1‑150 (MQ=11)
                                                                                                  cGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTtgggtgg                                                    <  2:5105364/150‑1 (MQ=32)
                                                                                                  cGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTtgggtgg                                                    <  2:3308710/150‑1 (MQ=32)
                                                                                                  cGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCCGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTAGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTtgggtgg                                                    >  2:5967648/1‑150 (MQ=17)
                                                                                                         cAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGa                                             <  1:996787/150‑1 (MQ=32)
                                                                                                          aGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGTTATCCGGTGCGGTGATTTGGGTGGGCAGCGAg                                            >  1:6872559/1‑150 (MQ=17)
                                                                                                          aGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAg                                            <  2:5177542/150‑1 (MQ=32)
                                                                                                            cacaGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCGCGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGa                                          <  1:2903898/150‑1 (MQ=17)
                                                                                                            cacaGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGa                                          <  2:4812079/150‑1 (MQ=32)
                                                                                                                gtggtgCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAaca                                      <  1:5140548/150‑1 (MQ=32)
                                                                                                                    tgCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACAccgc                                  <  2:3866424/150‑1 (MQ=32)
                                                                                                                     gCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACAccgcc                                 <  2:1156357/150‑1 (MQ=32)
                                                                                                                      ccGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACAccgcca                                <  2:4343475/150‑2 (MQ=17)
                                                                                                                        gTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGtc                              <  2:5635793/150‑1 (MQ=32)
                                                                                                                           ggcggcCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGtcatc                           >  1:6282556/1‑150 (MQ=32)
                                                                                                                             cgggCAGCACCGTACGCGTGCTGGCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCaa                         >  1:667581/1‑150 (MQ=11)
                                                                                                                              gggcAGCACCGTACGCGTGCTGGCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAAt                        >  1:2485543/4‑150 (MQ=11)
                                                                                                                               gcCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAATg                       <  1:4324517/150‑1 (MQ=32)
                                                                                                                                gcAGCACCGTACGCGTGCTGGCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAgg                                           >  2:3230190/2‑129 (MQ=11)
                                                                                                                                 cAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAATGGc                     >  2:3879611/1‑150 (MQ=32)
                                                                                                                                 cAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAATGGc                     >  1:4033645/1‑150 (MQ=32)
                                                                                                                                 cAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAATGGc                     >  2:2181625/1‑150 (MQ=32)
                                                                                                                                    cACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAATGGccac                  <  2:2855161/150‑1 (MQ=32)
                                                                                                                                       gcgTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAATGGCcaccac               >  1:2349115/1‑150 (MQ=32)
                                                                                                                                            cGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAATGGCCACCACGTCca          <  1:5566733/150‑1 (MQ=32)
                                                                                                                                             ggTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAATGGCCACCACGTCcag         <  1:476662/150‑1 (MQ=32)
                                                                                                                                             ggTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAATGGCCACCACGTCcag         >  1:6455311/1‑150 (MQ=32)
                                                                                                                                                tgctgcCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAATGGCCACCACGTCcagca       >  2:294168/1‑150 (MQ=32)
                                                                                                                                                   tgcCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAATGGCCACCACGTCCAGCAGtt    <  1:3154047/150‑1 (MQ=32)
                                                                                                                                                     cccGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAATGGCCACCACGTCCAGCAGTTCg  >  1:5054723/1‑150 (MQ=32)
                                                                                                                                                     cccGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAATGGCCACCACGTCCAGCAGTTCg  >  1:779777/1‑150 (MQ=32)
                                                                                                                                                     |                                                                                                                                                     
GTTGCCTGCCGCATCGGTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAATGGCCACCACGTCCAGCAGTTCG  >  NC_002947/939787‑940084

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