Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,107,730 C→A 26.9% E27* (GAG→TAG)  wcaN ← putative regulatory subunit for GalU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,107,7300CA26.9% 625.0 / 219.0 284E27* (GAG→TAG) wcaNputative regulatory subunit for GalU
Reads supporting (aligned to +/- strand):  ref base C (82/125);  new base A (28/49);  total (110/174)
Fisher's exact test for biased strand distribution p-value = 6.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 289 total aligned reads displayed.

AGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGTTCGATTTTTTTCGCCAGCATCAATTAC  >  CP009273/2107582‑2107877
                                                                                                                                                    |                                                                                                                                                   
 gTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCtt                                                                                                                                                   <  2:1013483/150‑1 (MQ=255)
  tGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCttt                                                                                                                                                  <  2:401225/150‑1 (MQ=255)
    aTTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTgg                                                                                                                                                <  1:461374/150‑1 (MQ=255)
    aTTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTgg                                                                                                                                                >  2:2212462/1‑150 (MQ=255)
       cTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTa                                                                                                                                             <  2:3046558/150‑1 (MQ=255)
       cTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTa                                                                                                                                             >  2:3190506/1‑150 (MQ=255)
         aaCTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATc                                                                                                                                           >  1:3805843/1‑150 (MQ=255)
         aaCTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATc                                                                                                                                           >  1:2100250/1‑150 (MQ=255)
          aCTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCg                                                                                                                                          <  1:4378978/150‑1 (MQ=255)
              aTAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGTGTATCGCCtt                                                                                                                                      >  2:4191087/1‑150 (MQ=255)
              aTAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTATTTGGGTCTCGCCtt                                                                                                                                      >  2:73971/1‑150 (MQ=255)
               tAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTATTTGGGTATCGCCTTa                                                                                                                                     >  2:2267024/1‑150 (MQ=255)
                  gagGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTg                                                                                                                                  <  2:270386/150‑1 (MQ=255)
                     gtgtCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTggca                                                                                                                               <  1:2880420/150‑1 (MQ=255)
                              tGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACata                                                                                                                      <  1:2239779/150‑1 (MQ=255)
                               ggTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACatat                                                                                                                     <  1:125048/150‑1 (MQ=255)
                               ggTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACatat                                                                                                                     >  1:3993991/1‑150 (MQ=255)
                               ggTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACatat                                                                                                                     >  1:3880621/1‑150 (MQ=255)
                               ggTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACatat                                                                                                                     <  2:1818201/150‑1 (MQ=255)
                                  tttCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCa                                                                                                                  <  1:1335391/150‑1 (MQ=255)
                                  tttCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTATTTGGGTATCGCCTTAGTGGCAGGCAACATATGCa                                                                                                                  <  1:3384561/150‑1 (MQ=255)
                                    tCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATc                                                                                                                <  1:3065354/150‑1 (MQ=255)
                                      gACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATccc                                                                                                              >  1:4697346/1‑150 (MQ=255)
                                      gACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATccc                                                                                                              <  1:3223623/150‑1 (MQ=255)
                                         cgcgTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCgag                                                                                                           >  1:3418405/1‑150 (MQ=255)
                                             ttCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTATTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGAccc                                                                                                       <  1:537000/150‑1 (MQ=255)
                                              tCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCg                                                                                                      >  1:2201283/1‑150 (MQ=255)
                                                 tGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTa                                                                                                   >  1:2333629/1‑150 (MQ=255)
                                                  ggACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTAc                                                                                                  >  2:2192669/1‑150 (MQ=255)
                                                  ggACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTAc                                                                                                  <  2:2423245/150‑1 (MQ=255)
                                                   gACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTATTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACa                                                                                                 <  1:1618101/150‑1 (MQ=255)
                                                     cgcgTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAgg                                                                                               <  1:4407317/150‑1 (MQ=255)
                                                      gtgTGAGTTACCATGAGTTTTTCTTTGATCACTGCAGCTACTTTCTCTTCATAAATGTTATGTATATTTGGCTTGTCGTCGTTTGTTTGCTTCTCTTTTGGTTTCGCCTTAGTGGCAGGAAAAATATGCATCCCGAGACCCGCTACAGGa                                                                                              <  2:1952083/148‑1 (MQ=255)
                                                       cgTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGaa                                                                                             <  2:3838244/150‑1 (MQ=255)
                                                        gTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAAt                                                                                            <  1:373101/150‑1 (MQ=255)
                                                          gAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACTATCTCGTCAACAATGTTCTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATaa                                                                                          <  2:2570877/150‑1 (MQ=255)
                                                           aGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAAc                                                                                         >  1:943301/1‑150 (MQ=255)
                                                             ttACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTg                                                                                       <  1:2009709/150‑1 (MQ=255)
                                                             ttACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTATTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTg                                                                                       <  2:1609196/150‑1 (MQ=255)
                                                                ccAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTATTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCtt                                                                                    <  2:390227/150‑1 (MQ=255)
                                                                   ggaggaTTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTa                                                                                 >  2:1224905/1‑150 (MQ=255)
                                                                    gaggaTTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTATTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTaa                                                                                >  1:4844257/1‑150 (MQ=255)
                                                                        aTTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTATTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTc                                                                            <  2:1943984/150‑1 (MQ=255)
                                                                         tttCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCg                                                                           <  2:889791/150‑1 (MQ=255)
                                                                         tttCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCg                                                                           >  1:1088420/1‑150 (MQ=255)
                                                                         tttCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTATTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCg                                                                           <  2:3976260/150‑1 (MQ=255)
                                                                          ttCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGt                                                                          >  2:4184460/1‑150 (MQ=255)
                                                                           tCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTATTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTc                                                                         <  1:3779330/150‑1 (MQ=255)
                                                                              ttGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCatt                                                                      <  1:2074416/150‑1 (MQ=255)
                                                                                gATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCattat                                                                    >  2:1494614/1‑150 (MQ=255)
                                                                                 aTCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCattatt                                                                   <  1:4723647/150‑1 (MQ=255)
                                                                                  tCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTAttt                                                                  <  1:4013425/150‑1 (MQ=255)
                                                                                  tCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTATTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATACCTGCTTTTAAATTCGTCATTAttt                                                                  <  2:2041549/150‑1 (MQ=255)
                                                                                    ccTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCa                                                                >  2:4848456/1‑150 (MQ=255)
                                                                                     cTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTATTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCAt                                                               >  1:2741442/1‑150 (MQ=255)
                                                                                     cTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTATTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCAt                                                               >  1:2741584/1‑150 (MQ=255)
                                                                                       gCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTATTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATcc                                                             <  2:4755844/150‑1 (MQ=255)
                                                                                       gCACCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATcc                                                             >  2:3861212/1‑150 (MQ=255)
                                                                                        cAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCa                                                            <  2:670727/150‑1 (MQ=255)
                                                                                           ccACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCt                                                         <  2:4380038/150‑1 (MQ=255)
                                                                                             acaATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTATTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGt                                                       >  2:4342409/1‑150 (MQ=255)
                                                                                              caATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTa                                                      <  2:335024/150‑1 (MQ=255)
                                                                                               aaTCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTaa                                                     <  1:3536456/150‑1 (MQ=255)
                                                                                               aaTCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTaa                                                     <  2:1546668/150‑1 (MQ=255)
                                                                                               aaTCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTATTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTaa                                                     <  1:3170017/150‑1 (MQ=255)
                                                                                                aTCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCGTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTaaa                                                    <  2:4435107/150‑1 (MQ=255)
                                                                                                aTCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTATTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTaaa                                                    <  2:4631214/150‑1 (MQ=255)
                                                                                                  ctcGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAAt                                                  <  2:3084543/150‑1 (MQ=255)
                                                                                                   tcGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATg                                                 <  1:782682/150‑1 (MQ=255)
                                                                                                   tcGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTATTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATg                                                 <  2:3402272/150‑1 (MQ=255)
                                                                                                       caacaaTGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTg                                             >  2:4735431/1‑150 (MQ=255)
                                                                                                       caacaaTGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTATTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTg                                             <  2:4710053/150‑1 (MQ=255)
                                                                                                         acaaTGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCt                                           <  2:3790714/150‑1 (MQ=255)
                                                                                                           aaTGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGa                                         >  2:2531254/1‑150 (MQ=255)
                                                                                                           aaTGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGa                                         >  2:1130718/1‑150 (MQ=255)
                                                                                                           aaTGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGa                                         <  2:1768179/150‑1 (MQ=255)
                                                                                                             tGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAAt                                       <  2:4551302/150‑1 (MQ=255)
                                                                                                               tACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATta                                     <  2:3063833/150‑1 (MQ=255)
                                                                                                                aCTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATtat                                    >  1:3433180/1‑150 (MQ=255)
                                                                                                                aCTGAATCATTGGCTTGTCGACGATTGGTAGCATCTATTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATtat                                    <  1:1792946/150‑1 (MQ=255)
                                                                                                                 cTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATtata                                   <  1:283740/150‑1 (MQ=255)
                                                                                                                 cTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATtata                                   >  2:496061/1‑150 (MQ=255)
                                                                                                                  tGAATCATTGGCTTGTCGACGATTGGTAGCATCTATTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAg                                  <  1:2445920/150‑1 (MQ=255)
                                                                                                                    aaTCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCt                                <  1:3697189/150‑1 (MQ=255)
                                                                                                                      tCATTGGCTTGTCGACGATTGGTAGCATCTATTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTg                              <  1:1683532/150‑1 (MQ=255)
                                                                                                                        aTTGGCTTGTCGACGATTGGTAGCATCTATTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGtt                            >  2:1217611/1‑150 (MQ=255)
                                                                                                                               tGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGTTCGAtttt                     >  1:3608823/1‑150 (MQ=255)
                                                                                                                                gTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGTTCGAttttt                    >  2:3780439/1‑150 (MQ=255)
                                                                                                                                gTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGTTCGAttttt                    >  2:4520359/1‑150 (MQ=255)
                                                                                                                                 tCGACGATTGGTAGCATCTATTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGTTCGAtttttt                   <  1:890537/150‑1 (MQ=255)
                                                                                                                                  cgacgaTTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGTTCGAttttttt                  >  1:2819318/1‑150 (MQ=255)
                                                                                                                                   gacgaTTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGTTCGATTTTTTTc                 >  1:4283781/1‑150 (MQ=255)
                                                                                                                                     cgaTTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGTTCGATTTTTTTCGc               >  2:2129015/1‑150 (MQ=255)
                                                                                                                                      gaTTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGTTCGATTTTTTTCGcc              >  1:1002790/1‑150 (MQ=255)
                                                                                                                                       aTTGGTAGCATCTATTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGTTCGATTTTTTTCGCCa             >  2:3397325/1‑150 (MQ=255)
                                                                                                                                             aGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGTTCGATTTTTTTCGCCAGCATCa       <  1:2549058/150‑1 (MQ=255)
                                                                                                                                               cATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGTTCGATTTTTTTCGCCAGCATCAAt     <  2:2918012/150‑1 (MQ=255)
                                                                                                                                               cATCTATTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGTTCGATTTTTTTCGCCAGCATCAAt     >  2:3061376/1‑150 (MQ=255)
                                                                                                                                                 tctcTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGTTCGATTTTTTTCGCCAGCATCAATTa   <  2:4800875/150‑1 (MQ=255)
                                                                                                                                                  ctctTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGTTCGATTTTTTTCGCCAGCATCAATTAc  >  1:1513237/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGTTCGATTTTTTTCGCCAGCATCAATTAC  >  CP009273/2107582‑2107877

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: