Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,653,365 C→A 15.7% S347* (TCG→TAG)  ynfE → putative selenate reductase, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,653,3650CA15.7% 539.8 / 57.2 178S347* (TCG→TAG) ynfEputative selenate reductase, periplasmic
Reads supporting (aligned to +/- strand):  ref base C (80/70);  new base A (19/9);  total (99/79)
Fisher's exact test for biased strand distribution p-value = 2.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 179 total aligned reads displayed.

TCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTAT  >  CP009273/1653219‑1653512
                                                                                                                                                  |                                                                                                                                                   
tctcGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCAAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCg                                                                                                                                                    <  2:22828/148‑1 (MQ=255)
 ctcGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGc                                                                                                                                                   >  2:869307/1‑148 (MQ=255)
 ctcGTTGATCAACCTTTTCTCGATAAATACTGAGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGTTCACTATAAAGCGTATATTCATGGTGAAGGTGACGATAAATAAGCGAAAACGCCGCAGTGGGCTTCGc                                                                                                                                                   >  2:547521/1‑148 (MQ=255)
   cGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCaa                                                                                                                                                 >  1:1449668/1‑148 (MQ=255)
    gTTGATCAACCTTTTCTAGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCTAAAACGCCGCAGTGGGCTTCGCaaa                                                                                                                                                <  2:2404359/148‑1 (MQ=255)
     ttGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAAt                                                                                                                                               >  2:1842431/1‑148 (MQ=255)
     ttGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAAt                                                                                                                                               <  2:2001538/148‑1 (MQ=255)
         tCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTAGCAAATTAcc                                                                                                                                           >  1:1171241/1‑148 (MQ=255)
          cAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCg                                                                                                                                          >  1:612064/1‑148 (MQ=255)
          cAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCg                                                                                                                                          >  1:612487/1‑148 (MQ=255)
             ccTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTa                                                                                                                                       >  2:2509886/1‑148 (MQ=255)
             ccTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTAGCAAATTACCGGTa                                                                                                                                       >  1:2539740/1‑148 (MQ=255)
             ccTTTTCTCGATAAATACTGCGTCGGATATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAACCATATATTCTTGGTGAAGGTGACTATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTa                                                                                                                                       <  1:3261864/148‑1 (MQ=255)
               ttttCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAACGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATc                                                                                                                                     <  2:1191098/148‑1 (MQ=255)
               ttttCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATc                                                                                                                                     <  2:2068426/148‑1 (MQ=255)
                     ctaTAAATACTGCGTCGGATATGACTAAAAAACCTTACCCGCAGATGCACCCAAAATTGGTCAATATAACGCCTATTTTCTTGGTGATGGTGACGATAAAACCGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCggag                                                                                                                               <  1:1784086/146‑1 (MQ=255)
                     cGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCggag                                                                                                                               >  2:357727/1‑148 (MQ=255)
                     cGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCggag                                                                                                                               >  1:2814627/1‑148 (MQ=255)
                       aTAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAATGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTAGCAAATTACCGGTATCCCggagga                                                                                                                             >  2:3312482/1‑148 (MQ=255)
                       aTAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCACATTACCGGTATCCCggagga                                                                                                                             >  2:896240/1‑148 (MQ=255)
                       aTAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCggagga                                                                                                                             >  2:1381537/1‑148 (MQ=255)
                       aTAAATACGGCGTCGGATATGACGAAAAAACCTAACCCGCAGATGCACCCAAAAATGGTCAATATAAAGCCTATATTCGTGGTTATGGTGACGATAAAACAGCCAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCggagga                                                                                                                             <  1:1258767/148‑1 (MQ=255)
                        tAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGAc                                                                                                                            >  2:1240418/1‑148 (MQ=255)
                        tAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTAGCAAATTACCGGTATCCCGGAGGAc                                                                                                                            >  2:2120700/1‑148 (MQ=255)
                          aaTACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTAGCAAATTACCGGTATCCCGGAGGACCg                                                                                                                          >  2:1350939/1‑148 (MQ=255)
                            tACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTa                                                                                                                        >  2:3105035/1‑148 (MQ=255)
                             aCTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTAGCAAATTACCGGTATCCCGGAGGACCGTat                                                                                                                       >  1:1553284/1‑148 (MQ=255)
                             aCTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTAGCAAATTACCGGTATCCCGGAGGACCGTat                                                                                                                       >  1:1554265/1‑148 (MQ=255)
                             aCTGCGTCGGCTATGACCAAAAAACCTTACCCCCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTAGCAAATTACCGGTATCCCGGAGGACCGTat                                                                                                                       <  2:721292/148‑1 (MQ=255)
                              cTGCGTCTGCTATTACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAATGCCTATATTGTTGGTGATGGTTACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTAGCCCGGAGGACCGTatc                                                                                                                      <  1:1039727/148‑1 (MQ=255)
                              cTGCGTCGGCTATGACTAAAAAACCTTACCCGCAGATGCACCCAAAAAGGGTCACTATAAAGCCTATATTCTGGGGGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTatc                                                                                                                      <  1:1039716/148‑1 (MQ=255)
                               tGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTGGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTatca                                                                                                                     <  1:525002/148‑1 (MQ=255)
                                   tCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCaa                                                                                                                 >  2:1522334/1‑148 (MQ=255)
                                    cGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCaaa                                                                                                                >  1:23826/1‑148 (MQ=255)
                                        tATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTgg                                                                                                            >  1:2026320/1‑148 (MQ=255)
                                        tATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTgg                                                                                                            >  1:137902/1‑148 (MQ=255)
                                           gACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGgcgc                                                                                                         <  2:2523935/148‑1 (MQ=255)
                                           gACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGgcgc                                                                                                         >  2:424155/1‑148 (MQ=255)
                                             cGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGt                                                                                                       >  2:2613055/1‑148 (MQ=255)
                                                 aaaaCCATACCCGCAGATGCACCCAAAAAAGGTCACTATAAAGACTATATTCTTGGAGACTGTGACGAGAAAACAGCCAAAACGACGCAGTGGGCTTCGCAAATTACCCTTATCCAGGAGGACCGTATCATCAAACTGGCGCCTCaaa                                                                                                   >  1:2332178/1‑148 (MQ=255)
                                                  aaaCCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAAt                                                                                                  >  2:510534/1‑148 (MQ=255)
                                                       ttACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGca                                                                                             <  1:204424/148‑1 (MQ=255)
                                                        tACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGcac                                                                                            >  1:1115145/1‑148 (MQ=255)
                                                         aCCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGcaca                                                                                           <  1:2462374/148‑1 (MQ=255)
                                                         aCCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGcaca                                                                                           >  2:839907/1‑148 (MQ=255)
                                                              cAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAACACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCaaa                                                                                      <  1:3220957/148‑1 (MQ=255)
                                                              cAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCaaa                                                                                      >  2:1135861/1‑148 (MQ=255)
                                                              cAGATGCACCCAAAAATGGTCACTATAAACCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCaaa                                                                                      <  1:1581228/148‑1 (MQ=255)
                                                                gATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAAc                                                                                    <  2:419213/148‑1 (MQ=255)
                                                                 aTGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAAcc                                                                                   <  2:2181034/148‑1 (MQ=255)
                                                                 aTGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAAcc                                                                                   >  1:2867680/1‑148 (MQ=255)
                                                                    cACCCAAAAATGGTCACTATAACGCATATTTTCTTGGTGATGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGc                                                                                <  1:1521239/148‑1 (MQ=255)
                                                                        cAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTtat                                                                            <  2:442750/148‑1 (MQ=255)
                                                                        aaaaaaTGGTCACTATAAAGCCTATATTCTGGGTGACGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTtat                                                                            <  1:2296135/147‑1 (MQ=255)
                                                                          aaaaTGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTatat                                                                          <  2:1569244/148‑1 (MQ=255)
                                                                             aTGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTg                                                                       <  2:1456988/148‑1 (MQ=255)
                                                                              tgttaaCTATAAACCCTATATTCTTGGTGAAGGTGACTATAAAACAGCCAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGc                                                                      <  1:1269899/143‑1 (MQ=255)
                                                                               ggTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGcc                                                                     >  1:1338148/1‑148 (MQ=255)
                                                                                gTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCa                                                                    <  1:493185/148‑1 (MQ=255)
                                                                                 tCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACACCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAg                                                                   <  1:691538/148‑1 (MQ=255)
                                                                                      ataAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTAGCAAATTACCGGTAGCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTg                                                              <  2:2851350/148‑1 (MQ=255)
                                                                                        aaaGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTggg                                                            >  2:1231638/1‑148 (MQ=255)
                                                                                           gCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGc                                                         >  1:1574022/1‑148 (MQ=255)
                                                                                           gCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTAGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGc                                                         >  1:300527/1‑148 (MQ=255)
                                                                                               atatTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCATTGGGCTTAGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCcaca                                                     >  1:871791/1‑148 (MQ=255)
                                                                                                  ttgttGGTGAAGGTGACGATAAAACCGCGAACACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGTGGGCCACAACg                                                  <  1:1957892/145‑1 (MQ=255)
                                                                                                  ttCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTAGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACg                                                  <  2:1340697/148‑1 (MQ=255)
                                                                                                  ttCTTGGTGAAGGTGACGATAAAACAGCCAAAACGCCGCAGTGGGCTTAGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACg                                                  <  2:6142/148‑1 (MQ=255)
                                                                                                      tGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAg                                              >  2:273812/1‑148 (MQ=255)
                                                                                                           aaGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCaaa                                         <  2:1115145/148‑1 (MQ=255)
                                                                                                           aaGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTAGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCaaa                                         >  2:2640220/1‑148 (MQ=255)
                                                                                                             ggTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACg                                       >  1:188359/1‑148 (MQ=255)
                                                                                                              gTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACgg                                      >  2:1666915/1‑148 (MQ=255)
                                                                                                                  cGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAATCTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGaa                                  <  2:191794/148‑1 (MQ=255)
                                                                                                                      aaaaCAGTGAAAACGCCGCAGTGGGCTTAGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAactga                              <  2:1923160/148‑1 (MQ=255)
                                                                                                                      aaaaCAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAactga                              >  1:1487675/1‑148 (MQ=255)
                                                                                                                      aaaaCAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAactga                              >  1:1487553/1‑148 (MQ=255)
                                                                                                                      aaaaCAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAactga                              >  1:558140/1‑148 (MQ=255)
                                                                                                                         aCAGCGAAAACGCCGCAGTGGGCTTAGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAactgactg                           <  2:3139511/148‑1 (MQ=255)
                                                                                                                             cGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCAcg                       >  1:2607452/1‑148 (MQ=255)
                                                                                                                             cGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCAcg                       >  1:2607758/1‑148 (MQ=255)
                                                                                                                             cGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCAcg                       >  1:914480/1‑148 (MQ=255)
                                                                                                                               aaaaCGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCAcgcg                     >  1:844836/1‑148 (MQ=255)
                                                                                                                               aaaaCGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCAcgcg                     >  1:844807/1‑148 (MQ=255)
                                                                                                                               aaaaCGCCGCAGTGGGCTTCGCAAATGACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCAcgcg                     <  2:2441280/148‑1 (MQ=255)
                                                                                                                                aaaCGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACgcgc                    >  2:1628493/1‑148 (MQ=255)
                                                                                                                                aaaCGCCGCAGTGGGCTTAGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACgcgc                    >  1:1969071/1‑148 (MQ=255)
                                                                                                                                        cAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGcgc            >  2:1688814/1‑146 (MQ=255)
                                                                                                                                        cAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTa            >  2:1799605/1‑148 (MQ=255)
                                                                                                                                        cAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTa            <  1:3254428/148‑1 (MQ=255)
                                                                                                                                          gTGGGCTTCGCTAATTACCGGTATCCCGGAGGACCGTATGATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATg          <  1:3027609/148‑1 (MQ=255)
                                                                                                                                          gTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATg          <  1:896391/148‑1 (MQ=255)
                                                                                                                                          gTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATg          <  1:356607/148‑1 (MQ=255)
                                                                                                                                          gTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATg          <  1:3132744/148‑1 (MQ=255)
                                                                                                                                          gTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATg          <  1:3132503/148‑1 (MQ=255)
                                                                                                                                           tGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGt         >  1:2731877/1‑148 (MQ=255)
                                                                                                                                                ttCGCAAATTACCTGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCt    <  2:612485/148‑1 (MQ=255)
                                                                                                                                                ttCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCt    <  2:2894440/148‑1 (MQ=255)
                                                                                                                                                ttCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCt    >  2:93795/1‑148 (MQ=255)
                                                                                                                                                  cGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTAt  <  1:1842431/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAAATTACCGGTATCCCGGAGGACCGTATCATCAAACTGGCGCGTGAAATTGGCACAGCAAAACCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTAT  >  CP009273/1653219‑1653512

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: