Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,951,334 C→A 100% T352T (ACG→ACT torY ← TMAO reductase III (TorYZ), cytochrome c‑type subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,951,3340CA99.2% 907.1 / ‑2.0 277T352T (ACG→ACTtorYTMAO reductase III (TorYZ), cytochrome c‑type subunit
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base A (130/144);  total (130/147)
Fisher's exact test for biased strand distribution p-value = 2.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 277 total aligned reads displayed.

CGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGCGTA  >  CP009273/1951187‑1951480
                                                                                                                                                   |                                                                                                                                                  
      tgcCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAAt                                                                                                                                            <  1:675105/150‑1 (MQ=255)
        ccGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGt                                                                                                                                          <  2:5134123/150‑1 (MQ=255)
               caccacCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTc                                                                                                                                   <  2:4823031/150‑1 (MQ=255)
               caccacCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTc                                                                                                                                   >  2:529221/1‑150 (MQ=255)
                  caccaGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGtt                                                                                                                                <  1:4058141/150‑1 (MQ=255)
                   accaGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGttt                                                                                                                               <  2:1675698/150‑1 (MQ=255)
                    ccaGCGCTCCAGCGGCTATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGTAAATACCGAGTTAAAATGTCCAGTTCGtttt                                                                                                                              <  1:4935824/150‑1 (MQ=255)
                     caGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTc                                                                                                                             <  1:3772107/150‑1 (MQ=255)
                      aGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCg                                                                                                                            <  2:3652728/150‑1 (MQ=255)
                      aGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCg                                                                                                                            <  2:4228302/150‑1 (MQ=255)
                          cTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCa                                                                                                                        <  2:744613/150‑1 (MQ=255)
                              aGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCt                                                                                                                    <  1:4013958/150‑1 (MQ=255)
                               gCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTg                                                                                                                   <  2:4903063/150‑1 (MQ=255)
                                    aaTACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCg                                                                                                              <  2:1888212/150‑1 (MQ=255)
                                     aTACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGt                                                                                                             >  2:4768863/1‑150 (MQ=255)
                                     aTACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGt                                                                                                             >  2:4393525/1‑150 (MQ=255)
                                       aCCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGc                                                                                                           >  1:5035219/1‑150 (MQ=255)
                                        ccGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGcc                                                                                                          <  2:1330309/150‑1 (MQ=255)
                                             gtgtTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCAt                                                                                                     >  1:445421/1‑150 (MQ=255)
                                                    ataaATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTg                                                                                              >  2:3599111/1‑150 (MQ=255)
                                                     taaATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTgg                                                                                             <  1:1199745/150‑1 (MQ=255)
                                                         ttCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAAt                                                                                         <  1:2244296/150‑1 (MQ=255)
                                                          tCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATg                                                                                        >  2:1423378/1‑150 (MQ=255)
                                                           cACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATgg                                                                                       <  1:648128/150‑1 (MQ=255)
                                                           cACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATgg                                                                                       <  1:647903/150‑1 (MQ=255)
                                                              gTCTTGTTAATGTCATAACTTCCTCCCTGATCAAAGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACg                                                                                    <  1:2643959/150‑1 (MQ=255)
                                                                cTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGc                                                                                  >  2:4654582/1‑150 (MQ=255)
                                                                 ttGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGcc                                                                                 >  1:522882/1‑150 (MQ=255)
                                                                   gTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCAc                                                                               <  2:3852984/150‑1 (MQ=255)
                                                                   gTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCAc                                                                               <  2:4685302/150‑1 (MQ=255)
                                                                     tAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGc                                                                             >  1:1036610/1‑150 (MQ=255)
                                                                     tAATGTCATAACTACCTCACAGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTACGTGCCCCCATTCCTTTGGCACTGGACGGCCACGc                                                                             >  1:1921949/1‑150 (MQ=255)
                                                                       aTGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCAt                                                                           <  1:2789491/150‑1 (MQ=255)
                                                                         gTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTg                                                                         <  2:4450038/150‑1 (MQ=255)
                                                                           cATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGAc                                                                       >  1:4158785/1‑150 (MQ=255)
                                                                            aTAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACg                                                                      >  1:3233309/1‑150 (MQ=255)
                                                                             tAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACgg                                                                     >  1:4099670/1‑150 (MQ=255)
                                                                             tAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACgg                                                                     >  2:3537431/1‑150 (MQ=255)
                                                                              aaCTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGgt                                                                    >  2:2135139/1‑150 (MQ=255)
                                                                               aCTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGgtg                                                                   <  1:2707323/150‑1 (MQ=255)
                                                                                 ttCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTa                                                                 <  2:4079250/150‑1 (MQ=255)
                                                                                 ttCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTa                                                                 <  2:3052969/150‑1 (MQ=255)
                                                                                   cctccCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGATGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAAt                                                               >  2:811584/1‑150 (MQ=255)
                                                                                     tccCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATgg                                                             >  2:930195/1‑150 (MQ=255)
                                                                                      cccTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGt                                                            <  2:1444111/150‑1 (MQ=255)
                                                                                        cTGATCACCGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTcg                                                          <  1:1834177/150‑1 (MQ=255)
                                                                                        cTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTcg                                                          >  2:2989617/1‑150 (MQ=255)
                                                                                        cTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTcg                                                          >  2:390180/1‑150 (MQ=255)
                                                                                          gATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTcggc                                                        <  1:291648/150‑1 (MQ=255)
                                                                                           aTCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTcggcg                                                       <  1:2169641/150‑1 (MQ=255)
                                                                                           aTCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTcggcg                                                       <  1:2379125/150‑1 (MQ=255)
                                                                                            tCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTcggcgg                                                      >  2:170305/1‑150 (MQ=255)
                                                                                             cAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCggcggc                                                     >  2:4527748/1‑150 (MQ=255)
                                                                                                cGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAAt                                                  <  2:3448512/150‑1 (MQ=255)
                                                                                                cGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAAt                                                  >  2:2377619/1‑150 (MQ=255)
                                                                                                cGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAAt                                                  >  1:2486372/1‑150 (MQ=255)
                                                                                                 gAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATa                                                 >  2:3409045/1‑150 (MQ=255)
                                                                                                 gAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATa                                                 >  2:3507230/1‑150 (MQ=255)
                                                                                                  aGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCCATAg                                                >  1:992131/1‑150 (MQ=255)
                                                                                                  aGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAg                                                <  1:2185741/150‑1 (MQ=255)
                                                                                                  aGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAg                                                >  1:3572125/1‑150 (MQ=255)
                                                                                                  aGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACCGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAg                                                <  2:1455280/150‑1 (MQ=255)
                                                                                                   ggATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAgg                                               >  1:3408350/1‑150 (MQ=255)
                                                                                                     aTCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAgggg                                             >  2:430395/1‑150 (MQ=255)
                                                                                                     aTCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAgggg                                             <  2:3319530/150‑1 (MQ=255)
                                                                                                      tCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGc                                            >  1:1457617/1‑150 (MQ=255)
                                                                                                       cACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCa                                           <  1:3198361/150‑1 (MQ=255)
                                                                                                        aCTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCAt                                          <  2:3205818/150‑1 (MQ=255)
                                                                                                             ttCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACaa                                     >  2:1828274/1‑150 (MQ=255)
                                                                                                                 cGGTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAg                                 >  2:22918/1‑150 (MQ=255)
                                                                                                                  ggTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAgc                                >  2:3865855/1‑150 (MQ=255)
                                                                                                                  ggTAATATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAgc                                >  1:1685654/1‑150 (MQ=255)
                                                                                                                     aaTATCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCa                             <  1:5001669/150‑1 (MQ=255)
                                                                                                                       tatCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAAt                           >  1:4871153/1‑150 (MQ=255)
                                                                                                                        atCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATa                          <  2:4234462/150‑1 (MQ=255)
                                                                                                                        atCTTTGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATa                          <  2:2041371/150‑1 (MQ=255)
                                                                                                                        atCTTTGGCGTTGTACTGGAAATAACGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGTCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCACTa                          >  2:5067839/1‑150 (MQ=255)
                                                                                                                          cTTTGGCGTTGTTCTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAgg                        <  1:3176870/150‑1 (MQ=255)
                                                                                                                           tttGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGt                       >  2:108038/1‑150 (MQ=255)
                                                                                                                             tGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTAt                     >  2:1785192/1‑150 (MQ=255)
                                                                                                                             tGGCGTTGTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTAt                     <  2:1413427/150‑1 (MQ=255)
                                                                                                                                    gTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAg              <  1:2759863/150‑1 (MQ=255)
                                                                                                                                    gTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAg              <  2:2800937/150‑1 (MQ=255)
                                                                                                                                    gTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAg              <  2:2556300/150‑1 (MQ=255)
                                                                                                                                    gTACTGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAg              <  2:4030499/150‑1 (MQ=255)
                                                                                                                                       atGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCAATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGttt           >  1:3774065/2‑150 (MQ=255)
                                                                                                                                        tGGAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTCTCGTCAAGtttt          >  1:2324499/1‑150 (MQ=255)
                                                                                                                                         ggAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGttttt         >  1:3056323/1‑150 (MQ=255)
                                                                                                                                         ggAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGttttt         <  1:3909964/150‑1 (MQ=255)
                                                                                                                                         ggAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGttttt         <  2:3274002/150‑1 (MQ=255)
                                                                                                                                          gAAATACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTc        <  1:1034386/150‑1 (MQ=255)
                                                                                                                                           acagaCCGAGTTAAACTCTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCa       <  1:2945886/146‑1 (MQ=255)
                                                                                                                                           aaaTACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCa       <  2:12634/150‑1 (MQ=255)
                                                                                                                                            aaTACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAg      <  2:4502141/150‑1 (MQ=255)
                                                                                                                                            aaTACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAg      >  1:4331396/1‑150 (MQ=255)
                                                                                                                                            aaTACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAg      <  1:5096665/150‑1 (MQ=255)
                                                                                                                                             aTACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGc     >  1:4488116/1‑150 (MQ=255)
                                                                                                                                              tACCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGCg    >  2:2249923/1‑150 (MQ=255)
                                                                                                                                               aCCGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGCGt   <  1:2895092/150‑1 (MQ=255)
                                                                                                                                                ccGAGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGCGTa  >  1:4767627/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTCGTGCCCCCATTCCTTTGGCAATGGACGGCCACGCATTGACGGTGTAATGGTCGGCGGCAATAGGGGCATGACAACCAGCGCAATAGGTATCGTCAAGTTTTTCAGCGTA  >  CP009273/1951187‑1951480

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: