Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,199,349 C→G 5.6% T66T (ACC→ACG gspF → protein secretion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,199,3490CG 5.6% 334.6 / 10.1 108T66T (ACC→ACGgspFprotein secretion protein
Reads supporting (aligned to +/- strand):  ref base C (48/54);  new base G (3/3);  total (51/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.75e-01

Only 100 of 115 total aligned reads displayed.

AGGCCAGCCTGCAGGCTGACAGTGAACGCCACGCCCGCCAGCTGTTGCGCGAGCAGGGCCTGTTCGCCCGCCAGTTGCAGCGTCATGACGCCGGCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCACCCTTGAACGGCAACTGCGTCAGCCTGCCCTGCACAGTGTGCTGGTGGCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCT  >  NC_002947/1199201‑1199498
                                                                                                                                                    |                                                                                                                                                      
aGGCCAGCCTGCAGGCTGACAGTGAACGCCACGCCCGCCAGCTGTTGCGCGAGCAGGGCCTGTTCGCCCGCCAGTTGCAGCGTCATGACGCCGGCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGAccc                                                                                                                                                       <  2:42827/150‑1 (MQ=255)
          gCAGGGTGACAGTGAACGCCACGCCCGCCAGCTGTTGCGCGAGCAGGGCCTGTTCGCCCGCCAGTTGCAGCGTCATGACGCCGGCAGCCGACAACCCCGACACCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGc                                                                                                                                             >  2:1561180/1‑150 (MQ=255)
             ggCTGACAGTGAACGCCACGCCCGCCAGCTGTTGTGCGAGCAGGGCCTGTTCGCCCGCCAGTTGCAGCGTCATGACGCCGGCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCAc                                                                                                                                          <  1:1631187/150‑1 (MQ=255)
                tGACAGTGAACGCCACGCCCGCCAGCTGTTGCGCGAGCAGGGCCTGTTCGCCCGCCAGTTGCAGCGTCATGACGCCGGCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCtt                                                                                                                                       >  2:744522/1‑150 (MQ=255)
                  aCAGTGAACGCCACGCCCGCAAGCTGTTGCGCGAGCAGGGCCTGTTCGCCCGCCAGTTGCAGCGTCATGACGCCGGCAGCCGACAACCACGACGCCAGCGCGTTAGCCGCGCCCAGCTGTGTGAGCTGACCAGTCAGCTGGCCACCTTGa                                                                                                                                     >  2:1019697/1‑150 (MQ=255)
                    aGTGAACGCCACGCCCGCCAGCTGTTGCGCGAGCAGGGCCTGTTCGCCCGCCAGTTGCAGCGTCATGACGCCGGCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGAcc                                                                                                                                   <  1:569702/150‑1 (MQ=255)
                     gTGAACGCCACGCCCGCCAGCTGTTGCGCGAGCAGGGCCTGTTCGCCCGCCAGTTGCAGCGTCATGACGCCGGCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCg                                                                                                                                  >  2:368456/1‑150 (MQ=255)
                        aaCGCCACGCCCGCCAGCTGTTGCGCGAGCAGGGCCTGTTCGCCCGCCAGTTGCAGCGTCATGACGCCGGCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCgggg                                                                                                                              >  1:2162484/1‑150 (MQ=255)
                        aaCGCCACGCCCGCCAGCTGTTGCGCGAGCAGGGCCTGTTCGCCCGCCAGTTGCAGCGTCATGACGCCGGCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGAGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCgggg                                                                                                                              <  1:409238/150‑1 (MQ=255)
                          cgccacgccCGCCAGCTGTTGCGCGAGCAGGGCCTGTTCGCCCGCCAGTTGCAGCGTCATGACGCCGGCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGcc                                                                                                                            <  2:378929/150‑1 (MQ=255)
                           gccacgccCGCCAGCTGTTGCGCGAGCAGGGCCTGTTCGCCCGCCAGTTTCAGCGTCATGACGCCGGCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCg                                                                                                                           >  2:1497000/1‑150 (MQ=255)
                           gccacgccCGCCAGCTGTTGCGCGAGCAGGGCCTGTTCGCCCGCCAGTTGCAGCGTCATGACGCCGGCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCg                                                                                                                           >  2:290534/1‑150 (MQ=255)
                           gccacgccCGCCAGCTGTTGCGCGAGCAGGGCCTGTTCGCCCGCCAGTTGCAGCGTCATGACGCCGGCAACCGACAACCCCGACGCCAGCGCCTGAGCCGGTCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCg                                                                                                                           >  2:848946/1‑150 (MQ=255)
                             cacgccCGCCAGCTGTTGCGCGAGCAGGGCCTGTTCGCCCGCCAGTTGCAGCGTCATGACGCCGGCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGc                                                                                                                         <  2:180534/150‑1 (MQ=255)
                             cacgccCGCCAGCTGTTGCGCGAGCAGGGCCTGTTCGCCCGCCAGTTGCAGCGTCATGACGCCGGCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGc                                                                                                                         >  1:1187186/1‑150 (MQ=255)
                                  ccgccAGCTGTTGCGCGAGCAGGGCCTGTTCGCCCGCCAGTTGCAGCGTCATGACGCCGGCAGCCGCCCACCGCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGCACTGTCAGCTGGCCACCTTGACCGGGGCCGGCATTcc                                                                                                                    <  1:1354889/150‑1 (MQ=255)
                                  ccgccAGCTGTTGCGCGAGCAGGGCCTGTTCGCCCGCCAGTTGCAGCGTCATGACGCCGGCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGc                                                                                                                             <  1:1435335/141‑1 (MQ=255)
                                  ccgccAGCTGTTGCGCGAGCAGGGCCTGTTCGCCCGCCAGTTGCAGCGTCATGACGCCGGCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGc                                                                                                                             >  2:1435335/1‑141 (MQ=255)
                                  ccgccAGCTGTTGCGCGAGCAGGGCCTGTTCGCCCGCCAGTTGCAGCGTCATGACGCCGGCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTcc                                                                                                                    >  1:923717/1‑150 (MQ=255)
                                     ccAGCTGTTGCGCGAGCAGGGCCTGTTCGCCCGCCAGTTGCAGCGTCATGACGCCGGCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGt                                                                                                                                                     <  2:694851/115‑1 (MQ=255)
                                     ccAGCTGTTGCGCGAGCAGGGCCTGTTCGCCCGCCAGTTGCAGCGTCATGACGCCGGCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGt                                                                                                                                                     >  1:694851/1‑115 (MQ=255)
                                      cAGCTGTTGCGCGAGCAGGGCCTGTTCGCCCGCCAGTTGCAGCGTCATGACGCCGGCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTg                                                                                                                <  1:972993/150‑1 (MQ=255)
                                      cAGCTGTTGCGCGAGCAGGGCCTGTTCGCCCGCCAGTTGCAGCGTCATGACGCCGGCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTg                                                                                                                <  1:244918/150‑1 (MQ=255)
                                         cTGTTGCGCGAGCAGGGCCTGTTCGCCCGCCAGTTGCAGCGTCATGACGCCGGCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTc                                                                                                             <  1:820328/150‑1 (MQ=255)
                                            ttGCGCGAGCAGGGCCTGTTCGCCCGCCAGTTGCAGCGTCATGACGCCGGCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGAt                                                                                                          <  1:1077643/150‑1 (MQ=255)
                                            ttGCGCGAGCAGGGCCTGTTCGCCAGCAAGTTGCAGCGTCAGGACGCCGGCAGCCGAATACCCCGCCAACATCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGAt                                                                                                          <  1:1797151/150‑1 (MQ=255)
                                               cgcgAGCAGGGCCTGTTCGCCCGCCAGTTGCAGCGTCCTGACGCCGGCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCGTGACCGGGGCCGGCATTCCCTTGGTCGATGcc                                                                                                       >  1:11230/1‑150 (MQ=255)
                                                gcgAGCAGGGCCTGTTCGCCCGCCAGTTGCAGCGTCATGACGCCGGCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGccc                                                                                                      >  2:1707529/1‑150 (MQ=255)
                                                gcgAGCAGGGCCTGTTCGCCCGCCAGTTGCAGCGTCATGACGCCGGCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGccc                                                                                                      >  2:905231/1‑150 (MQ=255)
                                                  gAGCAGGCCCTCTTCGCCCGCAAGTCGCCCCGTCATTACGCCGCCTGCCGACAACACCGACGCGAGCGCCTGACCCGCGCCCAGCTGTGTGAGCTCACCCGTCACCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCtt                                                                                                    <  2:923623/150‑1 (MQ=255)
                                                    ccAGGGCCTGTTCGCCCGCCAGTTGCAGCGTCATGACGCCGGCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGc                                                                                                  <  1:838095/149‑1 (MQ=255)
                                                        ggCCTGTTCGCCCGCCAGTTGCAGCGTCATGACGCCGGCAGCCGACAACCCCGGCGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCAcc                                                                                              >  2:2259142/1‑150 (MQ=255)
                                                        ggCCTGTTCGCCCGCCAGTTGCAGCGTCATGACGCCGGCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAg                                                                                                                                                  <  1:1837773/99‑1 (MQ=255)
                                                        ggCCTGTTCGCCCGCCAGTTGCAGCGTCATGACGCCGGCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCAcc                                                                                              <  1:1434566/150‑1 (MQ=255)
                                                            tGTTCGCCCGCCAGTTGCAGCGTCATGACGCCGGCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCACCCTTg                                                                                          >  2:859857/1‑150 (MQ=255)
                                                                 gcccgccAGTTGCAGCGTCATGACGCCGGCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCACCCTTGAACgg                                                                                     <  1:1882774/150‑1 (MQ=255)
                                                                 gcccgccAGTTGCAGCGTCATGACGCCGGCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCACCCTTGAACgg                                                                                     <  1:880117/150‑1 (MQ=255)
                                                                  cccgccAGTTGCAGCGTCATGACGCCGGCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCACCCTTGAACGGc                                                                                    <  1:740019/150‑1 (MQ=255)
                                                                   ccgccAGTTGCAGCGTCATGACGCCGGCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCACCCTTGAACGGCa                                                                                   >  1:628369/1‑150 (MQ=255)
                                                                    cgccAGTTGCAGCGTCATGACGCCGGCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACGGGTCAGCtcaca                                                                                                                                            >  1:499513/1‑89 (MQ=255)
                                                                    cgccAGTTGCAGCGTCATGACGCCGGCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACGGGTCAGCtcaca                                                                                                                                            <  2:499513/93‑5 (MQ=255)
                                                                    cgccAGTTGCAGCGTCATGACGCCGGCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACGGGTCAGCtcaca                                                                                                                                            >  1:1705486/1‑89 (MQ=255)
                                                                    cgccAGTTGCAGCGTCATGACGCCGGCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACGGGTCAGCtcaca                                                                                                                                            <  2:1705486/93‑5 (MQ=255)
                                                                     gccAGTTGCAGCGTCATGACGCCGGCAGCCGACAACCCCGGCGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGc                                                                                                                                             >  1:2231593/1‑91 (MQ=255)
                                                                     gccAGTTGCAGCGTCATGACGCCGGCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGc                                                                                                                                             <  2:2231593/91‑1 (MQ=255)
                                                                       cAGTTGCAGCGTCATGACGCCGGCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCACCCTTGAACGGCAACTg                                                                               >  2:1743125/1‑150 (MQ=255)
                                                                        aGTTGCAGCGTCATGACGCCGGCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCACCCTTGAACGGCAACTGc                                                                              <  2:1773264/150‑1 (MQ=255)
                                                                        aGTTGCAGCGTCATGACGCCGGCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCACCCTTGAACGGCAACTGc                                                                              <  2:2162484/150‑1 (MQ=255)
                                                                             cAGCGTCATGACGCCGGCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGAccc                                                                                                                                                       >  2:1043443/1‑73 (MQ=255)
                                                                             cAGCGTCATGACGCCGGCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGAccc                                                                                                                                                       <  1:1043443/73‑1 (MQ=255)
                                                                             cAGCGTCATGACGCCGGCAGCAGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCACCCTTGAACGGCAACTGCGTCAg                                                                         >  1:643227/1‑150 (MQ=255)
                                                                              aGCTTCATGACGCCGGCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCACCCTTGAACGGCAACTGCGTCAGc                                                                        >  1:1930356/1‑150 (MQ=255)
                                                                              aGCGTCATGACGCCGGCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCACCCTTGAACGGCAACTGCGTCAGc                                                                        <  1:368456/150‑1 (MQ=255)
                                                                               gCGTCATGACGCCGGCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCACCCTTGAACGGCAACTGCGTCAGcc                                                                       >  2:568027/1‑150 (MQ=255)
                                                                                cGTCATGACGCCGGCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCACCCTTGaa                                                                                        <  1:57772/132‑1 (MQ=255)
                                                                                   cATGACGCCGGCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCACCCTTGAACGGCAACTGCGTCAGcctgcc                                                                   >  2:1100373/1‑150 (MQ=255)
                                                                                        cGCCGGCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCACCCTTGAACGGCAACTGCGTCAGCCTGCCCTGca                                                              <  2:1930356/150‑1 (MQ=255)
                                                                                          ccGGCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCACCCTTGAACGGCAACTGCGTCAGCCTGCCCTGcaca                                                            <  2:633717/150‑1 (MQ=255)
                                                                                            tccAGCTGACAACCCCGACGTCAGCGCCTGAGCCGCGCCTAGCTTTGTGAGCTTATCCGTCAGCTGGCCCGCCTACCGGTCGCCACGATTCTCTTCGTCGATGTCCTTCCCACCCTTGCCCTGCTCCTGCGTCAGCCTGCCCTGCACAgt                                                          <  2:37903/148‑1 (MQ=255)
                                                                                            ggCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCtt                                                                                                                 <  2:1059815/95‑1 (MQ=255)
                                                                                            ggCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCtt                                                                                                                 >  1:1059815/1‑95 (MQ=255)
                                                                                            ggCAGCAGACAACCGCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCACCCTTGAACGGCAACTGCGTCAGCCTGCCCTGCACAgt                                                          >  2:144801/1‑150 (MQ=255)
                                                                                              cAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCACCCTTGAACGGCAACTGCGTCAGCCTGCCCTGCACAgtgt                                                        <  2:628369/150‑1 (MQ=255)
                                                                                               aGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGAt                                                                                                          >  1:338995/1‑99 (MQ=255)
                                                                                               aGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGAt                                                                                                          <  2:338995/99‑1 (MQ=255)
                                                                                                        agccGACGCCAGCCCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCACCCTTGAACGGCAACTGCGTCAGCCTGCCCTGCACAGTGTGCTGGTGGcc                                              >  2:1744010/3‑150 (MQ=255)
                                                                                                         cccGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCACCCTTGAACGGCAACTGCGTCAGCCTGCCCTGCACAGTGTGCTGGTGGCCt                                             >  1:262593/1‑150 (MQ=255)
                                                                                                          ccGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCACCCTTGAACGGCAACTGCGTCAGCCTGCCCTGCACAgtg                                                         >  1:669790/1‑137 (MQ=255)
                                                                                                          ccGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCACCCTTGAACGGCAACTGCGTCAGCCTGCCCTGCACAgtg                                                         <  2:669790/137‑1 (MQ=255)
                                                                                                                 cAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCATTGGTCAATGCCCTTGCAACACTTGAACGGCAACTGCGTCAGCCTGCCCTGCACAGTGTGCTGGTGGCCTTGCGCGGc                                     >  2:2302754/1‑150 (MQ=255)
                                                                                                                 cAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTACCGTGGGCGATGCCCTTGCCACCCTTGAACGGCAACTGCGTCAGCCTGCCCTGCACAGTGTGCTGGTGGCCTTGCGCGGc                                     >  2:449337/1‑150 (MQ=255)
                                                                                                                  aGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCACCCTTGAACGGCAACTGCGTCAGCCTGCCCTGCACAGTGTGCTGGTGGCCTTGCGCGGCt                                    <  1:2304410/150‑1 (MQ=255)
                                                                                                                       cTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCACCCTTGAACGGCAACTGCGTCAGCCTGCCCTGCACAGTGTGCTGGTGGCCTTGCGCGGCTCACTc                               >  1:832320/1‑150 (MQ=255)
                                                                                                                          aGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCACCCTTGAACGGCAACTGCGTCAGCCTGCCCTGCACAGTGTGCTGGTGGCCTTGCGCGGCTCACTCGcc                            <  1:1497000/150‑1 (MQ=255)
                                                                                                                             cgcgCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCACCCTTGAACGGCAACTGCGTCAGACTGCGCTGCACAGTGTGCTGGTGGCCTTGCGCGGCTCACTCGcaagg                         <  1:1707529/150‑5 (MQ=255)
                                                                                                                              gcgcCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCACCCTTGAACGGCAACTGCGTCAGCCTGCCCTGCACAGTGTGCTGGTGGCCTTGCGCGGCTCACTCGCCGAAg                        <  1:1150768/150‑1 (MQ=255)
                                                                                                                                gcCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCACCCTTGAACGGCAACTGCGTCAGCCTGCCCTGCACAg                                                           >  1:721386/1‑113 (MQ=255)
                                                                                                                                gcCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCACCCTTGAACGGCAACTGCGTCAGCCTGCCCTGCACAg                                                           <  2:721386/113‑1 (MQ=255)
                                                                                                                                gcCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCACCCTTGAACGGCAACTGCGTCAGCCTGCCCTGCACAGTGTGCTGGTGGCCTTGCGCGGCTCACTCGCCGAAggc                      <  1:2259142/150‑1 (MQ=255)
                                                                                                                                 cccAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCACCCTTGAACGGCAACTGCGTCAGCCTGCCCTGCACAGTGTGCTGGTGGCCTTGCGCGGCTCACTCGCCGCAggcc                     >  2:429913/1‑150 (MQ=255)
                                                                                                                                 cccAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCACCCTTGAACGGCAACTGCGTCAGCCTGCCCTGCACAGTGTGCTGGTGGCCTTGCGCGGCTCACTCGCCGAAggcc                     <  1:744522/150‑1 (MQ=255)
                                                                                                                                  ccAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCACCCTTGAACGGCAACTGCGTCAGCCTGCCCTGCACAGTGTGCTGGTGGCCt                                             >  2:271914/1‑125 (MQ=255)
                                                                                                                                   cAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCACCCTTGAACGGCAACTGCGTCAGCCTGCCCTGCACAGTGTGCTGGTGGCCTTgc                                          <  2:758238/127‑1 (MQ=255)
                                                                                                                                   cAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCACCCTTGAACGGCAACTGCGTCAGCCTGCCCTGCACAGTGTGCTGGTGGCCTTGCGCGGCTCACTCGCCGAAggcctg                   <  1:855882/150‑1 (MQ=255)
                                                                                                                                    aGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCACCCTTGAACGGCAACTGCGTCAGCCTGCCCTGCACAGTGTGCTGGTGGCCTTGCGCGGCTCACCCGCCGAAggcctgg                  <  2:1750047/150‑1 (MQ=255)
                                                                                                                                      cTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCACCCTTGAACGGCAACTGCGTCAGCCTGCCCTGCACAGTGTGCTGGTGGCCTTGCGCGGCTCACTCGCCGAAggcctgggc                <  1:1100373/150‑1 (MQ=255)
                                                                                                                                       tgtgtgAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCACCCTTGAACGGCAACTGCGTCAGCCTGCCCTGCACAGTGTGCTGGTGGCCTTGCGCGGCTCACTCGCCGAAggcctgggcc               <  2:262593/150‑1 (MQ=255)
                                                                                                                                          gtgAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGcc                                                                                                       <  1:2020553/59‑1 (MQ=255)
                                                                                                                                          gtgAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGcc                                                                                                       >  2:2020553/1‑59 (MQ=255)
                                                                                                                                          gtgAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCACCCTTGAACGGCAACTGCGTCAGCCTGCCCTGCACAGTGTGCTGGTGGCCTTGCGCGGCTCACTCGCCGAAGgcctgggcctgg            <  2:340290/150‑1 (MQ=255)
                                                                                                                                          gtgAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCACCCTTGAACGGCAACTGCGTCAGCCTGCCCTGCACAGTGTGCTGGTGGCCTTGCGCGGCTCACTCGCCGAAGgcctgggcctgg            <  1:848946/150‑1 (MQ=255)
                                                                                                                                             aGCTGACGGGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCACCCTTGAACGGCAACTGCGTCAGCCTGCCCTGCACAGTGTGCTGGTGGCCTTg                                           >  1:913845/1‑116 (MQ=255)
                                                                                                                                             aGCTGACGGGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCACCCTTGAACGGCAACTGCGTCAGCCTGCCCTGCACAGTGTGCTGGTGGCCTTg                                           <  2:913845/116‑1 (MQ=255)
                                                                                                                                                 gACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCACCCTTGAACGGCAACTGCGTCAGCCTGCCCTGCACAGTGTGCtggt                                                  <  1:901781/105‑1 (MQ=255)
                                                                                                                                                 gACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCACCCTTGAACGGCAACTGCGTCAGCCTGCCCTGCACAGTGTGCtggt                                                  >  2:901781/1‑105 (MQ=255)
                                                                                                                                                 gACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCACCCTTGAACGGCAACTGCGTCAGCCTGCCCTGCACAGTGTGCTGGTGGCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGccc         <  2:1993148/146‑1 (MQ=255)
                                                                                                                                                 gACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCACCCTTGAACGGCAACTGCGTCAGCCTGCCCTGCACAGTGTGCTGGTGGCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAg     >  2:1363851/1‑150 (MQ=255)
                                                                                                                                                    ccGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCACCCTTGAACGGCAACTGCGTCAGCCTGCCCTGCACAGTGTGCTGGTGGCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGcc          >  1:1314709/1‑142 (MQ=255)
                                                                                                                                                    ccGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCACCCTTGAACGGCAACTGCGTCAGCCTGCCCTGCACAGTGTGCTGGTGGCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGcc          <  2:1314709/142‑1 (MQ=255)
                                                                                                                                                    ccGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCACCCTTGAACGGCAACTGCGTCAGCCTGCCCTGCACAGTGTGCTGGTGGCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCt  >  1:1339981/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
AGGCCAGCCTGCAGGCTGACAGTGAACGCCACGCCCGCCAGCTGTTGCGCGAGCAGGGCCTGTTCGCCCGCCAGTTGCAGCGTCATGACGCCGGCAGCCGACAACCCCGACGCCAGCGCCTGAGCCGCGCCCAGCTGTGTGAGCTGACCCGTCAGCTGGCCACCTTGACCGGGGCCGGCATTCCCTTGGTCGATGCCCTTGCCACCCTTGAACGGCAACTGCGTCAGCCTGCCCTGCACAGTGTGCTGGTGGCCTTGCGCGGCTCACTCGCCGAAGGCCTGGGCCTGGCCCGCAGCCT  >  NC_002947/1199201‑1199498

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