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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,950,4310TA30.1% 436.4 / 270.8 302N605Y (AAT→TAT) ptrAprotease III
Reads supporting (aligned to +/- strand):  ref base T (104/107);  new base A (38/53);  total (142/160)
Fisher's exact test for biased strand distribution p-value = 2.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.10e-01

Only 100 of 306 total aligned reads displayed.

AAACGCTTTACCCTTTTCTGCGGAATCCATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTAC  >  CP009273/2950282‑2950580
                                                                                                                                                     |                                                                                                                                                     
aaaCGCTTTACCCTTTTCTGCGGAATCCATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCAtt                                                                                                                                                       >  2:4253504/1‑150 (MQ=255)
   cGCTTTACCCTTTTCTGCGGAATCCATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCGGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAAc                                                                                                                                                    >  2:1031928/1‑150 (MQ=255)
       ttACCCTTTTCTGCGGAATCCATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATAAACCATa                                                                                                                                                <  1:5028233/150‑1 (MQ=255)
            cTTTTCTGCGGAATCCATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGcc                                                                                                                                           >  2:4903674/1‑150 (MQ=255)
                  tGCGGAATCCATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCgttgtt                                                                                                                                     <  2:2495332/150‑1 (MQ=255)
                  tGCGGAATCCATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATAAACCATAAGGCCgttgtt                                                                                                                                     <  2:3118034/150‑1 (MQ=255)
                    cGGAATCCATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAg                                                                                                                                   >  1:3971542/1‑150 (MQ=255)
                    cGGAATCCATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAg                                                                                                                                   >  1:3356973/1‑150 (MQ=255)
                     ggAATCCATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGc                                                                                                                                  <  2:154879/150‑1 (MQ=255)
                        aTCCATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGtt                                                                                                                               >  2:2528282/1‑150 (MQ=255)
                        aTCCATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGtt                                                                                                                               <  2:399894/150‑1 (MQ=255)
                          ccATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTtgg                                                                                                                             >  1:2118519/1‑150 (MQ=255)
                            atcatcTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTtggtg                                                                                                                           >  1:3698630/1‑150 (MQ=255)
                            atcatcTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTtggtg                                                                                                                           >  1:2973264/1‑150 (MQ=255)
                             tcatcTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTtggtgg                                                                                                                          <  1:5295875/150‑1 (MQ=255)
                              catcTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATAAACCATAAGGCCGTTGTTAGCGTTGGTGGa                                                                                                                         >  1:938113/1‑150 (MQ=255)
                               atcTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATAAACCATAAGGCCGTTGTTAGCGTTGGTGGaa                                                                                                                        >  1:4611048/1‑150 (MQ=255)
                                 cTGGTTAGACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATAAACCATAAGGCCGTTGTTAGCGTTGGTGGaaaa                                                                                                                      <  1:4502887/150‑1 (MQ=255)
                                     ttATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACtt                                                                                                                  <  2:2346202/150‑1 (MQ=255)
                                       ataCCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTAt                                                                                                                >  2:912692/1‑150 (MQ=255)
                                       ataCCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATAAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTAt                                                                                                                <  1:3447565/150‑1 (MQ=255)
                                         aCCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGc                                                                                                              <  1:332692/150‑1 (MQ=255)
                                          ccAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGcc                                                                                                             <  1:1961552/150‑1 (MQ=255)
                                          ccAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGcc                                                                                                             <  1:3766732/150‑1 (MQ=255)
                                          ccAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATAAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATTcc                                                                                                             <  2:2598397/150‑1 (MQ=255)
                                          ccAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATAAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGcc                                                                                                             >  2:294349/1‑150 (MQ=255)
                                           cAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATAAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGcca                                                                                                            >  2:2762549/1‑150 (MQ=255)
                                            aGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGccac                                                                                                           >  2:345401/1‑150 (MQ=255)
                                            aGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATAAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGccac                                                                                                           <  1:4185867/150‑1 (MQ=255)
                                             ggACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATAAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGccacc                                                                                                          >  1:2541734/1‑150 (MQ=255)
                                             ggACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATAAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGccacc                                                                                                          <  2:4396736/150‑1 (MQ=255)
                                               aCTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATAAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCaa                                                                                                        >  1:5118248/1‑150 (MQ=255)
                                                 ttCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAAcc                                                                                                      <  2:4592027/150‑1 (MQ=255)
                                                 ttCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAAcc                                                                                                      <  2:954939/150‑1 (MQ=255)
                                                 ttCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATAAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAAcc                                                                                                      <  2:2248736/150‑1 (MQ=255)
                                                  tCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATAAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCg                                                                                                     >  2:4299475/1‑150 (MQ=255)
                                                        gCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATAAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCt                                                                                               <  1:3649508/150‑1 (MQ=255)
                                                         cTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTg                                                                                              <  2:2141212/150‑1 (MQ=255)
                                                         cTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTg                                                                                              >  1:1353493/1‑150 (MQ=255)
                                                         cTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGAAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTg                                                                                              <  1:115492/150‑1 (MQ=255)
                                                          tCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTgg                                                                                             <  2:3171845/150‑1 (MQ=255)
                                                           cAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGt                                                                                            >  2:1773784/1‑150 (MQ=255)
                                                           cAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGt                                                                                            <  2:139294/150‑1 (MQ=255)
                                                               cTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCt                                                                                        <  2:5301500/150‑1 (MQ=255)
                                                               cTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATAAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCt                                                                                        >  2:1899592/1‑150 (MQ=255)
                                                                tGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCtt                                                                                       >  2:2858118/1‑150 (MQ=255)
                                                                 gATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTa                                                                                      >  1:5299620/1‑150 (MQ=255)
                                                                 gATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTa                                                                                      >  1:4972530/1‑150 (MQ=255)
                                                                 gATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATAAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTa                                                                                      >  1:5067456/1‑150 (MQ=255)
                                                                    cTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACt                                                                                   <  1:1712981/150‑1 (MQ=255)
                                                                      tCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATAAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGa                                                                                 <  2:5126112/150‑1 (MQ=255)
                                                                        cGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATAAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATc                                                                               >  2:5224221/1‑150 (MQ=255)
                                                                            gCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAgc                                                                           <  1:4423105/150‑1 (MQ=255)
                                                                            gCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAgc                                                                           <  1:3738461/150‑1 (MQ=255)
                                                                             cGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATAAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAgcg                                                                          <  2:1323822/150‑1 (MQ=255)
                                                                                tataGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCa                                                                       >  1:1220611/1‑150 (MQ=255)
                                                                                tataGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATAAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCa                                                                       <  1:4886214/150‑1 (MQ=255)
                                                                                 ataGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAg                                                                      >  1:793933/1‑150 (MQ=255)
                                                                                  taGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGc                                                                     >  2:2713537/1‑150 (MQ=255)
                                                                                    gCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGccc                                                                   <  2:3174812/150‑1 (MQ=255)
                                                                                     cTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCt                                                                  <  1:3356462/150‑1 (MQ=255)
                                                                                     cTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATAAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCt                                                                  <  2:3846234/150‑1 (MQ=255)
                                                                                       aaaGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGc                                                                >  1:3021368/1‑150 (MQ=255)
                                                                                       aaaGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATAAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGc                                                                <  2:4778674/150‑1 (MQ=255)
                                                                                        aaGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCg                                                               <  1:4296875/150‑1 (MQ=255)
                                                                                            aCCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGAt                                                           <  2:642467/150‑1 (MQ=255)
                                                                                            aCCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATAAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGAt                                                           <  2:1678106/150‑1 (MQ=255)
                                                                                            aCCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATAAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGAt                                                           <  2:392804/150‑1 (MQ=255)
                                                                                               ccTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATAAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAAt                                                        <  2:114936/150‑1 (MQ=255)
                                                                                                   gAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCAtt                                                    >  1:4076833/1‑150 (MQ=255)
                                                                                                   gAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCAtt                                                    >  2:3376761/1‑150 (MQ=255)
                                                                                                       aaTGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATAAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAgc                                                >  2:4961352/1‑150 (MQ=255)
                                                                                                          gCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATAAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCa                                             <  1:3579920/150‑1 (MQ=255)
                                                                                                            cTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCaaa                                           <  1:1016067/150‑1 (MQ=255)
                                                                                                             tGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAAc                                          >  1:2862159/1‑150 (MQ=255)
                                                                                                              ggAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGTCCGTTGTTAGCGTTGGTGGAAAAACTTATGACACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACa                                         >  2:470337/1‑150 (MQ=255)
                                                                                                                 aCAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCa                                      >  1:2340793/1‑150 (MQ=255)
                                                                                                                    gCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCt                                   <  2:3747916/150‑1 (MQ=255)
                                                                                                                    gCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCt                                   <  2:553004/150‑1 (MQ=255)
                                                                                                                       gCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGAt                                <  1:1904782/150‑1 (MQ=255)
                                                                                                                        cGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGAtt                               <  2:4476335/150‑1 (MQ=255)
                                                                                                                         ggCAGACGCTGGGTGTAACCATTAGCATAAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTg                              <  1:1122423/150‑1 (MQ=255)
                                                                                                                         ggCAGACGCTGGGTGTAACCATTAGCATAAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTg                              <  1:1525131/150‑1 (MQ=255)
                                                                                                                         ggCAGACGCTGGGTGTAACCATTAGCATAAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTg                              <  2:679834/150‑1 (MQ=255)
                                                                                                                            aGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCgg                           >  2:3572240/1‑150 (MQ=255)
                                                                                                                             gACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCggg                          <  2:5145014/150‑1 (MQ=255)
                                                                                                                             gACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCggg                          >  1:1241140/1‑150 (MQ=255)
                                                                                                                                 cTGGGTGTAACCATTAGCATAAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCt                      >  1:2720894/1‑150 (MQ=255)
                                                                                                                                    ggTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTc                   >  2:5190299/1‑150 (MQ=255)
                                                                                                                                     gtgtAACCATTAGCATAAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTcc                  <  2:2820480/150‑1 (MQ=255)
                                                                                                                                      tgtAACCATTAGCATAAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCa                 >  1:2062191/1‑150 (MQ=255)
                                                                                                                                        tAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATg               >  2:2013746/1‑150 (MQ=255)
                                                                                                                                        tAACCATTAGCATAAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATg               >  2:1206096/1‑150 (MQ=255)
                                                                                                                                         aaCCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATgg              >  2:2456505/1‑150 (MQ=255)
                                                                                                                                          aCCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGc             >  1:1858004/1‑150 (MQ=255)
                                                                                                                                              ttAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTc         >  2:4638223/1‑150 (MQ=255)
                                                                                                                                              ttAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTc         <  2:4140599/150‑1 (MQ=255)
                                                                                                                                                  cATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGAt     >  1:3516721/1‑150 (MQ=255)
                                                                                                                                                   aTTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGAtt    >  1:589423/1‑150 (MQ=255)
                                                                                                                                                     tAACCATACGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTAc  <  1:4210690/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAACGCTTTACCCTTTTCTGCGGAATCCATCATCTGGTTATACCAGGACTTCGCCTGCTCAAGCTGATCTTCCGTAGCGGTATAGCTAAAGTACCCCTCGAGCAATGCCTGGAACAGCTGCGGCAGACGCTGGGTGTAACCATTAGCATTAACCATAAGGCCGTTGTTAGCGTTGGTGGAAAAACTTATGCCACCAACCGACGCCTGGTTGCTTAACTGATCAAGCGCCAGCCCTGCGAGATAATCATTGAGCGCAAACATCACCTGATTGCGGGCGCTGTCCATGGCTTTCGGATTAC  >  CP009273/2950282‑2950580

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