Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,950,431 T→A 7.7% N605Y (AAT→TAT)  ptrA ← protease III

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,950,4310TA 7.7% 508.8 / 26.2 169N605Y (AAT→TAT) ptrAprotease III
Reads supporting (aligned to +/- strand):  ref base T (89/67);  new base A (6/7);  total (95/74)
Fisher's exact test for biased strand distribution p-value = 5.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 170 total aligned reads displayed.

GCTTTACCCTTTTCTGCGGAATCCATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATT  >  CP009273/2950286‑2950578
                                                                                                                                                 |                                                                                                                                                   
gCTTTACCCTTTTCTGCGGAATCCATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTaa                                                                                                                                                   >  1:501646/1‑148 (MQ=255)
 cTTTACCCTTTTCTGCGGAATCCATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAAc                                                                                                                                                  <  2:2481251/148‑1 (MQ=255)
 cTTTACCCTTTTCTGCGGAATCCATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATAAAc                                                                                                                                                  >  2:2965706/1‑148 (MQ=255)
  tttACCCTTTTCTGCGGAATCCATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAAcc                                                                                                                                                 >  1:3514130/1‑148 (MQ=255)
     aCCCTTTTCTGCGGAATCCATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATa                                                                                                                                              >  2:79041/1‑148 (MQ=255)
        cTTTTCTGCGGAATCCATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAgg                                                                                                                                           >  2:730355/1‑148 (MQ=255)
         ttttCTGCGGAATCCATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGc                                                                                                                                          >  2:1435488/1‑148 (MQ=255)
          tttCTGCGGAATCCATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGcc                                                                                                                                         >  2:585682/1‑148 (MQ=255)
          tttCTGCGGAATCCATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGcc                                                                                                                                         >  2:1720398/1‑148 (MQ=255)
           ttCTGCGGAATCCATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCg                                                                                                                                        >  2:1365797/1‑148 (MQ=255)
            tCTGCGGAATCCATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCgt                                                                                                                                       <  2:3116067/148‑1 (MQ=255)
               gCGGAATCCATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCgttgt                                                                                                                                    >  1:2752600/1‑148 (MQ=255)
                cGGAATCCATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCgttgtt                                                                                                                                   <  2:2761468/148‑1 (MQ=255)
                cGGAATCCATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCgttgtt                                                                                                                                   >  2:2907144/1‑148 (MQ=255)
                     tCCATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAATGCCGTTGTTAGCGt                                                                                                                              >  2:1028986/1‑148 (MQ=255)
                        atcatcTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTtgg                                                                                                                           >  1:3260216/1‑148 (MQ=255)
                        atcatcTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTtgg                                                                                                                           <  2:3438413/148‑1 (MQ=255)
                        atcatcTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTtgg                                                                                                                           <  2:3447786/148‑1 (MQ=255)
                        atcatcTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTtgg                                                                                                                           <  2:3448643/148‑1 (MQ=255)
                         tcatcTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTtggt                                                                                                                          >  1:3119726/1‑148 (MQ=255)
                          catcTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTtggtg                                                                                                                         >  2:3206270/1‑148 (MQ=255)
                                gTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGaaaaa                                                                                                                   <  2:1664561/148‑1 (MQ=255)
                                  tataCCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACt                                                                                                                 >  1:2790632/1‑148 (MQ=255)
                                  tataCCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGATGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACt                                                                                                                 >  2:2804695/1‑148 (MQ=255)
                                    taCCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTa                                                                                                               >  2:348078/1‑148 (MQ=255)
                                     aCCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTAt                                                                                                              <  1:2832560/148‑1 (MQ=255)
                                       cAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGc                                                                                                            <  1:2619516/148‑1 (MQ=255)
                                            cTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGccacca                                                                                                       >  1:3528532/1‑148 (MQ=255)
                                             ttCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCaa                                                                                                      <  1:597959/148‑1 (MQ=255)
                                             ttCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCaa                                                                                                      <  1:2508957/148‑1 (MQ=255)
                                                 ccTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGa                                                                                                  >  1:3162349/1‑148 (MQ=255)
                                                   tGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACg                                                                                                >  2:2954415/1‑148 (MQ=255)
                                                    gCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGc                                                                                               <  2:62588/148‑1 (MQ=255)
                                                    gCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGc                                                                                               <  1:1696322/148‑1 (MQ=255)
                                                    gCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGc                                                                                               <  1:1748065/148‑1 (MQ=255)
                                                     cTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGcc                                                                                              <  1:779875/148‑1 (MQ=255)
                                                       cAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTg                                                                                            >  1:2368294/1‑148 (MQ=255)
                                                        aaGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATAAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTgg                                                                                           <  1:2965706/148‑1 (MQ=255)
                                                         aGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGt                                                                                          >  2:2671330/1‑148 (MQ=255)
                                                         aGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATAAACCATCAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGt                                                                                          <  2:967167/148‑1 (MQ=255)
                                                          gCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGtt                                                                                         >  1:326887/1‑148 (MQ=255)
                                                          gCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGtt                                                                                         >  1:326868/1‑148 (MQ=255)
                                                          gCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGtt                                                                                         >  2:872657/1‑148 (MQ=255)
                                                          gCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATAAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGtt                                                                                         >  1:1643885/1‑148 (MQ=255)
                                                           cTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTg                                                                                        >  1:2182/1‑148 (MQ=255)
                                                           cTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTg                                                                                        >  1:467070/1‑148 (MQ=255)
                                                             gATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCt                                                                                      <  1:103557/148‑1 (MQ=255)
                                                             gATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAAGGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCt                                                                                      <  1:2992963/148‑1 (MQ=255)
                                                              aTCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCtt                                                                                     >  1:2400762/1‑148 (MQ=255)
                                                              aTCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCtt                                                                                     >  2:547273/1‑148 (MQ=255)
                                                                 ttCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAAc                                                                                  >  1:2596019/1‑148 (MQ=255)
                                                                    cGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGa                                                                               <  2:2900310/148‑1 (MQ=255)
                                                                     gTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATAAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGATGCCTGGTTGCTTAACTGAt                                                                              <  2:448646/148‑1 (MQ=255)
                                                                      tAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATc                                                                             >  2:138287/1‑148 (MQ=255)
                                                                       aGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATAAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCa                                                                            <  2:2864838/148‑1 (MQ=255)
                                                                        gCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCaa                                                                           >  1:2174723/1‑148 (MQ=255)
                                                                        gCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCaa                                                                           <  2:1767192/148‑1 (MQ=255)
                                                                             ataGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGcc                                                                      <  2:1141397/148‑1 (MQ=255)
                                                                             ataGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGcc                                                                      <  1:894604/148‑1 (MQ=255)
                                                                             ataGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATAAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGcc                                                                      >  1:1383194/1‑148 (MQ=255)
                                                                                 cTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGTCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGcc                                                                  <  1:3513000/148‑1 (MQ=255)
                                                                                 cTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGcc                                                                  >  2:1243144/1‑148 (MQ=255)
                                                                                 cTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGcc                                                                  >  1:2671089/1‑148 (MQ=255)
                                                                                    aaGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTg                                                               <  2:3162349/148‑1 (MQ=255)
                                                                                         ccccTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCgaga                                                          <  2:2599342/148‑1 (MQ=255)
                                                                                          cccTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGAt                                                         >  2:1860428/1‑148 (MQ=255)
                                                                                          cccTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGAt                                                         >  1:1495594/1‑148 (MQ=255)
                                                                                           ccTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATa                                                        >  1:2323313/1‑148 (MQ=255)
                                                                                            cTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATaa                                                       >  1:675021/1‑148 (MQ=255)
                                                                                               gAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCa                                                    <  1:1025068/148‑1 (MQ=255)
                                                                                                  cAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTg                                                 >  1:2598620/1‑148 (MQ=255)
                                                                                                  cAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTg                                                 >  2:1575319/1‑148 (MQ=255)
                                                                                                  cAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTg                                                 <  2:1975404/148‑1 (MQ=255)
                                                                                                    aTGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAg                                               >  2:2801360/1‑148 (MQ=255)
                                                                                                     tGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAgc                                              >  1:3175396/1‑148 (MQ=255)
                                                                                                      gCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAgcg                                             <  1:2955723/148‑1 (MQ=255)
                                                                                                         tGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCaa                                          <  1:626114/148‑1 (MQ=255)
                                                                                                         tGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCaa                                          <  1:1928135/148‑1 (MQ=255)
                                                                                                         tGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCaa                                          >  2:2162004/1‑148 (MQ=255)
                                                                                                           gAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAAc                                        >  1:1241668/1‑148 (MQ=255)
                                                                                                           gAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAAc                                        >  1:3374002/1‑148 (MQ=255)
                                                                                                           gAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAAc                                        <  2:1438609/148‑1 (MQ=255)
                                                                                                            aaCAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACa                                       >  2:1861571/1‑148 (MQ=255)
                                                                                                            aaCAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACa                                       <  2:2573755/148‑1 (MQ=255)
                                                                                                             aCAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACAt                                      >  2:2955595/1‑148 (MQ=255)
                                                                                                              cAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATc                                     <  2:2872889/148‑1 (MQ=255)
                                                                                                              cAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATc                                     <  1:2538835/148‑1 (MQ=255)
                                                                                                                   gCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATACTCATTGAGCGCAAACATCACCTg                                >  1:2327471/1‑148 (MQ=255)
                                                                                                                   gCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTg                                >  1:2327597/1‑148 (MQ=255)
                                                                                                                    cGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGa                               <  1:1096238/148‑1 (MQ=255)
                                                                                                                     ggCAGACGCTGGGTGTAACCATTAGCATAAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGAt                              <  2:2361162/148‑1 (MQ=255)
                                                                                                                      gCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGAtt                             >  2:523000/1‑148 (MQ=255)
                                                                                                                         gACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCg                          <  1:1365797/148‑1 (MQ=255)
                                                                                                                           cGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCggg                        >  2:378362/1‑148 (MQ=255)
                                                                                                                                     aaCCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCAt              >  1:1761232/1‑148 (MQ=255)
                                                                                                                                       ccATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATgg            >  1:549995/1‑148 (MQ=255)
                                                                                                                                          ttAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCtt         <  2:10856/148‑1 (MQ=255)
                                                                                                                                          ttAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGAGAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCtt         >  1:395231/1‑148 (MQ=255)
                                                                                                                                               aTTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGa    <  2:515826/148‑1 (MQ=255)
                                                                                                                                                 tAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGAtt  >  1:126214/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GCTTTACCCTTTTCTGCGGAATCCATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATT  >  CP009273/2950286‑2950578

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: