Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,292,914 C→A 52.6% D121Y (GAT→TAT)  yhbP ← UPF0306 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,292,9140CA52.6% 50.5 / inf 298D121Y (GAT→TAT) yhbPUPF0306 family protein
Reads supporting (aligned to +/- strand):  ref base C (84/57);  new base A (87/70);  total (171/127)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 299 total aligned reads displayed.

TGATAAAGATCAGCTGGTTACCAGCCGGACACCGGACGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGAATTAACGCTACCGTTTTCG  >  CP009273/3292767‑3293063
                                                                                                                                                   |                                                                                                                                                     
  aTAAAGATCAGCTGGTTACCAGCCGGACACCGGACGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGc                                                                                                                                                   >  2:2195450/1‑150 (MQ=255)
       gATCAGCTGGTTACCAGCCGGACACCGGACGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGAt                                                                                                                                              >  1:440127/1‑150 (MQ=255)
        aTCAGCTGGTTACCAGCCGGACACCGGACGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATAGAGCCGGAtt                                                                                                                                             <  2:4145002/150‑1 (MQ=255)
        aTCAGCTGGTTACCAGCCGGACACCGGACGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATAGAGCCGGAtt                                                                                                                                             <  2:4301264/150‑1 (MQ=255)
         tCAGCTGGGTACCAGCCGGACACCGGACGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGAttt                                                                                                                                            <  1:4290676/150‑1 (MQ=255)
          cAGCTGGTTACCAGCCGGACACCGGACGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTc                                                                                                                                           >  2:704607/1‑150 (MQ=255)
            gCTGGTTACCAGCCGGACACCGGACGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTccc                                                                                                                                         <  1:3362732/150‑1 (MQ=255)
            gCTGGTTACCAGCCGGACACCGGACGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATAGAGCCGGATTTccc                                                                                                                                         <  1:659304/150‑1 (MQ=255)
                gTTACCAGCCGGACACCGGACGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATAGAGCCGGATTTCCCacac                                                                                                                                     <  1:3926073/150‑1 (MQ=255)
                 ttACCAGCCGGACACCGGACGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATAGAGCCGGATTTCCCACAcc                                                                                                                                    <  2:3792109/150‑1 (MQ=255)
                       gCCGGACACCGGACGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATAGAGCCGGATTTCCCACACCGGGGcc                                                                                                                              <  1:2701524/150‑1 (MQ=255)
                        ccGGACACCGGACGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCg                                                                                                                             >  1:3095902/1‑150 (MQ=255)
                        ccGGACACCGGACGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATAGAGCCGGATTTCCCACACCGGGGCCg                                                                                                                             >  2:2002435/1‑150 (MQ=255)
                        ccGGACACCGGACGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGAGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATAGAGCCGGATTTCCCACACCTGGGCCg                                                                                                                             >  2:4252180/1‑150 (MQ=255)
                         cGGACACCGGACGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGa                                                                                                                            >  2:1753338/1‑150 (MQ=255)
                          ggACACCGGACGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGAc                                                                                                                           >  2:2612914/1‑150 (MQ=255)
                            acacCGGACGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAg                                                                                                                         <  1:3416413/150‑1 (MQ=255)
                                 ggACGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTc                                                                                                                    <  1:1455042/150‑1 (MQ=255)
                                 ggACGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATAGAGCCGGATTTCCCACACCGGGGCCGACAGCATTc                                                                                                                    >  2:892112/1‑150 (MQ=255)
                                   aCGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTg                                                                                                                  <  1:4214666/150‑1 (MQ=255)
                                   aCGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTg                                                                                                                  <  1:2788718/150‑1 (MQ=255)
                                   aCGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATAGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTg                                                                                                                  >  2:3835038/1‑150 (MQ=255)
                                   aCGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATAGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTg                                                                                                                  >  2:1667137/1‑150 (MQ=255)
                                     gATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGc                                                                                                                >  1:2852094/1‑150 (MQ=255)
                                         tGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATAGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAAcc                                                                                                            <  1:519791/150‑1 (MQ=255)
                                           ccAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCgg                                                                                                          >  2:4492244/1‑150 (MQ=255)
                                             aGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATAGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGaa                                                                                                        >  1:632122/1‑150 (MQ=255)
                                             aGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATAGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGaa                                                                                                        >  1:632904/1‑150 (MQ=255)
                                             aGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATAGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGaa                                                                                                        >  1:693653/1‑150 (MQ=255)
                                               cGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAg                                                                                                      >  1:3834573/1‑150 (MQ=255)
                                               cGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAg                                                                                                      >  1:4473276/1‑150 (MQ=255)
                                                gTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATAGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGc                                                                                                     <  1:3242592/150‑1 (MQ=255)
                                                     aCCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATAGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGcgacga                                                                                                <  1:3739575/150‑1 (MQ=255)
                                                      ccGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATAGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGAt                                                                                               <  1:243796/150‑1 (MQ=255)
                                                           aGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATAGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTAc                                                                                          <  2:719767/150‑1 (MQ=255)
                                                                ttACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATAGAGCCGGATTTTCCACACCGGGGCCGACCGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCtt                                                                                     >  1:1239162/1‑150 (MQ=255)
                                                                ttACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATAGAGCCGGATTTTCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCtt                                                                                     >  1:3703268/1‑150 (MQ=255)
                                                                     ccTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATAGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTgcgcg                                                                                <  1:3109829/150‑1 (MQ=255)
                                                                       tGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGcgcgcg                                                                              >  1:475729/1‑150 (MQ=255)
                                                                       tGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATAGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGcgcgcg                                                                              >  1:24892/1‑150 (MQ=255)
                                                                        gCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATAGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGa                                                                             >  1:3120001/1‑150 (MQ=255)
                                                                         cTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAg                                                                            >  2:3815423/1‑150 (MQ=255)
                                                                         cTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATAGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAg                                                                            <  1:4112757/150‑1 (MQ=255)
                                                                          tCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAgg                                                                           >  2:2421800/1‑150 (MQ=255)
                                                                          tCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAgg                                                                           >  1:1765205/1‑150 (MQ=255)
                                                                             gTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATAGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCg                                                                        <  2:2142464/150‑1 (MQ=255)
                                                                              tGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATAGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGc                                                                       >  2:2792334/1‑150 (MQ=255)
                                                                              tGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATAGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGc                                                                       >  2:506880/1‑150 (MQ=255)
                                                                                   gAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTcc                                                                  >  1:1921929/1‑150 (MQ=255)
                                                                                    aGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCt                                                                 <  1:2131673/150‑1 (MQ=255)
                                                                                       cgcgTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCAc                                                              >  1:3799232/1‑150 (MQ=255)
                                                                                          gTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATAGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCtt                                                           <  1:1079677/150‑1 (MQ=255)
                                                                                          gTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATAGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCtt                                                           >  1:302633/1‑150 (MQ=255)
                                                                                          gTACCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATAGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCtt                                                           >  2:3190131/1‑150 (MQ=255)
                                                                                            aGCCAAATCATTTATTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTcc                                                         <  1:4413182/150‑1 (MQ=255)
                                                                                             cccAACTCAATTTCTTACCCCAGCCCAGCGTGCTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCa                                                        <  1:3279172/149‑1 (MQ=255)
                                                                                              ccAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAg                                                       >  1:1025290/1‑150 (MQ=255)
                                                                                              ccAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATAGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAg                                                       >  1:3323133/1‑150 (MQ=255)
                                                                                                aaaTCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGcc                                                     <  2:3445228/150‑1 (MQ=255)
                                                                                                aaaTCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATAGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGcc                                                     >  2:2974682/1‑150 (MQ=255)
                                                                                                aaaTCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATAGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGcc                                                     >  2:2630111/1‑150 (MQ=255)
                                                                                                  aTCATTTTTTGACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTg                                                   <  2:2796898/150‑1 (MQ=255)
                                                                                                    cATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCg                                                 >  2:3961659/1‑150 (MQ=255)
                                                                                                     aTTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATAGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCgg                                                >  2:878169/1‑150 (MQ=255)
                                                                                                        tttttACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGAtc                                             >  2:752430/1‑150 (MQ=255)
                                                                                                         ttttACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATAGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGAtct                                            <  2:4374731/150‑1 (MQ=255)
                                                                                                            tACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCAc                                         <  1:83990/150‑1 (MQ=255)
                                                                                                              ccAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCt                                       >  2:3829947/1‑150 (MQ=255)
                                                                                                               cAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCtt                                      <  1:2694824/150‑1 (MQ=255)
                                                                                                               cAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCtt                                      <  1:2042004/150‑1 (MQ=255)
                                                                                                               cAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATAGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCtt                                      >  2:1468813/1‑150 (MQ=255)
                                                                                                                caaGCCCAGCGTGTTGTCGGTGAATTTGATTTCATAGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCttt                                     <  2:3407494/149‑1 (MQ=255)
                                                                                                                 aaGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTa                                    >  1:3559023/1‑150 (MQ=255)
                                                                                                                 aaGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTa                                    <  1:2613259/150‑1 (MQ=255)
                                                                                                                  aGCCCAGCGTGTTGTCGGTGAATTTGATTTCATAGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTaa                                   <  2:2864239/150‑1 (MQ=255)
                                                                                                                   gCCCAGCGTGTTGTCGGTGAATTTGATTTCATAGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTaaa                                  <  1:2973787/150‑1 (MQ=255)
                                                                                                                      cAGCGTGTTGTCGGTGAATTTGATTTCATAGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTg                               <  1:2240910/150‑1 (MQ=255)
                                                                                                                      cAGCGTGTTGTCGGTGAATTTGATTTCATAGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTg                               >  2:1684201/1‑150 (MQ=255)
                                                                                                                       aGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGc                              >  1:121551/1‑150 (MQ=255)
                                                                                                                         cGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGcac                            <  1:2338380/150‑1 (MQ=255)
                                                                                                                          gtgtTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGcaca                           >  1:795453/1‑150 (MQ=255)
                                                                                                                          gtgtTGTCGGTGAATTTGATTTCATAGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGcaca                           >  2:3257443/1‑150 (MQ=255)
                                                                                                                          gtgtTGTCGGTGAATTTGATTTCATAGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGcaca                           >  2:2930113/1‑150 (MQ=255)
                                                                                                                           tgttgtCGGTGAATTTGATTTCATAGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCacac                          >  2:2120642/1‑150 (MQ=255)
                                                                                                                            gttgtCGGTGAATTTGATTTCATAGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCACAcc                         <  1:1957972/150‑1 (MQ=255)
                                                                                                                            gttgtCGGTGAATTTGATTTCATAGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACATTCCAGCCTGCGGATCTCACCTTTAAACTGCACAcc                         >  2:933944/1‑150 (MQ=255)
                                                                                                                                tCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGa                     <  1:3750971/150‑1 (MQ=255)
                                                                                                                                  ggTGAATTTGATTTCATAGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGAAt                   >  2:4241757/1‑150 (MQ=255)
                                                                                                                                  ggTGAATTTGATTTCATAGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGAAt                   <  1:892112/150‑1 (MQ=255)
                                                                                                                                   gTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGAAtt                  <  1:4374076/150‑1 (MQ=255)
                                                                                                                                    tGAATTTGATTTCATAGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGAATTa                 >  2:2785047/1‑150 (MQ=255)
                                                                                                                                      aaTTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGAATTAAc               >  1:1752835/1‑150 (MQ=255)
                                                                                                                                      aaTTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGAATTAAc               >  1:1191572/1‑150 (MQ=255)
                                                                                                                                      aaTTTGATTTCATAGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGAATTAAc               <  2:3639851/150‑1 (MQ=255)
                                                                                                                                       aTTTGATTTCATAGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGAATTAACg              >  2:3914552/1‑150 (MQ=255)
                                                                                                                                           gATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGAATTAACGCTAc          <  1:752430/150‑1 (MQ=255)
                                                                                                                                              ttCATAGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGAATTAACGCTACCGt       <  2:4016018/150‑1 (MQ=255)
                                                                                                                                              ttCATAGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGAATTAACGCTACCGt       <  1:1684357/150‑1 (MQ=255)
                                                                                                                                                cATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGAATTAACGCTACCGttt     >  2:2083644/1‑150 (MQ=255)
                                                                                                                                                cATCGAGCCGGAGTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGAATTAACGCTACCGttt     <  2:757716/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TGATAAAGATCAGCTGGTTACCAGCCGGACACCGGACGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGAATTAACGCTACCGTTTTCG  >  CP009273/3292767‑3293063

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: ‑