Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,879,449 G→T 100% S377* (TCA→TAA)  ygcB ← R‑loop helicase‑annealase Cas3 needed for Cascade anti‑viral activity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,879,4490GT100.0% 981.0 / ‑6.7 298S377* (TCA→TAA) ygcBR‑loop helicase‑annealase Cas3 needed for Cascade anti‑viral activity
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (148/149);  total (149/149)
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 298 total aligned reads displayed.

TGATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATAAGTTTCCA  >  CP009273/2879302‑2879594
                                                                                                                                                   |                                                                                                                                                 
tGATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTaa                                                                                                                                                 >  1:485940/1‑150 (MQ=255)
 gATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAAt                                                                                                                                                >  2:228492/1‑150 (MQ=255)
 gATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAAt                                                                                                                                                >  1:4705783/1‑150 (MQ=255)
 gATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAAt                                                                                                                                                <  2:1599333/150‑1 (MQ=255)
  aTCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTCATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAAtt                                                                                                                                               >  2:2149367/1‑150 (MQ=255)
  aTCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAAtt                                                                                                                                               <  2:1085757/150‑1 (MQ=255)
  aTCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAAtt                                                                                                                                               <  1:763392/150‑1 (MQ=255)
    cAATCGTGCAACCGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTTATTGc                                                                                                                                             <  1:3305329/150‑1 (MQ=255)
     aaTCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGcc                                                                                                                                            >  1:1238388/1‑150 (MQ=255)
     aaTCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGcc                                                                                                                                            >  1:2601009/1‑150 (MQ=255)
      aTCGTGTAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCa                                                                                                                                           <  1:5084657/150‑1 (MQ=255)
      aTCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCa                                                                                                                                           <  1:5314405/150‑1 (MQ=255)
      aTCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCa                                                                                                                                           >  1:3663939/1‑150 (MQ=255)
        cGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATg                                                                                                                                         <  1:2931053/150‑1 (MQ=255)
              aaCGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAg                                                                                                                                   >  2:3018941/1‑150 (MQ=255)
               aCGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGa                                                                                                                                  >  1:2436576/1‑150 (MQ=255)
                 gCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAAt                                                                                                                                <  2:3112936/150‑1 (MQ=255)
                  ccGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAAACGTTAATTGCCATGAGAAAGAATa                                                                                                                               >  1:1505892/1‑150 (MQ=255)
                     aTTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGa                                                                                                                            >  1:3838080/1‑150 (MQ=255)
                     aTTTGCCCAAGAAACACTTTCTTATTGCTTTGGGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTCATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGa                                                                                                                            >  1:968987/1‑150 (MQ=255)
                      tttGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGAt                                                                                                                           <  2:763003/150‑1 (MQ=255)
                       ttGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGAtt                                                                                                                          <  2:62527/150‑1 (MQ=255)
                            cAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGACTAAGATTTgggg                                                                                                                     >  1:4542340/1‑150 (MQ=255)
                              aGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGAt                                                                                                                   >  2:4557068/1‑150 (MQ=255)
                                aaaCACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGa                                                                                                                 <  2:3752962/150‑1 (MQ=255)
                                    acTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAAt                                                                                                             >  2:3440895/1‑150 (MQ=255)
                                    acTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAAt                                                                                                             <  2:478683/150‑1 (MQ=255)
                                        tCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGt                                                                                                         >  1:1302923/1‑150 (MQ=255)
                                        tCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGt                                                                                                         <  2:1498950/150‑1 (MQ=255)
                                           tATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGgc                                                                                                      <  2:4138547/150‑1 (MQ=255)
                                            aTTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGgct                                                                                                     >  1:934773/1‑150 (MQ=255)
                                            aTTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCGGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGgct                                                                                                     >  2:362644/1‑150 (MQ=255)
                                             ttGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGgctc                                                                                                    <  2:2667449/150‑1 (MQ=255)
                                             ttGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGgctc                                                                                                    <  1:3612030/150‑1 (MQ=255)
                                               gCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGgctcgc                                                                                                  <  2:375912/150‑1 (MQ=255)
                                                cTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCCAGAATAAGATTTGGGGATGAAAATAAGTGgctcgcg                                                                                                 >  2:3453683/1‑150 (MQ=255)
                                                cTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATNAGTGgctcgcg                                                                                                 <  1:2761700/150‑1 (MQ=255)
                                                cTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATCAGTGgctcgcg                                                                                                 >  1:4497610/1‑150 (MQ=255)
                                                     tgACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGgctcgcgctcgc                                                                                            >  1:1853517/1‑150 (MQ=255)
                                                      gACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGcgctcgct                                                                                           <  1:68087/150‑1 (MQ=255)
                                                            cACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCAtt                                                                                     <  2:4368225/150‑1 (MQ=255)
                                                              cTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCAGtct                                                                                   >  1:601738/1‑150 (MQ=255)
                                                                     cacTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAagc                                                                            <  2:3329614/150‑1 (MQ=255)
                                                                      acTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAagca                                                                           <  1:4718063/150‑1 (MQ=255)
                                                                      acTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAagca                                                                           >  1:4212737/1‑150 (MQ=255)
                                                                       cTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAagcat                                                                          <  2:1741389/150‑1 (MQ=255)
                                                                         gAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAagcatag                                                                        <  2:3579288/150‑1 (MQ=255)
                                                                          aaCCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAagcatagc                                                                       >  2:4488059/1‑150 (MQ=255)
                                                                          aaCCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAagcatagc                                                                       <  1:1911019/150‑1 (MQ=255)
                                                                                     cttcttGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAg                                                            >  1:3433020/1‑150 (MQ=255)
                                                                                      ttcttGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGc                                                           >  1:2201695/1‑150 (MQ=255)
                                                                                       tcttGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCt                                                          <  1:1172364/150‑1 (MQ=255)
                                                                                        cttGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCtt                                                         <  2:277208/150‑1 (MQ=255)
                                                                                           gCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTtgtg                                                      >  1:1964106/1‑150 (MQ=255)
                                                                                           gCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTtgtg                                                      <  1:3388486/150‑1 (MQ=255)
                                                                                                tGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTggg                                                 <  1:3783866/150‑1 (MQ=255)
                                                                                                 gTTCAGTACTCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGa                                                >  2:1826398/1‑150 (MQ=255)
                                                                                                       tAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCaa                                          >  2:3181204/1‑150 (MQ=255)
                                                                                                       tAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCaa                                          >  2:2015067/1‑150 (MQ=255)
                                                                                                         aTCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCaaaa                                        >  2:4290271/1‑150 (MQ=255)
                                                                                                           cgcgcgTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAAt                                      <  2:799868/150‑1 (MQ=255)
                                                                                                           cgcgcgTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAAt                                      <  2:393660/150‑1 (MQ=255)
                                                                                                           cgcgcgTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAAt                                      <  1:2250838/150‑1 (MQ=255)
                                                                                                             cgcgTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATaa                                    <  2:1439151/150‑1 (MQ=255)
                                                                                                              gcgTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAac                                   >  1:4376873/1‑150 (MQ=255)
                                                                                                              gcgTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAac                                   <  2:3309481/150‑1 (MQ=255)
                                                                                                               cgTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAaca                                  <  1:2323735/150‑1 (MQ=255)
                                                                                                                 tGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACt                                >  2:2376627/1‑150 (MQ=255)
                                                                                                                 tGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACt                                <  1:2882905/150‑1 (MQ=255)
                                                                                                                 tGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACt                                >  1:5070172/1‑150 (MQ=255)
                                                                                                                  gATTTTATTGATTGAAAGATTTGGTTAAACCGTTTATTGCCTTGTGCAAGAATAAGTTTTGGGTATGAACATAAGTGGCTCGATCTCGCTTCCATTCTCGTAAGCATAGCTTTCGCGGTAGCTTGTGTTGTGAGGTCAAAAATAacccga                               <  2:504877/150‑5 (MQ=255)
                                                                                                                  gATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTa                               >  1:102593/1‑150 (MQ=255)
                                                                                                                  gATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTa                               >  2:372455/1‑150 (MQ=255)
                                                                                                                  gATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTa                               <  1:843442/150‑1 (MQ=255)
                                                                                                                   aTTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTAt                              >  2:5116773/1‑150 (MQ=255)
                                                                                                                      ttATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCg                           >  1:4261152/1‑150 (MQ=255)
                                                                                                                           gattgaAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAAtt                      <  2:5243013/150‑1 (MQ=255)
                                                                                                                            attgaAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAAttt                     <  1:1088359/150‑1 (MQ=255)
                                                                                                                            attgaAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAAttt                     <  2:5305124/150‑1 (MQ=255)
                                                                                                                             ttgaAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATttg                    <  1:2097934/150‑1 (MQ=255)
                                                                                                                              tgaAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATttgt                   >  2:4998134/1‑150 (MQ=255)
                                                                                                                              tgaAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATttgt                   <  2:1171643/150‑1 (MQ=255)
                                                                                                                                aaaGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATttgttg                 <  1:3552267/150‑1 (MQ=255)
                                                                                                                                 aaGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGa                >  1:3357777/1‑150 (MQ=255)
                                                                                                                                  agagGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGAt               >  2:2398969/1‑150 (MQ=255)
                                                                                                                                   gagGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATc              <  2:2898985/150‑1 (MQ=255)
                                                                                                                                   gagGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATc              <  1:3893338/150‑1 (MQ=255)
                                                                                                                                   gagGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATGTGTTGATc              >  1:3162611/1‑150 (MQ=255)
                                                                                                                                    agGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCa             <  1:2484092/150‑1 (MQ=255)
                                                                                                                                    agGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCa             >  2:3517622/1‑150 (MQ=255)
                                                                                                                                      gtggtTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAAt           >  2:4649918/1‑150 (MQ=255)
                                                                                                                                      gtggtTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAAt           <  2:4874855/150‑1 (MQ=255)
                                                                                                                                       tggtTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATa          >  1:3952682/1‑150 (MQ=255)
                                                                                                                                         gtTAAACCTTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATAAg        <  1:4604354/150‑1 (MQ=255)
                                                                                                                                         gtTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATAAg        >  1:1737443/1‑150 (MQ=255)
                                                                                                                                          ttAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATAAGt       <  1:5118313/150‑1 (MQ=255)
                                                                                                                                          ttAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATAAGt       <  2:1710902/150‑1 (MQ=255)
                                                                                                                                            aaaCCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATAAGttt     <  2:3663939/150‑1 (MQ=255)
                                                                                                                                             aaCCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATAAGTTTc    >  2:2553087/1‑150 (MQ=255)
                                                                                                                                               ccGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATAAGTTTCCa  <  1:4672219/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TGATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATAAGTTTCCA  >  CP009273/2879302‑2879594

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: