Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,818,895 C→T 10.0% intergenic (+103/+26) PP_5092 → / ← pilT NLP/P60 protein/twitching motility protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,818,8950CT10.0% 279.5 / 16.3 100intergenic (+103/+26)PP_5092/pilTNLP/P60 protein/twitching motility protein
Reads supporting (aligned to +/- strand):  ref base C (48/42);  new base T (7/3);  total (55/45)
Fisher's exact test for biased strand distribution p-value = 5.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 101 total aligned reads displayed.

AGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGCCGCCCCGCCCTGCATCACTGAAGACAACTGCGCCAACCGCCCCTCCCGAAC  >  NC_002947/5818746‑5819043
                                                                                                                                                     |                                                                                                                                                    
aGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGgcgc                                                                                                                                                      <  2:3048279/150‑1 (MQ=255)
    ggTGCTGCCTAAAGAGCAGCTGCTGAGTCTGGCAAAGCTTTAAGCTTATGGGCGCTATCCGGTACCGGCCCCTTCGGGGGCACGCCCGCTCCCACTGGTAGCGGGTGGCTCTCGGGCGCAGCTCACATACCTGTGGGAGGGGGCGCGccc                                                                                                                                                  <  1:1438722/150‑1 (MQ=255)
      ggctgcCGAAAGAGCCGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGCGCGGGCGCGCCcgc                                                                                                                                                <  2:1834416/149‑1 (MQ=255)
          gcCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGaaa                                                                                                                                            >  2:3989073/1‑150 (MQ=255)
                      ggcggctgagggagcgcgggttgagccacggcttgtagtcccggcaacggccccttcgggggcaCGCCCGCTCCCACCGGTTGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAg                                                                                                                                >  1:315948‑M1/65‑150 (MQ=255)
                          cTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGtac                                                                                                                            <  2:2023365/150‑1 (MQ=255)
                            gAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCCCGCCCGCTCCCACAGGTAGCGCGTCGCTCACGCGCGCCGCTCACCTACCGGTGGGCGCGGGCCCGCCCGCGACAGGGCCTGTGCAGGtacag                                                                                                                          >  2:3435627/1‑148 (MQ=255)
                               tCTGGCCTAGCGTTAAGCTGATGGGCGCTAGCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACTGGTCGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGCGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACaa                                                                                                                       <  1:2902828/150‑1 (MQ=38)
                               tCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACCGGTAGCGCGGCGCACGCGCGCGCAGCACACCGACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCATTGCAGGTACTaaaa                                                                                                                       >  2:1176458/1‑150 (MQ=38)
                               tCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACaa                                                                                                                       <  2:3732982/150‑1 (MQ=255)
                               tCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACaa                                                                                                                       <  1:4878742/150‑1 (MQ=255)
                                    ccAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCtt                                                                                                                  <  2:5217521/150‑1 (MQ=255)
                                     cAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCGTTGCAGGTACTACAAGCCttt                                                                                                                 >  2:4034887/1‑150 (MQ=255)
                                      aaGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGGCGCTCTCGCGCGCAGCTCACCTACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACAACACGCCTTTg                                                                                                                >  2:2141499/1‑150 (MQ=255)
                                       aGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCCCAGGTAGCGCGTCGCTCTCGCGCGCAGCACACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGa                                                                                                               <  1:356592/150‑1 (MQ=255)
                                            tAAGCGGATCGGCGCTATCCGGTACCGGCCCCTCCGCGGGCCCGCCCGCGCCCCCTGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATCCCTGTGGGAGGGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAAc                                                                                                          <  1:4590653/150‑1 (MQ=255)
                                             aaGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCCCCCCCGCGCCCACCCGTAGCGCGTCGCGCGCGCGCGCCGCTCACCAACCGGGGGGGGCGGGCGCGCCCGCGCAAGGGCCTGTGGGGGTACAACAAGCCTTTGATCAAcc                                                                                                         >  2:4590653/1‑150 (MQ=255)
                                             aaGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAAcc                                                                                                         <  1:5432566/150‑1 (MQ=255)
                                                 tGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACCCCCGCTCCCACAGGTAGCCCGTCGCTCTCGCGCGCAGATAACCTAACTGTGGGAGCGGGCGCGCCCGCGAAACGGCCTTTTCCTGTACTACAACCCAGTGAACAACCCACg                                                                                                     >  2:2698499/1‑150 (MQ=37)
                                                               ctgcaaaggccctttcgcgggcgcgcccgctcccacaggtatgtgagctgcgcgcgagagcgacgcgcaaccggtgggagcgggcgTGCCCCCGAAGGGGCCGTTGCCGGGACTACAAGCCGTGGCTCAACCCGCGCTCCCTCAgccgcc                                                                                       >  1:315948‑M2/87‑150 (MQ=255)
                                                                    ccGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAAc                                                                                  <  1:4034887/150‑1 (MQ=255)
                                                                     cGGCCCCTTCGCGGGCACGCCCGCTCCCACACGTAGCGCGTCGCTCGCGCGCCCCGCGCACATACCGGGGGGCGCGGGCGCGCCCGCCAAAGGGCCGTTGCCGGGACTACAAGCCGTTGTTCAACCCGCGCGCCCTCAGCCACCGCAACg                                                                                 >  2:1179131/1‑150 (MQ=255)
                                                                       gCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGGGCGGGCGCGCCCCCGAAAGGGCCATTGCAGGTACTACAACCCGTTGATCAAACCGCGCCCCCACCGCCGCCGCACCGcc                                                                               >  2:5193370/1‑150 (MQ=37)
                                                                            ttCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAg                                                                                                                                           >  2:3481361/1‑85 (MQ=255)
                                                                            ttCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAg                                                                                                                                           <  1:3481361/85‑1 (MQ=255)
                                                                            ttCGCGGGCACGCCCGCTCCCACACGTATCGCGTCGCTCACGCGCGCAGCGCACAAAACGGTGGGGGCGGGCGCCCCCGCCAAAGGGCCGTTTCAGGGACAACAACCCGTTTATAAACCCACGCTCCCAAACCCACCGCAACGcccccgc                                                                          >  1:4374145/1‑148 (MQ=255)
                                                                               gcgGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCATTGTTCAACCCTCGCTCCCTCAGCCGCCGCAACGCCCCCGccct                                                                       >  1:4371719/1‑148 (MQ=255)
                                                                                   gCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTc                                                                   <  2:2343669/150‑1 (MQ=255)
                                                                                    cACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCa                                                                  <  1:3463920/150‑1 (MQ=255)
                                                                                    cACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCa                                                                  <  2:4574248/150‑1 (MQ=255)
                                                                                    cACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCa                                                                  <  1:5235491/150‑1 (MQ=255)
                                                                                           gctcccacaggtatgtgagctgcgcgcgagagcgacgcgctacctgtgggagcgggcgTGCCCGCGAAGGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATccc                                                           >  2:4086801‑M2/59‑150 (MQ=255)
                                                                                           GCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATccc                                                           <  1:832024‑M1/150‑1 (MQ=255)
                                                                                             TCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCt                                                         <  2:4393134‑M1/150‑1 (MQ=255)
                                                                                              CCCACCGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTc                                                        >  2:2405021‑M1/1‑150 (MQ=255)
                                                                                              CCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTc                                                        >  2:2405125‑M1/1‑150 (MQ=255)
                                                                                               ccacaggtatgtgagctgcgcgcgagagcgacgcgctacctgtgggagcgggcgTGCCCGCGAAGGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCAGCGTCAGCATCCCCTcc                                                       >  2:4469493‑M2/55‑150 (MQ=255)
                                                                                               ccacaggtatgtgagctgcgcgcgagagcgacgcgctacctgtgggagcgggcgCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTcc                                                       >  2:4116384‑M2/55‑150 (MQ=255)
                                                                                               ccacaggtatgtgagctgcgcgcgagagcgacgcgctacctgtgggagcgggcgCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTcc                                                       >  1:4925895‑M2/55‑150 (MQ=255)
                                                                                               ccacaggtatgtgagctgcgcgcgagagcgacgcgctacctgtgggagcgggcgCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCGcc                                                       >  1:3968134‑M2/55‑150 (MQ=255)
                                                                                               CCACTGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTcc                                                       >  1:3456577‑M1/1‑150 (MQ=255)
                                                                                               CCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGTGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTcc                                                       <  1:3545731‑M1/150‑1 (MQ=255)
                                                                                               CCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTcc                                                       <  2:4145320‑M1/150‑1 (MQ=255)
                                                                                               CCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCGATGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTcc                                                       <  2:3088728‑M1/150‑1 (MQ=255)
                                                                                                cacaggtatgtgagctgcgcgcgagagcgacgcgctacctgtgggagcgggcgTGCCCGCGAAGGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTccg                                                      >  2:3055536‑M2/54‑150 (MQ=255)
                                                                                                cacaggtatgtgagctgcgcgcgagagcgacgcgctacctgtgggagcgggcgTGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCtt                                                                                                                  >  1:3306795‑M2/54‑90 (MQ=255)
                                                                                                cacaggtatgtgagctgcgcgcgagagcgacgcgctacctgtgggagcgggcgTGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCtt                                                                                                                  <  2:3306795‑M2/37‑1 (MQ=255)
                                                                                                cacaggtatgtgagctgcgcgcgagagcgacgcgctacctgtgggagcgggcgCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTccg                                                      >  1:1722149‑M2/54‑150 (MQ=255)
                                                                                                CACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTccg                                                      <  2:4451147‑M1/150‑1 (MQ=255)
                                                                                                 acaggtatgtgagctgcgcgcgagagcgacgcgctacctgtgggagcgggcgTGCCCGCGAAGGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCAccgc                                                     >  1:2549342‑M2/53‑150 (MQ=255)
                                                                                                 acaggtatgtgagctgcgcgcgagagcgacgcgctacctgtgggagcgggcgCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTccgc                                                     >  2:5354043‑M2/53‑150 (MQ=255)
                                                                                                 acaggtatgtgagctgcgcgcgagagcgacgcgctacctgtgggagcgggcgCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTccgc                                                     >  2:5354088‑M2/53‑150 (MQ=255)
                                                                                                 ACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTccgc                                                     <  1:2859146‑M1/150‑1 (MQ=255)
                                                                                                  CAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTATGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTccgcc                                                    >  2:1992243‑M1/1‑150 (MQ=255)
                                                                                                   AGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGc                                                                <  1:3286894‑M1/137‑1 (MQ=255)
                                                                                                   AGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGc                                                                >  2:3286894‑M1/1‑137 (MQ=255)
                                                                                                     GTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCgccgcc                                                 <  2:4663399‑M1/150‑1 (MQ=255)
                                                                                                      gAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGCCGccc                                                >  1:4047364‑M1/2‑150 (MQ=255)
                                                                                                       atgtgagctgcgcgcgagagcgacgcgctacctgtgggagcgggcgCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGccgcccc                                               >  2:2047991‑M2/47‑150 (MQ=255)
                                                                                                        tgtgagctgcgcgcgagagcgacgcgctacctgtgggagcgggcgCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCACCGccgccccg                                              >  1:1654757‑M2/46‑150 (MQ=255)
                                                                                                         cgcgTCGCTCTCGCGCGCCGCTCCCATCCCTGTGGGCGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGccgccccgc                                             <  1:1409645/150‑1 (MQ=255)
                                                                                                         cgcgTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAAc                                                                                                          <  1:4606449/89‑1 (MQ=255)
                                                                                                         cgcgTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAAc                                                                                                          >  2:4606449/1‑89 (MQ=255)
                                                                                                            gTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGCCGCCCCGCCCt                                          >  2:815824/1‑150 (MQ=255)
                                                                                                                ctctCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGGACGACAAGCCGTTGATCAACCCTCGCTCCCTCCGCCTCCGCAACGCCCCCTCCAACGGCAGCAACCCCTCCCCCCCCCCGCCCGGcca                                      >  1:2118174/1‑148 (MQ=255)
                                                                                                                 tctcGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGCCGCCCCGCCCTGCATc                                     >  2:2454664/1‑150 (MQ=255)
                                                                                                                 tctcGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCAACGTCAGCATCCCCACCCCCGCCCCGCCCGGCaac                                     >  1:4899285/1‑148 (MQ=255)
                                                                                                                   tcGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCa                                                                                                                                 >  2:2648481/1‑56 (MQ=255)
                                                                                                                   tcGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCa                                                                                                                                 <  1:2648481/56‑1 (MQ=255)
                                                                                                                   tcGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGCCGCCCCGCCCTGCATCAc                                   >  1:3546158/1‑150 (MQ=255)
                                                                                                                   tcGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCAACGTCAGCATCCCCACCGCCGCCCCCCCCAGCAacac                                   >  1:2496672/1‑150 (MQ=255)
                                                                                                                   tcGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCACCATCGTCAGCATCCCCTCCGCCGCCCCCCCCAGCAacac                                   >  1:1059138/1‑150 (MQ=255)
                                                                                                                    cgcgcgCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGCCGCCCCGCCCTGCATCACt                                  <  1:1179131/150‑1 (MQ=255)
                                                                                                                     gcgcgcAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATc                                                             <  1:1425384/122‑1 (MQ=255)
                                                                                                                     gcgcgcAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATc                                                             >  2:1425384/1‑122 (MQ=255)
                                                                                                                      cgcgcAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGCCGCCCCGCCCTGCATCACTGa                                <  2:1231109/150‑1 (MQ=255)
                                                                                                                      cgcgcAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGCCGCCCCGCCCTGCATCACTGa                                <  1:1213618/150‑1 (MQ=255)
                                                                                                                       gcgcAGCTCACATACCTGTGGGTGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGCCGCCCCGCCCTGCATCACTGaa                               <  1:1176458/150‑1 (MQ=255)
                                                                                                                       gcgcAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGCCGCCCCGCCCTGCATCACTGaa                               <  2:3466220/150‑1 (MQ=255)
                                                                                                                       gcgcAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGCCGCCCCGCCCTGCATCACTGaa                               <  2:3687362/150‑1 (MQ=255)
                                                                                                                       gcgcAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGCCGCCCCGCCCTGCATCACTGaa                               <  2:4169491/150‑1 (MQ=255)
                                                                                                                       agagcgacgcgctacctgtgggagcgggcgTGCCCGCGAAGGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGCCGCCCCGCCCTGCATCACTGaa                               <  1:3055536‑M2/120‑1 (MQ=255)
                                                                                                                       agagcgacgcgctacctgtgggagcgggcgCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGCCGCCCCGCCCTGCATCACTGaa                               <  1:972347‑M2/120‑1 (MQ=255)
                                                                                                                          cAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGCCGCCCCGCCCTGCATCACTGAAGAc                            <  2:5321519/150‑1 (MQ=255)
                                                                                                                              tCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGCCGCCCCGCCCTGCATCACTGAAGACAACt                        >  2:947017/1‑150 (MQ=255)
                                                                                                                                acaTACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGCCGCCCCGCCCTGCATCACTGAAGACAACTgc                      <  1:1703554/150‑1 (MQ=255)
                                                                                                                                acaTACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGCCCCCCCCCCCAGCAACACTGAAGACCACAgc                      >  1:3461446/1‑150 (MQ=255)
                                                                                                                                    acctgtgggagcgggcgCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCGCCGCCGCCCCGCCCTGCATCACAGAAGACAACAGCgccc                  >  1:4370767‑M2/18‑149 (MQ=255)
                                                                                                                                      ctgtgggagcgggcgCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGCCGCCCCCCCCTGCATCACAGAAGACAACAGCGCCCAc                >  1:4804104‑M2/16‑150 (MQ=255)
                                                                                                                                      ctgtgggagcgggcgCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCACCGCCGCCCCGCCCTGCATCACTGAAGACAACAGCGCCAAc                >  1:4360777‑M2/16‑150 (MQ=255)
                                                                                                                                         tgggagcgggcgCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGCCGCCCCGCCCTGCATCACTGAAGACAACTGCGCCAACCGc             >  2:4734840‑M2/13‑150 (MQ=255)
                                                                                                                                         ctttcgcgggcgCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGCCGCCCCGCCCTGCATCACTGAAGACCACTGCGCCAACCGc             >  2:3363222‑M2/13‑150 (MQ=255)
                                                                                                                                          gggagcgggcgCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGCCGCCCCGCCCTGCATCACTGAAGACAACTGCGCCAACCGcc            <  1:1135351‑M2/139‑1 (MQ=255)
                                                                                                                                           ggagcgggcgCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGCCGCCCCGCCCTGCATCACTGAAGACAACTGCGCCAACCGccc           <  1:3989073‑M2/140‑1 (MQ=255)
                                                                                                                                           ggagcgggcgCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGCCGCCCCGCCCTGCATCACTGAAGACAACAGCGCCAACCGccc           >  2:4271895‑M2/11‑150 (MQ=255)
                                                                                                                                                  gcgTGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGCCGCCCCGCCCTGCATCACTGAAGACAACTGCGCCAACCGCCCCTCCCGa    >  1:2682405‑M2/4‑150 (MQ=255)
                                                                                                                                                  gcgCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGCCGCCCCGCCCTGCATCACTGAAGACAACTGCGCCAACCGccc           <  2:1597582‑M2/140‑1 (MQ=255)
                                                                                                                                                  gcgCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGCCGCCCCGCCCTGCATCACTGAAGACAACTGCGCCAACCGccc           >  1:1597582‑M2/4‑143 (MQ=255)
                                                                                                                                                   cgCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGCCGCCCCGCCCTGCATCACTGAAGACAACTGCGCCAACCGCCCCTCCCGaa   <  1:4113643‑M2/148‑1 (MQ=255)
                                                                                                                                                    gCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGCCGCCCCGCCCAGCAACACTGAAGACAACTGCGCCAACCGCCCCTCCCCAAc  >  1:4325995‑M2/2‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AGCGGGTGCTGCCGAAAGAGCAGCAGCTGAGTCTGGCCAAGCGTTAAGCTGATCGGCGCTATCCGGTACCGGCCCCTTCGCGGGCACGCCCGCTCCCACAGGTAGCGCGTCGCTCTCGCGCGCAGCTCACATACCTGTGGGAGCGGGCGCGCCCGCGAAAGGGCCTTTGCAGGTACTACAAGCCTTTGATCAACCCTCGCTCCCTCAGCCTCCGCAACGCCCCCTCCATCGTCAGCATCCCCTCCGCCGCCCCGCCCTGCATCACTGAAGACAACTGCGCCAACCGCCCCTCCCGAAC  >  NC_002947/5818746‑5819043

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: