Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,886,725 G→T 100% I256I (ATC→ATA yoaA ← putative ATP‑dependent helicase, DinG family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,886,7250GT99.5% 673.2 / ‑4.3 203I256I (ATC→ATAyoaAputative ATP‑dependent helicase, DinG family
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (86/116);  total (87/116)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 203 total aligned reads displayed.

GAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCT  >  CP009273/1886580‑1886869
                                                                                                                                                 |                                                                                                                                                
gAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTTACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTAt                                                                                                                                                <  2:2688539/148‑1 (MQ=255)
   tttGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGt                                                                                                                                             >  1:672094/1‑148 (MQ=255)
    ttGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGACAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTg                                                                                                                                            <  1:3417511/148‑1 (MQ=255)
    ttGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTg                                                                                                                                            <  1:1678596/148‑1 (MQ=255)
       cGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATg                                                                                                                                         <  2:2998020/148‑1 (MQ=255)
       cGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATg                                                                                                                                         <  2:2490067/148‑1 (MQ=255)
         gATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTc                                                                                                                                       <  1:908778/148‑1 (MQ=255)
          aTTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCt                                                                                                                                      >  2:308711/1‑148 (MQ=255)
           ttAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCtt                                                                                                                                     <  1:3740/148‑1 (MQ=255)
            tAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCttt                                                                                                                                    <  1:536574/148‑1 (MQ=255)
               cTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGc                                                                                                                                 <  1:2085112/148‑1 (MQ=255)
                      cTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCg                                                                                                                          <  1:1286591/148‑1 (MQ=255)
                      cTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCg                                                                                                                          >  2:2283767/1‑148 (MQ=255)
                        cACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAg                                                                                                                        >  1:1181877/1‑148 (MQ=255)
                         aCGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCATGCGCACTCTGGGCAATACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAAATAGGTGCGGTAGGCTATGGTGATGTCTTTTGACAGGTCGAGc                                                                                                                       >  2:545683/1‑148 (MQ=255)
                           gCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAg                                                                                                                     <  1:2999918/148‑1 (MQ=255)
                           gCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAg                                                                                                                     <  1:904369/148‑1 (MQ=255)
                           gCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAg                                                                                                                     <  2:158930/148‑1 (MQ=255)
                            cAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGt                                                                                                                    <  1:3487440/148‑1 (MQ=255)
                                ttACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTc                                                                                                                >  2:912194/1‑148 (MQ=255)
                                 tACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCg                                                                                                               <  1:1908951/148‑1 (MQ=255)
                                    cgcgATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACt                                                                                                            >  2:2956977/1‑148 (MQ=255)
                                      cgATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTgg                                                                                                          <  2:2256640/148‑1 (MQ=255)
                                       gATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGga                                                                                                         <  2:233417/148‑1 (MQ=255)
                                         tAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGgaga                                                                                                       >  2:796769/1‑148 (MQ=255)
                                         tAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGgaga                                                                                                       >  1:1149349/1‑148 (MQ=255)
                                          aaCCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGagag                                                                                                      <  1:2295661/148‑1 (MQ=255)
                                            ccTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGtg                                                                                                    >  2:2943523/1‑148 (MQ=255)
                                             cTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGtga                                                                                                   <  2:2599963/148‑1 (MQ=255)
                                             cTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGtga                                                                                                   <  2:6999/148‑1 (MQ=255)
                                             cTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGtga                                                                                                   >  1:499462/1‑148 (MQ=255)
                                              tGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGtgac                                                                                                  >  1:2543490/1‑148 (MQ=255)
                                              tGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGtgac                                                                                                  <  1:752984/148‑1 (MQ=255)
                                              tGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGtgac                                                                                                  >  1:718771/1‑148 (MQ=255)
                                                gCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTCATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGtgactg                                                                                                >  1:335397/1‑148 (MQ=255)
                                                gCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGtgactg                                                                                                >  1:2853840/1‑148 (MQ=255)
                                                     ccGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCaa                                                                                           >  2:731607/1‑148 (MQ=255)
                                                     ccGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCaa                                                                                           <  2:1158164/148‑1 (MQ=255)
                                                      cGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCaaa                                                                                          <  2:2341314/148‑1 (MQ=255)
                                                       gAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCaaaa                                                                                         <  1:2795249/148‑1 (MQ=255)
                                                       gAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCaaaa                                                                                         <  1:2793862/148‑1 (MQ=255)
                                                              aGACGAAAATCCTGCGCACTCTGGGCAAGGCGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATActggc                                                                                  <  2:399570/148‑1 (MQ=255)
                                                              aGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATActggc                                                                                  <  2:2805151/148‑1 (MQ=255)
                                                               gACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATActggct                                                                                 <  2:841903/148‑1 (MQ=255)
                                                                   aaaaTCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCa                                                                             >  1:3040264/1‑148 (MQ=255)
                                                                    aaaTCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCaa                                                                            <  2:320780/148‑1 (MQ=255)
                                                                    aaaTCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCaa                                                                            >  2:1177814/1‑148 (MQ=255)
                                                                    aaaTCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCaa                                                                            <  2:1282159/148‑1 (MQ=255)
                                                                      aTCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATg                                                                          >  1:3190454/1‑148 (MQ=255)
                                                                       tCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGt                                                                         <  1:628302/148‑1 (MQ=255)
                                                                        ccTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTc                                                                        >  2:120083/1‑148 (MQ=255)
                                                                        ccTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCACTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTc                                                                        >  2:1628011/1‑148 (MQ=255)
                                                                              cACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAg                                                                  >  2:3359371/1‑148 (MQ=255)
                                                                               aCTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGc                                                                 <  2:1041297/148‑1 (MQ=255)
                                                                                     ggCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTACGGTAGCtggtgg                                                           >  2:1985258/1‑148 (MQ=255)
                                                                                      gCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTggg                                                          >  2:1738628/1‑148 (MQ=255)
                                                                                       cAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGc                                                         <  1:1968314/148‑1 (MQ=255)
                                                                                       cAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGc                                                         <  1:1964850/148‑1 (MQ=255)
                                                                                         aGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCtt                                                       >  2:1478113/1‑148 (MQ=255)
                                                                                          gACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTtc                                                      >  2:488277/1‑148 (MQ=255)
                                                                                            cGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTtcgt                                                    >  1:2639448/1‑148 (MQ=255)
                                                                                            cGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTtcgt                                                    >  1:1633353/1‑148 (MQ=255)
                                                                                                  gcgcACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGa                                              <  1:2735802/148‑1 (MQ=255)
                                                                                                  gcgcACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGa                                              <  2:2412506/148‑1 (MQ=255)
                                                                                                     cACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCa                                           >  2:1223453/1‑148 (MQ=255)
                                                                                                       cTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATg                                         <  2:2125507/148‑1 (MQ=255)
                                                                                                       cTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATg                                         <  2:2853840/148‑1 (MQ=255)
                                                                                                           tGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGt                                     >  2:2824620/1‑148 (MQ=255)
                                                                                                           tGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGt                                     <  1:1063988/148‑1 (MQ=255)
                                                                                                               aCTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGc                                 >  2:1416773/1‑148 (MQ=255)
                                                                                                                ctgctgGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCt                                <  2:672094/148‑1 (MQ=255)
                                                                                                                ctgctgGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCt                                <  1:2880460/148‑1 (MQ=255)
                                                                                                                ctgctgGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCt                                <  1:1372969/148‑1 (MQ=255)
                                                                                                                 tgctgGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCtt                               <  2:748208/148‑1 (MQ=255)
                                                                                                                 tgctgGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCtt                               <  2:2506935/148‑1 (MQ=255)
                                                                                                                  gctgGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTc                              >  2:3294271/1‑148 (MQ=255)
                                                                                                                    tgGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCAGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCg                            >  1:237952/1‑148 (MQ=255)
                                                                                                                      ggTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCggg                          >  2:2642789/1‑148 (MQ=255)
                                                                                                                      ggTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCggg                          >  2:2365115/1‑148 (MQ=255)
                                                                                                                       gtgtCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGa                         <  1:203095/148‑1 (MQ=255)
                                                                                                                        tgtCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGAt                        <  2:2543757/148‑1 (MQ=255)
                                                                                                                         gtCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATc                       <  1:952418/148‑1 (MQ=255)
                                                                                                                           cTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAg                     >  1:445067/1‑148 (MQ=255)
                                                                                                                             tttAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGtt                   <  2:2986409/148‑1 (MQ=255)
                                                                                                                             tttAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGtt                   <  2:2124163/148‑1 (MQ=255)
                                                                                                                              ttAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTc                  <  2:3479343/148‑1 (MQ=255)
                                                                                                                                 aTTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGcc               >  2:2155798/1‑148 (MQ=255)
                                                                                                                                 aTTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGCTCGAGCAGTTGTCGACAGTAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGAGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGATCGcc               >  1:1535748/1‑148 (MQ=255)
                                                                                                                                  ttCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCa              >  1:2064137/1‑148 (MQ=255)
                                                                                                                                      gTGCGGTAGGCTATGGTTATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATc          >  2:560372/1‑148 (MQ=255)
                                                                                                                                      gTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATc          <  2:478690/148‑1 (MQ=255)
                                                                                                                                          ggTTGGCTATGTTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCGGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATAATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCCAct      <  1:308711/148‑1 (MQ=255)
                                                                                                                                          ggTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCCAct      <  1:1142799/148‑1 (MQ=255)
                                                                                                                                           gTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCCActc     >  2:1821713/1‑148 (MQ=255)
                                                                                                                                            tAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCCActct    <  1:2458146/148‑1 (MQ=255)
                                                                                                                                             aGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCCActctc   >  1:1096927/1‑148 (MQ=255)
                                                                                                                                              ggCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCCActctct  <  1:1419676/148‑1 (MQ=255)
                                                                                                                                              ggCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCCActctct  >  2:2780325/1‑148 (MQ=255)
                                                                                                                                              ggCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCCActctct  >  1:2867377/1‑148 (MQ=255)
                                                                                                                                              ggCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCCActctct  <  1:776915/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCT  >  CP009273/1886580‑1886869

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: