Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,571,233 G→A 100% S16F (TCT→TTT)  hsdS ← specificity determinant for hsdM and hsdR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,571,2330GA98.8% 887.5 / ‑0.1 258S16F (TCT→TTT) hsdSspecificity determinant for hsdM and hsdR
Reads supporting (aligned to +/- strand):  ref base G (3/0);  new base A (154/101);  total (157/101)
Fisher's exact test for biased strand distribution p-value = 2.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01
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 258 total aligned reads displayed.

TTTTAGGAACAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAG  >  CP009273/4571087‑4571377
                                                                                                                                                  |                                                                                                                                                
ttttAGGAACAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAAt                                                                                                                                                <  1:4072596/149‑1 (MQ=255)
ttttAGGAACAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAAt                                                                                                                                                <  1:710621/149‑1 (MQ=255)
 tttAGGAACAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATa                                                                                                                                               <  1:2760486/149‑1 (MQ=255)
  ttAGGAAC‑AAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACt                                                                                                                                             <  1:1763971/149‑1 (MQ=255)
    aGGAACAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATCATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTg                                                                                                                                            >  2:2107242/1‑149 (MQ=255)
    aGGAACAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTg                                                                                                                                            >  1:1315327/1‑149 (MQ=255)
    aGGAACAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTg                                                                                                                                            <  2:140902/149‑1 (MQ=255)
    aGGAACAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTg                                                                                                                                            <  1:235346/149‑1 (MQ=255)
     ggAACAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTgg                                                                                                                                           >  2:515677/1‑149 (MQ=255)
     ggAACAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTgg                                                                                                                                           <  2:4078090/149‑1 (MQ=255)
     ggAACAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTgg                                                                                                                                           <  2:1821551/149‑1 (MQ=255)
      gAACAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTggg                                                                                                                                          >  1:1136409/1‑149 (MQ=255)
       aaCAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATTCTgggg                                                                                                                                         >  2:312074/1‑149 (MQ=255)
       aaCAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATTCTgggg                                                                                                                                         >  2:898317/1‑149 (MQ=255)
       aaCAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTgggg                                                                                                                                         >  1:1965126/1‑149 (MQ=255)
        aCAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGc                                                                                                                                        <  2:3419227/149‑1 (MQ=255)
            aaaCCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATa                                                                                                                                    >  1:1892397/1‑149 (MQ=255)
             aaCCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATaa                                                                                                                                   >  2:3938840/1‑149 (MQ=255)
              aCCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAac                                                                                                                                  <  1:1667997/149‑1 (MQ=255)
                 aaGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAaccca                                                                                                                               >  1:4548919/1‑149 (MQ=255)
                  aGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAacccac                                                                                                                              <  2:2448268/149‑1 (MQ=255)
                   gTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAacccacc                                                                                                                             >  2:4365913/1‑149 (MQ=255)
                   gTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAacccacc                                                                                                                             >  1:4098950/1‑149 (MQ=255)
                     ccGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCAcccc                                                                                                                           <  2:3939533/149‑1 (MQ=255)
                       gtagtaTCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCAcaccac                                                                                                                         >  2:603210/1‑147 (MQ=255)
                        tagtaTCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACccaccccacc                                                                                                                        >  2:1595337/1‑149 (MQ=255)
                          gtaTCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCgg                                                                                                                      >  1:4271895/1‑149 (MQ=255)
                           taTCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGc                                                                                                                     >  1:2229323/1‑149 (MQ=255)
                            aTCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGTCAAATAATCATCTTTTAGATAATTTATTGCGTGCTCTTTTTTATACGTTACTCATAGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCa                                                                                                                    >  2:25031/1‑149 (MQ=255)
                                aaCTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAAttt                                                                                                                >  1:1118590/1‑149 (MQ=255)
                                  cTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTcc                                                                                                              >  1:225341/1‑149 (MQ=255)
                                  cTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTcc                                                                                                              >  1:277390/1‑149 (MQ=255)
                                      ccATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGc                                                                                                          >  1:1168167/1‑149 (MQ=255)
                                      ccATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCACCGGCAATTTCCCCGc                                                                                                          >  1:3556929/1‑149 (MQ=255)
                                            tGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCAt                                                                                                    >  2:2219327/1‑149 (MQ=255)
                                            tGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCAt                                                                                                    >  1:2953096/1‑149 (MQ=255)
                                             gAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCCtt                                                                                                   >  2:3141688/1‑149 (MQ=255)
                                               atatTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTcc                                                                                                 >  2:724789/1‑149 (MQ=255)
                                               atatTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTcc                                                                                                 >  1:461923/1‑149 (MQ=255)
                                               atatTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTcc                                                                                                 <  2:3058084/149‑1 (MQ=255)
                                               atatTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAAGTTCCCCGCACTCATTcc                                                                                                 >  1:2445848/1‑149 (MQ=255)
                                                 atTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCtt                                                                                               <  1:3504386/149‑1 (MQ=255)
                                                  ttGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTc                                                                                              >  2:3811819/1‑149 (MQ=255)
                                                          aCGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATTCTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCGTCACCccacc                                                                                      >  1:2304512/1‑149 (MQ=255)
                                                          aCGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCccacc                                                                                      >  1:2516759/1‑149 (MQ=255)
                                                           ctgATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCccacca                                                                                     <  2:3262754/147‑1 (MQ=255)
                                                             gATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCaaa                                                                                   <  1:3330174/149‑1 (MQ=255)
                                                                aaGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGc                                                                                >  1:1471926/1‑149 (MQ=255)
                                                                aaGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAAATTCCCCGCACACATTCCTTCACCCCACCAAACGc                                                                                >  1:3668188/1‑149 (MQ=255)
                                                                  gagGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGctt                                                                              >  1:1775448/1‑149 (MQ=255)
                                                                   agGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGcttc                                                                             <  2:2697120/149‑1 (MQ=255)
                                                                   agGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCGCCGGCAATTTCCCCGCACTCATTCCATCACCCCACCAAACGcgtc                                                                             >  1:3181790/1‑149 (MQ=255)
                                                                    ggCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGcttct                                                                            >  1:4412282/1‑149 (MQ=255)
                                                                    ggCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGcttct                                                                            >  1:4453239/1‑149 (MQ=255)
                                                                    ggCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGcttct                                                                            <  2:1744391/149‑1 (MQ=255)
                                                                      cAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCttcttc                                                                          <  1:2076307/149‑1 (MQ=255)
                                                                        aataatCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCa                                                                        <  2:2745315/149‑1 (MQ=255)
                                                                        aataatCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCa                                                                        >  2:4519601/1‑149 (MQ=255)
                                                                        aataatCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCa                                                                        >  1:2269927/1‑149 (MQ=255)
                                                                         ataatCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAg                                                                       <  2:2009627/149‑1 (MQ=255)
                                                                          taatCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTAACCCCACCAAACGCTTCATACAGc                                                                      >  1:2912136/1‑149 (MQ=255)
                                                                           aatCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCCTTCCTTCACCCCACCAAACGCTTCATCCAGCa                                                                     >  1:164882/1‑149 (MQ=255)
                                                                            atcatcTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCaa                                                                    <  2:4482969/149‑1 (MQ=255)
                                                                             tcatcTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAAc                                                                   <  1:2005305/149‑1 (MQ=255)
                                                                             tcatcTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAAc                                                                   <  1:2005276/149‑1 (MQ=255)
                                                                               atcTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTg                                                                 >  1:2565116/1‑149 (MQ=255)
                                                                                   tttAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGc                                                             >  2:4724975/1‑149 (MQ=255)
                                                                                   tttAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGc                                                             <  1:3964018/149‑1 (MQ=255)
                                                                                    ttAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCt                                                            <  2:1111455/149‑1 (MQ=255)
                                                                                           aTTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACAGACGCTGCAAATc                                                     >  1:2641265/1‑149 (MQ=255)
                                                                                              tATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCACATCGGc                                                  >  2:3046607/1‑149 (MQ=255)
                                                                                              tATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGc                                                  >  1:3707467/1‑149 (MQ=255)
                                                                                              tATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTATTCCAGCAACTGACGCTGCAAATCGGc                                                  >  1:935003/1‑149 (MQ=255)
                                                                                               aTTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGcc                                                 >  2:3629824/1‑149 (MQ=255)
                                                                                                ttGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCt                                                <  1:1196707/149‑1 (MQ=255)
                                                                                                  gCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCa                                              <  1:4667804/149‑1 (MQ=255)
                                                                                                  gCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCa                                              >  1:983702/1‑149 (MQ=255)
                                                                                                  gCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCa                                              >  1:634766/1‑149 (MQ=255)
                                                                                                      gCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCAtcgc                                          >  1:4498998/1‑149 (MQ=255)
                                                                                                      gCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCAtcgc                                          <  1:4200031/149‑1 (MQ=255)
                                                                                                          ttttttATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCAtcgctcgc                                      >  1:1087262/1‑149 (MQ=255)
                                                                                                               tataCGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAAATTCCCCGCACTCATTCCATACCCCCCCCAAACGCTTCTTAAAGCAACTAACGATGCAAATCGCCAAAATCGCTCGCCCCCa                                 >  2:2338190/1‑149 (MQ=255)
                                                                                                                 taCGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGt                               <  2:1133050/149‑1 (MQ=255)
                                                                                                                 taCGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAAACGGCCTCATCGCTCGCCCCCAGt                               >  1:4392776/1‑149 (MQ=255)
                                                                                                                 gaCGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGt                               <  1:3262698/148‑1 (MQ=255)
                                                                                                                        ctcctcGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCACATCGCGCGCCCCCAGTTCACGCa                        >  2:4768002/1‑149 (MQ=255)
                                                                                                                          cctcGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATc                      <  2:461923/149‑1 (MQ=255)
                                                                                                                            tcGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAg                    >  2:3773850/1‑149 (MQ=255)
                                                                                                                               gATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAgcgc                 >  2:1858253/1‑149 (MQ=255)
                                                                                                                                aTTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCa                >  2:1290085/1‑149 (MQ=255)
                                                                                                                                aTTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCATCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCa                >  2:1208122/1‑149 (MQ=255)
                                                                                                                                  tAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATc              <  1:4729960/149‑1 (MQ=255)
                                                                                                                                     agTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAg           >  2:4241662/1‑149 (MQ=255)
                                                                                                                                      gTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGt          >  1:844764/1‑149 (MQ=255)
                                                                                                                                         gtgACCGTAAATTCTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACAGACGCTGCAAATCGGCCTCATCGCTCGCCCCCCGTTCACGCATCAGCGCATCCAGTTCa       >  1:580866/1‑149 (MQ=255)
                                                                                                                                         gtgACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCa       <  2:1248772/149‑1 (MQ=255)
                                                                                                                                           gACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGa     <  1:1836142/149‑1 (MQ=255)
                                                                                                                                           gACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGa     <  2:2379122/149‑1 (MQ=255)
                                                                                                                                             ccGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACa   >  2:2814009/1‑149 (MQ=255)
                                                                                                                                              cGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAg  <  2:4596709/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
TTTTAGGAACAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAG  >  CP009273/4571087‑4571377

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑