Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,217,098 C→A 18.4% I369I (ATC→ATA ebgA → evolved beta‑D‑galactosidase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,217,0980CA18.4% 794.6 / 114.8 261I369I (ATC→ATAebgAevolved beta‑D‑galactosidase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (90/123);  new base A (27/21);  total (117/144)
Fisher's exact test for biased strand distribution p-value = 1.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 264 total aligned reads displayed.

TCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCA  >  CP009273/3216950‑3217247
                                                                                                                                                    |                                                                                                                                                     
  aaaGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATAAAc                                                                                                                                                    >  2:3050714/1‑150 (MQ=255)
     gTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTcc                                                                                                                                                 >  1:985521/1‑150 (MQ=255)
      tgagcgACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACTATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGAGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCg                                                                                                                                                <  2:3510631/147‑1 (MQ=255)
       gcgcGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGt                                                                                                                                               <  2:862495/150‑1 (MQ=255)
         gcgACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATAAACTCCGTGc                                                                                                                                             <  1:3713152/150‑1 (MQ=255)
          cgACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCg                                                                                                                                            <  2:4353385/150‑1 (MQ=255)
           gACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGt                                                                                                                                           >  2:448551/1‑150 (MQ=255)
            aCGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTa                                                                                                                                          >  1:3442939/1‑150 (MQ=255)
             cGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTAc                                                                                                                                         >  1:2326861/1‑150 (MQ=255)
               gTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCg                                                                                                                                       >  1:3598897/1‑150 (MQ=255)
                  tGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTc                                                                                                                                    >  1:1500849/1‑150 (MQ=255)
                   gTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCa                                                                                                                                   <  1:259172/150‑1 (MQ=255)
                       tGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTAc                                                                                                                               <  2:350112/150‑1 (MQ=255)
                          aTCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCg                                                                                                                            >  1:3399411/1‑150 (MQ=255)
                            cAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGaa                                                                                                                          <  2:2183693/150‑1 (MQ=255)
                             aaTAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAAc                                                                                                                         <  2:2741359/150‑1 (MQ=255)
                             aaTAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAAc                                                                                                                         >  1:3658073/1‑150 (MQ=255)
                             aaTAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAAc                                                                                                                         >  2:226992/1‑150 (MQ=255)
                                aaCCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGAt                                                                                                                      <  1:1898519/150‑1 (MQ=255)
                                    ggtATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCcgc                                                                                                                  >  2:3357518/3‑150 (MQ=255)
                                    gTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATAAACTCCGTGCGTACCGCTCACTACCCGAACGATCcgc                                                                                                                  <  1:1245921/150‑1 (MQ=255)
                                      tATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGt                                                                                                                >  2:761103/1‑150 (MQ=255)
                                        tgtgATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGttt                                                                                                              <  1:400147/150‑1 (MQ=255)
                                        tgtgATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGttt                                                                                                              <  1:1988924/150‑1 (MQ=255)
                                        tgtgATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGttt                                                                                                              <  1:1988913/150‑1 (MQ=255)
                                        tgtgATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGttt                                                                                                              >  2:2035733/1‑150 (MQ=255)
                                        tgtgATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGttt                                                                                                              <  1:639873/150‑1 (MQ=255)
                                          tgATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGttttt                                                                                                            >  2:2168801/1‑150 (MQ=255)
                                          tgATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGttttt                                                                                                            <  1:3309440/150‑1 (MQ=255)
                                           gATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTa                                                                                                           >  1:2607166/1‑150 (MQ=255)
                                            aTGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATAAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTAc                                                                                                          >  1:1575112/1‑150 (MQ=255)
                                                   aCGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACtgt                                                                                                   >  2:1357675/1‑150 (MQ=255)
                                                     ggCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACtgtgt                                                                                                 >  1:254402/1‑150 (MQ=255)
                                                     ggCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACtgtgt                                                                                                 <  1:1929921/150‑1 (MQ=255)
                                                     ggCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACtgtgt                                                                                                 >  2:4123246/1‑150 (MQ=255)
                                                      gCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATAAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACtgtgtg                                                                                                <  1:2959418/150‑1 (MQ=255)
                                                      gCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATAAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACtgtgtg                                                                                                <  1:312984/150‑1 (MQ=255)
                                                        gTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGat                                                                                              <  1:4235536/150‑1 (MQ=255)
                                                         tCAACCGTCACGACAACGATCAGCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGTTGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCTCTACCCGAACGATCCGCGTTTTTACGAACTGTGTGata                                                                                             <  1:4211918/150‑1 (MQ=255)
                                                          cAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGatat                                                                                            <  1:1085058/150‑1 (MQ=255)
                                                           aaCCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGTACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATc                                                                                           <  2:3133013/150‑1 (MQ=255)
                                                               gTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACg                                                                                       >  1:1373684/1‑150 (MQ=255)
                                                                 cACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATAAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGc                                                                                     >  2:4485259/1‑150 (MQ=255)
                                                                   cGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATAAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCt                                                                                   >  1:2367282/1‑150 (MQ=255)
                                                                     aCAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGt                                                                                 <  1:3972331/150‑1 (MQ=255)
                                                                      cAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGtt                                                                                <  2:4413746/150‑1 (MQ=255)
                                                                       aaCGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGACCGATCCGCGATTTTACGAACTGTGTGATATCTACGGCCTGttt                                                                               <  2:1221352/150‑1 (MQ=255)
                                                                       aaCGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGttt                                                                               <  2:1468936/150‑1 (MQ=255)
                                                                         cGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTtgt                                                                             <  1:4299973/150‑1 (MQ=255)
                                                                         cGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTtgt                                                                             <  1:2678162/150‑1 (MQ=255)
                                                                           atcatcGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATAAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGAGTTTTTACGAACTGTGTGATATCTAAGGCCTGTTTGTGa                                                                           >  2:2240408/1‑150 (MQ=255)
                                                                               tcGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGc                                                                       <  1:240376/150‑1 (MQ=255)
                                                                               tcGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGc                                                                       <  1:4382258/150‑1 (MQ=255)
                                                                                cGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATggag                                                                      >  2:1400127/1‑148 (MQ=255)
                                                                                 gCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCgg                                                                     >  1:749385/1‑150 (MQ=255)
                                                                                   aaaGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGaa                                                                   >  2:3939561/1‑150 (MQ=255)
                                                                                    aaGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGaaa                                                                  >  1:2165745/1‑150 (MQ=255)
                                                                                      ggCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAAcc                                                                >  1:3009754/1‑150 (MQ=255)
                                                                                      ggCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAAcc                                                                >  1:2497739/1‑150 (MQ=255)
                                                                                         cgcgCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGAc                                                             <  2:3806979/150‑1 (MQ=255)
                                                                                          gcgcCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACg                                                            <  2:338181/150‑1 (MQ=255)
                                                                                            gcCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATAAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTc                                                          <  2:2108001/150‑1 (MQ=255)
                                                                                             ccGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCg                                                         >  2:133696/1‑150 (MQ=255)
                                                                                               gTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGaa                                                       >  1:3231551/1‑150 (MQ=255)
                                                                                                  ggAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCg                                                    <  2:3405360/150‑1 (MQ=255)
                                                                                                  ggAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCg                                                    <  2:3405381/150‑1 (MQ=255)
                                                                                                  ggAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCg                                                    <  1:3325616/150‑1 (MQ=255)
                                                                                                    aaTGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCa                                                  <  1:1426948/150‑1 (MQ=255)
                                                                                                         aTCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCt                                             <  1:2843723/150‑1 (MQ=255)
                                                                                                          tCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATAAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCtt                                            <  2:972476/150‑1 (MQ=255)
                                                                                                              gTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCt                                        >  1:611524/1‑150 (MQ=255)
                                                                                                                 gagaAAGATCTCCAGTTGATGAAGCAGCACAATATAAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAAt                                     <  2:1852607/150‑1 (MQ=255)
                                                                                                                  agaAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATg                                    >  2:1417364/1‑150 (MQ=255)
                                                                                                                  agaAAGATCTCCAGTTGATGAAGCAGCACAATATAAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATg                                    >  2:1108784/1‑150 (MQ=255)
                                                                                                                      aGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCgg                                >  2:391197/1‑150 (MQ=255)
                                                                                                                          ctcCAGTTGATGAAGCAGCACAATATAAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGat                            >  2:2406925/1‑150 (MQ=255)
                                                                                                                           tcCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGata                           <  1:448541/150‑1 (MQ=255)
                                                                                                                            ccAGTTGATGAAGCAGCACAATATAAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGatat                          >  1:3257213/1‑150 (MQ=255)
                                                                                                                             cAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATAtt                         >  1:1602359/1‑150 (MQ=255)
                                                                                                                              aGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTa                        <  1:4328652/150‑1 (MQ=255)
                                                                                                                                 tgatgaAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGcc                     <  1:566441/150‑1 (MQ=255)
                                                                                                                                  gatgaAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGCTCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTTCTGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCg                    <  1:4509888/150‑1 (MQ=255)
                                                                                                                                  gatgaAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCg                    >  1:881688/1‑150 (MQ=255)
                                                                                                                                   atgaAGCAGCACAATATAAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGt                   >  1:2101440/1‑150 (MQ=255)
                                                                                                                                      aaGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTAtt                <  2:2581084/150‑1 (MQ=255)
                                                                                                                                      aaGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTAtt                <  1:1760326/150‑1 (MQ=255)
                                                                                                                                      aaGCAGCACAATATAAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTAtt                >  2:4413438/1‑150 (MQ=255)
                                                                                                                                         cagcaCAATATCAACTCCGTGCGTACCGCTCACTTCCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAATCCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTAcc             <  2:2832567/150‑1 (MQ=255)
                                                                                                                                         cagcaCAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTAcc             <  1:226999/150‑1 (MQ=255)
                                                                                                                                          agcaCAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACcg            <  2:13224/150‑1 (MQ=255)
                                                                                                                                            cacaATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACcgac          <  2:1563262/150‑1 (MQ=255)
                                                                                                                                            cacaATATAAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACcgac          <  2:356884/150‑1 (MQ=255)
                                                                                                                                             acaATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACcgacg         >  1:123479/1‑150 (MQ=255)
                                                                                                                                               aaTATAAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGAt       >  2:2651054/1‑150 (MQ=255)
                                                                                                                                                atatCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATc      >  1:3952832/1‑150 (MQ=255)
                                                                                                                                                atatCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATc      >  1:4574016/1‑150 (MQ=255)
                                                                                                                                                 tatCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATcc     <  2:633598/150‑1 (MQ=255)
                                                                                                                                                 tatCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATcc     <  2:3993936/150‑1 (MQ=255)
                                                                                                                                                    cAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCa  <  1:374763/150‑1 (MQ=255)
                                                                                                                                                    cAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCa  >  1:4196741/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TCAAAGTGCGCGACGGTCTGTTCTGGATCAATAACCGTTATGTGATGCTGCACGGCGTCAACCGTCACGACAACGATCATCGCAAAGGCCGCGCCGTTGGAATGGATCGCGTCGAGAAAGATCTCCAGTTGATGAAGCAGCACAATATCAACTCCGTGCGTACCGCTCACTACCCGAACGATCCGCGTTTTTACGAACTGTGTGATATCTACGGCCTGTTTGTGATGGCGGAAACCGACGTCGAATCGCACGGCTTTGCTAATGTCGGCGATATTAGCCGTATTACCGACGATCCGCA  >  CP009273/3216950‑3217247

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: