Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,879,449 G→T 86.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,4490GT86.6% 741.4 / 114.2 283S377* (TCA→TAA) ygcBR‑loop helicase‑annealase Cas3 needed for Cascade anti‑viral activity
Reads supporting (aligned to +/- strand):  ref base G (16/22);  new base T (131/114);  total (147/136)
Fisher's exact test for biased strand distribution p-value = 2.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 283 total aligned reads displayed.

TGATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATAAGTTTCCAAG  >  CP009273/2879302‑2879596
                                                                                                                                                   |                                                                                                                                                   
tGATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTaa                                                                                                                                                   >  2:155377/1‑150 (MQ=255)
 gATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTCATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAAt                                                                                                                                                  >  1:939147/1‑150 (MQ=255)
  aTCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAAtt                                                                                                                                                 >  1:1081435/1‑150 (MQ=255)
  aTCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAAtt                                                                                                                                                 <  1:1757616/150‑1 (MQ=255)
  aTCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTAGTGATTGAAAGAGGTGGTTAAACCGTTAAtt                                                                                                                                                 >  2:3927119/1‑150 (MQ=255)
   tCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTg                                                                                                                                                >  1:718675/1‑150 (MQ=255)
    cAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGc                                                                                                                                               <  2:2650621/150‑1 (MQ=255)
      aTCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCa                                                                                                                                             <  2:2904246/150‑1 (MQ=255)
       tCGTGCAAACGCCGATTTGCCCAAGAAACACTTTATTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCAt                                                                                                                                            <  2:4356040/150‑1 (MQ=255)
            cAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCa                                                                                                                                       >  2:1647790/1‑150 (MQ=255)
                cGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGaa                                                                                                                                   >  2:11533/1‑150 (MQ=255)
                  ccGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATa                                                                                                                                 >  2:3811657/1‑150 (MQ=255)
                   cGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATaa                                                                                                                                >  1:4017879/1‑150 (MQ=255)
                   cGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATaa                                                                                                                                >  1:1231139/1‑150 (MQ=255)
                    gATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAg                                                                                                                               <  2:2060794/150‑1 (MQ=255)
                    gATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAg                                                                                                                               >  2:2992434/1‑150 (MQ=255)
                    gATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAg                                                                                                                               <  1:2589403/150‑1 (MQ=255)
                         gCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTg                                                                                                                          >  1:4233749/1‑150 (MQ=255)
                          cccAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTCCTCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTgg                                                                                                                         >  1:1966110/1‑150 (MQ=255)
                          cccAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTgg                                                                                                                         >  1:1105446/1‑150 (MQ=255)
                            cAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTgggg                                                                                                                       >  2:4072650/1‑150 (MQ=255)
                            cAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTgggg                                                                                                                       <  1:3647374/150‑1 (MQ=255)
                             aaGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGa                                                                                                                      >  1:1574539/1‑150 (MQ=255)
                               gAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATg                                                                                                                    <  1:3120447/150‑1 (MQ=255)
                                 aaCACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGaa                                                                                                                  >  1:407099/1‑150 (MQ=255)
                                  acacTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGaaa                                                                                                                 >  1:2221422/1‑150 (MQ=255)
                                     cTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATa                                                                                                              >  1:362312/1‑150 (MQ=255)
                                     cTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATa                                                                                                              >  1:928834/1‑150 (MQ=255)
                                         cTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAATTg                                                                                                          >  1:1769483/1‑150 (MQ=255)
                                         cTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTg                                                                                                          <  2:332927/150‑1 (MQ=255)
                                            aTTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGAGTTGGGGATGAAAATAAGTGgct                                                                                                       >  2:1656873/1‑150 (MQ=255)
                                              tGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGgctcg                                                                                                     <  1:1605210/150‑1 (MQ=255)
                                               gCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGgctcgc                                                                                                    >  1:3154268/1‑150 (MQ=255)
                                                cTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGTTTTGGGGATGAAAATAAGTGgctcgcg                                                                                                   <  1:1208206/150‑1 (MQ=255)
                                                cTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGgctcgcg                                                                                                   <  1:3399539/150‑1 (MQ=255)
                                                cTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGgctcgcg                                                                                                   >  1:1724025/1‑150 (MQ=255)
                                                         aaCCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTcc                                                                                          <  1:1943491/150‑1 (MQ=255)
                                                         aaCCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTcc                                                                                          <  2:2728066/150‑1 (MQ=255)
                                                         aaCCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTcc                                                                                          <  1:2252314/150‑1 (MQ=255)
                                                         aaCCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTcc                                                                                          <  1:3205664/150‑1 (MQ=255)
                                                                gACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCg                                                                                   >  2:1072909/1‑150 (MQ=255)
                                                                gACAACACTGAACCCACGCTTCTTCTTGCCCCTGGGCAGTCCGCGCGCGTGATGTTAGTGATTGCAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCGCGCGTcccatcgcg                                                                                   >  1:33074/1‑143 (MQ=255)
                                                                   aacaCTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTACTCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCACGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTaa                                                                                >  2:248207/1‑150 (MQ=255)
                                                                      acTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAagca                                                                             >  1:85603/1‑150 (MQ=255)
                                                                        tGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAagcata                                                                           <  2:2254969/150‑1 (MQ=255)
                                                                          aaCCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAagcatagc                                                                         >  1:2812777/1‑150 (MQ=255)
                                                                          aaCCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAagcatagc                                                                         <  2:1448714/150‑1 (MQ=255)
                                                                            cccACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAagcatagcat                                                                       <  2:2321665/150‑1 (MQ=255)
                                                                             ccACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCAtt                                                                      <  2:2100807/150‑1 (MQ=255)
                                                                                 gCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTcgcg                                                                  >  2:1106371/1‑150 (MQ=255)
                                                                                   ttcttcttGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGt                                                                >  1:1499181/1‑150 (MQ=255)
                                                                                   ttcttcttGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGt                                                                <  2:3794427/150‑1 (MQ=255)
                                                                                   ttcttcttGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGt                                                                >  2:3934220/1‑150 (MQ=255)
                                                                                    tcttcttGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTa                                                               <  1:3366712/150‑1 (MQ=255)
                                                                                    tcttcttGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCTCGCTCGCTTCCATTCTCGTAAGCATAGCATTAGCGGTa                                                               >  2:3770608/1‑150 (MQ=255)
                                                                                          tGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTtgt                                                         >  2:1289086/1‑150 (MQ=255)
                                                                                           gCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTtgtg                                                        <  2:3048838/150‑1 (MQ=255)
                                                                                            ccccTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATATCATTCGCGGTAGCTTgtgt                                                       >  2:4285085/1‑150 (MQ=255)
                                                                                            ccccTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTgtgt                                                       <  1:539939/150‑1 (MQ=255)
                                                                                             cccTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGtt                                                      <  2:2554399/150‑1 (MQ=255)
                                                                                             cccTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGtt                                                      <  2:4099769/150‑1 (MQ=255)
                                                                                                  ttCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAg                                                 <  2:562846/150‑1 (MQ=255)
                                                                                                   tCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAgg                                                <  2:4043040/150‑1 (MQ=255)
                                                                                                   tCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAgg                                                <  2:858954/150‑1 (MQ=255)
                                                                                                    cAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAggg                                               >  2:1347542/1‑150 (MQ=255)
                                                                                                       tAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAACATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCaa                                            >  2:2385808/1‑150 (MQ=255)
                                                                                                       tAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCaa                                            <  2:1396261/150‑1 (MQ=255)
                                                                                                           cgcgcgTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAAt                                        <  2:3211285/150‑1 (MQ=255)
                                                                                                             cgcgTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTTGGAGGGCAAAAATaa                                      >  1:204345/1‑150 (MQ=255)
                                                                                                             cgcgTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATaa                                      <  2:1966110/150‑1 (MQ=255)
                                                                                                              gcgTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAac                                     <  2:3278597/150‑1 (MQ=255)
                                                                                                               cgTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAaca                                    <  1:4173455/150‑1 (MQ=255)
                                                                                                               cgTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAaca                                    >  2:3309640/1‑150 (MQ=255)
                                                                                                                gTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAacac                                   <  2:936080/150‑1 (MQ=255)
                                                                                                                 tGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACt                                  >  2:1076816/1‑150 (MQ=255)
                                                                                                                 tGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACt                                  >  2:576306/1‑150 (MQ=255)
                                                                                                                  gATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTa                                 <  2:4233749/150‑1 (MQ=255)
                                                                                                                     tttATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATcc                              <  2:3753286/150‑1 (MQ=255)
                                                                                                                      ttATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCg                             <  1:1388808/150‑1 (MQ=255)
                                                                                                                       tATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGc                            >  2:2002245/1‑150 (MQ=255)
                                                                                                                          tgattgaAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAAt                         >  2:1026218/1‑150 (MQ=255)
                                                                                                                           gattgaAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAAtt                        <  2:3346523/150‑1 (MQ=255)
                                                                                                                            attgaAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAAttt                       <  2:809148/150‑1 (MQ=255)
                                                                                                                            attgaAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAAttt                       >  1:1464160/1‑150 (MQ=255)
                                                                                                                              tgaAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATttgt                     >  2:4303061/1‑150 (MQ=255)
                                                                                                                              tgaAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATttgt                     >  1:4113891/1‑150 (MQ=255)
                                                                                                                                aaaGAGGTGGTTAAACCATGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATttgttg                   <  2:2354995/150‑1 (MQ=255)
                                                                                                                                   gagGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATc                >  1:3899266/1‑150 (MQ=255)
                                                                                                                                   gagGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATc                >  1:1278327/1‑150 (MQ=255)
                                                                                                                                   gagGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATc                >  2:2958355/1‑150 (MQ=255)
                                                                                                                                    agGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCa               <  2:2983245/150‑1 (MQ=255)
                                                                                                                                       tggtTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAGAAATAACACTATCCGCAATTTGTTGATCAAtt            >  2:659799/1‑149 (MQ=255)
                                                                                                                                       tggtTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATa            >  2:3102250/1‑150 (MQ=255)
                                                                                                                                            aaaCCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATAAGttt       <  1:2994450/150‑1 (MQ=255)
                                                                                                                                            aaaCCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCAGTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATAAGttt       <  2:2357341/150‑1 (MQ=255)
                                                                                                                                               ccGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATAAGTTTCCa    >  2:397708/1‑150 (MQ=255)
                                                                                                                                               ccGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATAAGTTTCCa    >  2:1854921/1‑150 (MQ=255)
                                                                                                                                               ccGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATAAGTTTCCa    >  2:3542881/1‑150 (MQ=255)
                                                                                                                                                cGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATAAGTTTCCaa   >  1:3511022/1‑150 (MQ=255)
                                                                                                                                                cGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATAAGTTTCCaa   >  1:186567/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TGATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATAAGTTTCCAAG  >  CP009273/2879302‑2879596

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: