Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,204,284 C→A 45.2% D98Y (GAT→TAT)  tfaE ← e14 prophage; putative tail fiber assembly protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,204,2840CA45.2% 52.1 / 90.7 73D98Y (GAT→TAT) tfaEe14 prophage; putative tail fiber assembly protein
Reads supporting (aligned to +/- strand):  ref base C (33/7);  new base A (15/18);  total (48/25)
Fisher's exact test for biased strand distribution p-value = 1.25e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATTGCAGCATCAATCAGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGG  >  CP009273/1204135‑1204433
                                                                                                                                                     |                                                                                                                                                     
cATTGCAGCATCAATCAGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATc                                                                                                                                                       <  2:3444364/150‑1 (MQ=255)
    gcagcaTCAATCAGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCAccc                                                                                                                                                   >  2:3939838/1‑150 (MQ=255)
    gcagcaTCAATCAGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCAccc                                                                                                                                                   >  2:3939836/1‑150 (MQ=255)
      agcaTCAATCAGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCgg                                                                                                                                                 >  2:3185997/1‑150 (MQ=255)
        caTCAATCAGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAg                                                                                                                                               <  1:213426/150‑1 (MQ=255)
         atcaatcaGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGc                                                                                                                                              >  1:3517787/1‑150 (MQ=255)
           caatcaGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCtt                                                                                                                                            <  1:724801/150‑1 (MQ=255)
             atcaGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTg                                                                                                                                          >  1:1185524/1‑150 (MQ=255)
             atcaGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTg                                                                                                                                          >  1:1185554/1‑150 (MQ=255)
                   gACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCt                                                                                                                                    >  1:2704863/1‑150 (MQ=255)
                   gACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCt                                                                                                                                    >  1:2704867/1‑150 (MQ=255)
                   gACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCt                                                                                                                                    >  1:2704967/1‑150 (MQ=255)
                   gACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCt                                                                                                                                    >  1:2704992/1‑150 (MQ=255)
                              tgCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATtctca                                                                                                                         >  2:1236578/1‑149 (MQ=255)
                                  tCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGtt                                                                                                                     >  2:1690498/1‑150 (MQ=255)
                                  tCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTGGATTTCTTTTGATTCTCCAGtt                                                                                                                     >  2:3949341/1‑150 (MQ=255)
                                   cTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGttt                                                                                                                    <  2:3755233/150‑1 (MQ=255)
                                   cTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGttt                                                                                                                    <  2:3755268/150‑1 (MQ=255)
                                     gCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTg                                                                                                                  <  2:3738396/150‑1 (MQ=255)
                                     gCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTg                                                                                                                  <  2:3738376/150‑1 (MQ=255)
                                          gTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCt                                                                                                             <  1:1265962/150‑1 (MQ=255)
                                          gTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCt                                                                                                             <  1:1265958/150‑1 (MQ=255)
                                          gTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCt                                                                                                             <  1:1265957/150‑1 (MQ=255)
                                           tCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCta                                                                                                            <  1:4082540/150‑1 (MQ=255)
                                           tCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCta                                                                                                            <  1:4080816/150‑1 (MQ=255)
                                              aCTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAg                                                                                                         >  2:3509620/1‑150 (MQ=255)
                                                 gcggcgCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGAcc                                                                                                      <  2:993691/150‑1 (MQ=255)
                                                               tgCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATg                                                                                        >  1:392940/1‑150 (MQ=255)
                                                                  ttCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGtt                                                                                     >  1:2810462/1‑149 (MQ=255)
                                                                  ttCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTc                                                                                     >  1:2810238/1‑150 (MQ=255)
                                                                  ttCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTc                                                                                     >  1:2810444/1‑150 (MQ=255)
                                                                  ttCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTc                                                                                     >  1:2810453/1‑150 (MQ=255)
                                                                     aGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTgg                                                                                  <  2:2435510/150‑1 (MQ=255)
                                                                        aTCGGTCACCCATTTCTCTCCATCCCATATATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGcac                                                                               >  1:4433652/1‑150 (MQ=255)
                                                                          cgctCAGCCATTTCTCTCCATGCCATTTATCGTTTTGAGTTATAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCTCCCGGAGCTTTGATTTTTATTGATTCTCCAGTTTTGGTGATATAGACCGTTTCACCGCTATGGTCTGGcactt                                                                             <  2:4370769/147‑3 (MQ=255)
                                                                           ggTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACata                                                                            >  1:3152490/1‑150 (MQ=255)
                                                                                     ttCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAg                                                                  >  1:3903435/1‑150 (MQ=255)
                                                                                     ttCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAg                                                                  >  1:3903434/1‑150 (MQ=255)
                                                                                     ttCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAg                                                                  >  1:3902646/1‑150 (MQ=255)
                                                                                        tcACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTa                                                               <  1:2445475/150‑1 (MQ=255)
                                                                                        tcACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTa                                                               <  1:2445521/150‑1 (MQ=255)
                                                                                           ccATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTaaaa                                                            <  1:4205553/150‑1 (MQ=255)
                                                                                            cATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAAt                                                           >  2:3919318/1‑150 (MQ=255)
                                                                                            cATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGAACGTTTCACCGCGAGGGTCTGGCACATATTCCCATGAGTTAAAAt                                                           >  2:3919317/1‑150 (MQ=255)
                                                                                            cATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTAGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCGGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAAt                                                           >  2:3921912/1‑150 (MQ=255)
                                                                                              tCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCt                                                         <  2:1180912/150‑1 (MQ=255)
                                                                                                 cATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCa                                                      >  2:3502891/1‑150 (MQ=255)
                                                                                                 cATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCa                                                      >  2:3502879/1‑150 (MQ=255)
                                                                                                 cATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCa                                                      >  2:3502632/1‑150 (MQ=255)
                                                                                                    ttATCGTATGGAGATAAAGGGCCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGaa                                                   >  1:1017738/1‑150 (MQ=255)
                                                                                                     tATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAAc                                                  >  2:404595/1‑150 (MQ=255)
                                                                                                     tATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAAc                                                  >  2:404617/1‑150 (MQ=255)
                                                                                                      aTCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACg                                                 >  2:1595712/1‑150 (MQ=255)
                                                                                                                  tAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATa                                     <  1:1114188/150‑1 (MQ=255)
                                                                                                                    aaGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAAc                                   >  1:521608/1‑150 (MQ=255)
                                                                                                                        ggCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATATACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGc                               >  2:3882243/1‑150 (MQ=255)
                                                                                                                        ggCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGc                               >  2:4481082/1‑150 (MQ=255)
                                                                                                                                 ggTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATgtgt                      >  2:4268313/1‑150 (MQ=255)
                                                                                                                                 ggTTGGATTTTCAGGGTAATAACCCGGAGATTTGATTTATTTTGATTCTCCAGTTTTGGTTCTATAGAGCGTTTCTCAGCGATGGTCTGGCACATAGTCCCATGAGTTAAAATCTGCAGAACGGACGATTGCATAACCAGCCGTAtgtgg                      >  2:4266324/1‑149 (MQ=255)
                                                                                                                                  gTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTa                     <  1:2313544/150‑1 (MQ=255)
                                                                                                                                      tATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAg                 <  2:2840149/150‑1 (MQ=255)
                                                                                                                                         tttCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAgggg              <  2:1185524/150‑1 (MQ=255)
                                                                                                                                         tttCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAgggg              <  2:1185554/150‑1 (MQ=255)
                                                                                                                                          ttCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGc             >  2:4250552/1‑150 (MQ=255)
                                                                                                                                          ttCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGc             >  2:3581505/1‑150 (MQ=255)
                                                                                                                                          ttCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGc             >  2:3442986/1‑150 (MQ=255)
                                                                                                                                          ttCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGc             >  2:3442872/1‑150 (MQ=255)
                                                                                                                                           tCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCa            <  2:1789852/150‑1 (MQ=255)
                                                                                                                                            cAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCAt           >  2:1678737/1‑150 (MQ=255)
                                                                                                                                            cAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCAt           >  2:1678617/1‑150 (MQ=255)
                                                                                                                                            cAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCAt           >  2:1678118/1‑150 (MQ=255)
                                                                                                                                                  aaTCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAAc     <  1:4406673/150‑1 (MQ=255)
                                                                                                                                                  aaTAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTGCCAGGGGCATCTAAAc     >  2:2144177/1‑150 (MQ=255)
                                                                                                                                                  aaTAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAAc     >  2:2144226/1‑150 (MQ=255)
                                                                                                                                                     cACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAgg  <  1:1856877/150‑1 (MQ=255)
                                                                                                                                                     cACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAgg  <  1:383350/150‑1 (MQ=255)
                                                                                                                                                     cACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAgg  <  1:383355/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CATTGCAGCATCAATCAGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGG  >  CP009273/1204135‑1204433

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