Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,952,931 A→C 8.7% P37P (CCT→CCG PP_3483 ← type II secretion system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,952,9310AC 8.7% 802.8 / 5.7 296P37P (CCT→CCGPP_3483type II secretion system protein
Reads supporting (aligned to +/- strand):  ref base A (141/111);  new base C (11/24);  total (154/142)
Fisher's exact test for biased strand distribution p-value = 6.95e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 315 total aligned reads displayed.

CGGGTGAAGGTGTCGGGGGTGTCGCAGGCCAGGCGTTCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCT  >  NC_002947/3952782‑3953080
                                                                                                                                                     |                                                                                                                                                     
   gTGAAGGTGTCGGGGGTGTCGCAGGCCAGGCGTTCCAGGCAGTCGAGCAGGCGCTCGTCTGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGt                                                                                                                                                    <  1:2355202/150‑1 (MQ=255)
   gTGAAGGTGTCGGGGGTGTCGCAGGCCAGGCGTTCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGg                                                                                                                                                    >  1:5433982‑M1/1‑149 (MQ=255)
         gtgtCGGGGGTGTCGCAGGCAAGGCGTTCCAGTCAGTCGAGCAGGCGCTCGTCAGTCTCTTGCAGACGGGCCTGGGCAAGGAGGTCGCGGGGGTAGGCGGAAATTATTCCTGGGTTCATCGCGGGCATCCCCGCTCCCACACCAACATca                                                                                                                                              >  2:3793770/1‑150 (MQ=255)
           gtCGGGGGTGTCGCAGGCCAGGCGTTCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTGGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATcaca                                                                                                                                            <  2:8048332/150‑1 (MQ=255)
           gtCGGGGGTGTCGCAGGCCAGGCGTTCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATcaca                                                                                                                                            <  1:1391936/150‑1 (MQ=255)
           gtCGGGGGTGTCGCAGGCCAGGCGTTCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATcaca                                                                                                                                            <  1:2038187/150‑1 (MQ=255)
           gtCGGGGGTGTCGCAGGCCAGGCGTTCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATcaca                                                                                                                                            <  2:7042330/150‑1 (MQ=255)
            tCGGGGGTGTCGCAGGCCAGGCGTTCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAAGTGTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGGCCATCACAg                                                                                                                                           >  2:5275538/1‑150 (MQ=255)
             cGGGGGTGTCGCAGGCCAGGCGTTCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGatgtatgcagc                                                                                                                                          <  2:8764377‑M1/150‑12 (MQ=255)
               ggggTGTCGCAGGCCAGGCGTTCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTc                                                                                                                                        <  1:6038479/150‑1 (MQ=255)
                  gtgtCGCAGGCCAGGCGTTCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCtgt                                                                                                                                     <  2:5002466/150‑1 (MQ=255)
                     tCGCAGGCCAGGCGTTCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGatgtatgcagcccacagac                                                                                                                                  >  2:2305239‑M1/1‑131 (MQ=255)
                     tCGCAGGCCAGGCGTTCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTggg                                                                                                                                  <  2:5485204/150‑1 (MQ=255)
                      cGCAGGCCAGGCGTTCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGc                                                                                                                                 <  1:5768946/150‑1 (MQ=255)
                       gCAGGCCAGGCGTTCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACCGGTCTGTGGGCt                                                                                                                                >  2:2642057/1‑150 (MQ=255)
                       gCAGGCCAGGCGTTCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCGCCCACCGGTACATCACAGGTGTGTGGGCt                                                                                                                                >  1:3983898/1‑150 (MQ=255)
                         aggccaggcGTTCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGc                                                                                                                              <  1:1687146/150‑1 (MQ=255)
                         aggccaggcGTTCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCGGTGGGCTGc                                                                                                                              >  2:8227334/1‑150 (MQ=255)
                          ggccaggcGTTCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCa                                                                                                                             <  1:1257509/150‑1 (MQ=255)
                            ccaggcGTTCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATa                                                                                                                           <  2:6032446/150‑1 (MQ=255)
                             caggcGTTCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATAc                                                                                                                          >  2:3527744/1‑150 (MQ=255)
                             caggcGTTCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATAc                                                                                                                          <  2:1342670/150‑1 (MQ=255)
                              aggcGTTCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACa                                                                                                                         <  2:4846231/150‑1 (MQ=255)
                                gcGTTCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCCCCGGGGCCTCACAGGTCGGTGGGCTGCATACagc                                                                                                                       >  1:6854812/1‑148 (MQ=255)
                                  gTTCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCCCCGGGCCATCCCCGGTCGGTGGGCTGCATACCTccg                                                                                                                     <  1:6403757‑M2/150‑4 (MQ=255)
                                    tCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATcctgt                                                                                                                   <  2:1396190‑M2/150‑6 (MQ=255)
                                     ccAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATcctgtg                                                                                                                  >  1:7771973‑M2/1‑144 (MQ=255)
                                     ccAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATcctgtg                                                                                                                  >  2:3067257‑M2/1‑144 (MQ=255)
                                     ccAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATcctgtg                                                                                                                  >  1:2570390‑M2/1‑144 (MQ=255)
                                      cAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGatgtatgcagcccacagacctgtgatgtacctgtgg                                                                                                                 <  1:4835513‑M1/150‑37 (MQ=255)
                                      cAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATcctgtgg                                                                                                                 >  2:3392367‑M2/1‑143 (MQ=255)
                                      cAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATcctgtgg                                                                                                                 >  1:2920989‑M2/1‑143 (MQ=255)
                                      cAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATcctgtgg                                                                                                                 >  2:5046444‑M2/1‑143 (MQ=255)
                                      cAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGACCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCCTGCCCGCGCCCACCGGGCCCTCACCGGTGGGTGGGCGGCATACAGcctgtgg                                                                                                                 <  1:3886973‑M2/150‑8 (MQ=255)
                                       aGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTCCATCACAGGTCTGTGGGCTGCATACATcctgtggg                                                                                                                <  2:6570840‑M2/150‑9 (MQ=255)
                                        ggCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGatgtatgcagcccacagacctgtgatgtacctgtggga                                                                                                               <  2:3932763‑M1/150‑39 (MQ=255)
                                        ggCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGatgtatgcagcccacagacctgtgatgtacctgtggga                                                                                                               <  2:2825477‑M1/150‑39 (MQ=255)
                                        ggCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGatgtatgcagcccacagacctgtgatgtacctgtggga                                                                                                               <  2:1074900‑M1/150‑39 (MQ=255)
                                        ggCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGatgtatgcagcccacagacctgtgatgtacctgtggga                                                                                                               <  2:5901747‑M1/150‑39 (MQ=255)
                                        ggCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGatgtatgcagcccacagacctgtgatgtacctgtggga                                                                                                               <  1:5266611‑M1/150‑39 (MQ=255)
                                        ggCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGatgtatgcagcccacagacctgtgatgtacctgtggga                                                                                                               <  1:6225578‑M1/150‑39 (MQ=255)
                                          cAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGatgtatgcagcccacagacctgtgatgtacctgtgggagc                                                                                                             <  1:3754374‑M1/150‑41 (MQ=255)
                                          cAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATcctgtgggagc                                                                                                             >  1:500945‑M2/1‑139 (MQ=255)
                                           aGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACCGGTACATCACAGGTCTGTGGGCTGCATACATccggtgggggcg                                                                                                            <  1:7750611‑M2/150‑13 (MQ=255)
                                            gTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGatgtatgcagcccacagacctgtgatgtacctgtgggagcgg                                                                                                           <  1:8815559‑M1/150‑43 (MQ=255)
                                            gTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATcctgtgggagcgg                                                                                                           >  1:8121137‑M2/1‑137 (MQ=255)
                                             tCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGatgtatgcagcccacagacctgtgatgtacctgtgggagcggg                                                                                                          <  1:4600486‑M1/150‑44 (MQ=255)
                                                aGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAACTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATcctgtgggagcgggcat                                                                                                       >  2:6812535‑M2/1‑133 (MQ=255)
                                                    ggCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATcctgtgggagcgggcatgccc                                                                                                   >  1:4240349‑M2/1‑129 (MQ=255)
                                                     gcgcTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATcctgtgggagcgggcatgcccg                                                                                                  >  2:567024‑M2/1‑128 (MQ=255)
                                                         tCGTCTGTGGCTTGCATAGGGGCCTGGGCCTGGTGGTCGGGGGGGTTGGCGGAAATTTTTGCTGGGTTCTTCGGGGGCATGCCCGCTCCCACCGGTACATCACCGGTCTGTGGGCTGCATACATcctgtgggagcgggcatgcccgcgaa                                                                                              >  2:852198‑M2/1‑124 (MQ=255)
                                                         tCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATcctgtgggagcgggcatgcccgcgaa                                                                                              <  2:4878750‑M2/150‑27 (MQ=255)
                                                         tCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCGGTGGGCGGCATACATactgtgggagcgggcatggccgcgaa                                                                                              <  1:5069326‑M2/150‑27 (MQ=255)
                                                           gTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGGACATCACAGGTCTGTGGGCTGCATACATactgtgggagcgggcctgcccgcgacgc                                                                                            <  1:6154027‑M2/150‑29 (MQ=255)
                                                                    ttGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGAAGTCGGTGc                                                                                   >  2:3384979/1‑150 (MQ=255)
                                                                     ggctgcggggCCGGGGCCCGGCGGTCGGGGGGGTTGGCGTAATTTTTTGCGGGGTCCTTCGGGGGCAGGCCCGCCCCCACAGGTACCTCACCGGTCGGTGGGCTGCATACAGCCTGTGGGAGGGGGAATGCCCGCGAAGACGTCGGTGCa                                                                                  <  2:8373452/143‑1 (MQ=255)
                                                                       cAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACGTCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGAAGGCGGTGCAgg                                                                                >  1:2793259/1‑150 (MQ=255)
                                                                                      cAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCATGTGGGAGCGGGCATGCCCGCGACGACGTCGGTGCAGGCTTTGCAGATTCCgg                                                                 >  1:8049400/1‑150 (MQ=255)
                                                                                      cAGGAGGTCGCGGGGGGAGGCGGAAATTTTTGCGGGGTTCTTCGCGGGCATGCCCGCTCCCACGGGTACATCACAGGTCTGTGGGCAGCATACATCCTGTGGGAGCGGGCCTGCCCGCGACCACGTCGGTGCAGGCGTTTCAGATTCAgg                                                                 >  2:4765834/1‑150 (MQ=38)
                                                                                      aaggaggTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTAAAAAAAAGGTAAGTGGGCTGCATAAATCCTTTTGTATCGGGCCAGCCCGCGAAAACGTCGTTGAAGGCTATGCACATTCCgg                                                                 >  2:1557530/2‑150 (MQ=255)
                                                                                         gaggTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGCAGACGTCGGTGCAGGCGTTGCGGATTCCGGGGt                                                              >  2:6382090/1‑150 (MQ=255)
                                                                                               gggggggTCGGGGGAAATTTGTGCGGGGGGCTCCGCGGGCCGGCCCGCGCCCCCCGGTCCCGCCCCGGTCTGGGGGCTGCATCCAGCCGGGGGGAGGGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATc                                                        <  2:2374309/148‑1 (MQ=255)
                                                                                               gcgGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCGCCGGTCGGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGACGACGTCGGTGCAGGCTTTGCGGATGCCGGGGTTGCATc                                                        >  2:7230381/1‑150 (MQ=39)
                                                                                                          cgaatttttttGCTGGGTTGTTCGGGGGCATGCCCGCGCCCACCGGTACAGCACCGGTCGGTGGGCTGCATACAGCCGGTGGGCGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTg                                             <  2:3085213/144‑1 (MQ=255)
                                                                                                           ggAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCGGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTgg                                            >  2:6449322/1‑150 (MQ=255)
                                                                                                                     cccgtgttcgtCGGGGGCTTGCCCGCTCCCACCGGatgtatgcagcccacagacctgtgatgtacctgtgggagcgggcatgcccgcgaagaacccagcaaaaatttccgcctacccccgcgacctcctggcccaggcccgtctgcaagc                                  <  1:2863680‑M1/145‑116 (MQ=255)
                                                                                                                              ttcGCGGGCATGCCCGCTCCCACCGGatgtatgcagcccacagaccggtgatgtacctgtgggagcgggcctgcccgcgaagaacccagcaaaaatttccgcctacccccgcgacctcctggcccaggcccgtctgccagccacagacgg                         >  2:6921931‑M1/1‑26 (MQ=255)
                                                                                                                              gccgggggcagccCCGCGCCCACCGGatgtatgcagcccaccgacctgtgatggacctgtgggcgcgggcatgcccgcgaagaacccagcaaaaatttccgcctacccccgcgacctcctggcccaggcccgtctgcaagccactgacga                         <  1:3198826‑M1/145‑125 (MQ=255)
                                                                                                                                  ggggCTTGCCCGCTCCCACCGTTACATCACTGGTCTGTGGGCTGCATACATCCTGTGGGCGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAgg                     <  2:2143765/149‑1 (MQ=255)
                                                                                                                                   gggCATGCCCGCTCCCACAGGatgtatgcagcccacagacctgtgatgtacctgtgggagcgggcatgcccgcgaagaacccagccaaaatttccgcctacccccgcgccctcctggccccggcccgtctgcacgccactgacgagcgcc                    >  2:6670674‑M1/1‑21 (MQ=255)
                                                                                                                                     gcctGCCCGCTCCCCCCGGTCCCTCACCGGTCTGTGGGCTGCATACCTCCGGGGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGgc                  <  2:8626171/147‑1 (MQ=255)
                                                                                                                                     gCATGCCCGCTCCCACTGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGgc                  <  2:7864014/150‑1 (MQ=255)
                                                                                                                                     gCATGCCCGCTCCCACAGGatgtatgcagcccacagacctgtgatgtacctgtgggagcgggcatgcccgcgaagaacccagcaaaaatttccgcctacccccgcgacctcctggcccaggcccgtctgcaagccactgacgagcgcctg                  >  2:7418312‑M1/1‑19 (MQ=255)
                                                                                                                                       aTGCCCGCTCCCACCGGatgtatgcagcccaccgacctgtgatgtacctgtgggagcgggcatgcccgcgacgaacccagcaaaaatttccgcctacccccgcgacctcctggcccaggcccgtctgcaagccactgacgagcgcctgct                >  2:6217179‑M1/1‑17 (MQ=255)
                                                                                                                                       aTGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGgcgc                >  2:7506962/1‑150 (MQ=255)
                                                                                                                                        tGCCCGCTCCCACAGGatgtatgcagcccacagacctgtgatgtacctgtgggagcgggcatgcccgcgaagaacccagcaaaaatttccgcctacccccgcgacctcctggcccaggcccgtctgcaagccactgacgagcgcctgctc               >  1:8638776‑M1/1‑16 (MQ=255)
                                                                                                                                         gCCCGCTCCCACAGGatgtatgcagcccacagacctgtgatgtacctgtgggagcgggcatgcccgcgaagaacccagcaaaaatttccgcctacccccgcgacctcctggcccaggcccgtctgcaagccactgacgagcgcctgctcg              >  1:8687450‑M1/1‑15 (MQ=255)
                                                                                                                                         gCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATcctgtgggagcgggcatgcccgcgaagaacccagcaaaaatttccgcctacccccgcgacctcctggcccaggcccgtctgcaagccactgacgagcgcctgctcg              <  1:2305046‑M2/150‑107 (MQ=255)
                                                                                                                                          cccGCTCCCACAGGatgtatgcagcccacagacctgtgatgtacctgtgggagcgggcatgcccgcgaagaacccagcaaaaatttccgcctacccccgcgacctcctggcccaggcccgtctgcaagccactgacgagcgcctgctcga             >  1:8530443‑M1/1‑14 (MQ=255)
                                                                                                                                          cccGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATcctgtgggagcgggcatgcccgcgaagaacccagcaaaaatttccgcctacccccgcgacctcctggcccaggcccgtctgcaagccactgacgagcgcctgctcga             <  2:25860‑M2/150‑108 (MQ=255)
                                                                                                                                           tcggccccACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCg            <  2:2537668/145‑1 (MQ=255)
                                                                                                                                           ccGCTCCCACAGGatgtatgcagcccacagacctgtgatgtacctgtgggagcgggcatgcccgcgaagaacccagcaaaaatttccgcctacccccgcgacctcctggcccaggcccgtctgcaagccactgacgagcgcctgctcgac            >  2:5209254‑M1/1‑13 (MQ=255)
                                                                                                                                           ccGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATcctgtgggagcgggcatgcccgcgaagaacccagcaaaaatttccgcctacccccgcgacctcctggcccaggcccgtctgcaagccactgacgagcgcctgctcgac            <  2:894102‑M2/150‑109 (MQ=255)
                                                                                                                                           ccGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCg            >  1:329990/1‑150 (MQ=255)
                                                                                                                                            cGCTCCCACAGGatgtatgcagcccacagacctgtgatgtacctgtgggagcgggcatgcccgcgaagaacccagcaaaaatttccgcctacccccgcgacctcctggcccaggcccgtctgcaagccactgacgagcgcctgctcgact           >  2:2977551‑M1/1‑12 (MQ=255)
                                                                                                                                            cGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATcctgtgggagcgggcatgcccgcgaagaacccagcaaaaatttccgcctacccccgcgacctcctggcccaggcccgtctgcaagccactgacgagcgcctgctcgact           <  2:4336786‑M2/150‑110 (MQ=255)
                                                                                                                                             gCTCCAACTGTTACATCAGTGGTGTGTGGGCTGAATACTGCGTGTGGGCGGGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCa          <  1:1141145/150‑1 (MQ=255)
                                                                                                                                              cccccacAGGatgtacgcagcccacacacctgtgatgtacctgtgggagcgggcatgcccgcgaagaacccagcaaaaatttccgcctacccccgcgaccccctggcccaggcccgtctgcaagccactgacgagcgcctgctcgactgc         <  1:6112467‑M1/148‑141 (MQ=255)
                                                                                                                                              cTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAg         >  1:7233279/1‑150 (MQ=255)
                                                                                                                                               tCCCACCGGatgtatgcagcccacagacctgtgatgtacctgtgggagcgggcatgcccgcgaagaacccagcaaaaatttccgcctacccccgcgacctcctggcccaggcccgtctgcaagccactgacgagcgcctgctcgactgcc        <  1:6950293‑M1/150‑142 (MQ=255)
                                                                                                                                               tCCCACAGGatgtatgcagcccacagacgtgtgatgtacctgtgggagcgggcatgcccgcgaagaacccagcaaaaatttccgcctacccccgcgacctcctggcccaggcccgtctgcaagccactgacgagcgcctgctcgactgcc        >  2:4088668‑M1/1‑9 (MQ=255)
                                                                                                                                               tCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATcctgtgggagcgggcatgcccgcgaagaacccagcaaaaatttccgcctacccccgcgacctcctggcccaggcccgtctgcaagccactgacgagcgcctgctcgactgcc        <  1:2086573‑M2/150‑113 (MQ=255)
                                                                                                                                                 ccACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAggct      >  1:7511103/1‑150 (MQ=255)
                                                                                                                                                  cacaGGatgtatgcagcccacagacctgtgatgtacctgtgggagcgggcatgcccgcgaagaacccagcaaaaatttccgcctacccccgcgacctcctggcccaggcccgtctgcaagccactgacgagcgcctgctcgactgcctgg     >  1:1930839‑M1/1‑6 (MQ=255)
                                                                                                                                                  cacaGGatgtatgcagcccacagacctgtgatgtacctgtgggagcgggcatgcccgcgaagaacccagcaaaaatttccgcctacccccgcgacctcctggcccaggcccgtctgcaagccactgacgagcgcctgctcgactgcctgg     >  1:8806330‑M1/1‑6 (MQ=255)
                                                                                                                                                  cacaGGatgtatgcagcccacagacctgtgatgtacctgtgggagcgggcatgcccgcgaagaacccagcaaaaatttccgcctacccccgcgacctcctggcccaggcccgtctgcaagccactgacgagcgcctgctcgactgcctgg     >  1:4375364‑M1/1‑6 (MQ=255)
                                                                                                                                                  cacaGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAggctg     >  2:7735096/1‑150 (MQ=255)
                                                                                                                                                   acaGGatgtatgcagcccacagacctgtgatgtacctgtgggagcgggcatgcccgcgaagaacccagcaaaaatttccgcctacccccgcgacctcctggcccaggcccgtctgcaagccactgacgagcgcctgctcgactgcctgga    >  2:5989437‑M1/1‑5 (MQ=255)
                                                                                                                                                   acaGGTACATCACCGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGTAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAggctgg    <  2:3886973/150‑1 (MQ=255)
                                                                                                                                                     aGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAggctggct  >  1:7171368/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGGGTGAAGGTGTCGGGGGTGTCGCAGGCCAGGCGTTCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCT  >  NC_002947/3952782‑3953080

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: