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,6680GT99.6% 737.2 / ‑4.2 231V637V (GTC→GTA) ygcBR‑loop helicase‑annealase Cas3 needed for Cascade anti‑viral activity
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (101/129);  total (101/130)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-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 231 total aligned reads displayed.

AAATGGTGGCAACAGGAATCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACGACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCGTGTTATTTA  >  CP009273/2878523‑2878810
                                                                                                                                                 |                                                                                                                                              
aaaTGGTGGCAACAGGAATCTCAAAACCAGCGGGACGATATGTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTAc                                                                                                                                              <  2:1847233/148‑1 (MQ=255)
     gtggCAACAGGAATCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTggg                                                                                                                                         <  1:3232843/148‑1 (MQ=255)
      tggCAACAGGAATCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGt                                                                                                                                        <  2:3480169/148‑1 (MQ=255)
      tggCAACAGGAATCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGt                                                                                                                                        <  1:3102056/148‑1 (MQ=255)
      tggCAACAGGAATCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCCGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGt                                                                                                                                        <  2:3414584/148‑1 (MQ=255)
           aCAGGAATCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGa                                                                                                                                   <  2:1157703/148‑1 (MQ=255)
                 aTCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTa                                                                                                                             <  1:3286641/148‑1 (MQ=255)
                 aTCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTa                                                                                                                             >  2:3406703/1‑148 (MQ=255)
                   ctcAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATc                                                                                                                           <  2:2103897/148‑1 (MQ=255)
                     cAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCg                                                                                                                         >  2:2664062/1‑148 (MQ=255)
                       aaaCCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAACTCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTc                                                                                                                       <  1:1657666/148‑1 (MQ=255)
                       aaaCCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTc                                                                                                                       >  1:3555057/1‑148 (MQ=255)
                       aaaCCAGCGGGACGATATATGCGATGATGTCGATGTAAACGGACCAATCGTTGTAAAAGCAAATCTTCAGGACAATGATGAGTACTAAAACAATCAAAATCAAAGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTc                                                                                                                       >  2:1758883/1‑148 (MQ=255)
                           cAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGTGACAAGTATCCGTCCAAc                                                                                                                   <  2:2973408/148‑1 (MQ=255)
                            aGCGTGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACa                                                                                                                  >  1:757085/1‑148 (MQ=255)
                             gCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACAt                                                                                                                 >  2:306194/1‑148 (MQ=255)
                               gggACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACAttt                                                                                                               <  2:1549619/148‑1 (MQ=255)
                                 gACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCg                                                                                                             >  2:2354442/1‑148 (MQ=255)
                                 gACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCg                                                                                                             >  2:967198/1‑148 (MQ=255)
                                 gACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCACGACAATGCTGAGTAATTAACCAATCAACATCAACGTCGAGTGATTGTTCCACTACCTGTGTTGCGACAAGTATCAGTCCAACATTTCg                                                                                                             >  2:2281050/1‑148 (MQ=255)
                                  aCGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGc                                                                                                            >  1:2584669/1‑148 (MQ=255)
                                    gATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCtt                                                                                                          <  1:1910593/148‑1 (MQ=255)
                                    gATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCtt                                                                                                          <  1:3225011/148‑1 (MQ=255)
                                     atatTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTc                                                                                                         <  1:3054407/148‑1 (MQ=255)
                                      tatTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTcc                                                                                                        >  1:2641948/1‑148 (MQ=255)
                                       atTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTccc                                                                                                       <  1:1443377/148‑1 (MQ=255)
                                        tttGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCa                                                                                                      <  2:1115240/148‑1 (MQ=255)
                                        tttGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCa                                                                                                      >  2:134692/1‑148 (MQ=255)
                                        tttGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCa                                                                                                      >  1:2896262/1‑148 (MQ=255)
                                         ttGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCAt                                                                                                     >  1:2022344/1‑148 (MQ=255)
                                            cGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCAtgat                                                                                                  >  1:2998859/1‑145 (MQ=255)
                                               tgatgGCGATTTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGc                                                                                               <  1:176735/148‑1 (MQ=255)
                                                       atcgccccGGCCCACTCGTGCCAACAGCACCTCTGCAGGACAATGCTGAGAAATTAACCCAACAAACCCAACCTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTg                                                                                       <  2:3167854/141‑1 (MQ=255)
                                                       aTGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTg                                                                                       >  2:694161/1‑148 (MQ=255)
                                                       aTGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTg                                                                                       <  2:2480945/148‑1 (MQ=255)
                                                        tGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGc                                                                                      >  2:1153182/1‑148 (MQ=255)
                                                        tGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGc                                                                                      >  1:3226874/1‑148 (MQ=255)
                                                         gTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCt                                                                                     >  1:2813822/1‑148 (MQ=255)
                                                          tAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCta                                                                                    <  2:1617147/148‑1 (MQ=255)
                                                          tAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCta                                                                                    <  2:3279890/148‑1 (MQ=255)
                                                           aaaCGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCtaa                                                                                   >  1:627019/1‑148 (MQ=255)
                                                           aaaCGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCtaa                                                                                   >  1:3187074/1‑148 (MQ=255)
                                                            aaCGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCtaat                                                                                  <  2:2895611/148‑1 (MQ=255)
                                                            aaCGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCtaat                                                                                  >  1:2459778/1‑148 (MQ=255)
                                                             aCGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCtaata                                                                                 <  1:3212670/148‑1 (MQ=255)
                                                              cGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCtaataa                                                                                >  1:84725/1‑148 (MQ=255)
                                                                gCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACt                                                                              <  1:1950506/148‑1 (MQ=255)
                                                                gCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACt                                                                              <  2:1434375/148‑1 (MQ=255)
                                                                gCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACt                                                                              <  1:145848/148‑1 (MQ=255)
                                                                   cAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGa                                                                           <  2:865066/148‑1 (MQ=255)
                                                                    aaTCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGAt                                                                          <  1:714645/148‑1 (MQ=255)
                                                                     aTCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGAtt                                                                         >  2:3058951/1‑148 (MQ=255)
                                                                     aTCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGAtt                                                                         <  1:2072752/148‑1 (MQ=255)
                                                                         ttGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTctct                                                                     >  2:2557449/1‑148 (MQ=255)
                                                                         ttGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTctct                                                                     >  2:2731591/1‑148 (MQ=255)
                                                                         ttGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTctct                                                                     <  2:914172/148‑1 (MQ=255)
                                                                         ttGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATGTTTGCCGAAATTGCTAATAACTCGATTctct                                                                     >  2:78323/1‑148 (MQ=255)
                                                                             aaaaGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCttttt                                                                 >  2:1171998/1‑148 (MQ=255)
                                                                              aaaGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTc                                                                <  1:1673428/148‑1 (MQ=255)
                                                                              aaaGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTc                                                                <  1:580559/148‑1 (MQ=255)
                                                                               aaGCAAATCTGCAGTACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCa                                                               <  1:3281889/148‑1 (MQ=255)
                                                                               aaGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGTGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCa                                                               <  2:2184646/148‑1 (MQ=255)
                                                                               aaGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCa                                                               >  2:854766/1‑148 (MQ=255)
                                                                                aGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTAGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATGCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCac                                                              <  1:2666096/148‑1 (MQ=255)
                                                                                  cAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCacga                                                            <  1:3406703/148‑1 (MQ=255)
                                                                                  cAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCacga                                                            <  1:2961749/148‑1 (MQ=255)
                                                                                      tCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGAt                                                        <  2:965045/148‑1 (MQ=255)
                                                                                       cTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATc                                                       >  1:1547535/1‑148 (MQ=255)
                                                                                         gCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGt                                                     >  2:2412251/1‑148 (MQ=255)
                                                                                           aGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTc                                                   <  2:890724/148‑1 (MQ=255)
                                                                                              aCAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGc                                                >  1:2280586/1‑148 (MQ=255)
                                                                                                aaTGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGt                                              >  2:1804433/1‑148 (MQ=255)
                                                                                                  tGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTaa                                            >  1:1862721/1‑148 (MQ=255)
                                                                                                  tGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTaa                                            >  2:2616759/1‑148 (MQ=255)
                                                                                                   gCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTaaa                                           >  2:981239/1‑148 (MQ=255)
                                                                                                    cTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAg                                          <  2:2584669/148‑1 (MQ=255)
                                                                                                        gTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAgcgcg                                      >  2:2137115/1‑148 (MQ=255)
                                                                                                          aaTTAAACAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCa                                    <  1:2545719/148‑1 (MQ=255)
                                                                                                            ttAACCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGTCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATg                                  >  1:933098/1‑148 (MQ=255)
                                                                                                               aCCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGaaa                               <  1:1243074/148‑1 (MQ=255)
                                                                                                               aCCAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGaaa                               >  2:1037553/1‑148 (MQ=255)
                                                                                                                ccAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAAc                              >  2:1984062/1‑148 (MQ=255)
                                                                                                                ccAATCAAAATCAACGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAAc                              >  2:1671808/1‑148 (MQ=255)
                                                                                                                   atcaaaatcaaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACaaa                           >  1:3024389/1‑148 (MQ=255)
                                                                                                                   atcaaaatcaaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACaaa                           >  1:465532/1‑148 (MQ=255)
                                                                                                                    tcaaaatcaaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAAt                          <  2:1260230/148‑1 (MQ=255)
                                                                                                                       aaatcaaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCta                       >  2:2382389/1‑148 (MQ=255)
                                                                                                                        aatcaaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCtat                      >  2:894288/1‑148 (MQ=255)
                                                                                                                        aatcaaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCtat                      <  2:1047740/148‑1 (MQ=255)
                                                                                                                         atcaaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCtata                     <  1:2575587/148‑1 (MQ=255)
                                                                                                                          tcaaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTatat                    <  2:3566388/148‑1 (MQ=255)
                                                                                                                          tcaaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTatat                    <  1:2731591/148‑1 (MQ=255)
                                                                                                                          tcaaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTatat                    >  1:1208165/1‑148 (MQ=255)
                                                                                                                           caaCGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAATCGCGCATGAAACAAATCTATATc                   >  1:217350/1‑148 (MQ=255)
                                                                                                                              cGTCGAGTGATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCTCGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTAc                >  2:2171430/1‑148 (MQ=255)
                                                                                                                                      gATTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCgtgt        <  1:694161/148‑1 (MQ=255)
                                                                                                                                       aTTGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCGTGtt       <  1:1887203/148‑1 (MQ=255)
                                                                                                                                         tGTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCGTGTTAt     <  2:2230100/148‑1 (MQ=255)
                                                                                                                                          gTTCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCGTGTTAtt    <  1:1101985/148‑1 (MQ=255)
                                                                                                                                           ttCCACTACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCGTGTTAttt   <  1:1034375/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
AAATGGTGGCAACAGGAATCTCAAAACCAGCGGGACGATATTTGCGATGATGGCGATGTAAACGGCCCAATCGTTGGAAAAGCAAATCTGCAGGACAATGCTGAGTAATTAACCAATCAAAATCAACGTCGAGTGATTGTTCCACGACCTGGGTTGCGACAAGTATCCGTCCAACATTTCGCTTCCCATTTTTGCCGAAATTGCTAATAACTCGATTCTCTTTTTCACGACGATCGTTCAGCGTAAAGCGCGCATGAAACAAATCTATATCTACTTGCGTGTTATTTA  >  CP009273/2878523‑2878810

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑