Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,320,020 G→T 9.1% P92Q (CCG→CAG)  ftsH ← protease, ATP‑dependent zinc‑metallo

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,320,0200GT 9.1% 965.2 / 45.8 290P92Q (CCG→CAG) ftsHprotease, ATP‑dependent zinc‑metallo
Reads supporting (aligned to +/- strand):  ref base G (123/140);  new base T (8/19);  total (131/159)
Fisher's exact test for biased strand distribution p-value = 1.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 293 total aligned reads displayed.

CATGCGCGCTTTGCTCTTACCAAACGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTA  >  CP009273/3319872‑3320169
                                                                                                                                                    |                                                                                                                                                     
cATGCGCGCTTTGCTCTTACCAAACGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCgg                                                                                                                                                      <  1:4539443/150‑1 (MQ=255)
    cgcgcTTTGCTCTTACCAAACGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCa                                                                                                                                                  <  2:2922879/150‑1 (MQ=255)
    cgcgcTTTGCTCTTACCAAACGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCa                                                                                                                                                  >  1:2421049/1‑150 (MQ=255)
           tGCTCTTACCAAACGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACaa                                                                                                                                           >  2:3883295/1‑150 (MQ=255)
            gCTCTTACCAAACGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAAc                                                                                                                                          >  1:3391034/1‑150 (MQ=255)
             ctctTACCAAACGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCTGTTCACCGACAAcc                                                                                                                                         >  2:4236888/1‑150 (MQ=255)
               ctTACCAAACGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCtt                                                                                                                                       <  2:1436539/150‑1 (MQ=255)
                ttACCAAACGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCTGTTCACCGACAACCTTg                                                                                                                                      <  2:1370109/150‑1 (MQ=255)
                    cAAACGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGt                                                                                                                                  >  1:5031422/1‑150 (MQ=255)
                        cGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCtt                                                                                                                              <  2:2350317/150‑1 (MQ=255)
                        cGACATGGCACCTTTGCCACCGCCGCCATGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCtt                                                                                                                              >  2:366430/1‑150 (MQ=255)
                        cGACATGGCACCTGTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCtt                                                                                                                              <  2:4339062/150‑1 (MQ=255)
                            aTGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTc                                                                                                                          >  2:2941477/1‑150 (MQ=255)
                            aTGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTc                                                                                                                          >  2:3549738/1‑150 (MQ=255)
                             tGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTTGTTCTTCAGGCTGTTCACCGACAACCATGACGTGCTTGTTCa                                                                                                                         >  1:2682905/1‑150 (MQ=255)
                             tGGCACCTGTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCa                                                                                                                         <  2:4098013/150‑1 (MQ=255)
                              ggCACCTTTGCAACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGAAGTTCTTGGTCaa                                                                                                                        >  2:4417800/1‑150 (MQ=255)
                                cACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACa                                                                                                                      >  2:2491844/1‑150 (MQ=255)
                                     ttGCCACCGCCGCCCTGCATTTGACGCATGAAGATGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTa                                                                                                                 <  2:4553805/150‑1 (MQ=255)
                                      tGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCTGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTAt                                                                                                                <  2:1874606/150‑1 (MQ=255)
                                        ccACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATcc                                                                                                              >  1:2286997/1‑150 (MQ=255)
                                        ccACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATcc                                                                                                              >  1:2996947/1‑150 (MQ=255)
                                        ccACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATcc                                                                                                              >  1:2163236/1‑150 (MQ=255)
                                         cACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCa                                                                                                             >  2:4633664/1‑150 (MQ=255)
                                             gccgccCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTaa                                                                                                         <  1:665187/150‑1 (MQ=255)
                                               cgccCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAAtt                                                                                                       >  2:1230621/1‑150 (MQ=255)
                                                gccCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAAttt                                                                                                      <  2:1103286/150‑1 (MQ=255)
                                                gccCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAAttt                                                                                                      <  2:235288/150‑1 (MQ=255)
                                                  ccTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCg                                                                                                    >  1:1445120/1‑150 (MQ=255)
                                                   cTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCgg                                                                                                   >  2:3968485/1‑150 (MQ=255)
                                                     gCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGAt                                                                                                 <  1:4007643/150‑1 (MQ=255)
                                                         ttGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTg                                                                                             <  1:4328203/150‑1 (MQ=255)
                                                         ttGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTg                                                                                             >  2:2351554/1‑150 (MQ=255)
                                                         ttGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTg                                                                                             <  2:4596788/150‑1 (MQ=255)
                                                            aCGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAAc                                                                                          <  2:4370822/150‑1 (MQ=255)
                                                             cGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAAcc                                                                                         <  2:2991933/150‑1 (MQ=255)
                                                              gCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCg                                                                                        <  1:4134526/150‑1 (MQ=255)
                                                              gCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCg                                                                                        <  1:3623880/150‑1 (MQ=255)
                                                              gCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCg                                                                                        <  2:3819629/150‑1 (MQ=255)
                                                               cATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCTGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCgg                                                                                       >  1:1024443/1‑150 (MQ=255)
                                                                aTGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGa                                                                                      <  1:351217/150‑1 (MQ=255)
                                                                aTGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGa                                                                                      >  2:674701/1‑150 (MQ=255)
                                                                 tGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGaa                                                                                     <  2:936524/150‑1 (MQ=255)
                                                                  gaagaaGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAAt                                                                                    <  1:4827151/150‑1 (MQ=255)
                                                                  gaagaaGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAAt                                                                                    >  2:2115873/1‑150 (MQ=255)
                                                                   aagaagATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATg                                                                                   >  1:4391104/1‑150 (MQ=255)
                                                                        gaTCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGt                                                                              >  2:824364/1‑150 (MQ=255)
                                                                        gaTCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGt                                                                              >  2:2809156/1‑150 (MQ=255)
                                                                        gaTCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGt                                                                              >  2:4965852/1‑150 (MQ=255)
                                                                        gaTCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGt                                                                              <  2:5085611/150‑1 (MQ=255)
                                                                          tCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAAGCGGAATGTAAGTgg                                                                            >  2:4685143/1‑150 (MQ=255)
                                                                            cAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGta                                                                          >  1:832531/1‑150 (MQ=255)
                                                                            cAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGta                                                                          >  1:3934839/1‑150 (MQ=255)
                                                                                    aaTCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGtt                                                                  >  2:1161464/1‑150 (MQ=255)
                                                                                        aGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTa                                                              >  1:1199345/1‑150 (MQ=255)
                                                                                         gCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTAt                                                             >  1:4765475/1‑150 (MQ=255)
                                                                                            aCAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTAtctt                                                          >  1:2534517/1‑150 (MQ=255)
                                                                                            aCAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTAtctt                                                          <  1:4217950/150‑1 (MQ=255)
                                                                                             cAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCTGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTAtcttt                                                         <  1:3854521/150‑1 (MQ=255)
                                                                                             cAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTAtcttt                                                         <  1:3049371/150‑1 (MQ=255)
                                                                                              aGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTAtctttc                                                        >  2:240911/1‑150 (MQ=255)
                                                                                                cATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTAtctttctt                                                      >  2:2695002/1‑150 (MQ=255)
                                                                                                  tCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGACTGTAAGTGGTATAACGGTTACTATCTTTCTTgg                                                    >  2:3024817/1‑150 (MQ=255)
                                                                                                  tCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTgg                                                    <  1:2872076/150‑1 (MQ=255)
                                                                                                    gggAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTa                                                  <  1:4916063/150‑1 (MQ=255)
                                                                                                       aaCCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACg                                               >  2:2650398/1‑150 (MQ=255)
                                                                                                       aaCCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACg                                               <  2:743350/150‑1 (MQ=255)
                                                                                                              agaTGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCTGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGAttt                                        <  2:1969589/150‑1 (MQ=255)
                                                                                                              agaTGAAGATAGAAGCCAGCAGGCTTGGTTATTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGAttt                                        <  1:928533/150‑1 (MQ=255)
                                                                                                               gaTGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTc                                       <  2:5031422/150‑1 (MQ=255)
                                                                                                                aTGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCTGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCa                                      <  2:1065063/150‑1 (MQ=255)
                                                                                                                aTGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCa                                      >  1:4639889/1‑150 (MQ=255)
                                                                                                                 tGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCAc                                     >  2:3202748/1‑150 (MQ=255)
                                                                                                                 tGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCAc                                     <  2:4845080/150‑1 (MQ=255)
                                                                                                                 tGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCAc                                     >  2:531083/1‑150 (MQ=255)
                                                                                                                    aGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTc                                  >  1:443741/1‑150 (MQ=255)
                                                                                                                     gATAGAAGCCAGCAGGCTTGGTTCTTCAGGCTGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTcc                                 <  1:1107085/150‑1 (MQ=255)
                                                                                                                     gATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTcc                                 <  1:2446277/150‑1 (MQ=255)
                                                                                                                      aTAGAAGCCAGCAGGCTTGGTTCTTCAGGCTGTTCACCGACAACCTTGACGTTCTTGTTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCg                                <  1:4701130/150‑1 (MQ=255)
                                                                                                                      aTAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCg                                <  1:2783200/150‑1 (MQ=255)
                                                                                                                      aTAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCg                                <  2:3430520/150‑1 (MQ=255)
                                                                                                                      aTAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCg                                <  1:1364533/150‑1 (MQ=255)
                                                                                                                       tAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGt                               >  1:1813068/1‑150 (MQ=255)
                                                                                                                       tAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGt                               >  1:1759522/1‑150 (MQ=255)
                                                                                                                        aGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGtt                              >  1:3764803/1‑150 (MQ=255)
                                                                                                                         gAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTg                             <  1:1938185/150‑1 (MQ=255)
                                                                                                                          aaGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGa                            <  2:1107722/150‑1 (MQ=255)
                                                                                                                           aGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGAt                           <  2:4726264/150‑1 (MQ=255)
                                                                                                                             ccAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATAc                         <  1:3297445/150‑1 (MQ=255)
                                                                                                                                    gCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTc                  <  2:488933/150‑1 (MQ=255)
                                                                                                                                       tGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACg               >  2:3119219/1‑150 (MQ=255)
                                                                                                                                       tGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACg               >  1:1661509/1‑150 (MQ=255)
                                                                                                                                             ttcAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTg         >  2:2734843/1‑150 (MQ=255)
                                                                                                                                               cAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGt       <  2:4759256/150‑1 (MQ=255)
                                                                                                                                               cAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGt       >  2:1830308/1‑150 (MQ=255)
                                                                                                                                               cAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGt       <  2:1659001/150‑1 (MQ=255)
                                                                                                                                                aGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTc      <  1:1095421/150‑1 (MQ=255)
                                                                                                                                                 ggCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCg     <  2:219335/150‑1 (MQ=255)
                                                                                                                                                   cGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGtt   <  1:2077059/150‑1 (MQ=255)
                                                                                                                                                   cGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGtt   <  2:1537997/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CATGCGCGCTTTGCTCTTACCAAACGACATGGCACCTTTGCCACCGCCGCCCTGCATTTGACGCATGAAGAAGATCCAGACACCAATCAGCAACAGCATCGGGAACCAGGAGATGAAGATAGAAGCCAGCAGGCTTGGTTCTTCAGGCGGTTCACCGACAACCTTGACGTTCTTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTA  >  CP009273/3319872‑3320169

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