Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,879,449 G→T 91.4% 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,4490GT91.4% 799.9 / 64.4 267S377* (TCA→TAA) ygcBR‑loop helicase‑annealase Cas3 needed for Cascade anti‑viral activity
Reads supporting (aligned to +/- strand):  ref base G (10/13);  new base T (117/127);  total (127/140)
Fisher's exact test for biased strand distribution p-value = 8.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

Only 100 of 267 total aligned reads displayed.

GATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATAAGTTTCCAAGC  >  CP009273/2879303‑2879597
                                                                                                                                                  |                                                                                                                                                    
  tCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTg                                                                                                                                                 <  1:2662945/150‑1 (MQ=255)
  tCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCGGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTg                                                                                                                                                 >  1:2839187/1‑150 (MQ=255)
   cAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGc                                                                                                                                                <  2:3106066/150‑1 (MQ=255)
        gTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGa                                                                                                                                           >  1:4131374/1‑150 (MQ=255)
           cAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTACTCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCa                                                                                                                                        >  2:2828976/1‑150 (MQ=255)
           cAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCa                                                                                                                                        <  1:679396/150‑1 (MQ=255)
            aaaCGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCaa                                                                                                                                       >  2:2192909/1‑150 (MQ=255)
               cGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGaa                                                                                                                                    <  2:2044782/150‑1 (MQ=255)
               cGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGaa                                                                                                                                    >  1:725389/1‑150 (MQ=255)
                  cGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATaa                                                                                                                                 <  1:2760561/150‑1 (MQ=255)
                   gATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAg                                                                                                                                >  2:4698139/1‑150 (MQ=255)
                            aaGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGa                                                                                                                       <  2:450380/150‑1 (MQ=255)
                             aGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGAt                                                                                                                      <  1:147173/150‑1 (MQ=255)
                             aGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGAt                                                                                                                      >  1:3196414/1‑150 (MQ=255)
                              gAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTCATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATg                                                                                                                     >  1:3799559/1‑150 (MQ=255)
                              gAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATg                                                                                                                     <  2:530466/150‑1 (MQ=255)
                                aaCACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGaa                                                                                                                   >  1:363091/1‑150 (MQ=255)
                                 acacTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGaaa                                                                                                                  <  2:2051860/150‑1 (MQ=255)
                                  cacTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGaaaa                                                                                                                 <  1:2383592/150‑1 (MQ=255)
                                    cTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATa                                                                                                               >  1:44841/1‑150 (MQ=255)
                                      ttCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAg                                                                                                             >  2:3266522/1‑150 (MQ=255)
                                        cTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTg                                                                                                           <  2:3680015/150‑1 (MQ=255)
                                         ttATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTgg                                                                                                          >  1:3853350/1‑150 (MQ=255)
                                         ttATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTgg                                                                                                          <  1:1325250/150‑1 (MQ=255)
                                         ttATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTgg                                                                                                          <  2:605187/150‑1 (MQ=255)
                                         ttATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGAGTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTgg                                                                                                          <  1:3735146/150‑1 (MQ=255)
                                          tATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGgc                                                                                                         >  2:3026273/1‑150 (MQ=255)
                                            ttGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGgctc                                                                                                       >  2:268247/1‑150 (MQ=255)
                                            ttGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGgctc                                                                                                       <  1:1627609/150‑1 (MQ=255)
                                             tGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGgctcg                                                                                                      >  2:3977610/1‑150 (MQ=255)
                                              gCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGgctcgc                                                                                                     >  1:737601/1‑150 (MQ=255)
                                               cTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTACTCGCGCGTGATTTTCTTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGgctcgcg                                                                                                    >  2:634360/1‑150 (MQ=255)
                                               cTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGgctcgcg                                                                                                    >  2:262858/1‑150 (MQ=255)
                                               cTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGgctcgcg                                                                                                    >  2:3647179/1‑150 (MQ=255)
                                               cTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGgctcgcg                                                                                                    >  2:1235617/1‑150 (MQ=255)
                                                  tgtgACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGgctcgcgctc                                                                                                 >  1:1271121/1‑150 (MQ=255)
                                                  tgtgACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGgctcgcgctc                                                                                                 >  1:1131259/1‑150 (MQ=255)
                                                    tgACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGgctcgcgctcgc                                                                                               >  2:827016/1‑150 (MQ=255)
                                                     gACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGcgctcgct                                                                                              <  2:3501083/150‑1 (MQ=255)
                                                      aCAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCtt                                                                                             <  1:4257460/150‑1 (MQ=255)
                                                      aCAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCtt                                                                                             <  1:2998203/150‑1 (MQ=255)
                                                           cACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCAtt                                                                                        <  2:4174981/150‑1 (MQ=255)
                                                             cTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATtct                                                                                      <  2:868215/150‑1 (MQ=255)
                                                             cTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATtct                                                                                      >  1:3967221/1‑150 (MQ=255)
                                                              tGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATtctc                                                                                     >  2:4592897/1‑150 (MQ=255)
                                                              tGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATtctc                                                                                     >  2:3934569/1‑150 (MQ=255)
                                                               gACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCg                                                                                    <  2:2132765/150‑1 (MQ=255)
                                                                 caacaCTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTa                                                                                  <  1:1169069/150‑1 (MQ=255)
                                                                  aacaCTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTACTCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTaa                                                                                 >  2:3985059/1‑150 (MQ=255)
                                                                       tGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAagcata                                                                            >  2:4773590/1‑150 (MQ=255)
                                                                         aaCCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAagcatagc                                                                          >  1:3998269/1‑150 (MQ=255)
                                                                           cccACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAagcatagcat                                                                        >  1:2535534/1‑150 (MQ=255)
                                                                           cccACGCTTCTTCTTGCCCCGGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAagcatagcat                                                                        >  1:1572424/1‑150 (MQ=255)
                                                                            ccACGCTTCTTCTTGCCCCTGTTCAGTCATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCAtt                                                                       >  2:3464174/1‑150 (MQ=255)
                                                                                gCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTcgcg                                                                   >  2:4573473/1‑150 (MQ=255)
                                                                                  ttcttcttGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGt                                                                 <  2:1984709/150‑1 (MQ=255)
                                                                                      tcttGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCt                                                             >  1:1960690/1‑150 (MQ=255)
                                                                                         tGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTtgt                                                          <  1:2930324/150‑1 (MQ=255)
                                                                                             ccTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTg                                                      >  1:3100103/1‑150 (MQ=255)
                                                                                              cTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTgg                                                     >  1:2686779/1‑150 (MQ=255)
                                                                                              cTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTgg                                                     <  1:4788368/150‑1 (MQ=255)
                                                                                                gTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGa                                                   >  1:4652950/1‑150 (MQ=255)
                                                                                                gTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGa                                                   <  1:4380971/150‑1 (MQ=255)
                                                                                                gTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGa                                                   <  2:1918363/150‑1 (MQ=255)
                                                                                                  tCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAgg                                                 <  1:3583034/150‑1 (MQ=255)
                                                                                                   cAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAggg                                                <  2:3527831/150‑1 (MQ=255)
                                                                                                    aGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGc                                               <  2:1491606/150‑1 (MQ=255)
                                                                                                     gTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCa                                              <  2:2763393/150‑1 (MQ=255)
                                                                                                      tAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCaa                                             >  1:4702206/1‑150 (MQ=255)
                                                                                                      tAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCaa                                             <  1:2192909/150‑1 (MQ=255)
                                                                                                       aaTCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCaaa                                            <  1:634360/150‑1 (MQ=255)
                                                                                                        aTCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCaaaa                                           <  2:1050028/150‑1 (MQ=255)
                                                                                                        aTCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCaaaa                                           <  2:1409593/150‑1 (MQ=255)
                                                                                                         tCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCaaaaa                                          <  2:1048290/150‑1 (MQ=255)
                                                                                                            cgcgTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATaa                                       >  2:324204/1‑150 (MQ=255)
                                                                                                             gcgGGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAac                                      <  1:2756779/150‑1 (MQ=255)
                                                                                                               gTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAacac                                    >  1:987799/1‑150 (MQ=255)
                                                                                                                tGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACt                                   <  1:2120855/150‑1 (MQ=255)
                                                                                                                 gATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTa                                  <  1:3395429/150‑1 (MQ=255)
                                                                                                                  aTTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTAt                                 <  2:3749723/150‑1 (MQ=255)
                                                                                                                   ttttATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATc                                <  2:3607919/150‑1 (MQ=255)
                                                                                                                   ttttATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATc                                >  1:524635/1‑150 (MQ=255)
                                                                                                                    tttATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATcc                               <  2:1727601/150‑1 (MQ=255)
                                                                                                                    tttATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATcc                               <  2:4131374/150‑1 (MQ=255)
                                                                                                                     ttATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCg                              >  2:1252115/1‑150 (MQ=255)
                                                                                                                        ttgattgaAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCaa                           <  2:11016/150‑1 (MQ=255)
                                                                                                                        ttgattgaAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCaa                           <  2:3574789/150‑1 (MQ=255)
                                                                                                                            ttgaAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATttg                       >  2:1136179/1‑150 (MQ=255)
                                                                                                                             tgaAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATttgt                      >  1:4345213/1‑150 (MQ=255)
                                                                                                                              gaAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATttgtt                     <  2:2618445/150‑1 (MQ=255)
                                                                                                                                aaGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGa                   >  2:2516961/1‑150 (MQ=255)
                                                                                                                                 agagGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGAt                  >  1:1976963/1‑150 (MQ=255)
                                                                                                                                  gagGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATc                 >  1:1219790/1‑150 (MQ=255)
                                                                                                                                      tggtTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATa             <  2:2350415/150‑1 (MQ=255)
                                                                                                                                      tggtTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATa             <  2:3197103/150‑1 (MQ=255)
                                                                                                                                        gtTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATAAg           >  2:1058637/1‑150 (MQ=255)
                                                                                                                                         ttAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATAAGt          <  2:725526/150‑1 (MQ=255)
                                                                                                                                          tAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATAAGtt         >  1:3273750/1‑150 (MQ=255)
                                                                                                                                             aCCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATAAGTTTcc      >  1:1747651/1‑150 (MQ=255)
                                                                                                                                                 tGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATAAGTTTCCAAGc  <  1:3465874/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATAAGTTTCCAAGC  >  CP009273/2879303‑2879597

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: