Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,503,471 A→T 6.9% intergenic (+112/‑107) PP_4834 → / → PP_4835 SPFH domain‑containing protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,503,4710AT 6.9% 325.6 / 10.7 102intergenic (+112/‑107)PP_4834/PP_4835SPFH domain‑containing protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (46/49);  new base T (3/4);  total (49/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

Only 100 of 111 total aligned reads displayed.

AGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTCCGCCGTACAATAGCCACCACGATTCCTGACCGGATCTCTCCATGCCCCCACGTCGGCACATTGCCTGGATAGCCTGCCTTGCAG  >  NC_002947/5503322‑5503620
                                                                                                                                                     |                                                                                                                                                     
agaCCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGca                                                                                                                                                       <  1:3910563/150‑1 (MQ=255)
  aCCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGcaca                                                                                                                                                     <  2:4582092/150‑1 (MQ=255)
     gTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCAcagca                                                                                                                                                  <  1:1905079/150‑1 (MQ=255)
     gTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCAcagca                                                                                                                                                  <  1:1905326/150‑1 (MQ=255)
         gcgGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTTGGCCCCAATAGCGACACAAAGCCCCACAGGACACCGCCGTAATCGAGGCTGCCGCAAATATGTTGGATCGGCCGTTTTTCGCGATGGGGTGCACAGCAGccc                                                                                                                                              >  1:1268869/1‑150 (MQ=255)
             gTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCTGTGCAGcccatcg                                                                                                                                          <  2:4610942/150‑5 (MQ=255)
              ttGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCTGTGCAGcccatcgc                                                                                                                                         <  2:2949631/150‑6 (MQ=255)
              ttGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGAGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCAtt                                                                                                                                         <  1:400741/150‑1 (MQ=255)
                   gggCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTTCGCCCCTATAGCGACAAAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAAGTCTGTTGGGGCGGCCTTGTTTGGCGATTGGGTGCACAGCAGCCCCCATTACCAg                                                                                                                                    >  2:78976/1‑150 (MQ=255)
                   gggCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTTCGCCCCTATAGCCACACAAAGCCCCACAGGAAAGCGCAGTAATCGAGGCTGCCCCAAATCGGTGGGAGCGGCCGTGTGTAGCGATTGGCTGCACAGCAGCCCCCAATaaaac                                                                                                                                    >  2:3456991/1‑146 (MQ=255)
                   gggCCGCAACAAGATTTGAGAGTGTCGGGGGCGCGTTTCGCCCCAATAGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAAATCAGTGGGGGCGGCCTAGTGTCGCGATGGGCTGCACAGCAGCCCCCAAAACCAg                                                                                                                                    >  2:3456956/1‑150 (MQ=255)
                     gCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGac                                                                                                                                  <  2:3243910/150‑1 (MQ=255)
                     gCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGac                                                                                                                                  <  2:3243819/150‑1 (MQ=255)
                      ccGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGacc                                                                                                                                 <  2:423093/150‑1 (MQ=255)
                          aacaaGATTTGAGAGTGTCGGGGGCGCTTTGGGCCCCAATAGCGGCACAAGGCCCCACAGGGCAGCGCAGTATTAGAGGCTGCCGCAAATCTTTTGGTGCGGCCGGGTGTCGCGATGGGCTGCACAGCAGCCCCCAATAACACaccaccg                                                                                                                             >  2:4144949/1‑149 (MQ=255)
                            caaGATTTTAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACTAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGATGCACAGCAGCCCCCATTACCAGACCACCTGc                                                                                                                           <  2:4434705/150‑1 (MQ=255)
                            caaGATTTGAGAGTGTCGGGGGGGGTTTTGGCCCCTATGGCGTGACAATGGCCCACAGGTCATCGCAGTTTTGGTGGGTGGGGGAATTGTGTGGGTGCGGGCTTGTGTGGGGTTTGGCTGGACAGCGGGCCCCATTAACAGtacaccagg                                                                                                                           >  2:404001/1‑147 (MQ=255)
                            caaGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGGGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGc                                                                                                                           <  2:3015060/150‑1 (MQ=255)
                                 tttGTGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACACGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCAc                                                                                                                      <  1:1560220/150‑1 (MQ=255)
                                           tCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCAc                                                                                                            <  2:1789675/150‑1 (MQ=255)
                                           tCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCAc                                                                                                            <  2:1789858/150‑1 (MQ=255)
                                              ggggCGCTTTGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCAcccc                                                                                                         <  2:1751599/150‑1 (MQ=255)
                                                     tttGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCAttt                                                                                                  >  1:4238368/1‑150 (MQ=255)
                                                       tGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCAttttt                                                                                                <  2:4238368/150‑1 (MQ=255)
                                                        gcgcCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTTGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTa                                                                                               >  1:2193541/1‑150 (MQ=255)
                                                         cgcCCCTATCGCTACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGAAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTAc                                                                                              <  1:2429648/150‑1 (MQ=255)
                                                         cgcCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGTGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTAc                                                                                              <  1:9612/150‑1 (MQ=255)
                                                         cgcCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTAc                                                                                              <  1:1429529/150‑1 (MQ=255)
                                                          gcCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTAcc                                                                                             >  1:2263205/1‑150 (MQ=255)
                                                          gcCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTAcc                                                                                             >  1:2263211/1‑150 (MQ=255)
                                                          gcCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTAcc                                                                                             >  1:2932764/1‑150 (MQ=255)
                                                          gcCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTAcc                                                                                             <  1:4317642/150‑1 (MQ=255)
                                                          gcCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTAcc                                                                                             <  1:4317705/150‑1 (MQ=255)
                                                            cccTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTTGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCta                                                                                           <  2:2193541/150‑1 (MQ=255)
                                                            cccTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCt                                                                                                                             <  2:130396/116‑1 (MQ=255)
                                                            cccTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCt                                                                                                                             >  1:130396/1‑116 (MQ=255)
                                                             ccTCTCGCGCCCCACGGCCCCACAGGACAGCGCCGTATTCGAGGCCGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACCGCACCCCCCACGACCAGCCCACCAGCGTCACGCCCCCGCACCCCGGCATTATTACCtat                                                                                          <  2:302161/150‑1 (MQ=255)
                                                             ccTCTCGCGACACAAGGCCCCACTGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGAGATGGGCTGCACAGAAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCtat                                                                                          <  2:302238/150‑1 (MQ=255)
                                                              cTATCGCGACACAAGGCCCCACAGGACAGCGCAGTTTTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCtata                                                                                         <  1:2305011/150‑1 (MQ=255)
                                                              cTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCtata                                                                                         <  1:2396673/150‑1 (MQ=255)
                                                              cTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCtata                                                                                         >  1:3958400/1‑150 (MQ=255)
                                                                 tCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTc                                                                                      <  1:3269279/150‑1 (MQ=255)
                                                                 tCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGACCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTAACTATACTc                                                                                      >  2:3504159/1‑150 (MQ=255)
                                                                  cgcgACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTcc                                                                                     >  2:497671/1‑150 (MQ=255)
                                                                   gcgACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTccg                                                                                    >  1:2618798/1‑150 (MQ=255)
                                                                   gcgACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTccg                                                                                    >  1:5039790/1‑150 (MQ=255)
                                                                          aaGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTCCGCCGTACa                                                                             >  2:1116751/1‑150 (MQ=255)
                                                                          aaGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTCCGCCGTACa                                                                             >  2:383819/1‑150 (MQ=255)
                                                                          aaGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTCCGCCGTACa                                                                             >  2:383865/1‑150 (MQ=255)
                                                                          aaGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTCCGCCGTACa                                                                             >  1:1219544/1‑150 (MQ=255)
                                                                          aaGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTCCGCCGTACa                                                                             >  1:2235212/1‑150 (MQ=255)
                                                                           aGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTCCGCCGTACaa                                                                            <  2:2741508/150‑1 (MQ=255)
                                                                            ggCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCTGTGCAGcccat                                                                                                                                            <  1:48723/85‑3 (MQ=25)
                                                                            ggCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCTGTGCAGcccat                                                                                                                                            >  2:48723/1‑83 (MQ=25)
                                                                            ggCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCTGTGCAGcccat                                                                                                                                            <  1:48880/85‑3 (MQ=25)
                                                                            ggCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCTGTGCAGcccat                                                                                                                                            >  2:48880/1‑83 (MQ=25)
                                                                               cccACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTCCGCCGTACAATAGc                                                                        <  2:4551287/150‑1 (MQ=255)
                                                                                  acaGGACAGCGCTGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTCCGCCGTACAATAGccac                                                                     <  1:563463/150‑1 (MQ=255)
                                                                                  acaGGACAGCGCATTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTCCGCCGTACAATAGccac                                                                     <  1:1070510/150‑1 (MQ=255)
                                                                                    aGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTCCGCCGTACAATAGccacca                                                                   <  1:4765792/150‑1 (MQ=255)
                                                                                       aCAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTCCGCCGTACAATAGCCACCACGa                                                                >  1:4391237/1‑150 (MQ=255)
                                                                                        cAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTCCGCCGTACAATAGCCACCACGAt                                                               <  2:4005021/150‑1 (MQ=255)
                                                                                        cAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTCCGCCGTACAATAGCCACCACGAt                                                               <  1:3089694/150‑1 (MQ=255)
                                                                                          gcgcAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTCCGCCGTACAATAGCCACCACGATTc                                                             >  1:4770112/1‑150 (MQ=255)
                                                                                          gcgcAGTATTCGAGGCTGCCGCAATACTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACAGCACCCCGGCATTTTTACCTATACTCCCCCGTACAATAGCCACCACGATTc                                                             <  2:2063232/150‑1 (MQ=255)
                                                                                            gcAGTATTCGAGGCTGCCTCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTCCGCCGTACAATAGCCACCACGATTCCt                                                           <  2:3298789/150‑1 (MQ=255)
                                                                                            gcAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTCCGCCGTACAATAGCCACCACGATTCCt                                                           <  1:2249122/150‑1 (MQ=255)
                                                                                              aGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTCCGCCGTACAATAGCCACCACGATTCCTGa                                                         <  2:3686362/150‑1 (MQ=255)
                                                                                                   tCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTCCGCCGTACAATAGCCACCACGATTCCTGACCGGa                                                    <  2:2100986/150‑1 (MQ=255)
                                                                                                   tCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTCCGCCGTACAATAGCCACCACGATTCCTGACCGGa                                                    <  1:754880/150‑1 (MQ=255)
                                                                                                    cGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTCCGCCGTACAATAGCCACCACGATTCCTGACCGGAt                                                   <  1:5226618/150‑1 (MQ=255)
                                                                                                        gCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTCCGCCGTACAATAGCCACCACGATTCCTGACCGGAtctct                                               >  2:4069690/1‑150 (MQ=255)
                                                                                                           gCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTCCGCCGTACAATAGCCACCACGATTCCTGACCGGATCTCTCCa                                            >  1:3964570/1‑150 (MQ=255)
                                                                                                           gCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGAAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTCCGCCGTACAATATCCACCACGATTCCTGACCGGATCTCTCCa                                            >  1:878815/1‑150 (MQ=255)
                                                                                                           gCCGCAAATCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTCCGCCGTACAATAGCCACCACGATTCCAGACCGGAGCTCGCCa                                            >  1:1915850/1‑150 (MQ=255)
                                                                                                             cGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTCCGCCGTACAATAGCCACCACGATTCCTGACCGGATCTCTCCATg                                          >  2:3475746/1‑150 (MQ=255)
                                                                                                             cGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTCCGCCGTACAATAGCCACCACGATTCCTGACCGGATCTCTCCATg                                          >  1:5299610/1‑150 (MQ=255)
                                                                                                               cAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTCCGCCGTACAATAGCCACCACGATTCCTGACCGGATCTCTCCATGcc                                        >  1:4756070/1‑150 (MQ=255)
                                                                                                                      gtgGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTCCGCCGTACAATAGCCACCACGATTCCTGACCGGATCTCTCCATGCCCCCACGt                                 >  2:2925508/1‑150 (MQ=255)
                                                                                                                         ggAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTCCGCCGTACAATAGCCACCACGATTCCTGACCGGATCTCTCCATGCCCCCACGTCgg                              >  2:4989830/1‑150 (MQ=255)
                                                                                                                         ggAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTCCGCCGTACAATAGCCACCACGATTCCTGACCGGATCTCTCCATGCCCCCACGTCgg                              >  2:4989853/1‑150 (MQ=255)
                                                                                                                         ggAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTCCGCCGTACAATAGCCACCACGATTCCTGACCGGATCTCTCCATGCCCCCACGTCgg                              >  2:4990104/1‑150 (MQ=255)
                                                                                                                              ggCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTCCGCCGTACAATAGCCACCACGATTCCTGACCGGATCTCTCCATGCCCCCACGTCGGCACAt                         <  2:1268869/150‑1 (MQ=255)
                                                                                                                                  ttGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTCCGCCGTACAATAGCCACCACGATTCCTGACCGGATCTCTCCATGCCCCCACGTCGGCACATTGcc                     >  1:2806062/1‑150 (MQ=255)
                                                                                                                                   tgtgTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTCCGCCGTACAATAGCCACCACGATTCCTGACCGGATCTCTCCATGCCCCCACGTCGGCACATTGCCt                    <  2:392341/150‑1 (MQ=255)
                                                                                                                                         aatgggggCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTCCGCCGTACAATAGCCACCACGATTCCTGACCGGATCTCTCCATGCCCCCACGTCGGCACATTGCCTGGATAg              >  2:544186/6‑150 (MQ=255)
                                                                                                                                          cgATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTCCGCCGTACAATAGCCACCACGATTCCTGACCGGATCTCTCCATGCCCCCACGTCGGCACATTGCCTGGATAgc             >  2:1108434/1‑150 (MQ=255)
                                                                                                                                          cgATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTCCGCCGTACAATAGCCACCACGATTCCTGACCGGATCTCTCCATGCCCCCACGTCGGCACATTGCCTGGATAgc             <  1:4775868/150‑1 (MQ=255)
                                                                                                                                          cgATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTCCGCCGTACAATAGCCACCACGATTCCTGACCGGATCTCTCCATGCCCCCACGTCGGCACATTGCCTGGATAgc             >  2:1108415/1‑150 (MQ=255)
                                                                                                                                             ggggCTGCTGTGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTCCGCCGTACAATAGCCACCACGATTCCTGACCGGATCTCTCCATGCCCCCACGTCGGCACATTGCCTGGATAgcctg          >  2:1496718/2‑150 (MQ=255)
                                                                                                                                                 cTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTCCGCCGTACAATAGCCACCACGATTCCTGACCGGATCTCTCCATGCCCCCACGTCGGCACATTGCCTGGATAGCCTGCCtt      >  1:3421399/1‑150 (MQ=255)
                                                                                                                                                  tgaacaGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTCCGCCGTACAATAGCCACCACGATTCCTGACCGGATCTCTCCATGCCCCCACGTCGGCACATTGCCTGGATAGCCTGCCTTg     <  1:2522560/147‑1 (MQ=255)
                                                                                                                                                  tGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTCCGCCGTACAATAGCCACCACGATTCCTGACCGGATCTCTCCATGCCCCCACGTCGGCACATTGCCTGGATAgcc            <  1:220322/143‑1 (MQ=255)
                                                                                                                                                  tGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTCCGCCGTACAATAGCCACCACGATTCCTGACCGGATCTCTCCATGCCCCCACGTCGGCACATTGCCTGGATAGCCTGCCTTg     >  1:1163247/1‑150 (MQ=255)
                                                                                                                                                  tGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTCCGCCGTACAATAGCCACCACGATTCCTGACCGGATCTCTCCATGCCACCACGTCGGCACATTGCCTGGATAGCCTGCCTTg     >  1:427268/1‑150 (MQ=255)
                                                                                                                                                    cacaGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGTCATTTTTACCTATACTCCGCCGTACAATAGCCACCACGATTCCTGACCGGATCTCTCCATGCCCCCACGTCGGCACATTGCCTGGATAGCCTGCCTTGCa   <  1:404001/150‑1 (MQ=255)
                                                                                                                                                    cacaGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTCCGCCGTACAATAGCCACCACGATTCCTGACCGGATCTCTCCATGCCCCCACGTCGGCACATTGCCTGGATAGCCTGCCTTGCa   <  1:3456956/150‑1 (MQ=255)
                                                                                                                                                    cacaGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTCCGCCGTACAATAGCCACCACGATTCCTGACCGGATCTCTCCATGCCCCCACGTCGGCACATTGCCTGGATAGCCTGCCTTGCa   <  1:3456991/150‑1 (MQ=255)
                                                                                                                                                    cacaGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTCCGCCGTACAATAGCCACCACGATTCCTGACCGGATCTCTCCATGCCCCCACGTCGGCACATTGCCTGGATAGCCTGCCTTGCa   <  2:599327/150‑1 (MQ=255)
                                                                                                                                                     acaGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTCCGCCGTACAATAGCCACCACGATTCCTGACCGGATCTCTCCATGCCCCCACGTCGGCACATTGCCTGGATAGCCTGCCTTGCAg  >  1:502361/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGACCGTGCGCGGGTTGATGGGCCGCAACAAGATTTGAGAGTGTCGGGGGCGCTTTGCGCCCCTATCGCGACACAAGGCCCCACAGGACAGCGCAGTATTCGAGGCTGCCGCAATTCTGTGGGAGCGGCCTTGTGTCGCGATGGGCTGCACAGCAGCCCCCATTACCAGACCACCTGCGTCACTTCACCGCACCCCGGCATTTTTACCTATACTCCGCCGTACAATAGCCACCACGATTCCTGACCGGATCTCTCCATGCCCCCACGTCGGCACATTGCCTGGATAGCCTGCCTTGCAG  >  NC_002947/5503322‑5503620

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: