Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,952,942 A→C 7.8% C34G (TGT→GGT)  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,9420AC 7.8% 722.9 / 4.5 249C34G (TGT→GGT) PP_3483type II secretion system protein
Reads supporting (aligned to +/- strand):  ref base A (110/107);  new base C (13/12);  total (125/124)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 271 total aligned reads displayed.

GTCGGGGGTGTCGCAGGCCAGGCGTTCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAGGGGT  >  NC_002947/3952793‑3953091
                                                                                                                                                     |                                                                                                                                                      
gtCGGGGGTGTCGCAGGCCAGGCGTTCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATcaca                                                                                                                                                        <  2:1704578/150‑1 (MQ=255)
gtCGGGGGTGTCGCAGGCCAGGCGTTCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATcaca                                                                                                                                                        <  2:5475681/150‑1 (MQ=255)
   gggggTGTCGCAGGCCAGGCGTTCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGt                                                                                                                                                     <  2:3103988/150‑1 (MQ=255)
   ggcgcccgcacccctggccgagccagcctgcgacgcgcccctgagcggagtcgcaggccatgccagaagcattcgatgcaaccccggaatctgcaaagcctgcaccgacgTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGt                                                                                                                                                     >  2:7440832‑M1/111‑150 (MQ=255)
       gtgtCGCAGGCCAGGCGTTCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCtgt                                                                                                                                                 <  1:1187974/150‑1 (MQ=255)
       gtgtCGCAGGCCAGGCGTTCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCtgt                                                                                                                                                 <  1:6942678/150‑1 (MQ=255)
       gtgtCGCAGGCCAGGCGTTCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCtgt                                                                                                                                                 <  2:1870181/150‑1 (MQ=255)
          tCGCAGGCCAGGCGTTCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACCGGGCTGTggg                                                                                                                                              >  1:1717720/1‑150 (MQ=255)
          tCGCAGGCCAGGCGTTCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTggg                                                                                                                                              >  2:3656581/1‑150 (MQ=255)
               ggccaggcGTTCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCa                                                                                                                                        <  1:3013887/150‑1 (MQ=255)
                  ggccgagccagcctgcgacgcgcccctgagcggagtcgcaggccatgccagaagcattcgatgcaaccccggaatctgcaaagcctgcaccgacgTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATAc                                                                                                                                     >  2:4649714‑M1/96‑150 (MQ=255)
                    ggcGTTCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACAt                                                                                                                                   <  1:3714625/150‑1 (MQ=255)
                        ttCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCtg                                                                                                                               >  2:5912493/1‑150 (MQ=255)
                         ccagcctgcgacgcgcccctgagcggagtcgcaggccatgccagaagcattcgatgcaaccccggaatctgcaaagcctgcaccgacgTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCtgt                                                                                                                              >  1:611573‑M1/89‑150 (MQ=255)
                          ccAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCtgtg                                                                                                                             <  1:662688/150‑1 (MQ=255)
                          ccAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCtgtg                                                                                                                             <  1:1670745/150‑1 (MQ=255)
                          cagcctgcgacgcgcccctgagcggagtcgcaggccatgccagaagcattcgatgcaaccccggaatctgcaaagcctgcaccgacgTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCtgtg                                                                                                                             >  1:5311516‑M1/88‑150 (MQ=255)
                          cagcctgcgacgcgcccctgagcggagtcgcaggccatgccagaagcattcgatgcaaccccggaatctgcaaagcctgcaccgacgTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCtgtg                                                                                                                             >  2:7309374‑M1/88‑150 (MQ=255)
                          cagcctgcgacgcgcccctgagcggagtcgcaggccatgccagaagcattcgatgcaaccccggaatctgcaaagcctgcaccgacgTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCtgtg                                                                                                                             >  1:1169560‑M1/88‑150 (MQ=255)
                           agcctgcgacgcgcccctgagcggagtcgcaggccatgccagaagcattcgatgcaaccccggaatctgcaaagcctgcaccgacgTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTgg                                                                                                                            >  1:4470519‑M1/87‑150 (MQ=255)
                            gcctgcgacgcgcccctgagcggagtcgcaggccatgccagaagcattcgatgcaaccccggaatctgcaaagcctgcaccgacgTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTggg                                                                                                                           >  2:3857439‑M1/86‑150 (MQ=255)
                            aGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCGGTggg                                                                                                                           >  2:2920384/1‑150 (MQ=255)
                             ggCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGa                                                                                                                          <  1:1880851/150‑1 (MQ=255)
                             ggCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGGGGGCTGCATACATCCGGTGGGa                                                                                                                          >  2:5176008/1‑150 (MQ=255)
                             cctgcgacgcgcccctgagcggagtcgcaggccatgccagaagcattcgatgcaaccccggaatctgcaaagcctgcaccgacgTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGa                                                                                                                          >  2:4249319‑M1/85‑150 (MQ=255)
                             cctgcgacgcgcccctgagcggagtcgcaggccatgccagaagcattcgatgcaaccccggaatctgcaaagcctgcaccgacgTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACACGTCTGTGGGCTGCATACATCCTGTGGGa                                                                                                                          >  2:390397‑M1/85‑150 (MQ=255)
                              gCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATcaca                                                                                                                                                        <  1:2205042/120‑1 (MQ=255)
                              gCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCACCCACAGGTACAGAACAGGTGTGTGGGCTGCGTACATACTGTGGGAg                                                                                                                         >  1:3118924/1‑150 (MQ=255)
                              gCAGTCAAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAg                                                                                                                         <  1:1038720/150‑1 (MQ=255)
                               cAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGc                                                                                                                        <  2:4065876/150‑1 (MQ=255)
                               cAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGc                                                                                                                        <  1:811270/150‑1 (MQ=255)
                                 gTCGAGCAGGCGCTCGTCTGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCgg                                                                                                                      <  1:7337413/150‑1 (MQ=255)
                                 gTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCgg                                                                                                                      <  2:4332974/150‑1 (MQ=255)
                                 gTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCgg                                                                                                                      <  2:5011481/150‑1 (MQ=255)
                                 gTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAAATTTTGCTGGGTACTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCgg                                                                                                                      <  1:1524793/150‑1 (MQ=255)
                                 cgacgcgcccctgagcggagtcgcaggccatgccagaagcattcgatgcaaccccggaatctgcaaagcctgcaccgacgTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCgg                                                                                                                      >  1:1578888‑M1/81‑150 (MQ=255)
                                 cgacgcgcccctgagcggagtcgcaggccatgccagaagcattcgatgcaaccccggaatctgcaaagcctgcaccgacgTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCgg                                                                                                                      >  1:6249095‑M1/81‑150 (MQ=255)
                                 cgacgcgcccctgagcggagtcgcaggccatgccagaagcattcgatgcaaccccggaatctgcaaagcctgcaccgacgTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCgg                                                                                                                      >  1:1402678‑M1/81‑150 (MQ=255)
                                 cgacgcgcccctgagcggagtcgcaggccatgccagaagcattcgatgcaaccccggaatctgcaaagcctgcaccgacgTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCgg                                                                                                                      >  1:6899194‑M1/81‑150 (MQ=255)
                                  tCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCggg                                                                                                                     <  1:3283784/150‑1 (MQ=255)
                                  tCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGGGGGCTGCATACATCCTGTGGGAGCggg                                                                                                                     >  2:6148561/1‑150 (MQ=255)
                                  gacgcgcccctgagcggagtcgcaggccatgccagaagcattcgatgcaaccccggaatctgcaaagcctgcaccgacgTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCggg                                                                                                                     >  1:5953477‑M1/80‑150 (MQ=255)
                                    gAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACCGGTCTGGGGGCTGCATACATCCGGTGGGAGCGGGCa                                                                                                                   >  1:2977036/1‑150 (MQ=255)
                                         ccctgagcggagtcgcaggccatgccagacgcattcgatgcaaccccggaatcggcaaagcctgcaccgccgTCTTCGCGGGCAGGCCCGCTCCCACCGGTCCATCACCGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGccc                                                                                                              >  1:111495‑M1/73‑150 (MQ=255)
                                         ccctgagcggagtcgcaggccatgccagaagcattcgatgcaaccccggaatctgcaaagcctgcaccgacgTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGccc                                                                                                              >  1:3922426‑M1/73‑150 (MQ=255)
                                              tCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCCCAGGTCCATCACCGGGCAGTGGGCGGCAAACAACCGGGGGGAGCGGGCGGGCCCCCGaa                                                                                                         >  1:7140951/1‑150 (MQ=255)
                                                ggccggggCTTGCAGCCGGGCCTGGGCCTGGTGGGCGCGGGGGTCGGGGGAAATTTTTGCTGGGTTCTTCGCGGGCAGGCCCGCGCCCACCGGTCCATCCCCGGTCGGTGGGCTGCATACCTCGGGTGGGAGCGGGCATGCCCGCGAAGa                                                                                                       <  2:7140951/144‑1 (MQ=255)
                                                 tCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCCTGCCCGCgacgac                                                                                                      >  2:6591068/1‑150 (MQ=255)
                                                     tGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTCGGCGGAAATGTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCCGGCCCGCGACGACGGCg                                                                                                  >  2:6972542/1‑150 (MQ=39)
                                                     tGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTTCATCCCAGGTCTGTGGGCTGCGTAAATACAGTGGGGTCGGGGCTGCCCGCGCAGACGGCg                                                                                                  >  1:5644058/1‑150 (MQ=255)
                                                     tGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATa                                                                                                                                      <  2:5739265/114‑1 (MQ=255)
                                                     tGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGGCg                                                                                                  >  2:785323/1‑150 (MQ=255)
                                                         ttGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCtgt                                                                                                                              <  2:399609/118‑1 (MQ=255)
                                                         ttGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCGCCCACAGGTACATCACAGGTCTGTGGGCTGGCTTACTCCTGGGGGAGCGGGCCGGCCCGGGGCGACGGCGGTgg                                                                                              >  1:5839170/1‑149 (MQ=255)
                                                             aGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGc                                                                                          >  2:1240478/1‑150 (MQ=255)
                                                             aGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCGGTGGGCTGCATACAACCGGTGGGAGCGGGCCTGCCCGCGAAGACGACGGTGCAGGc                                                                                          >  1:1541698/1‑150 (MQ=39)
                                                                 gggCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTGCGTCCCCGGTCTGTGGGCTGCATACAGCCGGGGGGAGCGGGCCTGCCCGCGCCGCCGACGGTGCAGGCtggg                                                                                      >  1:7417370/1‑147 (MQ=255)
                                                                 gggCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATACCAGGTCGGTGGGCAGCCTACATCCGGTGGGGGCGGGCGTGGCCGCGCCGGCGGCCGTGCAGGcggtg                                                                                      >  1:6553958/1‑146 (MQ=255)
                                                                  ggCGTGGGCCTGGCGGTCGGGGGTGTAGGCGTAAATTTTTGCGGGGTTCTTGGCGGGCTGGCCCGCTCCCACCGGTACATCACCGGTCTGGGGGCGGAATACACCGTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGc                                                                                     <  2:1950597/150‑1 (MQ=255)
                                                                     cTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATAACAGGTCGGTGGGCTGCTTACATCGAGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCGGTGCAgc                                                                                  >  2:3741971/1‑149 (MQ=255)
                                                                     cTGGGACTGGCGGGCGCGGGGGTCGGCGCATTTTTGTGCTGGGTTCTTCGCGGCCCGGCCCGCCCCCCCGGGTACAGCACCGGTGGGTGGGCTGCATACATCCGGGGGGAGGGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGa                                                                                  <  2:298231/150‑1 (MQ=255)
                                                                      tGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCCTGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGAt                                                                                 >  2:467454/1‑150 (MQ=255)
                                                                          ttcgatgcaaccccggaatctgcaaagcctgcaccgacgTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCggg                                                                                                                     >  1:5628120‑M1/40‑110 (MQ=255)
                                                                          ttcgatgcaaccccggaatctgcaaagcctgcaccgacgTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCggg                                                                                                                     <  2:5628120‑M1/71‑1 (MQ=255)
                                                                          cctggaggCCGCGGGGGTCGGGGCAAATTTTTGCGGGGTTCTGCGCGGGCAGGCCCGCTCCCACTGGTACATCACCGGTCGGTGGGCGGCATCCAGCCGGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCg                                                                             <  2:5229401/147‑1 (MQ=255)
                                                                             ggaggTCGCGGGGGTAGGCGGAACTTTTGGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCgggg                                                                          <  2:6785522/150‑1 (MQ=255)
                                                                                       gggtttggCTGAATTTTTTGCTGGGTTGTTGGCGGGCAGGCCCGCTCCCACTGTTACCTCACCGGTGTGTGGGCTTAATACATCCTGGGGGTGGGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGaa                                                                <  2:3317760/144‑1 (MQ=255)
                                                                                        gggTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCgg                                                                                                                      <  1:2517240/95‑1 (MQ=255)
                                                                                        gggTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCgg                                                                                                                      >  2:2517240/1‑95 (MQ=255)
                                                                                                ggAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCCCAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTgg                                                       >  2:4838514/1‑150 (MQ=255)
                                                                                                ggAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATAAAGTGGGAGCGGGCATGCCCCCGAACACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTgg                                                       >  2:226049/1‑150 (MQ=39)
                                                                                                      ttttGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACAACACAGGTCGGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGc                                                 >  2:2186718/1‑150 (MQ=255)
                                                                                                            tGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCtgtg                                                                                                                             <  2:1597563/68‑1 (MQ=255)
                                                                                                                          gcttGCCCGCTCCCACTGGTACATCCCTGGTCTGTGGGCTGCATACATCCGTTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGgc                             <  2:3316956/147‑1 (MQ=255)
                                                                                                                          GCAGGCCCGCTCCCCCAGGTCCACCACCGGTCTGTGGGCTGCATCCATCCTGGGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGgc                             >  2:865406‑M1/1‑150 (MQ=255)
                                                                                                                                CCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGTATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCg                       <  1:2723741‑M1/150‑1 (MQ=255)
                                                                                                                                CCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCg                       <  1:4008224‑M1/150‑1 (MQ=255)
                                                                                                                                CCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCg                       <  1:2581608‑M1/150‑1 (MQ=255)
                                                                                                                                 CGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGc                      <  1:6061146‑M1/150‑1 (MQ=255)
                                                                                                                                 CGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGc                      <  1:493986‑M1/150‑1 (MQ=255)
                                                                                                                                  GCTCCCACGGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCa                     <  2:4481730‑M1/150‑1 (MQ=255)
                                                                                                                                    TCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAgg                   <  2:5754571‑M1/150‑1 (MQ=255)
                                                                                                                                    TCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAgg                   <  1:6715977‑M1/150‑1 (MQ=255)
                                                                                                                                      CCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGt                                                                                                    >  2:1170100‑M1/1‑67 (MQ=255)
                                                                                                                                      CCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAggct                 <  1:5488561‑M1/150‑1 (MQ=255)
                                                                                                                                       CACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAggctg                <  2:3564410‑M1/150‑1 (MQ=255)
                                                                                                                                        cccGGTACATCGCGGGTCTGTGGTAGGCATACCTCCTGTGGGCGCGGGAATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAggctgg               >  1:3656480‑M1/4‑150 (MQ=255)
                                                                                                                                            GTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCCGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAggctg                <  1:414256‑M1/145‑1 (MQ=255)
                                                                                                                                              aCAGCACCGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGc         <  2:5947185/150‑1 (MQ=255)
                                                                                                                                                aTCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCa       >  1:3179884/1‑150 (MQ=255)
                                                                                                                                                aTCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCa       >  2:2779510/1‑150 (MQ=255)
                                                                                                                                                 tCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGtt                                                                        <  2:1453544/84‑1 (MQ=255)
                                                                                                                                                 tCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAg      >  1:5366366/1‑150 (MQ=255)
                                                                                                                                                 tCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAg      >  1:6977060/1‑150 (MQ=255)
                                                                                                                                                  caccGGTCGGTGGGCTGCAGCCATCCGGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAgg     <  2:6220456/150‑1 (MQ=255)
                                                                                                                                                  cacaGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAgg     >  1:3452164/1‑150 (MQ=255)
                                                                                                                                                  cacaGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAgg     >  1:1252079/1‑150 (MQ=255)
                                                                                                                                                  cacaGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAgg     >  1:3451801/1‑150 (MQ=255)
                                                                                                                                                    caGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAgggg   >  2:4788367/1‑150 (MQ=255)
                                                                                                                                                     aGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAGGGGt  >  1:2385279/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GTCGGGGGTGTCGCAGGCCAGGCGTTCCAGGCAGTCGAGCAGGCGCTCGTCAGTGGCTTGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAGGGGT  >  NC_002947/3952793‑3953091

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: