Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,785,915 G→T 63.8% intergenic (‑4/+271) kbl ← / ← yibB glycine C‑acetyltransferase/YibB family protein, function unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,785,9150GT63.8% 305.9 / inf 362intergenic (‑4/+271)kbl/yibBglycine C‑acetyltransferase/YibB family protein, function unknown
Reads supporting (aligned to +/- strand):  ref base G (61/70);  new base T (124/106);  total (186/176)
Fisher's exact test for biased strand distribution p-value = 1.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

Only 100 of 362 total aligned reads displayed.

ACAAAAGTTAATGACGTGGCTTCCATCAGCCACAGTGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCGATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTAC  >  CP009273/3785768‑3786061
                                                                                                                                                   |                                                                                                                                                  
acaAAAGTTAATGACGTGGCTTCCATCAGCCACAGTGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCGAt                                                                                                                                                  >  1:5112323/1‑150 (MQ=255)
acaAAAGTTAATGACGTGGCTTCCATCAGCCACAGTGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCGAt                                                                                                                                                  <  2:4434248/150‑1 (MQ=255)
acaAAAGTTAATGACGTGGCTTCCATCAGCCACAGTGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCGAt                                                                                                                                                  >  2:5003877/1‑150 (MQ=255)
acaAAAGTTAATGACGTGGCTTCCATCAGCCACAGTGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCGAt                                                                                                                                                  >  2:1485261/1‑150 (MQ=255)
 caAAAGTTAATGACGTGGCTTCCATCAGCCACAGTGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTAtt                                                                                                                                                 <  2:3476355/150‑1 (MQ=255)
  aaaaGTTAATGACGTGGCTTCCATCAGCCACAGTGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATtc                                                                                                                                                >  2:4038802/1‑150 (MQ=255)
  aaaaGTTAATGACGTGGCTTCCATCAGCCACAGTGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATtc                                                                                                                                                >  1:1303851/1‑150 (MQ=255)
    aaGTTAATGACGTGGCTTCCATCAGCCACAGTGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATtctc                                                                                                                                              >  2:3958194/1‑150 (MQ=255)
    aaGTTAATGACGTGGCTTCCATCAGCCACAGTGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCGATtctc                                                                                                                                              >  1:339880/1‑150 (MQ=255)
     aGTTAATGACGTGGCTTCCATCAGCCACAGTGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTcc                                                                                                                                             <  2:414530/150‑1 (MQ=255)
             aCGTGGCTTCCATCAGCCACAGTGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAggg                                                                                                                                     <  1:134484/150‑1 (MQ=255)
                  gCTTCCATCAGCCACAGTGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCGATTCTCCAGACAGGGCAAAt                                                                                                                                <  1:1568597/150‑1 (MQ=255)
                  gCTTCCATCAGCCACAGTGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCGATTCTCCAGACAGGGCAAAt                                                                                                                                <  1:4176312/150‑1 (MQ=255)
                     tCCATCAGCCACAGTGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCATTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCGATTCTCCAGACAGGGCAAATTcc                                                                                                                             >  1:3541154/1‑150 (MQ=255)
                      ccATCAGCCACAGTGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCa                                                                                                                            <  2:1162933/150‑1 (MQ=255)
                      ccATCAGCCACAGTGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCGATTCTCCAGACAGGGCAAATTCCa                                                                                                                            >  2:4532941/1‑150 (MQ=255)
                       cATCAGCCACAGTGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAg                                                                                                                           <  2:796607/150‑1 (MQ=255)
                       cATCAGCCACAGTGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCGATTCTCCAGACAGGGCAAATTCCAg                                                                                                                           <  1:768089/150‑1 (MQ=255)
                        aTCAGCCACAGTGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGc                                                                                                                          <  1:1498196/150‑1 (MQ=255)
                        aTCAGCCACAGTGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGc                                                                                                                          >  1:124392/1‑150 (MQ=255)
                        aTCAGCCACAGTGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCGATTCTCCAGACAGGGCAAATTCCAGc                                                                                                                          >  2:922008/1‑150 (MQ=255)
                         tCAGCCACAGTGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCGATTCTCCAGACAGGGCAAATTCCAGca                                                                                                                         <  2:1827000/150‑1 (MQ=255)
                             ccACAGTGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCGATTCTCCAGACAGGGCAAATTCCAGCACata                                                                                                                     <  2:4772105/150‑1 (MQ=255)
                               acaGTGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATAtt                                                                                                                   <  2:5105176/150‑1 (MQ=255)
                                caGTGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTa                                                                                                                  <  1:559794/150‑1 (MQ=255)
                                caGTGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTa                                                                                                                  <  2:2392034/150‑1 (MQ=255)
                                 aGTGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTAc                                                                                                                 >  1:5149609/1‑150 (MQ=255)
                                 aGTGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTAc                                                                                                                 <  1:1761998/150‑1 (MQ=255)
                                   tGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTAccc                                                                                                               >  1:2196130/1‑150 (MQ=255)
                                    gATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCGATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCa                                                                                                              >  2:2263684/1‑150 (MQ=255)
                                          tgcttgctGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCGATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTta                                                                                                        >  1:4609202/1‑150 (MQ=255)
                                          tgcttgctGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCGATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTta                                                                                                        <  2:3404122/150‑1 (MQ=255)
                                           gcttgctGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTtat                                                                                                       >  1:1915511/1‑150 (MQ=255)
                                            cttgctGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTtata                                                                                                      >  2:1454445/1‑150 (MQ=255)
                                             ttgctGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAg                                                                                                     >  2:4662653/1‑150 (MQ=255)
                                              tgctgcGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAgg                                                                                                    <  1:5194841/150‑1 (MQ=255)
                                               gctgcGCAGACGTAATAATGCGCTCTTCTTTAAACAACCATTCCGCCCGTGCGGTTTCCAGATAGTTGGTTAACTGCTGATAAAATTCTACACGCATTGCTTTTATCCAGACAGGGCAAATTCCAGCACATATTACCCATTCTTATATGt                                                                                                   <  2:4261504/150‑1 (MQ=255)
                                                ctgcGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTc                                                                                                  <  1:2223684/150‑1 (MQ=255)
                                                   cgcAGACGTAATTATGAGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTGTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTcgcg                                                                                               <  2:3587532/150‑1 (MQ=255)
                                                     cAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGcgac                                                                                             <  1:3615971/150‑1 (MQ=255)
                                                      aGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGcgacg                                                                                            >  1:3339056/1‑150 (MQ=255)
                                                        aCGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCGATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACgag                                                                                          <  2:2550137/150‑1 (MQ=255)
                                                          gTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGAt                                                                                        >  1:3759933/1‑150 (MQ=255)
                                                            aataatGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCGATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATaa                                                                                      >  2:218871/1‑150 (MQ=255)
                                                             ataatGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAAc                                                                                     <  1:3064086/150‑1 (MQ=255)
                                                             ataatGCGCTCTTATTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAAc                                                                                     <  2:4860254/150‑1 (MQ=255)
                                                              taatGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAAcg                                                                                    <  2:502722/150‑1 (MQ=255)
                                                                atGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAAcgcg                                                                                  >  2:3759735/1‑150 (MQ=255)
                                                                atGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAAcgcg                                                                                  >  1:4860652/1‑150 (MQ=255)
                                                                atGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAAcgcg                                                                                  >  1:4386139/1‑150 (MQ=255)
                                                                 tGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGt                                                                                 <  2:1964227/150‑1 (MQ=255)
                                                                  gcgcTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGtt                                                                                <  2:4487697/150‑1 (MQ=255)
                                                                   cgcTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTa                                                                               >  2:2050360/1‑150 (MQ=255)
                                                                      tcttctTTAAACAACCCTTCCGACAGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCTCGCATTGCGATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAaca                                                                            <  2:3541154/150‑1 (MQ=255)
                                                                      tcttctGTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAaca                                                                            <  1:3261988/150‑1 (MQ=255)
                                                                       cttcttTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAacac                                                                           >  1:4686806/1‑150 (MQ=255)
                                                                        ttcttTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCGATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACt                                                                          <  2:187413/150‑1 (MQ=255)
                                                                         tcttTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCGATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACtt                                                                         <  1:1068742/150‑1 (MQ=255)
                                                                         tcttTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCGATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACtt                                                                         >  2:1647057/1‑150 (MQ=255)
                                                                         tcttTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCGATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACtt                                                                         >  1:1042212/1‑150 (MQ=255)
                                                                            ttAAACAACCCTTCCGCCCGTGCTGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACACAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTg                                                                      >  2:2918482/1‑150 (MQ=255)
                                                                            ttAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTGTAGGTCGCGACGAGATAACGCGTTAACACTTCTg                                                                      >  2:5170340/1‑150 (MQ=255)
                                                                            ttAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTg                                                                      >  2:281205/1‑150 (MQ=255)
                                                                            ttAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCAAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTATg                                                                      >  1:4992710/1‑150 (MQ=255)
                                                                               aacaacCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTACTGCaa                                                                   >  2:642790/1‑150 (MQ=255)
                                                                                  aacCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCGATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAAt                                                                >  2:1091790/1‑150 (MQ=255)
                                                                                   acCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAAtt                                                               <  1:185054/150‑1 (MQ=255)
                                                                                      cTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAg                                                            <  1:1951091/150‑1 (MQ=255)
                                                                                         ccGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCGATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGat                                                         <  1:689269/150‑1 (MQ=255)
                                                                                           gCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGataa                                                       >  1:4240562/1‑150 (MQ=255)
                                                                                           gCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCGATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGagaa                                                       >  2:1729882/1‑150 (MQ=255)
                                                                                                  cGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCACATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACgt                                                >  2:724597/1‑150 (MQ=255)
                                                                                                   ggTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAATCTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACgtg                                               <  1:2431077/150‑1 (MQ=255)
                                                                                                   ggTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACgtg                                               >  1:3878926/1‑150 (MQ=255)
                                                                                                   ggTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCGATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACgtg                                               >  1:2619222/1‑150 (MQ=255)
                                                                                                   ggTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCGATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACgtg                                               >  1:45506/1‑150 (MQ=255)
                                                                                                     tttCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTa                                             <  1:4769947/150‑1 (MQ=255)
                                                                                                       tCCAGATCGTTGGTTAACAGCTGATAAAATTCTCCACTCATTGCGATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAAc                                           <  2:3098430/150‑1 (MQ=255)
                                                                                                         cAGATCGTGGGTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACtt                                         <  2:4733868/150‑1 (MQ=255)
                                                                                                             tCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCGATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATg                                     <  2:1573415/150‑1 (MQ=255)
                                                                                                              cGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCGATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGt                                    >  2:782299/1‑150 (MQ=255)
                                                                                                                   gTTAACTGCTGATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGtctctc                               <  1:2733765/150‑1 (MQ=255)
                                                                                                                        ctgctgATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAacac                          <  1:3725106/150‑1 (MQ=255)
                                                                                                                        ctgctgATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAacac                          >  1:1327054/1‑150 (MQ=255)
                                                                                                                        ctgctgATAAAATTCTCCACGCATTGCGATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAacac                          <  1:1208724/150‑1 (MQ=255)
                                                                                                                          gctgATAAAATTCTCCACGCATTGCGATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCt                        <  1:2263668/150‑1 (MQ=255)
                                                                                                                            tgATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTAc                      >  1:2014422/1‑150 (MQ=255)
                                                                                                                            tgATAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTAc                      <  2:3914376/150‑1 (MQ=255)
                                                                                                                              aTAAAATTCTCCTCGCATTGCTATTCTCCAGAAAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAg                    <  2:5200566/150‑1 (MQ=255)
                                                                                                                              aTAAAATTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAg                    >  2:3313900/1‑150 (MQ=255)
                                                                                                                                aaaaTTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGAt                  >  2:2946252/1‑150 (MQ=255)
                                                                                                                                 aaaTTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGAta                 >  2:4653007/1‑150 (MQ=255)
                                                                                                                                  aaTTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGAtag                <  2:3215555/150‑1 (MQ=255)
                                                                                                                                  aaTTCTCCACGCATTGCGATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGAtag                <  1:1679602/150‑1 (MQ=255)
                                                                                                                                   aTTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGAtagt               >  1:2961247/1‑150 (MQ=255)
                                                                                                                                   aTTCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGAtagt               >  2:800596/1‑150 (MQ=255)
                                                                                                                                    ttCTCCACGCATTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGAtagta              >  1:2509197/1‑150 (MQ=255)
                                                                                                                                     tctcCACGCATTGCGATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGAtagtag             <  1:2326398/150‑1 (MQ=255)
                                                                                                                                       tcCACGCATTGCGATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGaa           >  1:2482936/1‑150 (MQ=255)
                                                                                                                                              aTTGCTATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGt    <  1:4357492/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ACAAAAGTTAATGACGTGGCTTCCATCAGCCACAGTGATATCTGCTTGCTGCGCAGACGTAATAATGCGCTCTTCTTTAAACAACCCTTCCGCCCGTGCGGTTTCCAGATCGTTGGTTAACTGCTGATAAAATTCTCCACGCATTGCGATTCTCCAGACAGGGCAAATTCCAGCACATATTACCCAAACTTATAGGTCGCGACGAGATAACGCGTTAACACTTCTGCAAAATTCAGGATAAATAACGTGTAACTTGATGTCTCTCAACACCTACAGATAGTAGAAAAACTGTAC  >  CP009273/3785768‑3786061

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