Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,885,611 G→T 100% R628S (CGT→AGT)  yoaA ← putative ATP‑dependent helicase, DinG family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,885,6110GT96.1% 850.2 / 19.7 282R628S (CGT→AGT) yoaAputative ATP‑dependent helicase, DinG family
Reads supporting (aligned to +/- strand):  ref base G (5/6);  new base T (171/99);  total (177/105)
Fisher's exact test for biased strand distribution p-value = 3.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 282 total aligned reads displayed.

CGTCGTTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACGAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATCAGTCGCCCTACCCCCT  >  CP009273/1885464‑1885757
                                                                                                                                                   |                                                                                                                                                  
   cgtTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCg                                                                                                                                               >  1:1004166/1‑150 (MQ=255)
    gtTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGc                                                                                                                                              >  2:3876195/1‑150 (MQ=255)
     ttCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCa                                                                                                                                             >  1:1917662/1‑150 (MQ=255)
     ttCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCa                                                                                                                                             <  1:220133/150‑1 (MQ=255)
      tCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCAc                                                                                                                                            >  2:52018/1‑150 (MQ=255)
        cacaGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACgg                                                                                                                                          >  2:3310427/1‑150 (MQ=255)
        cacaGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACgg                                                                                                                                          <  2:3563866/150‑1 (MQ=255)
         acaGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACggg                                                                                                                                         <  2:280094/150‑1 (MQ=255)
          caGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGc                                                                                                                                        >  1:1070879/1‑150 (MQ=255)
          caGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGc                                                                                                                                        <  2:1885079/150‑1 (MQ=255)
                      gCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGt                                                                                                                            <  1:3743523/150‑1 (MQ=255)
                        aGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGc                                                                                                                          >  1:2335255/1‑150 (MQ=255)
                        aGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGc                                                                                                                          >  1:3086130/1‑150 (MQ=255)
                        aGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCAAACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGACACTAACCTCACGGGCAATGTCACGGGTGc                                                                                                                          >  1:2808810/1‑150 (MQ=255)
                         ggCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACGAACCGCACGGGCAATGTCACGGGTGCg                                                                                                                         >  2:3131747/1‑150 (MQ=255)
                             tctGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGgc                                                                                                                     <  1:5114011/150‑1 (MQ=255)
                               tGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACGAACCGCACGGGCAATGTCACGGGTGCGTGgcgc                                                                                                                   >  1:2362254/1‑150 (MQ=255)
                                gTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGcgcg                                                                                                                  >  2:1441653/1‑150 (MQ=255)
                                gTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCAGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGcgcg                                                                                                                  <  2:1078232/150‑1 (MQ=255)
                                    cgGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCgcgggcg                                                                                                              <  2:2942867/150‑1 (MQ=255)
                                     ggTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCgcgggcgg                                                                                                             >  2:119715/1‑150 (MQ=255)
                                        aTCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAg                                                                                                          <  1:4412871/150‑1 (MQ=255)
                                        aTCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAg                                                                                                          >  1:3971933/1‑150 (MQ=255)
                                        aTCGAGAGCCAGAATTCGCATGGGTCTTCGTGCTGAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAg                                                                                                          <  1:3524918/150‑1 (MQ=255)
                                          cGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGAc                                                                                                        >  2:3165592/1‑150 (MQ=255)
                                             tAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTgg                                                                                                     <  1:619872/150‑1 (MQ=255)
                                              aGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGc                                                                                                    >  1:2621128/1‑150 (MQ=255)
                                                ccAGAATTCGCATGGGTCTTCGTGCTTAGATCTATAAAAAGGCGCGCATATTACCTTTCTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCga                                                                                                  <  1:3888679/150‑1 (MQ=255)
                                                ccAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCACTGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCga                                                                                                  >  2:2459432/1‑150 (MQ=255)
                                                ccAGAAGTCGCATGGGTCTTCGTGCTTAGATCAATACAAAGGCGCGCATCATACCATACTCCGTAACAAAGTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCga                                                                                                  <  1:174272/150‑1 (MQ=255)
                                                  aGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATAATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCgaga                                                                                                <  1:5046886/150‑1 (MQ=255)
                                                    aaTTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGaaa                                                                                              <  2:4524179/150‑1 (MQ=255)
                                                     aTTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAAc                                                                                             <  1:640914/150‑1 (MQ=255)
                                                     aTTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAAGTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAAc                                                                                             <  1:4932231/150‑1 (MQ=255)
                                                       tCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGt                                                                                           >  1:4138329/1‑150 (MQ=255)
                                                          cATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTcgc                                                                                        <  2:739313/150‑1 (MQ=255)
                                                           aTGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTcgcg                                                                                       >  1:1218352/1‑150 (MQ=255)
                                                             gggTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGcc                                                                                     <  1:2342483/150‑1 (MQ=255)
                                                             gggTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGcc                                                                                     <  1:2892059/150‑1 (MQ=255)
                                                                tCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTa                                                                                  >  1:2992816/1‑150 (MQ=255)
                                                                 cTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTaa                                                                                 <  2:1729089/150‑1 (MQ=255)
                                                                  ttCGTGCTTAGATCAATAAAAAGGCGCTCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAg                                                                                <  1:972785/150‑1 (MQ=255)
                                                                  ttCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAg                                                                                <  2:181378/150‑1 (MQ=255)
                                                                    cGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGa                                                                              >  2:4213946/1‑150 (MQ=255)
                                                                        cTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCa                                                                          >  2:3976346/1‑150 (MQ=255)
                                                                          tAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATc                                                                        >  1:4039634/1‑150 (MQ=255)
                                                                          tAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATc                                                                        >  1:2928480/1‑150 (MQ=255)
                                                                           aGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCa                                                                       >  2:3846589/1‑150 (MQ=255)
                                                                           aGATCAAT‑AAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCAc                                                                      <  1:531044/150‑1 (MQ=255)
                                                                            gATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCAc                                                                      >  1:1489398/1‑150 (MQ=255)
                                                                            gATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCAc                                                                      <  2:240105/150‑1 (MQ=255)
                                                                               cAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAg                                                                   >  2:3246618/1‑150 (MQ=255)
                                                                                aaTAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGc                                                                  <  2:1632379/150‑1 (MQ=255)
                                                                                    aaaaGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGAt                                                              >  1:3203008/1‑150 (MQ=255)
                                                                                     aaaGGCGCGCATCATACCATACTCCATAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGAtt                                                             <  1:383432/150‑1 (MQ=255)
                                                                                      aaGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTg                                                            >  2:4961948/1‑150 (MQ=255)
                                                                                       aGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGt                                                           >  1:4156573/1‑150 (MQ=255)
                                                                                       aGGCGCGCATCATACCAAACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACAAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGt                                                           >  2:4573873/1‑150 (MQ=255)
                                                                                         gcgcgcATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTca                                                         <  2:4379489/150‑1 (MQ=255)
                                                                                           gcgcATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTcaca                                                       <  2:2131118/150‑1 (MQ=255)
                                                                                             gcATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACaaa                                                     >  1:2807997/1‑150 (MQ=255)
                                                                                                 cataccatacTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCAc                                                 >  2:2181974/1‑150 (MQ=255)
                                                                                                 cataccatacTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCAc                                                 >  2:1986937/1‑150 (MQ=255)
                                                                                                 cataccatacTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCAc                                                 >  1:4458497/1‑150 (MQ=255)
                                                                                                 cataccatacTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCAc                                                 >  2:1227325/1‑150 (MQ=255)
                                                                                                  ataccatacTCCGTAACAAATGACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCAcc                                                <  2:4007147/150‑1 (MQ=255)
                                                                                                    accatacTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAg                                              >  1:3582288/1‑150 (MQ=255)
                                                                                                    accatacTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAg                                              <  1:2608009/150‑1 (MQ=255)
                                                                                                    accatacTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAg                                              >  1:4857977/1‑150 (MQ=255)
                                                                                                     ccatacTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGc                                             >  1:4651233/1‑150 (MQ=255)
                                                                                                     ccatacTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGc                                             >  1:2081321/1‑150 (MQ=255)
                                                                                                     ccatacTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGc                                             >  2:3725363/1‑150 (MQ=255)
                                                                                                      catacTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCa                                            >  1:3667995/1‑150 (MQ=255)
                                                                                                      catacTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCa                                            >  2:118414/1‑150 (MQ=255)
                                                                                                      catacTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCa                                            >  2:4936632/1‑150 (MQ=255)
                                                                                                         acTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGc                                         >  2:4806835/1‑150 (MQ=255)
                                                                                                            ccGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCAc                                      <  2:2560394/150‑1 (MQ=255)
                                                                                                                 aCAAATTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGAtcg                                 >  2:13323/1‑150 (MQ=255)
                                                                                                                   aaaTTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGAtcgtc                               >  1:5118527/1‑150 (MQ=255)
                                                                                                                   aaaGTACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGAtcgtc                               <  2:1935379/150‑1 (MQ=255)
                                                                                                                     aTTACCTGGAGGATGGTATCGCAAGGAAACTACCCGCACGTGAAATGTCACGGGTGCGTGGCTCGGTCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACTAATCACCAGCACGCCACGATCGTCgg                             <  2:5077781/150‑1 (MQ=255)
                                                                                                                      ttACCTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGc                            <  2:2521713/150‑1 (MQ=255)
                                                                                                                          cTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACCAATCACCAGCACGCCACGATCGTCGGCGTcg                        >  2:1364884/1‑150 (MQ=255)
                                                                                                                          cTGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTca                        <  1:2372671/150‑2 (MQ=255)
                                                                                                                           tGGAGGATGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCTGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTcgc                       >  2:2387756/1‑150 (MQ=255)
                                                                                                                              aggaTGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGaa                    >  2:4853080/1‑150 (MQ=255)
                                                                                                                              aggaTGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGaa                    >  2:3090677/1‑150 (MQ=255)
                                                                                                                               ggaTGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAAt                   <  1:4765195/150‑1 (MQ=255)
                                                                                                                               ggaTGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAAt                   >  2:1124226/1‑150 (MQ=255)
                                                                                                                                 aTGGTATCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATCa                 <  2:1004166/150‑1 (MQ=255)
                                                                                                                                      aTCGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATCAGTCGc            >  1:3323806/1‑150 (MQ=255)
                                                                                                                                        cGCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATCAGTCGccc          <  2:3635686/150‑1 (MQ=255)
                                                                                                                                         gCAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATCAGTCGCCCt         <  1:4616427/150‑1 (MQ=255)
                                                                                                                                          cAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATCAGTCGCCCTa        >  1:1721142/1‑150 (MQ=255)
                                                                                                                                          cAAGGAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATCAGTCGCCCTa        >  1:1518504/1‑150 (MQ=255)
                                                                                                                                           aaGGAAACGAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATCAGTCGCCCTAc       <  1:748911/150‑1 (MQ=255)
                                                                                                                                              gAAACTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATCAGTCGCCCTAcccc    <  1:3009736/150‑1 (MQ=255)
                                                                                                                                               aaaCTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATCAGTCGCCCTAccccc   >  1:4532623/1‑150 (MQ=255)
                                                                                                                                                aaCTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGTCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATCAGTCGCCCTACCCCCt  <  2:640893/150‑1 (MQ=255)
                                                                                                                                                aaCTAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATCAGTCGCCCTACCCCCt  >  2:770752/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGTCGTTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACGAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGCATCACCAGCCGATTGTCACAAATCACCAGCACGCCACGATCGTCGGCGTCGCGAATCAGTCGCCCTACCCCCT  >  CP009273/1885464‑1885757

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