Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,879,449 G→T 93.6% 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,4490GT93.6% 903.4 / 51.6 296S377* (TCA→TAA) ygcBR‑loop helicase‑annealase Cas3 needed for Cascade anti‑viral activity
Reads supporting (aligned to +/- strand):  ref base G (11/8);  new base T (129/146);  total (141/155)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

Only 100 of 296 total aligned reads displayed.

GATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATAAGTTTCCAA  >  CP009273/2879303‑2879595
                                                                                                                                                  |                                                                                                                                                  
 aTCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAAtt                                                                                                                                                <  2:3206513/150‑1 (MQ=255)
  tCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTg                                                                                                                                               <  2:1683781/150‑1 (MQ=255)
       cGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATg                                                                                                                                          <  2:2040307/150‑1 (MQ=255)
       cGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATg                                                                                                                                          >  2:3413465/1‑150 (MQ=255)
       cGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATg                                                                                                                                          <  1:2591682/150‑1 (MQ=255)
       cGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATg                                                                                                                                          <  1:2651465/150‑1 (MQ=255)
          gCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTTCCATGAGc                                                                                                                                       >  1:3735962/1‑150 (MQ=255)
          gCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAAAACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGc                                                                                                                                       >  2:4633892/1‑150 (MQ=255)
           cAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCa                                                                                                                                      >  2:4335309/1‑150 (MQ=255)
              aCGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTACTCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGa                                                                                                                                   >  1:3637181/1‑150 (MQ=255)
                gCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTCATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAAt                                                                                                                                 >  1:1300346/1‑150 (MQ=255)
                gCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAAt                                                                                                                                 <  2:451316/150‑1 (MQ=255)
                 ccGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATa                                                                                                                                <  1:692153/150‑1 (MQ=255)
                  cGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATaa                                                                                                                               <  1:836344/150‑1 (MQ=255)
                    aTTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGa                                                                                                                             >  2:1943441/1‑150 (MQ=255)
                     tttGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGAt                                                                                                                            <  1:1251092/150‑1 (MQ=255)
                      ttGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGAtt                                                                                                                           <  2:3903410/150‑1 (MQ=255)
                         cccAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTgg                                                                                                                        >  2:4120672/1‑150 (MQ=255)
                          ccAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTggg                                                                                                                       >  2:2922490/1‑150 (MQ=255)
                          ccAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTggg                                                                                                                       >  2:1674627/1‑150 (MQ=255)
                          ccAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTggg                                                                                                                       <  2:598775/150‑1 (MQ=255)
                          ccAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTggg                                                                                                                       >  1:4022942/1‑150 (MQ=255)
                             aGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGAt                                                                                                                    <  2:246111/150‑1 (MQ=255)
                               aaaCACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGa                                                                                                                  <  1:715103/150‑1 (MQ=255)
                                 acacTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGaaa                                                                                                                >  1:1565117/1‑150 (MQ=255)
                                        cTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTg                                                                                                         <  2:4006841/150‑1 (MQ=255)
                                        cTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTg                                                                                                         <  2:522802/150‑1 (MQ=255)
                                        cTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTg                                                                                                         <  1:3758864/150‑1 (MQ=255)
                                              gCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGgctcgc                                                                                                   >  2:4269965/1‑150 (MQ=255)
                                              gCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGgctcgc                                                                                                   >  2:2180902/1‑150 (MQ=255)
                                              gCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGgctcgc                                                                                                   >  1:911205/1‑150 (MQ=255)
                                               cTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGgctcgcg                                                                                                  >  1:3412422/1‑150 (MQ=255)
                                                 ttGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGgctcgcgct                                                                                                <  1:4180608/150‑1 (MQ=255)
                                                 ttGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGgctcgcgct                                                                                                <  1:218338/150‑1 (MQ=255)
                                                  tgtgACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGgctcgcgctc                                                                                               <  2:2579788/150‑1 (MQ=255)
                                                  tgtgACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGgctcgcgctc                                                                                               <  2:2824624/150‑1 (MQ=255)
                                                   gtgACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGgctcgcgctcg                                                                                              >  1:2746784/1‑150 (MQ=255)
                                                      aCAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCtt                                                                                           <  2:366988/150‑1 (MQ=255)
                                                      aCAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCtt                                                                                           <  2:755433/150‑1 (MQ=255)
                                                       cAACCACAGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACAGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTc                                                                                          >  2:2380679/1‑150 (MQ=255)
                                                           cACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCAtt                                                                                      <  2:3889704/150‑1 (MQ=255)
                                                              tGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATtctc                                                                                   >  2:2942465/1‑150 (MQ=255)
                                                                 caacaCTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTa                                                                                <  2:2066779/150‑1 (MQ=255)
                                                                 caacaCTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTa                                                                                <  2:1080684/150‑1 (MQ=255)
                                                                  aacaCTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTaa                                                                               <  2:1484891/150‑1 (MQ=255)
                                                                  aacaCTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTaa                                                                               <  2:1424548/150‑1 (MQ=255)
                                                                   acacTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAag                                                                              >  1:4344207/1‑150 (MQ=255)
                                                                   acacTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAag                                                                              >  1:3386331/1‑150 (MQ=255)
                                                                      cTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAagcat                                                                           <  1:633770/150‑1 (MQ=255)
                                                                       tGAACCCACGCTTCTTCTTGCCCCTGTTCAGTACTCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAagcata                                                                          >  1:1572524/1‑150 (MQ=255)
                                                                          aCCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAagcatagca                                                                       <  1:2301446/150‑1 (MQ=255)
                                                                             cACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTc                                                                    <  2:646698/150‑1 (MQ=255)
                                                                               cGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTcgc                                                                  >  2:116269/1‑150 (MQ=255)
                                                                                gCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTcgcg                                                                 >  2:3716307/1‑150 (MQ=255)
                                                                                gCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTcgcg                                                                 <  2:2349508/150‑1 (MQ=255)
                                                                                  ttcttcttGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGt                                                               >  1:2915647/1‑150 (MQ=255)
                                                                                   tcttcttGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTa                                                              <  2:782150/150‑1 (MQ=255)
                                                                                   tcttcttGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTa                                                              <  1:1917984/150‑1 (MQ=255)
                                                                                   tcttcttGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTa                                                              >  1:3343920/1‑150 (MQ=255)
                                                                                   tcttcttGCCCCTGTTCAGTAATCGCGCGTGATTTTATCGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTa                                                              <  2:853247/150‑1 (MQ=255)
                                                                                     ttcttGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGc                                                            >  2:322875/1‑150 (MQ=255)
                                                                                      tcttGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCt                                                           >  1:3195204/1‑150 (MQ=255)
                                                                                      tcttGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCt                                                           >  2:2243965/1‑150 (MQ=255)
                                                                                           ccccTGTTCAGTAATCGAGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTgtgt                                                      >  2:3739855/1‑150 (MQ=255)
                                                                                            cccTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGtt                                                     >  2:587252/1‑150 (MQ=255)
                                                                                             ccTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTg                                                    <  1:3819705/150‑1 (MQ=255)
                                                                                                gTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGa                                                 >  1:3960476/1‑150 (MQ=255)
                                                                                                  tCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAgg                                               >  1:2458850/1‑150 (MQ=255)
                                                                                                  tCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAgg                                               >  2:2970379/1‑150 (MQ=255)
                                                                                                     gTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCa                                            >  2:836095/1‑150 (MQ=255)
                                                                                                     gTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCa                                            >  2:2607849/1‑150 (MQ=255)
                                                                                                      tAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCaa                                           <  2:1122667/150‑1 (MQ=255)
                                                                                                       aaTCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCaaa                                          >  1:1553561/1‑150 (MQ=255)
                                                                                                          cgcgcgTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAAt                                       >  2:2275455/1‑150 (MQ=255)
                                                                                                            cgcgTGATTGTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATaa                                     >  1:4135717/1‑150 (MQ=255)
                                                                                                             gcgTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAac                                    <  1:471606/150‑1 (MQ=255)
                                                                                                              cgTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAaca                                   >  2:1041235/1‑150 (MQ=255)
                                                                                                              cgTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAaca                                   <  2:3320697/150‑1 (MQ=255)
                                                                                                              cgTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAaca                                   <  1:2580240/150‑1 (MQ=255)
                                                                                                                         tgattgaAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAAt                        >  2:1541289/1‑150 (MQ=255)
                                                                                                                         tgattgaAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAAt                        <  2:1708492/150‑1 (MQ=255)
                                                                                                                         tgattgaAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAAt                        >  1:101720/1‑150 (MQ=255)
                                                                                                                          cattgaCAGAGCTGGTCAAACCCTAAATTGCCATGAGCCAGAATAAGATTAGCGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATCACACTATCCGCAAtt                       <  2:2611075/149‑1 (MQ=255)
                                                                                                                           attgaAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAAttt                      <  1:2966813/150‑1 (MQ=255)
                                                                                                                            ttgaAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATttg                     >  2:2366387/1‑150 (MQ=255)
                                                                                                                            ttgaAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATttg                     <  1:405972/150‑1 (MQ=255)
                                                                                                                               aaaGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATttgttg                  <  1:369590/150‑1 (MQ=255)
                                                                                                                               aaaGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATttgttg                  <  1:534952/150‑1 (MQ=255)
                                                                                                                                     gtggtTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAAt            <  1:4269965/150‑1 (MQ=255)
                                                                                                                                      tggtTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAACATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATa           <  1:2609864/150‑1 (MQ=255)
                                                                                                                                         ttAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATAAGt        >  2:444113/1‑150 (MQ=255)
                                                                                                                                         ttAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATAAGt        >  1:4546587/1‑150 (MQ=255)
                                                                                                                                          tAAAGCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATATCACTATCCGCATTTTGTTGATCAATAAGtt       <  2:1847289/150‑1 (MQ=255)
                                                                                                                                          tAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATAAGtt       <  2:4206512/150‑1 (MQ=255)
                                                                                                                                           aaaCCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATAAGttt      <  2:3489541/150‑1 (MQ=255)
                                                                                                                                            aaCCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATAAGTTTc     <  1:4317283/150‑1 (MQ=255)
                                                                                                                                            aaCCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATAAGTTTc     >  1:2352022/1‑150 (MQ=255)
                                                                                                                                             aCCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATAAGTTTcc    >  1:3248839/1‑150 (MQ=255)
                                                                                                                                              ccGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATAAGTTTCCa   >  1:1760924/1‑150 (MQ=255)
                                                                                                                                              ccGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATAAGTTTCCa   <  2:3339398/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATAAGTTTCCAA  >  CP009273/2879303‑2879595

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: