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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,881,6450CA94.5% 624.6 / 26.4 220D202Y (GAT→TAT) rndribonuclease D
Reads supporting (aligned to +/- strand):  ref base C (4/8);  new base A (104/104);  total (108/112)
Fisher's exact test for biased strand distribution p-value = 3.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.11e-01

Only 100 of 220 total aligned reads displayed.

ATAACGCGCTACCGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATCACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATACCAGACATCCGCC  >  CP009273/1881499‑1881790
                                                                                                                                                  |                                                                                                                                                 
    cgcgCTACCGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGcc                                                                                                                                              >  1:3064431/1‑148 (MQ=255)
    cgcgCTACCGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGcc                                                                                                                                              <  2:1187602/148‑1 (MQ=255)
          aCCGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCAt                                                                                                                                        >  1:2878565/1‑148 (MQ=255)
          aCCGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCAt                                                                                                                                        <  2:267047/148‑1 (MQ=255)
          aCCGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCAt                                                                                                                                        <  2:266978/148‑1 (MQ=255)
           ccGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATc                                                                                                                                       <  2:3003829/148‑1 (MQ=255)
               ccACAAATGCTCTTCACGCACGACAAAGTTCACCTCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTcc                                                                                                                                   <  1:938149/148‑1 (MQ=255)
                cacaAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCg                                                                                                                                  >  2:487976/1‑148 (MQ=255)
                   aaaTGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGgcg                                                                                                                               <  1:1999062/148‑1 (MQ=255)
                     aTGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCa                                                                                                                             <  2:281069/148‑1 (MQ=255)
                      tGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCaa                                                                                                                            >  1:1340997/1‑148 (MQ=255)
                       gCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATCTAACGCCAGGCATCTTCCGGCGCAAc                                                                                                                           >  2:1733875/1‑148 (MQ=255)
                       gCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAAc                                                                                                                           >  2:3496190/1‑148 (MQ=255)
                       gCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAAc                                                                                                                           >  1:1816363/1‑148 (MQ=255)
                         tctTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGa                                                                                                                         <  2:1443125/148‑1 (MQ=255)
                         tctTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGa                                                                                                                         <  1:2640201/148‑1 (MQ=255)
                          ctTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCCTCTTCCGGCGCAACGAc                                                                                                                        >  2:2431795/1‑148 (MQ=255)
                          ctTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGAc                                                                                                                        >  2:2431638/1‑148 (MQ=255)
                            tCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACtt                                                                                                                      >  2:1158438/1‑148 (MQ=255)
                             cacgcacgACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTc                                                                                                                     <  1:1166596/148‑1 (MQ=255)
                             cacgcacgACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTc                                                                                                                     <  1:806271/148‑1 (MQ=255)
                              acgcacgACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTcc                                                                                                                    >  2:2680658/1‑148 (MQ=255)
                              acgcacgACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTcc                                                                                                                    >  2:3335507/1‑148 (MQ=255)
                              acgcacgACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTcc                                                                                                                    >  2:2680882/1‑148 (MQ=255)
                              acgcacgACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTcc                                                                                                                    >  2:3487666/1‑148 (MQ=255)
                                gcacgACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTg                                                                                                                  <  2:1536614/148‑1 (MQ=255)
                                 cacgACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTga                                                                                                                 >  1:511084/1‑148 (MQ=255)
                                   cgACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATCACGCCAGGCATCTTCCGGCGCAACGACTTCCTgacg                                                                                                               <  2:3364935/148‑1 (MQ=255)
                                   cgACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTgacg                                                                                                               >  2:3320174/1‑148 (MQ=255)
                                     aCAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTgacgac                                                                                                             >  1:1556955/1‑148 (MQ=255)
                                      cAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGacgacg                                                                                                            >  1:1185818/1‑148 (MQ=255)
                                       aaaGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGc                                                                                                           <  2:498795/148‑1 (MQ=255)
                                        aaGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCa                                                                                                          <  2:1566809/148‑1 (MQ=255)
                                        aaGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCa                                                                                                          <  1:1293000/148‑1 (MQ=255)
                                           ttCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCAttt                                                                                                       >  2:2194053/1‑148 (MQ=255)
                                            tCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTg                                                                                                      >  1:2237542/1‑148 (MQ=255)
                                            tCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTg                                                                                                      >  2:2327723/1‑148 (MQ=255)
                                            tCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTg                                                                                                      <  1:902555/148‑1 (MQ=255)
                                                cGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATc                                                                                                  <  1:3546622/148‑1 (MQ=255)
                                                cGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATc                                                                                                  <  2:909243/148‑1 (MQ=255)
                                                 gCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCa                                                                                                 >  2:2226065/1‑148 (MQ=255)
                                                   cAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAgg                                                                                               >  2:2141727/1‑148 (MQ=255)
                                                      aTCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAggcgg                                                                                            >  2:1951233/1‑148 (MQ=255)
                                                      aTCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAggcgg                                                                                            >  2:2994386/1‑148 (MQ=255)
                                                       tCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCAGTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAggcggc                                                                                           >  2:2544899/1‑148 (MQ=255)
                                                        cgcgCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATCACGCCAGGCATCTTCCGGCGCAACGACTTCCTGAAGACGCATTTGCATCAGGCGGCa                                                                                          <  1:1492538/148‑1 (MQ=255)
                                                            ctctCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCa                                                                                      >  1:1600046/1‑148 (MQ=255)
                                                                 gcgcCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAg                                                                                 <  2:2423402/148‑1 (MQ=255)
                                                                    ccTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAgcgc                                                                              <  2:400419/148‑1 (MQ=255)
                                                                    ccTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAgcgc                                                                              <  2:470876/148‑1 (MQ=255)
                                                                    ccTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAgcgc                                                                              <  1:1346222/148‑1 (MQ=255)
                                                                      ttGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGcgccg                                                                            >  2:486781/1‑148 (MQ=255)
                                                                      ttGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGcgccg                                                                            <  1:2563677/148‑1 (MQ=255)
                                                                      ttGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGcgccg                                                                            >  2:2029924/1‑148 (MQ=255)
                                                                              tCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTagc                                                                    <  2:787953/148‑1 (MQ=255)
                                                                               cGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTagcc                                                                   <  2:29980/148‑1 (MQ=255)
                                                                                 ccAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTagccag                                                                 >  2:3154789/1‑148 (MQ=255)
                                                                                 ccAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTagccag                                                                 >  2:843346/1‑148 (MQ=255)
                                                                                     tCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCgg                                                             <  1:1531040/148‑1 (MQ=255)
                                                                                     tCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCgg                                                             <  1:2690449/148‑1 (MQ=255)
                                                                                          tAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGcc                                                        >  2:548072/1‑148 (MQ=255)
                                                                                            aCAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATCACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTc                                                      <  1:732161/148‑1 (MQ=255)
                                                                                            aCAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTc                                                      <  1:2083097/148‑1 (MQ=255)
                                                                                            aCAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTc                                                      <  1:1006010/148‑1 (MQ=255)
                                                                                             cAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTcc                                                     >  2:262014/1‑148 (MQ=255)
                                                                                             cAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTcc                                                     >  2:3311051/1‑148 (MQ=255)
                                                                                              aGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCg                                                    <  2:1384570/148‑1 (MQ=255)
                                                                                              aGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCg                                                    <  2:3175046/148‑1 (MQ=255)
                                                                                                ttGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGtt                                                  >  2:1498415/1‑148 (MQ=255)
                                                                                                 tGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTCGCCAGCCGGAGGCCTCCGttt                                                 >  2:2746079/1‑148 (MQ=255)
                                                                                                  gCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTc                                                <  1:540416/148‑1 (MQ=255)
                                                                                                    aGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTa                                              <  1:2773619/148‑1 (MQ=255)
                                                                                                      aCAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTAcc                                            >  1:1125153/1‑148 (MQ=255)
                                                                                                         ggCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATa                                         <  2:2638679/148‑1 (MQ=255)
                                                                                                         ggCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATa                                         >  2:2541562/1‑148 (MQ=255)
                                                                                                           ccAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATCACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAg                                       >  1:27049/1‑148 (MQ=255)
                                                                                                           ccAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAg                                       <  2:930582/148‑1 (MQ=255)
                                                                                                               ttGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTg                                   >  1:616126/1‑148 (MQ=255)
                                                                                                               ttGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTg                                   >  2:872471/1‑148 (MQ=255)
                                                                                                                  gCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCg                                >  1:1396657/1‑148 (MQ=255)
                                                                                                                     tgtgCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTg                             >  1:180873/1‑148 (MQ=255)
                                                                                                                      gtgCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGa                            <  2:3084054/148‑1 (MQ=255)
                                                                                                                      gtgCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGa                            >  1:460740/1‑148 (MQ=255)
                                                                                                                        gCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATc                          <  1:311638/148‑1 (MQ=255)
                                                                                                                         cGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCg                         <  1:3362136/148‑1 (MQ=255)
                                                                                                                          gTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGATTCTACCATAAGCTTGGCGGTGATCgg                        <  1:2321643/148‑1 (MQ=255)
                                                                                                                           tAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGt                       <  1:2939312/148‑1 (MQ=255)
                                                                                                                            aaTTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTa                      >  1:915784/1‑148 (MQ=255)
                                                                                                                              gtGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAAc                    <  2:2120503/147‑1 (MQ=255)
                                                                                                                                  cAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGAt                <  1:3318904/148‑1 (MQ=255)
                                                                                                                                  cAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGAt                >  1:3346767/1‑148 (MQ=255)
                                                                                                                                      cATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATACCa            >  1:3230864/1‑148 (MQ=255)
                                                                                                                                       aTTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATACCAg           <  1:2979604/148‑1 (MQ=255)
                                                                                                                                        ttGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATACCAGa          >  2:2514044/1‑148 (MQ=255)
                                                                                                                                         tGGTGATATCACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATACCAGAc         >  2:2935653/1‑148 (MQ=255)
                                                                                                                                          ggTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATACCAGACa        <  1:2116735/148‑1 (MQ=255)
                                                                                                                                           gTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCTTTTCTTCCATTAGCTTGGCGGTGATCGGTAACAGATACCAGACAt       <  2:2760748/148‑1 (MQ=255)
                                                                                                                                           gTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATACCAGACAt       <  1:2788159/148‑1 (MQ=255)
                                                                                                                                            tGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATACCAGACATc      <  2:450090/148‑1 (MQ=255)
                                                                                                                                                atCACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATACCAGACATccgcc  >  1:3390757/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
ATAACGCGCTACCGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATCACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATACCAGACATCCGCC  >  CP009273/1881499‑1881790

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