Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,881,645 C→A 100% D202Y (GAT→TAT)  rnd ← ribonuclease D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,881,6450CA96.4% 758.2 / 16.9 253D202Y (GAT→TAT) rndribonuclease D
Reads supporting (aligned to +/- strand):  ref base C (4/5);  new base A (98/146);  total (102/151)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.69e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 253 total aligned reads displayed.

ACGCGCTACCGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATCACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATACCAGACATCCG  >  CP009273/1881502‑1881788
                                                                                                                                               |                                                                                                                                                
aCGCGCTACCGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGc                                                                                                                                              <  2:2603354/148‑1 (MQ=255)
aCGCGCTACCGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGc                                                                                                                                              <  2:822790/148‑1 (MQ=255)
aCGCGCTACCGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGc                                                                                                                                              <  2:583545/148‑1 (MQ=255)
  gcgcTACCGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCTGTCGCCTGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCa                                                                                                                                            <  1:708012/148‑1 (MQ=255)
    gcTACCGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAgg                                                                                                                                          <  1:481258/148‑1 (MQ=255)
     cTACCGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGc                                                                                                                                         <  2:891218/148‑1 (MQ=255)
      tACCGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCa                                                                                                                                        >  1:2808177/1‑148 (MQ=255)
       aCCGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCAt                                                                                                                                       <  2:1041424/148‑1 (MQ=255)
       aCCGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCAt                                                                                                                                       <  2:3262806/148‑1 (MQ=255)
           aCCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATCACGCCAGGCATCTTc                                                                                                                                   >  2:142647/1‑148 (MQ=255)
           aCCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTc                                                                                                                                   >  2:3308164/1‑148 (MQ=255)
              acaAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCgg                                                                                                                                <  1:2989254/148‑1 (MQ=255)
              acaAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCgg                                                                                                                                <  1:2999587/148‑1 (MQ=255)
                aaaTGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGgcg                                                                                                                              >  2:2143644/1‑148 (MQ=255)
                 aaTGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGgcgc                                                                                                                             >  1:1247117/1‑148 (MQ=255)
                   tGCTCGTCACGCACGACAAAGGTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCaa                                                                                                                          <  2:3138354/148‑1 (MQ=255)
                     ctctTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACg                                                                                                                        <  1:1363521/148‑1 (MQ=255)
                     ctctTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACg                                                                                                                        <  1:1494174/148‑1 (MQ=255)
                      tctTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGa                                                                                                                       <  1:2451606/148‑1 (MQ=255)
                        ttCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACt                                                                                                                     >  1:2217558/1‑148 (MQ=255)
                         tCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACtt                                                                                                                    >  2:1906800/1‑148 (MQ=255)
                         tCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACtt                                                                                                                    >  2:108514/1‑148 (MQ=255)
                          cacgcacgACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTc                                                                                                                   <  2:1456969/148‑1 (MQ=255)
                          cacgcacgACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTc                                                                                                                   >  1:2902156/1‑148 (MQ=255)
                           acgcacgACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTcc                                                                                                                  >  1:738889/1‑148 (MQ=255)
                             gcacgACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTg                                                                                                                <  1:1876693/148‑1 (MQ=255)
                                cgACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTgacg                                                                                                             <  1:939507/148‑1 (MQ=255)
                                     aaGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCACCGACTTCCTGACGACGCa                                                                                                        >  2:3085599/1‑148 (MQ=255)
                                     aaGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCa                                                                                                        <  1:2583438/148‑1 (MQ=255)
                                      aGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCAt                                                                                                       >  2:1686724/1‑148 (MQ=255)
                                      aGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCAt                                                                                                       <  2:3544982/148‑1 (MQ=255)
                                       gTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGACAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCAtt                                                                                                      <  1:1083888/148‑1 (MQ=255)
                                        ttCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCAttt                                                                                                     >  2:1211009/1‑148 (MQ=255)
                                            ccGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCAt                                                                                                 >  1:1221029/1‑148 (MQ=255)
                                             cGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATc                                                                                                <  2:3404558/148‑1 (MQ=255)
                                              gCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAAACGACTTCCTGACGACGCATTTGCATc                                                                                                <  2:2179359/148‑1 (MQ=255)
                                              gCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCa                                                                                               <  1:2381036/148‑1 (MQ=255)
                                              gCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCa                                                                                               <  2:374751/148‑1 (MQ=255)
                                                      gcgcTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATCACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCAt                                                                                       >  1:1392557/1‑148 (MQ=255)
                                                      gcgcTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATCACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCAt                                                                                       >  1:3161825/1‑148 (MQ=255)
                                                       cgcTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCAtt                                                                                      <  1:2688839/148‑1 (MQ=255)
                                                       cgcTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCAtt                                                                                      <  2:231687/148‑1 (MQ=255)
                                                       cgcTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCAtt                                                                                      <  2:2948287/148‑1 (MQ=255)
                                                        gcTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTACAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTc                                                                                     <  2:1233904/148‑1 (MQ=255)
                                                          tctcGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCAt                                                                                   >  2:3018684/1‑148 (MQ=255)
                                                            tcGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATcc                                                                                 <  1:1520541/148‑1 (MQ=255)
                                                               cgcCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAgc                                                                              <  1:3354834/148‑1 (MQ=255)
                                                                    tGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGcgccgc                                                                         >  1:2378141/1‑148 (MQ=255)
                                                                     gcgcAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCa                                                                        <  2:2684962/148‑1 (MQ=255)
                                                                     gcgcAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCa                                                                        <  2:990905/148‑1 (MQ=255)
                                                                      cgcAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTTCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAg                                                                       >  1:2472462/1‑148 (MQ=255)
                                                                         aGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTa                                                                    >  1:2106509/1‑148 (MQ=255)
                                                                           tCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTagc                                                                  <  1:2828556/148‑1 (MQ=255)
                                                                              ccAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTagccag                                                               <  1:827685/148‑1 (MQ=255)
                                                                                aGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTagccagcc                                                             >  2:2649937/1‑148 (MQ=255)
                                                                                  tCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCgg                                                           >  2:3319706/1‑148 (MQ=255)
                                                                                   cGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGa                                                          <  1:1711441/148‑1 (MQ=255)
                                                                                     gCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAgg                                                        >  2:295796/1‑148 (MQ=255)
                                                                                     gCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAgg                                                        <  2:1174941/148‑1 (MQ=255)
                                                                                     gCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAgg                                                        >  2:2381182/1‑148 (MQ=255)
                                                                                      cTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGc                                                       <  1:1644953/148‑1 (MQ=255)
                                                                                      cTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGc                                                       >  2:2132099/1‑148 (MQ=255)
                                                                                        aaCAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCt                                                     <  1:3247815/148‑1 (MQ=255)
                                                                                        aaCAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCt                                                     >  1:2683513/1‑148 (MQ=255)
                                                                                         aCAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCATGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTc                                                    >  1:870798/1‑148 (MQ=255)
                                                                                          cAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTcc                                                   >  1:810575/1‑148 (MQ=255)
                                                                                          cAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTcc                                                   <  1:244321/148‑1 (MQ=255)
                                                                                           aGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCg                                                  <  2:2214553/148‑1 (MQ=255)
                                                                                             ttGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGtt                                                <  1:1900387/148‑1 (MQ=255)
                                                                                               gCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTc                                              <  2:2837467/148‑1 (MQ=255)
                                                                                                   aCAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTAcc                                          <  2:1619891/148‑1 (MQ=255)
                                                                                                   aCAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTAcc                                          <  1:2691262/148‑1 (MQ=255)
                                                                                                    cAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCa                                         >  2:1662959/1‑148 (MQ=255)
                                                                                                    cAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCa                                         >  2:625527/1‑148 (MQ=255)
                                                                                                          aGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCt                                   <  2:1099874/148‑1 (MQ=255)
                                                                                                           gTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCtt                                  >  1:401681/1‑148 (MQ=255)
                                                                                                           gTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCtt                                  <  1:3433054/148‑1 (MQ=255)
                                                                                                              ggCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGc                               >  1:2728789/1‑148 (MQ=255)
                                                                                                               gCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCg                              >  2:1353642/1‑148 (MQ=255)
                                                                                                                    tgCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGAt                         <  2:3472536/148‑1 (MQ=255)
                                                                                                                     gCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATc                        >  1:124670/1‑148 (MQ=255)
                                                                                                                      cGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCg                       <  1:3298892/148‑1 (MQ=255)
                                                                                                                        tAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGt                     <  2:48527/148‑1 (MQ=255)
                                                                                                                          aTTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTaa                   <  1:1902033/148‑1 (MQ=255)
                                                                                                                           ttGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAAc                  >  1:2276398/1‑148 (MQ=255)
                                                                                                                               cAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGAt              <  1:2579176/148‑1 (MQ=255)
                                                                                                                               cAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGAt              <  1:1206943/148‑1 (MQ=255)
                                                                                                                                aGGCATTGGTGATATCACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATa             <  2:1773225/148‑1 (MQ=255)
                                                                                                                                aGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATa             <  2:1608778/148‑1 (MQ=255)
                                                                                                                                 ggCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATAc            >  1:29469/1‑148 (MQ=255)
                                                                                                                                 ggCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATAc            <  2:1949764/148‑1 (MQ=255)
                                                                                                                                  gCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATAcc           <  2:740416/148‑1 (MQ=255)
                                                                                                                                   cATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATACCa          >  1:1810083/1‑148 (MQ=255)
                                                                                                                                      tGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATACCAGAc       <  1:2816793/148‑1 (MQ=255)
                                                                                                                                      tGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATACCAGAc       >  2:2149924/1‑148 (MQ=255)
                                                                                                                                       ggTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATACCAGACa      >  2:3487117/1‑148 (MQ=255)
                                                                                                                                        gTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATACCAGACAt     >  1:2069423/1‑148 (MQ=255)
                                                                                                                                        gTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATACCAGACAt     >  1:554666/1‑148 (MQ=255)
                                                                                                                                         tGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATACCAGACATc    >  1:1107616/1‑148 (MQ=255)
                                                                                                                                           atatAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATACCAGACATccg  >  1:971305/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                                
ACGCGCTACCGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATCACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATACCAGACATCCG  >  CP009273/1881502‑1881788

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: