Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,689,180 T→G 100% E383D (GAA→GAC uidA ← beta‑D‑glucuronidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,689,1800TG96.7% 835.5 / 5.1 210E383D (GAA→GACuidAbeta‑D‑glucuronidase
Reads supporting (aligned to +/- strand):  ref base T (4/3);  new base G (103/100);  total (107/103)
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.95e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 210 total aligned reads displayed.

TTGCTTCCGCCAGTGGCGCGAAATATTCCCGTGCACCTTGCGGACGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTTTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACGATGCCATGTTCATCTGCCCAGTCGAGCATCTCT  >  CP009273/1689035‑1689324
                                                                                                                                                 |                                                                                                                                                
ttGCTTCCGCCAGTGGCGCGAAATATTCCCGTGCACCTTGCGGACGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTc                                                                                                                                                >  2:2228762/1‑148 (MQ=255)
 tGCTTCCGCCAGTGGCGCGAAATATTCCCGTGCACCTTGCGGACGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTcc                                                                                                                                               <  1:2940897/148‑1 (MQ=255)
     tCCGCCAGTGGCGCGAAATATTCCCGTGCACCTTGCGGACGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGt                                                                                                                                           <  1:244784/148‑1 (MQ=255)
      ccGCCAGTGGCGCGAAATATTCCCGTGCACCTTGCGGACGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGtt                                                                                                                                          >  1:533444/1‑148 (MQ=255)
         ccAGTGGCGCGAAATATTCCCGTGCACCTTGCGGACGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTTTCCCCGTTGAc                                                                                                                                       >  1:1375683/1‑148 (MQ=255)
         ccAGTGGCGCGAAATATTCCCGTGCACCTTGCGGACGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGAc                                                                                                                                       <  1:1504018/148‑1 (MQ=255)
          cAGTGGCGCGAAATATTCCCGTGCACCTTGCGGACGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACt                                                                                                                                      >  2:859567/1‑148 (MQ=255)
           aGTGGCGCGAAATATTCCCGTGCACCTTGCGGACGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTg                                                                                                                                     >  1:2566434/1‑148 (MQ=255)
             tGGCGCGAAATATTCCCGTGCACCTTGCGGACGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGcc                                                                                                                                   <  2:1333698/148‑1 (MQ=255)
              ggCGCGAAATATTCCCGTGCACCTTGCGGACGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCct                                                                                                                                  <  2:2163213/148‑1 (MQ=255)
              ggCGCGAAATATTCCCGTGCACCTTGCGGACGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCct                                                                                                                                  <  2:3033745/148‑1 (MQ=255)
                cgcgAAATATTCCCGTGCACCTTGCGGACGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCctct                                                                                                                                >  2:3305360/1‑148 (MQ=255)
                 gcgAAATATTCCCGTGCACCTTGCGGACGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCtt                                                                                                                               <  1:1453334/148‑1 (MQ=255)
                     aaTATTCCCGTGCACCTTGCGGACGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCt                                                                                                                           >  1:164502/1‑148 (MQ=255)
                      atatTCCCGTGCACCTTGCGGACGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTg                                                                                                                          >  1:1799396/1‑148 (MQ=255)
                       tatTCCCGTGCACCTTGCGGACGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGt                                                                                                                         >  1:297099/1‑148 (MQ=255)
                         ttCCCGTGCACCTTGCGGACGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTAc                                                                                                                       >  1:2022053/1‑148 (MQ=255)
                          tCCCGTGCACCTTGCGGACGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACa                                                                                                                      >  1:3001388/1‑148 (MQ=255)
                             cGTGCACCTTGCGGACGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGtt                                                                                                                   <  2:1108871/148‑1 (MQ=255)
                              gTGCACCTTGCGGACGGGTATCCGGTTCGTTGGCAATACTCCGCATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTc                                                                                                                  <  1:1577543/148‑1 (MQ=255)
                              gTGCACCTTGCGGACGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTc                                                                                                                  <  1:2045610/148‑1 (MQ=255)
                               tGCACCTTGCGGACGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCt                                                                                                                 <  2:1314279/148‑1 (MQ=255)
                                 cACCTTGCGGACGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCttt                                                                                                               <  2:656391/148‑1 (MQ=255)
                                   ccTTGCGGACGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCg                                                                                                             >  1:377991/1‑148 (MQ=255)
                                       gCGGACGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCtt                                                                                                         <  1:2633782/148‑1 (MQ=255)
                                         ggACGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTTTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCttgt                                                                                                       <  2:559899/148‑1 (MQ=255)
                                           aCGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCttgttg                                                                                                     >  1:1147613/1‑148 (MQ=255)
                                           aCGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCttgttg                                                                                                     <  1:1221038/148‑1 (MQ=255)
                                           aCGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCttgttg                                                                                                     <  1:80880/148‑1 (MQ=255)
                                           aCGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCttgttg                                                                                                     >  1:3518044/1‑148 (MQ=255)
                                           aCGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCttgttg                                                                                                     <  1:2805654/148‑1 (MQ=255)
                                           aCGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCttgttg                                                                                                     >  1:1760707/1‑148 (MQ=255)
                                               gTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCg                                                                                                 >  2:2616990/1‑148 (MQ=255)
                                               gTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCg                                                                                                 <  2:1281878/148‑1 (MQ=255)
                                                     ggTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGa                                                                                           <  1:752528/148‑1 (MQ=255)
                                                       ttCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGaaa                                                                                         >  1:2736551/1‑148 (MQ=255)
                                                       ttCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGaaa                                                                                         <  2:2022053/148‑1 (MQ=255)
                                                       ttCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGaaa                                                                                         >  2:227142/1‑148 (MQ=255)
                                                       ttCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGaaa                                                                                         >  2:3232650/1‑148 (MQ=255)
                                                       ttCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGaaa                                                                                         <  2:3090974/148‑1 (MQ=255)
                                                              gCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGc                                                                                  <  2:1929587/148‑1 (MQ=255)
                                                               cAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGcc                                                                                 <  1:1681404/148‑1 (MQ=255)
                                                                 aTACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTTTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTa                                                                               >  1:1978709/1‑148 (MQ=255)
                                                                 aTACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTa                                                                               >  2:408473/1‑148 (MQ=255)
                                                                  tACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTaa                                                                              <  2:3193707/148‑1 (MQ=255)
                                                                      ccACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAagag                                                                          >  1:2744321/1‑148 (MQ=255)
                                                                        acaTCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAagagag                                                                        <  1:2374511/148‑1 (MQ=255)
                                                                        acaTCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAagagag                                                                        >  1:538629/1‑148 (MQ=255)
                                                                        acaTCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAagagag                                                                        <  2:766991/148‑1 (MQ=255)
                                                                         caTCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAgg                                                                       <  2:337117/148‑1 (MQ=255)
                                                                         caTCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGAGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACAAATGCCTAAAGAGAgg                                                                       >  2:3323478/1‑148 (MQ=255)
                                                                             accacGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTaa                                                                   <  1:1083801/148‑1 (MQ=255)
                                                                             accacGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTaa                                                                   >  2:1454379/1‑148 (MQ=255)
                                                                                 cGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGcc                                                               >  1:2447986/1‑148 (MQ=255)
                                                                                  gCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCg                                                              <  2:1760707/148‑1 (MQ=255)
                                                                                  gCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCg                                                              <  2:3518044/148‑1 (MQ=255)
                                                                                  gCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCg                                                              <  2:1147663/148‑1 (MQ=255)
                                                                                    ttGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTACTCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGAc                                                            <  1:607669/148‑1 (MQ=255)
                                                                                       ggtggtTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGAcagc                                                         >  2:666356/1‑148 (MQ=255)
                                                                                           gtTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGAcagcagca                                                     <  1:3435228/148‑1 (MQ=255)
                                                                                            tttttGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTTTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGAcagcagcag                                                    >  2:2535966/1‑148 (MQ=255)
                                                                                             ttttGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGt                                                   >  1:3524958/1‑148 (MQ=255)
                                                                                              tttGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGtt                                                  >  2:1543774/1‑148 (MQ=255)
                                                                                                tGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTtc                                                <  1:2874537/148‑1 (MQ=255)
                                                                                                   cACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTtcatc                                             >  2:2520454/1‑148 (MQ=255)
                                                                                                    aCGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCTAAACCAATGCCTAAAGAGAGGTTAAAGCCGAAAGCAGCAGTTtcatca                                            >  2:1058338/1‑148 (MQ=255)
                                                                                                     cgcgCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCatcaa                                           >  2:10260/1‑148 (MQ=255)
                                                                                                     cgcgCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCatcaa                                           >  2:3344796/1‑148 (MQ=255)
                                                                                                     cgcgCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCatcaa                                           <  1:888911/148‑1 (MQ=255)
                                                                                                           aTCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATcacca                                     <  2:592320/148‑1 (MQ=255)
                                                                                                            tCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATcaccac                                    >  2:18508/1‑148 (MQ=255)
                                                                                                             cAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACg                                   <  2:2456587/148‑1 (MQ=255)
                                                                                                             cAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACg                                   <  2:896122/148‑1 (MQ=255)
                                                                                                             cAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACg                                   >  2:1076234/1‑148 (MQ=255)
                                                                                                             cAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCGTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACg                                   <  2:472402/148‑1 (MQ=255)
                                                                                                              aGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACGa                                  >  1:481372/1‑148 (MQ=255)
                                                                                                               gCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACGAt                                 >  1:853607/1‑148 (MQ=255)
                                                                                                                ctctTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACGATg                                >  1:3113810/1‑148 (MQ=255)
                                                                                                                ctctTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACGATg                                >  1:1343003/1‑148 (MQ=255)
                                                                                                                ctctTTAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTACCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACGATg                                >  1:2125348/1‑148 (MQ=255)
                                                                                                                    ttAATCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACGATGCCAt                            >  2:3376077/1‑148 (MQ=255)
                                                                                                                      aaTCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACGATGCCATGt                          <  2:1514943/148‑1 (MQ=255)
                                                                                                                       aTCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACGATGCCATGtt                         <  2:2736551/148‑1 (MQ=255)
                                                                                                                       aTCGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACGATGCCATGtt                         <  1:325323/148‑1 (MQ=255)
                                                                                                                         cGCCTGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACGATGCCATGTTCa                       <  1:2769312/148‑1 (MQ=255)
                                                                                                                             tGTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACGATGCCATGTTCATCTg                   <  1:2663447/148‑1 (MQ=255)
                                                                                                                              gTAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACGATGCCATGTTCATCTGc                  <  1:705770/148‑1 (MQ=255)
                                                                                                                               tAAGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACGATGCCATGTTCATCTGcc                 >  2:2705314/1‑148 (MQ=255)
                                                                                                                                aaGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACGATGCCATGTTCATCTGccc                >  1:2284368/1‑148 (MQ=255)
                                                                                                                                 aGTGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACGATGCCATGTTCATCTGCCCa               >  1:1167649/1‑148 (MQ=255)
                                                                                                                                   tGCGCTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACGATGCCATGTTCATCTGCCCAGt             <  1:401345/148‑1 (MQ=255)
                                                                                                                                     cgcTTGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACGATGCCATGTTCATCTGCCCAGTCg           <  1:2769067/148‑1 (MQ=255)
                                                                                                                                         tGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACGATGCCATGTTCATCTGCCCAGTCGAGCa       <  1:2616990/148‑1 (MQ=255)
                                                                                                                                         tGCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACGATGCCATGTTCATCTGCCCAGTCGAGCa       >  1:132885/1‑148 (MQ=255)
                                                                                                                                          gCTGAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACGATGCCATGTTCATCTGCCCAGTCGAGCAt      <  2:1390521/148‑1 (MQ=255)
                                                                                                                                             gAGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACGATGCCATGTTCATCTGCCCAGTCGAGCAtctc   >  1:231484/1‑148 (MQ=255)
                                                                                                                                              aGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACGATGCCATGTTCATCTGCCCAGTCGAGCATctct  >  1:1916394/1‑148 (MQ=255)
                                                                                                                                              aGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACGATGCCATGTTCATCTGCCCAGTCGAGCATctct  >  1:1778304/1‑148 (MQ=255)
                                                                                                                                              aGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACGATGCCATGTTCATCTGCCCAGTCGAGCATctct  <  2:2925253/148‑1 (MQ=255)
                                                                                                                                              aGTGTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACGATGCCATGTTCATCTGCCCAGTCGAGCATctct  >  1:117075/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TTGCTTCCGCCAGTGGCGCGAAATATTCCCGTGCACCTTGCGGACGGGTATCCGGTTCGTTGGCAATACTCCACATCACCACGCTTGGGTGGTTTTTGTCACGCGCTATCAGCTCTTTAATCGCCTGTAAGTGCGCTTGCTGAGTTTCCCCGTTGACTGCCTCTTCGCTGTACAGTTCTTTCGGCTTGTTGCCCGCTTCGAAACCAATGCCTAAAGAGAGGTTAAAGCCGACAGCAGCAGTTTCATCAATCACCACGATGCCATGTTCATCTGCCCAGTCGAGCATCTCT  >  CP009273/1689035‑1689324

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: