Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 566,710 G→T 46.9% 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,7100GT46.9% 177.3 / inf 266G121V (GGC→GTC) ybcLinactive polymorphonuclear leukocyte migration suppressor; DLP12 prophage; secreted protein, UPF0098 family
Reads supporting (aligned to +/- strand):  ref base G (65/76);  new base T (71/54);  total (136/130)
Fisher's exact test for biased strand distribution p-value = 8.68e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 270 total aligned reads displayed.

ATATGATCCAGATGCACCTACAGGCAGTGGTTGGTGGCATTGGACTGTTGTTAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATAAAA  >  CP009273/566561‑566859
                                                                                                                                                     |                                                                                                                                                     
atatGATCCAGATGCACCTACAGGCAGTGGTTGGTGGCATTGGACTGTTGTTAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATGTTgg                                                                                                                                                       >  2:2383622/1‑150 (MQ=255)
    gATCCAGATGCACCTACAGGCAGTGGTTGGTGGCATTGGACTGTTGTTAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTAt                                                                                                                                                   <  1:363979/150‑1 (MQ=255)
     aTCCAGATGCACCTACAGGCAGTGGTTGGTGGCATTGGACTGTTGTTAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGGCGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATg                                                                                                                                                  <  1:3247173/150‑1 (MQ=255)
     aTCCAGATGCACCTACAGGCAGTGGTTGGTGGCATTGGACTGTTGTTAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCTAGGCCGAAATGATTTTGGCTATg                                                                                                                                                  <  2:564842/150‑1 (MQ=255)
     aTCCAGATGCACCTACAGGCAGTGGTTGGTGGCATTGGACTGTTGTTAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATg                                                                                                                                                  >  2:725049/1‑150 (MQ=255)
     aTCCAGATGCACCTACAGGCAGTGGTTGGTGGCATTGGACTGTTGTTAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATg                                                                                                                                                  >  2:724910/1‑150 (MQ=255)
     aTCCAGATGCACCTACAGGCAGGGGTTGGTGGCATTGGACTGGTGTTAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGTCTATg                                                                                                                                                  >  1:327316/1‑150 (MQ=255)
       ccAGATGCACCTACAGGCAGTGGTTGGTGGCATTGGACTGTTGTTAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCt                                                                                                                                                <  1:259824/150‑1 (MQ=255)
            tGCACCTACAGGCAGTGGTTGGTGGCATTGGACTGTTGTTAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGTCTATGCTGGGtt                                                                                                                                           >  2:1897245/1‑150 (MQ=255)
                      ggCAGTGGTTGGTGGCATTGGACTGTTGTTAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGTCTATGCTGGGTTTGGTGGCGCa                                                                                                                                 <  2:3753714/150‑1 (MQ=255)
                       gCAGTGGTTGGTGGCATTGGACTGTTGTTAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGTCTATGCTGGGTTTGGTGGCGCAt                                                                                                                                <  2:1468064/150‑1 (MQ=255)
                         aGTGGTTGGTGGCATTGGACTGTTGTTAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGTCTATGCTGGGTTTGGTGGCGCATGt                                                                                                                              <  1:1279153/150‑1 (MQ=255)
                         aGTGGTTGGTGGCATTGGACTGTTGTTAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGTCTATGCTGGGTTTGGTGGCGCATGt                                                                                                                              <  1:2823650/150‑1 (MQ=255)
                           tggttggtGGCATTGGACTGTTGTTAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGTCTATGCTGGGTTTGGTGGCGCATGtcc                                                                                                                            <  2:3023297/150‑1 (MQ=255)
                            ggttggtGGCATTGGACTGTTGTTAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGtcct                                                                                                                           >  1:1134843/1‑150 (MQ=255)
                              ttggtGGCATTGGACTGTTGTTAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGtcctcc                                                                                                                         <  2:2103215/150‑1 (MQ=255)
                              ttggtGGCATTGGACTGTTGTTAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGtcctcc                                                                                                                         <  1:686536/150‑1 (MQ=255)
                               tggtggCATTGGACTGTTGTTAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGTCTATGCTGGGTTTGGTGGCGCATGTcctcct                                                                                                                        <  2:1310554/150‑1 (MQ=255)
                               tggtggCATTGGACTGTTGTTAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTcctcct                                                                                                                        >  2:640415/1‑150 (MQ=255)
                                  tggCATTGGACTGTTGTTAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGTCTATGCTGGGTTTGGTGGCGCATGTCCTCCTaaa                                                                                                                     >  1:2447142/1‑150 (MQ=255)
                                  tggCATTGGACTGTTGTTAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTaaa                                                                                                                     <  1:4481084/150‑1 (MQ=255)
                                    gCATTGGACTGTTGTTAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAgg                                                                                                                   <  2:3894354/150‑1 (MQ=255)
                                     cATTGGACTGTTGTTAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGga                                                                                                                  <  2:1436053/150‑1 (MQ=255)
                                     aaatGGACTGATGTTAATATTCCAGCAACAGTAACATATTTGCCCGTAGATGCAGGGAGACGTGATGGAACAAAACAGCCGAATGGTGATGATCAAGGCCGAAATAATATTGTCTATGCTGTGTTTGGTGGCGCAGGTCCTCCTAAAGga                                                                                                                  >  2:4624814/4‑150 (MQ=255)
                                         ggACTGTTGTTAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATa                                                                                                              <  2:2321556/150‑1 (MQ=255)
                                          gACTGTTGTTAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATaa                                                                                                             <  2:4594906/150‑1 (MQ=255)
                                              gttgttAATTTTCCAGCTACAGTTACATTTTGGCCCGTTGATGCAGGGAGACGTGTTGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAATTGATTTTTGCTATGGTGGGTTTGGTGGCGCATGTCCTCCTAATGGAGATAAACca                                                                                                         <  1:4101043/150‑1 (MQ=255)
                                              gttgttAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACca                                                                                                         >  1:1614441/1‑150 (MQ=255)
                                               ttgttAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGTCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACcac                                                                                                        >  2:3186145/1‑150 (MQ=255)
                                                 gttAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGTCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCAcat                                                                                                      <  1:3776268/150‑1 (MQ=255)
                                                  ttAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGTCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCAcatc                                                                                                     <  2:3478469/150‑1 (MQ=255)
                                                  ttAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGTCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCAcatc                                                                                                     <  1:1322608/150‑1 (MQ=255)
                                                   tAATATTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGTCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCAcatca                                                                                                    >  2:3110071/1‑150 (MQ=255)
                                                      tatTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTa                                                                                                 <  1:930822/150‑1 (MQ=255)
                                                       atTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGTCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTAc                                                                                                <  2:3571850/150‑1 (MQ=255)
                                                       atTCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTAc                                                                                                >  2:871328/1‑150 (MQ=255)
                                                        ttCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTAcc                                                                                               <  2:1551421/150‑1 (MQ=255)
                                                         tCCAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCa                                                                                              >  2:4167761/1‑150 (MQ=255)
                                                           cAGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGTCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGt                                                                                            >  2:3392146/1‑150 (MQ=255)
                                                            aGCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGTCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGtt                                                                                           >  2:1320413/1‑150 (MQ=255)
                                                             gCAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGTCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGttt                                                                                          <  2:1084767/150‑1 (MQ=255)
                                                              cAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGTCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTa                                                                                         >  2:2946200/1‑150 (MQ=255)
                                                              cAACAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGTCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTa                                                                                         >  2:1653497/1‑150 (MQ=255)
                                                               aaCAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGTCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTaa                                                                                        >  1:457738/1‑150 (MQ=255)
                                                                 cAGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAg                                                                                      >  1:2053584/1‑150 (MQ=255)
                                                                  aGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGTCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGt                                                                                     >  1:2194354/1‑150 (MQ=255)
                                                                  aGTAACATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGt                                                                                     >  1:1220721/1‑150 (MQ=255)
                                                                     aaCATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATg                                                                                  >  2:1065410/1‑150 (MQ=255)
                                                                       cATATTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGTCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATggg                                                                                <  2:4301340/150‑1 (MQ=255)
                                                                        atatTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGTCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGc                                                                               <  1:1853202/150‑1 (MQ=255)
                                                                         tatTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGTCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGct                                                                              >  1:135271/1‑150 (MQ=255)
                                                                         tatTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGct                                                                              >  2:3788724/1‑150 (MQ=255)
                                                                          atTTGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGctc                                                                             <  2:2762271/150‑1 (MQ=255)
                                                                           tttGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGctct                                                                            <  2:3464373/150‑1 (MQ=255)
                                                                            ttGCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGTCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTa                                                                           >  2:4328742/1‑150 (MQ=255)
                                                                              gCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGTCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAATTATGGGCTCTaaa                                                                         >  2:2469279/1‑150 (MQ=255)
                                                                              gCCCGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGTCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTaaa                                                                         >  1:3482438/1‑150 (MQ=255)
                                                                                 cGTTGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGTCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAAc                                                                      >  1:3354369/1‑150 (MQ=255)
                                                                                    tGATGCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGa                                                                   >  1:3582325/1‑150 (MQ=255)
                                                                                        gCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGTCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAg                                                               <  2:3660771/150‑1 (MQ=255)
                                                                                        gCAGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGTCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAg                                                               <  2:3660702/150‑1 (MQ=255)
                                                                                          aGGGAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGAt                                                             >  2:2435736/1‑150 (MQ=255)
                                                                                           gggAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGTCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGAtt                                                            >  2:2667313/1‑150 (MQ=255)
                                                                                           gggAGACGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGAtt                                                            <  1:1203064/150‑1 (MQ=255)
                                                                                            ggAGACGTGATGGAACAAAACTGCCTACTGGTGCAGTTCAAGGCCGAAATGATTTTCGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTAACAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTc                                                           >  1:811839/1‑150 (MQ=255)
                                                                                             gagaCGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTcc                                                          <  1:3863715/150‑1 (MQ=255)
                                                                                             gagaCGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTcc                                                          <  1:1025196/150‑1 (MQ=255)
                                                                                              agaCGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCt                                                         >  2:2909248/1‑150 (MQ=255)
                                                                                                aCGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGTCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGt                                                       >  1:1594897/1‑150 (MQ=255)
                                                                                                aCGTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGt                                                       >  2:800168/1‑150 (MQ=255)
                                                                                                  gTGATGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGTCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAg                                                     >  1:2514909/1‑150 (MQ=255)
                                                                                                      tGGAACAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTc                                                 >  1:2992432/1‑150 (MQ=255)
                                                                                                         aaCAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGTCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTaa                                              <  1:3573858/150‑1 (MQ=255)
                                                                                                         aaCAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTaa                                              >  1:464075/1‑150 (MQ=255)
                                                                                                           cAAAACTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACt                                            <  1:1449820/150‑1 (MQ=255)
                                                                                                            aaaaCTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGTCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTc                                           >  2:2382736/1‑150 (MQ=255)
                                                                                                               aCTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGTCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAg                                        <  1:4213119/150‑1 (MQ=255)
                                                                                                                cTGCCGACTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGc                                       >  2:2728830/1‑150 (MQ=255)
                                                                                                                     gACTGGTGCTGTTCCAGGCCGAAATGATTTTGTCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGc                                  <  2:4031253/150‑1 (MQ=255)
                                                                                                                      aCTGGTGCTGTTCAAGGCCGAAATGATTTTGTCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCg                                 <  1:335931/150‑1 (MQ=255)
                                                                                                                      aCTGGTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCg                                 >  2:4391242/1‑150 (MQ=255)
                                                                                                                        tGGTGCTGTTCAAGGCCGAAATGATTTTGTCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGtt                               >  2:1203147/1‑150 (MQ=255)
                                                                                                                         ggTGCTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTa                              >  1:449541/1‑150 (MQ=255)
                                                                                                                             cTGTTCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTg                          <  1:3389962/150‑1 (MQ=255)
                                                                                                                               gTTCAAGGCCGAAATGATTTTGTCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGt                        >  2:4185618/1‑150 (MQ=255)
                                                                                                                                ttCAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGtt                       >  1:3557440/1‑150 (MQ=255)
                                                                                                                                 tCAAGGCCGAAATGATTTTGTCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGAATCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTta                      >  1:371142/1‑150 (MQ=255)
                                                                                                                                  cAAGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTtat                     <  2:3645445/150‑1 (MQ=255)
                                                                                                                                   aaGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTtata                    >  1:30295/1‑150 (MQ=255)
                                                                                                                                    aGGCCGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTatat                   <  2:3315897/150‑1 (MQ=255)
                                                                                                                                       ccGAAATGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCt                <  2:2542155/150‑1 (MQ=255)
                                                                                                                                         gAAATGATTTTGTCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTa              >  2:1099898/1‑150 (MQ=255)
                                                                                                                                          aaaTGATTTTGTCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTaa             <  2:2072249/150‑1 (MQ=255)
                                                                                                                                          aaaTGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTaa             <  2:1174122/150‑1 (MQ=255)
                                                                                                                                            aTGATTTTGTCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATg           >  2:4478474/1‑150 (MQ=255)
                                                                                                                                            aTGATTTTGTCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATg           >  2:2152180/1‑150 (MQ=255)
                                                                                                                                             tGATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGc          <  1:2882221/150‑1 (MQ=255)
                                                                                                                                              gATTTTGGCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCt         <  2:1523373/150‑1 (MQ=255)
                                                                                                                                                  ttGTCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCtaata     <  2:2169491/150‑1 (MQ=255)
                                                                                                                                                     gCTATGCTGGGTTTGGTGGCGCATGTCCTCCTAAAGGAGATAAACCACATCATTACCAGTTTAAAGTATGGGCTCTAAAAACTGAAAAGATTCCTGTAGATTCTAACTCCAGCGGAGCGTTAGTTGGTTATATGCTTAATGCTAATaaaa  <  1:703516/150‑1 (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: