Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 566,710 G→T 17.2% G121V (GGC→GTC)  ybcL → inactive polymorphonuclear leukocyte migration suppressor; DLP12 prophage; secreted protein, UPF0098 family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273566,7100GT17.2% 736.6 / 93.2 236G121V (GGC→GTC) ybcLinactive polymorphonuclear leukocyte migration suppressor; DLP12 prophage; secreted protein, UPF0098 family
Reads supporting (aligned to +/- strand):  ref base G (87/108);  new base T (15/26);  total (102/134)
Fisher's exact test for biased strand distribution p-value = 3.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 240 total aligned reads displayed.

ATATGATCCAGATGCACCTACAGGCAGTGGTTGGTGGCATTGGACTGTTGTTAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAA  >  CP009273/566561‑566859
                                                                                                                                                     |                                                                                                                                                     
 tatGATCCAGATGCACCTACAGGCAGTGGTTGGTGGCATTGGACTGTTGTTAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGc                                                                                                                                                      <  1:2624710/150‑1 (MQ=255)
   tGATCCAGATGCACCTACAGGCAGTGGTTGGTGGCATTGGACTGTTGTTAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGTCTa                                                                                                                                                    >  1:354374/1‑150 (MQ=255)
     aTCCAGATGCACCTACAGGCAGTGGTTGGTGGCATTGGACTGTTGTTAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGTCTATg                                                                                                                                                  >  1:2362671/1‑150 (MQ=255)
     aTCCAGATGCACCTACAGGCAGTGGTTGGTGGCATTGGACTGTTGTTAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGTCTATg                                                                                                                                                  <  1:2089807/150‑1 (MQ=255)
         aGATGCACCTACAGGCAGTGGTTGGTGGCATTGGACTGTTGTTAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTgg                                                                                                                                              <  2:3801219/150‑1 (MQ=255)
           aTGCACCTACAGGCAGTGGTTGGTGGCATTGGACTGTTGTTAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGt                                                                                                                                            <  2:1428833/150‑1 (MQ=255)
            tGCACCTACAGGCAGTGGTTGGTGGCATTGGACTGTTGTTAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGtt                                                                                                                                           <  2:73681/150‑1 (MQ=255)
            tGCACCTACAGGCAGTGGTTGGTGGCATTGGACTGTTGTTAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGtt                                                                                                                                           >  1:2931012/1‑150 (MQ=255)
             gCACCTACAGGCAGTGGTTGGTGGCATTGGACTGTTGTTAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGTCTATGCTGGGttt                                                                                                                                          <  2:4335958/150‑1 (MQ=255)
             gCACCTACAGGCAGTGGTTGGTGGCATTGGACTGTTGTTAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGttt                                                                                                                                          <  1:791808/150‑1 (MQ=255)
             gCACCTACAGGCAGTGGTTGGTGGCATTGGACTGTTGTTAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGttt                                                                                                                                          <  2:1369353/150‑1 (MQ=255)
             gCACCTACAGGCAGTGGTTGGTGGCATTGGACTGTTGTTAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGttt                                                                                                                                          <  2:1369323/150‑1 (MQ=255)
              cACCTACAGGCAGTGGTTGGTGGCATTGGACTGTTGTTAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTtg                                                                                                                                         <  1:706850/150‑1 (MQ=255)
                    cAGGCAGTGGTTGGTGGCATTGGACTGTTGTTAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGgcg                                                                                                                                   >  1:4505011/1‑150 (MQ=255)
                         aGTGGTTGGTGGCATTGGACTGTTGTTAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGt                                                                                                                              >  1:2763951/1‑150 (MQ=255)
                         aGTGGTTGGTGGCATTGGACTGTTGTTAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGt                                                                                                                              >  1:3939930/1‑150 (MQ=255)
                              ttggtGGCATTGGACTGTTGTTAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGTCTATGCTGGGTTTGGTGGCGCATGtcctcc                                                                                                                         <  1:3721297/150‑1 (MQ=255)
                              ttggtGGCATTGGACTGTTGTTAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGtcctcc                                                                                                                         >  1:1789707/1‑150 (MQ=255)
                               tggtggCATTGGACTGTTGTTAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTcctcct                                                                                                                        <  1:4588470/150‑1 (MQ=255)
                                    gCATTGGACTGTTGTTAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAgg                                                                                                                   <  2:4233959/150‑1 (MQ=255)
                                         ggACTGTTGTTAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGTCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATa                                                                                                              <  1:4349065/150‑1 (MQ=255)
                                           aCTGTTGTTAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATaaa                                                                                                            >  1:2253053/1‑150 (MQ=255)
                                              gttgttAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACca                                                                                                         <  2:3896670/150‑1 (MQ=255)
                                              gttgttAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACca                                                                                                         <  2:1179619/150‑1 (MQ=255)
                                              gttgttAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACca                                                                                                         <  2:354551/150‑1 (MQ=255)
                                                tgttAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACcaca                                                                                                       >  1:193736/1‑150 (MQ=255)
                                                 gttAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCAcat                                                                                                      <  1:1285003/150‑1 (MQ=255)
                                                 gttAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCAcat                                                                                                      >  2:1251220/1‑150 (MQ=255)
                                                  ttAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCAGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCAcatc                                                                                                     >  2:3509875/1‑150 (MQ=255)
                                                     atatTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGTCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCAtt                                                                                                  >  1:3840607/1‑150 (MQ=255)
                                                      tatTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGTCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTa                                                                                                 <  1:1526213/150‑1 (MQ=255)
                                                      tatTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGTCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTa                                                                                                 <  2:674288/150‑1 (MQ=255)
                                                       atTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGTCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTAc                                                                                                >  2:1600657/1‑150 (MQ=255)
                                                       atTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTAc                                                                                                <  1:1437087/150‑1 (MQ=255)
                                                         tCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCa                                                                                              >  2:1850779/1‑150 (MQ=255)
                                                         tCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCa                                                                                              <  1:1188740/150‑1 (MQ=255)
                                                          ccAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAg                                                                                             <  1:3389584/150‑1 (MQ=255)
                                                           cAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGt                                                                                            <  2:2993768/150‑1 (MQ=255)
                                                             gCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGttt                                                                                          >  2:1599817/1‑150 (MQ=255)
                                                             gCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGttt                                                                                          <  2:3926438/150‑1 (MQ=255)
                                                                   gTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTa                                                                                    >  1:4293261/1‑150 (MQ=255)
                                                                    tAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTAt                                                                                   >  1:2568851/1‑150 (MQ=255)
                                                                      aCATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATgg                                                                                 <  2:312476/150‑1 (MQ=255)
                                                                      aCATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATgg                                                                                 <  1:178776/150‑1 (MQ=255)
                                                                        atatTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGc                                                                               >  2:3219031/1‑150 (MQ=255)
                                                                          atTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGctc                                                                             <  2:3364078/150‑1 (MQ=255)
                                                                          atTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGctc                                                                             >  1:3118348/1‑150 (MQ=255)
                                                                            ttGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTa                                                                           <  1:1448780/150‑1 (MQ=255)
                                                                               cccGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTaaaa                                                                        >  1:1880110/1‑150 (MQ=255)
                                                                               cccGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTaaaa                                                                        >  1:85855/1‑150 (MQ=255)
                                                                                 cGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAAc                                                                      >  2:2355007/1‑150 (MQ=255)
                                                                                    tGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGa                                                                   >  2:2089740/1‑150 (MQ=255)
                                                                                     gATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGaa                                                                  <  1:1554084/150‑1 (MQ=255)
                                                                                     gATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGaa                                                                  >  2:199784/1‑150 (MQ=255)
                                                                                     gATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGaa                                                                  <  1:698182/150‑1 (MQ=255)
                                                                                     gATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGaa                                                                  <  1:3378808/150‑1 (MQ=255)
                                                                                     gATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGaa                                                                  >  2:4613300/1‑150 (MQ=255)
                                                                                          aGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGAt                                                             >  2:505168/1‑150 (MQ=255)
                                                                                          aGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGAt                                                             >  1:2134731/1‑150 (MQ=255)
                                                                                            ggAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTc                                                           <  2:165798/150‑1 (MQ=255)
                                                                                             gagaCGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATGTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTcc                                                          >  2:856963/1‑150 (MQ=255)
                                                                                                aCGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGt                                                       <  1:1158894/150‑1 (MQ=255)
                                                                                                   tGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGTCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGa                                                    <  1:3962072/150‑1 (MQ=255)
                                                                                                       ggAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCt                                                <  2:4313583/150‑1 (MQ=255)
                                                                                                        gAACAAAACTGCCGACTGGTGCTGTTCCAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTa                                               <  2:3431768/150‑1 (MQ=255)
                                                                                                        gAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTa                                               <  1:2179322/150‑1 (MQ=255)
                                                                                                        gAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTa                                               <  2:1523601/150‑1 (MQ=255)
                                                                                                           cAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACt                                            <  2:2171985/150‑1 (MQ=255)
                                                                                                           cAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACt                                            >  1:1950981/1‑150 (MQ=255)
                                                                                                             aaaCTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTcc                                          >  1:3752121/1‑150 (MQ=255)
                                                                                                              aaCTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCa                                         <  2:439305/150‑1 (MQ=255)
                                                                                                              aaCTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCa                                         >  1:2966308/1‑150 (MQ=255)
                                                                                                                cTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGTCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGc                                       <  2:1787151/150‑1 (MQ=255)
                                                                                                                cTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGc                                       >  2:1376158/1‑150 (MQ=255)
                                                                                                                    cGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAg                                   >  1:3024366/1‑150 (MQ=255)
                                                                                                                    cGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAg                                   >  1:3604693/1‑150 (MQ=255)
                                                                                                                        tGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCggt                               >  2:4160091/1‑148 (MQ=255)
                                                                                                                           tGCTGTTCAAGGCCGAAATGATTTTGTCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGt                            >  2:1417921/1‑150 (MQ=255)
                                                                                                                           tGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGt                            >  2:603784/1‑150 (MQ=255)
                                                                                                                            gCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGtt                           >  2:3627046/1‑150 (MQ=255)
                                                                                                                              tGTTCAAGGCCGAAATGATTTTGTCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAACAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTTGTTgg                         <  2:1109901/150‑1 (MQ=255)
                                                                                                                                ttCAAGGCCGAAATGATTTTGTCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGtt                       <  1:2418512/150‑1 (MQ=255)
                                                                                                                                 tCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTta                      >  1:404441/1‑150 (MQ=255)
                                                                                                                                 tCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTta                      <  2:2433858/150‑1 (MQ=255)
                                                                                                                                  cAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTtat                     <  2:2253053/150‑1 (MQ=255)
                                                                                                                                  cAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTtat                     >  1:946740/1‑150 (MQ=255)
                                                                                                                                       ccGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCt                >  2:4033337/1‑150 (MQ=255)
                                                                                                                                        cGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCtt               <  1:1413991/150‑1 (MQ=255)
                                                                                                                                        cGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCtt               <  1:4115783/150‑1 (MQ=255)
                                                                                                                                        cGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCtt               >  2:1073813/1‑150 (MQ=255)
                                                                                                                                         gAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTa              <  1:3509768/150‑1 (MQ=255)
                                                                                                                                          aaaTGATTTTGTCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTaa             >  1:249198/1‑150 (MQ=255)
                                                                                                                                          aaaTGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTaa             >  1:1306384/1‑150 (MQ=255)
                                                                                                                                           aaTGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAAt            >  2:3802184/1‑150 (MQ=255)
                                                                                                                                             tGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGc          >  1:4387134/1‑150 (MQ=255)
                                                                                                                                              gATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCt         >  2:2800715/1‑150 (MQ=255)
                                                                                                                                              gATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCt         <  2:778889/150‑1 (MQ=255)
                                                                                                                                               aTTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCta        >  1:1115366/1‑150 (MQ=255)
                                                                                                                                                 tttGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCtaat      <  2:1126482/150‑1 (MQ=255)
                                                                                                                                                     gCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATaaaa  >  1:2507623/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATATGATCCAGATGCACCTACAGGCAGTGGTTGGTGGCATTGGACTGTTGTTAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAA  >  CP009273/566561‑566859

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: