Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,879,449 G→T 46.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,4490GT46.6% 121.4 / inf 221S377* (TCA→TAA) ygcBR‑loop helicase‑annealase Cas3 needed for Cascade anti‑viral activity
Reads supporting (aligned to +/- strand):  ref base G (55/63);  new base T (43/60);  total (98/123)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 221 total aligned reads displayed.

CTGATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATAAGTTTCCAAGC  >  CP009273/2879301‑2879597
                                                                                                                                                    |                                                                                                                                                    
cTGATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGa                                                                                                                                                     >  1:188588/1‑150 (MQ=255)
 tGATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTaa                                                                                                                                                    >  2:3776291/1‑150 (MQ=255)
  gATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAAt                                                                                                                                                   <  2:3473463/150‑1 (MQ=255)
    tCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTg                                                                                                                                                 <  1:3396833/150‑1 (MQ=255)
       aTCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCa                                                                                                                                              >  1:1945140/1‑150 (MQ=255)
          gTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGa                                                                                                                                           >  2:2996134/1‑150 (MQ=255)
               aaCGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAg                                                                                                                                      <  1:1312468/150‑1 (MQ=255)
                  gCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAAt                                                                                                                                   <  1:1744477/150‑1 (MQ=255)
                   ccGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATa                                                                                                                                  <  1:4488296/150‑1 (MQ=255)
                    cGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATaa                                                                                                                                 >  1:4509161/1‑150 (MQ=255)
                     gATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAg                                                                                                                                <  1:1102157/150‑1 (MQ=255)
                     gATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAg                                                                                                                                >  1:1593850/1‑150 (MQ=255)
                      aTTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGa                                                                                                                               >  1:556164/1‑150 (MQ=255)
                          gCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTg                                                                                                                           <  2:3150432/150‑1 (MQ=255)
                            ccAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTggg                                                                                                                         <  2:909375/150‑1 (MQ=255)
                                   acacTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGaaa                                                                                                                  <  2:4030900/150‑1 (MQ=255)
                                   acacTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGaaa                                                                                                                  <  1:261554/150‑1 (MQ=255)
                                    cacTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGaaaa                                                                                                                 <  1:3348262/150‑1 (MQ=255)
                                    cacTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGaaaa                                                                                                                 >  2:2920767/1‑150 (MQ=255)
                                     acTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAAt                                                                                                                <  2:1445611/150‑1 (MQ=255)
                                      cTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTACTCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATa                                                                                                               >  1:2089353/1‑150 (MQ=255)
                                       tttCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATaa                                                                                                              <  1:948884/150‑1 (MQ=255)
                                         tCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGt                                                                                                            <  1:431362/150‑1 (MQ=255)
                                          cTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTg                                                                                                           >  2:169905/1‑150 (MQ=255)
                                           ttATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTgg                                                                                                          >  1:4253926/1‑150 (MQ=255)
                                              ttGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCCTGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGgctc                                                                                                       >  2:214550/1‑150 (MQ=255)
                                              ttGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGgctc                                                                                                       >  1:2676010/1‑150 (MQ=255)
                                              ttGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGgctc                                                                                                       >  1:302371/1‑150 (MQ=255)
                                               tGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGgctcg                                                                                                      <  1:4234638/150‑1 (MQ=255)
                                                   ttGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGgctcgcgct                                                                                                  <  1:2514176/150‑1 (MQ=255)
                                                   ttGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGgctcgcgct                                                                                                  <  2:4320507/150‑1 (MQ=255)
                                                     gtgACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGgctcgcgctcg                                                                                                >  1:207172/1‑150 (MQ=255)
                                                      tgACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGgctcgcgctcgc                                                                                               >  1:83574/1‑150 (MQ=255)
                                                        aCAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCtt                                                                                             >  2:344437/1‑150 (MQ=255)
                                                        aCAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCtt                                                                                             >  2:4633582/1‑150 (MQ=255)
                                                          aaCCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTcc                                                                                           <  1:3330178/150‑1 (MQ=255)
                                                          aaCCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTcc                                                                                           >  2:1489184/1‑150 (MQ=255)
                                                              aCTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATtc                                                                                       >  2:824834/1‑150 (MQ=255)
                                                               cTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATtct                                                                                      <  1:3129735/150‑1 (MQ=255)
                                                               cTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATtct                                                                                      <  2:2618506/150‑1 (MQ=255)
                                                                  acaacaCTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGt                                                                                   >  1:3109724/1‑150 (MQ=255)
                                                                         tGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGGTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAagcata                                                                            >  2:3356958/1‑150 (MQ=255)
                                                                           aaCCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAagcatagc                                                                          >  1:1043244/1‑150 (MQ=255)
                                                                           aaCCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAagcatagc                                                                          >  2:302185/1‑150 (MQ=255)
                                                                             cccACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAagcatagcat                                                                        <  2:3929092/150‑1 (MQ=255)
                                                                                aCGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTcg                                                                     <  1:3239521/150‑1 (MQ=255)
                                                                                aCGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTcg                                                                     <  1:2939764/150‑1 (MQ=255)
                                                                                   cttcttcttGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCgg                                                                  <  2:1567569/150‑1 (MQ=255)
                                                                                      cttcttGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAg                                                               >  2:2702554/1‑150 (MQ=255)
                                                                                        tcttGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCt                                                             <  1:3824346/150‑1 (MQ=255)
                                                                                        tcttGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCt                                                             <  1:1497519/150‑1 (MQ=255)
                                                                                          ttGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTtg                                                           <  2:4185226/150‑1 (MQ=255)
                                                                                            gCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTtgtg                                                         >  2:4567959/1‑150 (MQ=255)
                                                                                            gCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTtgtg                                                         >  2:1037357/1‑150 (MQ=255)
                                                                                            gCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTtgtg                                                         <  1:919408/150‑1 (MQ=255)
                                                                                              cccTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGtt                                                       <  2:3556350/150‑1 (MQ=255)
                                                                                                  gTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGa                                                   >  1:3481481/1‑150 (MQ=255)
                                                                                                  gTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGa                                                   >  1:2540767/1‑150 (MQ=255)
                                                                                                   ttCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAg                                                  <  2:1241737/150‑1 (MQ=255)
                                                                                                       gTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCa                                              <  2:188588/150‑1 (MQ=255)
                                                                                                        tAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCaa                                             >  2:3460072/1‑150 (MQ=255)
                                                                                                         aaTCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCaaa                                            >  1:3698340/1‑150 (MQ=255)
                                                                                                         aaTCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAACATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCaaa                                            >  2:2465108/1‑150 (MQ=255)
                                                                                                          aTCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCaaaa                                           <  1:4480425/150‑1 (MQ=255)
                                                                                                            cgcgcgTGATTTTCTTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCACGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCCTTCTCGTAAGCATAGCATTCGCGGTCGCTTGTGTTGGGAGGGCACAAAt                                         >  2:549199/1‑150 (MQ=255)
                                                                                                            cgcgcgTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAAt                                         >  2:549098/1‑150 (MQ=255)
                                                                                                              cgcgTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATaa                                       <  1:2898104/150‑1 (MQ=255)
                                                                                                               gcgTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAac                                      <  1:636319/150‑1 (MQ=255)
                                                                                                                  tGATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACt                                   >  1:2190927/1‑150 (MQ=255)
                                                                                                                   gATTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTa                                  >  2:3895141/1‑150 (MQ=255)
                                                                                                                   gATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTa                                  >  1:4179496/1‑150 (MQ=255)
                                                                                                                    aTTTTATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTAt                                 >  2:4197149/1‑150 (MQ=255)
                                                                                                                    aTTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTAt                                 >  2:3996492/1‑150 (MQ=255)
                                                                                                                    aTTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTAt                                 <  1:3614575/150‑1 (MQ=255)
                                                                                                                    aTTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTAt                                 <  1:1222023/150‑1 (MQ=255)
                                                                                                                      tttATTGATTGAAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATcc                               <  2:3517638/150‑1 (MQ=255)
                                                                                                                       ttATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCg                              >  2:3677085/1‑150 (MQ=255)
                                                                                                                       ttATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCg                              >  2:836331/1‑150 (MQ=255)
                                                                                                                         attgattgAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCa                            >  2:2578098/1‑150 (MQ=255)
                                                                                                                          ttgattgaAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCaa                           <  1:3498538/150‑1 (MQ=255)
                                                                                                                            gattgaAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAAtt                         >  1:2478008/1‑150 (MQ=255)
                                                                                                                             attgaAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAAttt                        <  1:1484521/150‑1 (MQ=255)
                                                                                                                              ttgaAAGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATttg                       >  1:2201853/1‑150 (MQ=255)
                                                                                                                                 aaaGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATttgttg                    <  1:3154510/150‑1 (MQ=255)
                                                                                                                                 aaaGAGGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCAGTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATttgttg                    <  1:3877929/150‑1 (MQ=255)
                                                                                                                                   agagGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGAt                  >  1:988027/1‑150 (MQ=255)
                                                                                                                                   agagGTGGTTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGAt                  <  2:2943580/150‑1 (MQ=255)
                                                                                                                                      ggtggtTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCaa               <  2:798022/150‑1 (MQ=255)
                                                                                                                                        tggtTAAACCGTTAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATa             <  2:1072349/150‑1 (MQ=255)
                                                                                                                                         ggtTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATaa            >  2:4407864/1‑150 (MQ=255)
                                                                                                                                           ttAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATAAGt          <  2:3798313/150‑1 (MQ=255)
                                                                                                                                           ttAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATAAGt          >  2:2146901/1‑150 (MQ=255)
                                                                                                                                            tAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATAAGtt         <  1:1033224/150‑1 (MQ=255)
                                                                                                                                            tAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCGTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATAAGtt         <  2:2089353/150‑1 (MQ=255)
                                                                                                                                             aaaCCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATAAGttt        <  2:2676951/150‑1 (MQ=255)
                                                                                                                                               aCCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATAAGTTTcc      <  1:4736935/150‑1 (MQ=255)
                                                                                                                                                ccGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATAAGTTTCCa     <  2:3934258/150‑1 (MQ=255)
                                                                                                                                                  gTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTTGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATAAGTTTCCAAg   >  1:782978/1‑150 (MQ=255)
                                                                                                                                                  gTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATAAGTTTCCAAg   <  2:4713958/150‑1 (MQ=255)
                                                                                                                                                   tGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATAAGTTTCCAAGc  >  2:697002/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTGATCAATCGTGCAAACGCCGATTTGCCCAAGAAACACTTTCTTATTGCTTTGTGACAACCACTGACAACACTGAACCCACGCTTCTTCTTGCCCCTGTTCAGTAATCGCGCGTGATTTTATTGATTGAAAGAGGTGGTTAAACCGTGAATTGCCATGAGCAAGAATAAGATTTGGGGATGAAAATAAGTGGCTCGCGCTCGCTTCCATTCTCGTAAGCATAGCATTCGCGGTAGCTTGTGTTGGGAGGGCAAAAATAACACTATCCGCAATTTGTTGATCAATAAGTTTCCAAGC  >  CP009273/2879301‑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: