Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,878,668 G→T 100% V637V (GTC→GTA)  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,878,6680GT98.9% 887.2 / 0.0 281V637V (GTC→GTA) ygcBR‑loop helicase‑annealase Cas3 needed for Cascade anti‑viral activity
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base T (134/143);  total (137/144)
Fisher's exact test for biased strand distribution p-value = 6.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.14e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 281 total aligned reads displayed.

CAAAATGGTGGCAACAGGAATCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACGACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCGTGTTATTTAGCTC  >  CP009273/2878521‑2878814
                                                                                                                                                   |                                                                                                                                                  
  aaaTGGTGGCAACATGAATCTCAAAACCAGCGGGACGATATTTGCGATGAAGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAAAATGCAGGACAATGCTGAGTAATTAACAAATAAAAATAAAAGACGAGTCATTGATCCACTAcaa                                                                                                                                                >  2:3406297/1‑148 (MQ=255)
    aTGGTGGCAACAGGAATCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAAACTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTgg                                                                                                                                              >  1:2131422/1‑150 (MQ=255)
     tggtggCAACAGGAATCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTggg                                                                                                                                             >  2:3229962/1‑150 (MQ=255)
      ggtggCAACAGGAATCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGt                                                                                                                                            >  1:1940616/1‑150 (MQ=255)
              cAGGAATCTCAAAACCAGCGGTACGATTTTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACaa                                                                                                                                    <  1:2621844/150‑1 (MQ=255)
               aGGAATCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAg                                                                                                                                   <  1:4805075/150‑1 (MQ=255)
                 gAATCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTa                                                                                                                                 >  2:4396050/1‑150 (MQ=255)
                  aaTCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTAt                                                                                                                                <  1:4231474/150‑1 (MQ=255)
                   aTCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATc                                                                                                                               <  1:1568068/150‑1 (MQ=255)
                    tctcAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATcc                                                                                                                              <  2:4234949/150‑1 (MQ=255)
                    tctcAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATcc                                                                                                                              <  1:3230758/150‑1 (MQ=255)
                      tcAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGt                                                                                                                            >  2:3823362/1‑150 (MQ=255)
                       cAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTc                                                                                                                           >  2:611152/1‑150 (MQ=255)
                         aaaCCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCa                                                                                                                         <  2:2352872/150‑1 (MQ=255)
                         aaaCCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCa                                                                                                                         <  1:966722/150‑1 (MQ=255)
                          aaCCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCaa                                                                                                                        >  2:3874511/1‑150 (MQ=255)
                          aaCCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCaa                                                                                                                        >  1:380550/1‑150 (MQ=255)
                           aCCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAAc                                                                                                                       <  1:4481816/150‑1 (MQ=255)
                               gCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACAttt                                                                                                                   <  1:2547993/150‑1 (MQ=255)
                                cGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTc                                                                                                                  <  2:2092167/150‑1 (MQ=255)
                                     cGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTc                                                                                                             >  1:1203930/1‑150 (MQ=255)
                                       atatTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTccc                                                                                                           <  1:4835841/150‑1 (MQ=255)
                                       atatTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTccc                                                                                                           <  1:318144/150‑1 (MQ=255)
                                        tatTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCa                                                                                                          <  1:3325241/150‑1 (MQ=255)
                                          tttGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCAtt                                                                                                        <  1:4980573/150‑1 (MQ=255)
                                              cGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTg                                                                                                    >  1:1658384/1‑150 (MQ=255)
                                               gatgatGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGc                                                                                                   >  2:1999331/1‑150 (MQ=255)
                                                atgatgGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGcc                                                                                                  <  1:1858605/150‑1 (MQ=255)
                                                 tgatgGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCg                                                                                                 >  2:1571823/1‑150 (MQ=255)
                                                  gatgGCGATGTAAACGGCACAATCGTTGGAAAAGCAATTCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGa                                                                                                <  1:2419598/150‑1 (MQ=255)
                                                     ggCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAAt                                                                                             <  2:848620/150‑1 (MQ=255)
                                                     ggCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATGTGTGCCGAAAt                                                                                             >  2:898852/1‑150 (MQ=255)
                                                       cGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTg                                                                                           <  1:1111780/150‑1 (MQ=255)
                                                       cGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTg                                                                                           <  2:4361370/150‑1 (MQ=255)
                                                         aTGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCt                                                                                         >  1:3458395/1‑150 (MQ=255)
                                                           gTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCtaa                                                                                       <  1:4629691/150‑1 (MQ=255)
                                                            tAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCtaat                                                                                      >  1:878497/1‑150 (MQ=255)
                                                              aaCGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCACCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCtaataa                                                                                    <  1:431036/150‑1 (MQ=255)
                                                              aaCGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCtaataa                                                                                    >  1:1151539/1‑150 (MQ=255)
                                                               aCGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAAc                                                                                   <  2:75735/150‑1 (MQ=255)
                                                                     cAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGAtt                                                                             >  1:1283594/1‑150 (MQ=255)
                                                                      aaTCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATtc                                                                            <  2:392090/150‑1 (MQ=255)
                                                                       aTCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATtct                                                                           <  2:270184/150‑1 (MQ=255)
                                                                       aTCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATtct                                                                           <  2:4552302/150‑1 (MQ=255)
                                                                          gTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCtt                                                                        >  1:2689619/1‑150 (MQ=255)
                                                                           ttGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCttt                                                                       >  1:5065895/1‑150 (MQ=255)
                                                                              gAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATACCTCGATTCTCTTTTTc                                                                    >  2:4846445/1‑150 (MQ=255)
                                                                              gAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTc                                                                    >  1:771499/1‑150 (MQ=255)
                                                                              gAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTc                                                                    >  1:1433275/1‑150 (MQ=255)
                                                                              gAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTc                                                                    <  1:4084100/150‑1 (MQ=255)
                                                                               aaaaGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCa                                                                   >  2:4266855/1‑150 (MQ=255)
                                                                               aaaaGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCa                                                                   <  2:1891653/150‑1 (MQ=255)
                                                                                 aaGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCacg                                                                 >  2:3495485/1‑150 (MQ=255)
                                                                                 aaGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCacg                                                                 >  2:1689233/1‑150 (MQ=255)
                                                                                   gCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCacgac                                                               <  2:412060/150‑1 (MQ=255)
                                                                                    cAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCacgacg                                                              <  1:4519608/150‑1 (MQ=255)
                                                                                     aaaTCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCAcgacga                                                             >  1:1230448/1‑150 (MQ=255)
                                                                                     aaaTCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCAATTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCAcgacga                                                             >  2:452726/1‑150 (MQ=255)
                                                                                       aTCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATc                                                           <  2:3091990/150‑1 (MQ=255)
                                                                                       aTCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATc                                                           <  1:1015476/150‑1 (MQ=255)
                                                                                        tCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCg                                                          <  1:305117/150‑1 (MQ=255)
                                                                                         cTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGt                                                         <  2:2872543/150‑1 (MQ=255)
                                                                                            cAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCa                                                      >  2:2715396/1‑150 (MQ=255)
                                                                                             aGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAg                                                     >  1:1795876/1‑150 (MQ=255)
                                                                                               gACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCg                                                   >  1:4436826/1‑150 (MQ=255)
                                                                                               gACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCg                                                   >  1:4691396/1‑150 (MQ=255)
                                                                                                   aTGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTaaa                                               >  1:1717666/1‑150 (MQ=255)
                                                                                                      cTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAgcg                                            <  2:2949815/150‑1 (MQ=255)
                                                                                                       tGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAgcgc                                           >  2:1525347/1‑150 (MQ=255)
                                                                                                        gAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAgcgcg                                          <  1:2279249/150‑1 (MQ=255)
                                                                                                        gAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAgcgcg                                          <  1:1755890/150‑1 (MQ=255)
                                                                                                          gTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCa                                        <  1:1296851/150‑1 (MQ=255)
                                                                                                           tAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCAt                                       >  2:4025195/1‑150 (MQ=255)
                                                                                                              ttAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGaa                                    <  2:5087951/150‑1 (MQ=255)
                                                                                                              ttAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGaa                                    <  2:3812234/150‑1 (MQ=255)
                                                                                                                aaCCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAAc                                  >  1:3131335/1‑150 (MQ=255)
                                                                                                                   cAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACaaa                               <  1:4397870/150‑1 (MQ=255)
                                                                                                                   cAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACaaa                               >  2:3089776/1‑150 (MQ=255)
                                                                                                                     atcaaaatcaaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTCAAGCGCGCATGAAACAAATc                             >  1:714861/1‑150 (MQ=255)
                                                                                                                     atcaaaatcaaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATc                             >  2:516928/1‑150 (MQ=255)
                                                                                                                        aaaatcaaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCtat                          <  2:4716576/150‑1 (MQ=255)
                                                                                                                         aaatcaaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCtata                         <  1:28758/150‑1 (MQ=255)
                                                                                                                          aatcaaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTatat                        >  2:255520/1‑150 (MQ=255)
                                                                                                                          aatcaaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTatat                        <  1:3495485/150‑1 (MQ=255)
                                                                                                                            tcaaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCt                      <  2:380550/150‑1 (MQ=255)
                                                                                                                             caaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTa                     >  1:295355/1‑150 (MQ=255)
                                                                                                                             caaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTa                     >  2:1145544/1‑150 (MQ=255)
                                                                                                                             caaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTa                     >  2:2888494/1‑150 (MQ=255)
                                                                                                                              aaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTAc                    <  2:3373552/150‑1 (MQ=255)
                                                                                                                               aCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACt                   <  2:1265382/150‑1 (MQ=255)
                                                                                                                                cGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACtt                  <  2:4436826/150‑1 (MQ=255)
                                                                                                                                    gAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCgt              >  2:1363579/1‑150 (MQ=255)
                                                                                                                                      gTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCgtgt            <  1:3905652/150‑1 (MQ=255)
                                                                                                                                       tGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCGTGtt           <  2:2139921/150‑1 (MQ=255)
                                                                                                                                       tGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCGTGtt           <  2:3351262/150‑1 (MQ=255)
                                                                                                                                           tGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCGTGTTAttt       >  2:3148146/1‑150 (MQ=255)
                                                                                                                                             ttCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTCGCGTGTTATTTAg     <  2:5027751/150‑1 (MQ=255)
                                                                                                                                             ttCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATATATATCTACTTGCGTGTTATTTAg     <  1:2467619/150‑1 (MQ=255)
                                                                                                                                               ccACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCGTGTTATTTAGCt   >  2:2285378/1‑150 (MQ=255)
                                                                                                                                                cACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCGTGTTATTTAGCTc  >  2:302018/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CAAAATGGTGGCAACAGGAATCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACGACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCGTGTTATTTAGCTC  >  CP009273/2878521‑2878814

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: ‑