Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,204,284 C→A 47.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,2840CA47.2% 96.4 / 239.0 172D98Y (GAT→TAT) tfaEe14 prophage; putative tail fiber assembly protein
Reads supporting (aligned to +/- strand):  ref base C (51/39);  new base A (37/44);  total (89/83)
Fisher's exact test for biased strand distribution p-value = 1.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 174 total aligned reads displayed.

ATTGCAGCATCAATCAGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGG  >  CP009273/1204136‑1204433
                                                                                                                                                    |                                                                                                                                                     
aTTGCAGCATCAATCAGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCa                                                                                                                                                      >  2:3255442/1‑150 (MQ=255)
  tGCAGCATCAATCAGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCAcc                                                                                                                                                    <  2:2215287/150‑1 (MQ=255)
  tGCAGCATCAATCAGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCAcc                                                                                                                                                    <  2:3866694/150‑1 (MQ=255)
  tGCAGCATCAATCAGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAAcc                                                                                                                                                    >  2:2269818/1‑150 (MQ=255)
       caTCAATCAGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAg                                                                                                                                               <  1:1116744/150‑1 (MQ=255)
       caTCAATCAGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAg                                                                                                                                               <  1:3373628/150‑1 (MQ=255)
         tcaatcaGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCt                                                                                                                                             <  2:507505/150‑1 (MQ=255)
           aatcaGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCttt                                                                                                                                           <  2:4675632/150‑1 (MQ=255)
           aatcaGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCttt                                                                                                                                           >  2:1156665/1‑150 (MQ=255)
           aatcaGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCttt                                                                                                                                           >  1:2422628/1‑150 (MQ=255)
                gTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGAttt                                                                                                                                      >  1:890537/1‑150 (MQ=255)
                gTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGAttt                                                                                                                                      >  1:3531740/1‑150 (MQ=255)
                     tGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCtttt                                                                                                                                 >  1:1701049/1‑150 (MQ=255)
                         gcTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGAtt                                                                                                                             >  2:3120999/1‑150 (MQ=255)
                           tgtgCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATtct                                                                                                                           >  2:1070631/1‑150 (MQ=255)
                            gtgCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATtctc                                                                                                                          >  1:1954283/1‑150 (MQ=255)
                            gtgCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATtctc                                                                                                                          >  1:762694/1‑150 (MQ=255)
                               cTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAg                                                                                                                       >  2:18930/1‑150 (MQ=255)
                                 tCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGtt                                                                                                                     <  1:2927464/150‑1 (MQ=255)
                                  cTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGttt                                                                                                                    >  2:3054799/1‑150 (MQ=255)
                                   tGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGtttt                                                                                                                   <  2:926292/150‑1 (MQ=255)
                                     ccGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTgg                                                                                                                 <  2:750020/150‑1 (MQ=255)
                                      cgcgTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGt                                                                                                                <  2:3348718/150‑1 (MQ=255)
                                       gcgTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTg                                                                                                               >  2:3430291/1‑150 (MQ=255)
                                         gTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCt                                                                                                             >  2:3868477/1‑150 (MQ=255)
                                         gTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGAGAGTGGTTGTATGTTCAGGGTAATAACCCGGAGCTTTGATTTCTGTTGATTCTCCAGTTTTGGTGCt                                                                                                             >  2:3868731/1‑150 (MQ=255)
                                          tCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCta                                                                                                            >  1:2336689/1‑150 (MQ=255)
                                            tACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGAGAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATGTCTTTTGATTCTCCAGTTTTGGTGCtata                                                                                                          >  2:1484579/1‑150 (MQ=255)
                                              cTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGa                                                                                                        <  2:4266846/150‑1 (MQ=255)
                                                 cggcgCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCg                                                                                                     >  2:4495384/1‑150 (MQ=255)
                                                 cggcgCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCg                                                                                                     <  2:2473174/150‑1 (MQ=255)
                                                  ggcgCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATCGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGt                                                                                                    >  1:4483943/1‑150 (MQ=255)
                                                  ggcgCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGt                                                                                                    <  2:613208/150‑1 (MQ=255)
                                                  ggcgCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGt                                                                                                    >  1:4678572/1‑150 (MQ=255)
                                                  ggcgCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGt                                                                                                    >  1:4678573/1‑150 (MQ=255)
                                                  ggcgCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGt                                                                                                    <  2:4431501/150‑1 (MQ=255)
                                                     gcTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTc                                                                                                 >  1:1460512/1‑150 (MQ=255)
                                                     gcTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTc                                                                                                 >  2:1386057/1‑150 (MQ=255)
                                                     gcTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTc                                                                                                 >  1:1460347/1‑150 (MQ=255)
                                                      cTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCa                                                                                                <  1:951283/150‑1 (MQ=255)
                                                       tATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCAc                                                                                               >  2:3797761/1‑150 (MQ=255)
                                                       tATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCAc                                                                                               <  1:2491192/150‑1 (MQ=255)
                                                       tATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCAc                                                                                               <  2:4157574/150‑1 (MQ=255)
                                                        aTGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCAcc                                                                                              <  1:1572733/150‑1 (MQ=255)
                                                        aTGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCAcc                                                                                              >  1:3337561/1‑150 (MQ=255)
                                                         tGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACcg                                                                                             <  2:1228871/150‑1 (MQ=255)
                                                         tGCTGTGCTTCAGTATCGGTCACCCATGTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACcg                                                                                             <  2:1228812/150‑1 (MQ=255)
                                                           cTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACcgcg                                                                                           >  1:963076/1‑150 (MQ=255)
                                                               gCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATgg                                                                                       >  2:1725580/1‑150 (MQ=255)
                                                                 ttCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTc                                                                                     >  1:3901345/1‑150 (MQ=255)
                                                                 ttCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTc                                                                                     >  1:4041836/1‑150 (MQ=255)
                                                                   cAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTg                                                                                   <  1:3588906/150‑1 (MQ=255)
                                                                   cAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTg                                                                                   <  1:2148651/150‑1 (MQ=255)
                                                                    aGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTgg                                                                                  >  2:636218/1‑150 (MQ=255)
                                                                         cGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACat                                                                             <  2:3632432/150‑1 (MQ=255)
                                                                           gTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACatat                                                                           <  1:1013556/150‑1 (MQ=255)
                                                                            tCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATAtt                                                                          <  2:2818354/150‑1 (MQ=255)
                                                                              aCCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTcc                                                                        >  1:5119432/1‑150 (MQ=255)
                                                                                ccATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCa                                                                      <  2:1811138/150‑1 (MQ=255)
                                                                                ccATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCa                                                                      <  2:1811182/150‑1 (MQ=255)
                                                                                 cATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCAt                                                                     <  1:4917645/150‑1 (MQ=255)
                                                                                  aTTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATg                                                                    >  2:4158822/1‑150 (MQ=255)
                                                                                    ttCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAg                                                                  >  1:3426946/1‑150 (MQ=255)
                                                                                       tcACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTa                                                               <  2:4155054/150‑1 (MQ=255)
                                                                                       tcACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTa                                                               <  2:5026282/150‑1 (MQ=255)
                                                                                            aTCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATc                                                          <  1:4340659/150‑1 (MQ=255)
                                                                                              cccATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTg                                                        <  1:1830791/150‑1 (MQ=255)
                                                                                              cccATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTg                                                        <  1:4384243/150‑1 (MQ=255)
                                                                                               ccATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGc                                                       <  1:4751293/150‑1 (MQ=255)
                                                                                                 aTTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAg                                                     <  1:347004/150‑1 (MQ=255)
                                                                                                    tATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAAc                                                  <  2:2193433/150‑1 (MQ=255)
                                                                                                      tCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACgg                                                <  1:5061377/150‑1 (MQ=255)
                                                                                                       cGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGc                                               >  2:2052093/1‑150 (MQ=255)
                                                                                                         tATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAg                                             <  2:1433458/150‑1 (MQ=255)
                                                                                                          aTGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGa                                            <  1:5099295/150‑1 (MQ=255)
                                                                                                              agaTAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGc                                        >  1:791559/1‑150 (MQ=255)
                                                                                                               gaTAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCa                                       >  2:3755491/1‑150 (MQ=255)
                                                                                                                aTAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCAt                                      >  1:2176065/1‑150 (MQ=255)
                                                                                                                 tAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGAAGGTCTGGCACATATTCCCATGAGTTAAAAACTGCAGAACGGCAGATTGCATa                                     >  1:2680792/1‑150 (MQ=255)
                                                                                                                  aaaGGGGCGATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATaa                                    >  2:1042045/1‑150 (MQ=255)
                                                                                                                   aaGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAAc                                   >  2:88758/1‑150 (MQ=255)
                                                                                                                      gggCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAg                                >  2:3118222/1‑150 (MQ=255)
                                                                                                                          gATAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCtt                            >  1:171445/1‑150 (MQ=255)
                                                                                                                           aTAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAAACAGCCTTa                           >  1:4778569/1‑150 (MQ=255)
                                                                                                                           aTAGTGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGAGTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTa                           <  2:351708/150‑1 (MQ=255)
                                                                                                                             aGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTAtg                         >  1:3658218/1‑150 (MQ=255)
                                                                                                                             aGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACAGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTAtg                         >  2:4311163/1‑150 (MQ=255)
                                                                                                                               tGGTTGTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTAtgtg                       >  1:4143568/1‑150 (MQ=255)
                                                                                                                                 gTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTa                     <  2:3619798/150‑1 (MQ=255)
                                                                                                                                  ttGTATTTTCAGGGTAATTACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTAc                    >  1:4292291/1‑150 (MQ=255)
                                                                                                                                   tGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTAcc                   <  1:158892/150‑1 (MQ=255)
                                                                                                                                    gTATTTTCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCa                  <  2:4844677/150‑1 (MQ=255)
                                                                                                                                     tATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAg                 <  2:890537/150‑1 (MQ=255)
                                                                                                                                          tCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCa            <  1:382040/150‑1 (MQ=255)
                                                                                                                                          tCAGGGTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGAGTGCATAACCAGCCTTATGTGTACCAGGGGCa            >  2:2690233/1‑150 (MQ=255)
                                                                                                                                             gggTAATAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCt         <  2:3040639/150‑1 (MQ=255)
                                                                                                                                                taataaCCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTaaa      >  2:1629436/1‑150 (MQ=255)
                                                                                                                                                tAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTaaa      <  2:3723890/150‑1 (MQ=255)
                                                                                                                                                 aaTAACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAAc     <  1:2620663/150‑1 (MQ=255)
                                                                                                                                                   tCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGGTGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAg   <  1:1325553/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ATTGCAGCATCAATCAGTGACTGGCGCTGTGCTTCTGCCGCGTCTACTGCGGCGCTATGCTGTGCTTCAGTATCGGTCACCCATTTCTCACCATCCCATTTATCGTATGGAGATAAAGGGGCGATAGTGGTTGTATTTTCAGGGTAATCACCCGGAGCTTTGATTTCTTTTGATTCTCCAGTTTTGGTGCTATAGACCGTTTCACCGCGATGGTCTGGCACATATTCCCATGAGTTAAAATCTGCAGAACGGCAGATTGCATAACCAGCCTTATGTGTACCAGGGGCATCTAAACAGG  >  CP009273/1204136‑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: