Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,630,770 G→T 91.5% T4T (ACG→ACT yjtD → putative methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,630,7700GT91.5% 556.9 / 46.3 211T4T (ACG→ACTyjtDputative methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (6/12);  new base T (98/95);  total (104/107)
Fisher's exact test for biased strand distribution p-value = 2.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.50e-01

Only 100 of 217 total aligned reads displayed.

AATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGAT  >  CP009273/4630625‑4630917
                                                                                                                                                 |                                                                                                                                                   
aaTAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGTTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTat                                                                                                                                                   <  1:104188/148‑1 (MQ=255)
aaTAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTat                                                                                                                                                   <  2:2295656/148‑1 (MQ=255)
 aTAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTatt                                                                                                                                                  >  2:2246276/1‑148 (MQ=255)
   aGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACtattat                                                                                                                                                >  1:2154670/1‑148 (MQ=255)
   aGTTTAGCCTTTTTGAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACGattat                                                                                                                                                <  2:2963721/148‑1 (MQ=255)
     tttAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTc                                                                                                                                              >  1:3141277/1‑148 (MQ=255)
     tttAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTc                                                                                                                                              >  2:310067/1‑148 (MQ=255)
      ttAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCt                                                                                                                                             >  1:2234366/1‑148 (MQ=255)
           cTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCg                                                                                                                                        >  2:2283928/1‑148 (MQ=255)
           cTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCg                                                                                                                                        >  2:1047100/1‑148 (MQ=255)
                  aaGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGc                                                                                                                                 <  1:145159/148‑1 (MQ=255)
                    gCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCa                                                                                                                               >  2:2445794/1‑148 (MQ=255)
                     cTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGAATTTGTACATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTGATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCag                                                                                                                              >  1:1511344/1‑148 (MQ=255)
                     cTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCag                                                                                                                              <  1:1897659/148‑1 (MQ=255)
                       aaGTAAAGGGCTTTTTCTGCGACTTACGTTAAGAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCagag                                                                                                                            >  2:380125/1‑148 (MQ=255)
                         gTAAAGGGCTTTTTCTGCGACTTACGTTAAGAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCa                                                                                                                          <  1:915962/148‑1 (MQ=255)
                           aaaGGGCTTTTTCTGCGACTTACGTTAAGAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGa                                                                                                                        >  2:3031995/1‑148 (MQ=255)
                             cgggCTTTTTCTGCGACTTACGTTAAGAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGaaa                                                                                                                      <  2:3119176/147‑1 (MQ=255)
                             aGGGCTTTTTCTGCGACTTACGTTAAGAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGaaa                                                                                                                      <  2:1207052/148‑1 (MQ=255)
                             aGGGCTTTTTCTGCGACTTACGTTAAGAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGaaa                                                                                                                      <  1:199776/148‑1 (MQ=255)
                              gggCTTTTTCTGCGACTTACGTTAAGAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGaaaa                                                                                                                     >  2:1538776/1‑148 (MQ=255)
                              gggCTTTTTCTGCGACTTACGTTAAGAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGaaaa                                                                                                                     >  2:3257423/1‑148 (MQ=255)
                              gggCTTTTTCTGCGACTTACGTTAAGAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGaaaa                                                                                                                     >  2:3179694/1‑148 (MQ=255)
                                gCTTTTTCTGCGACTTACGTTAAGAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAata                                                                                                                   <  2:2687138/148‑1 (MQ=255)
                                gCTTTTTCTGCGACTTACGTTAAGAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAata                                                                                                                   <  1:1180154/148‑1 (MQ=255)
                                gCTTTTTCTGCGACTTACGTTAAGAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAata                                                                                                                   >  2:705722/1‑148 (MQ=255)
                                 cTTTTTCTGCGACTTACGTTAAGAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAatat                                                                                                                  <  1:117152/148‑1 (MQ=255)
                                 cTTTTTCTGCGACTTACGTTAAGAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAatat                                                                                                                  <  1:3296467/148‑1 (MQ=255)
                                   ttttCTGCGACTTACGTTAAGAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTg                                                                                                                >  1:2810966/1‑148 (MQ=255)
                                     ttCTGCGACTTACGTTAAGAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTggg                                                                                                              <  1:451909/148‑1 (MQ=255)
                                      tCTGCGACTTACGTTAAGAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTgggg                                                                                                             >  2:1422586/1‑148 (MQ=255)
                                      tCTGCGACTTACGTTAAGAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTgggg                                                                                                             >  2:849719/1‑148 (MQ=255)
                                       cTGCGACTTACGTTAAGAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGc                                                                                                            <  1:1088691/148‑1 (MQ=255)
                                       cTGCGACTTACGTTAAGAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGc                                                                                                            <  2:2237994/148‑1 (MQ=255)
                                         gCGACTTACGTTAAGAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAg                                                                                                          >  1:345863/1‑148 (MQ=255)
                                            aCTTACGTTAAGAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAgcgg                                                                                                       <  2:2521518/148‑1 (MQ=255)
                                             cTTACGTTAAGAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAgcggc                                                                                                      >  1:1121394/1‑148 (MQ=255)
                                              ttACGTTAAGAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAgcggcg                                                                                                     <  2:101802/148‑1 (MQ=255)
                                                 cGTTAAGAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCgg                                                                                                  <  1:2441293/148‑1 (MQ=255)
                                                    tAAGAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCa                                                                                               >  2:48941/1‑148 (MQ=255)
                                                       gAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATg                                                                                            <  1:2988911/148‑1 (MQ=255)
                                                       gAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATg                                                                                            <  1:2989069/148‑1 (MQ=255)
                                                        aaTTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGa                                                                                           <  1:1852524/148‑1 (MQ=255)
                                                         aTTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGaa                                                                                          >  1:2728118/1‑148 (MQ=255)
                                                          tttGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGaaa                                                                                         <  2:3111188/148‑1 (MQ=255)
                                                            tGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAAc                                                                                       >  1:305499/1‑148 (MQ=255)
                                                             gTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACg                                                                                      <  2:3010974/148‑1 (MQ=255)
                                                             gTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACg                                                                                      >  1:2038734/1‑148 (MQ=255)
                                                             gTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACg                                                                                      <  2:941105/148‑1 (MQ=255)
                                                               aaaTTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGAt                                                                                    >  2:462954/1‑148 (MQ=255)
                                                                 aTTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATgg                                                                                  <  1:1551449/148‑1 (MQ=255)
                                                                  ttCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATggg                                                                                 >  2:1826472/1‑148 (MQ=255)
                                                                      cACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGttt                                                                             <  1:1005421/148‑1 (MQ=255)
                                                                      cACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGttt                                                                             >  1:3115863/1‑148 (MQ=255)
                                                                           cgTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGa                                                                        >  1:138528/1‑148 (MQ=255)
                                                                             taataaGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATc                                                                      <  2:1092693/148‑1 (MQ=255)
                                                                             taataaGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATc                                                                      <  2:3287364/148‑1 (MQ=255)
                                                                             taataaGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATc                                                                      >  1:1239301/1‑148 (MQ=255)
                                                                              aataaGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCt                                                                     >  2:598184/1‑148 (MQ=255)
                                                                               ataaGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTg                                                                    <  2:2648185/148‑1 (MQ=255)
                                                                                   gTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGa                                                                <  2:811014/148‑1 (MQ=255)
                                                                                    ttGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGAt                                                               >  2:2018808/1‑148 (MQ=255)
                                                                                    ttGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGAt                                                               >  2:1359978/1‑148 (MQ=255)
                                                                                         aGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCg                                                          <  1:527976/148‑1 (MQ=255)
                                                                                           tGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGAt                                                        <  1:2748399/148‑1 (MQ=255)
                                                                                             aTCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAg                                                      >  2:408410/1‑148 (MQ=255)
                                                                                              tCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGt                                                     <  2:1775671/148‑1 (MQ=255)
                                                                                              tCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGt                                                     <  2:2270642/148‑1 (MQ=255)
                                                                                                aCCCGGTTCGCGGTTATTTGATCAAGACGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCa                                                   <  1:339238/148‑1 (MQ=255)
                                                                                                  ccGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCAgg                                                 <  1:60636/148‑1 (MQ=255)
                                                                                                  ccGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCAgg                                                 >  2:871153/1‑148 (MQ=255)
                                                                                                    ggTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCAGGca                                               <  2:669559/148‑1 (MQ=255)
                                                                                                       tCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCAGGCacac                                            <  1:642896/148‑1 (MQ=255)
                                                                                                       tCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCAGGCacac                                            >  1:2292819/1‑148 (MQ=255)
                                                                                                          cgGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTg                                         >  2:1970439/1‑148 (MQ=255)
                                                                                                             ttATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGag                                      >  1:1002168/1‑148 (MQ=255)
                                                                                                              tATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGagc                                     >  2:80402/1‑148 (MQ=255)
                                                                                                              tATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGagc                                     <  2:1862263/148‑1 (MQ=255)
                                                                                                              tATTTGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGagc                                     <  1:2494030/148‑1 (MQ=255)
                                                                                                                 ttGATCAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGagccag                                  >  1:2711753/1‑148 (MQ=255)
                                                                                                                      cAAGAAGAGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGAGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCAcc                             >  1:2783922/1‑148 (MQ=255)
                                                                                                                         gaagaGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGc                          >  2:1316323/1‑148 (MQ=255)
                                                                                                                         gaagaGTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGc                          >  2:1799607/1‑148 (MQ=255)
                                                                                                                             agTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTggg                      >  2:1147763/1‑148 (MQ=255)
                                                                                                                              gTGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGt                     <  2:3141277/148‑1 (MQ=255)
                                                                                                                              gTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGt                     <  2:1761294/148‑1 (MQ=255)
                                                                                                                               tGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTc                    >  2:3120209/1‑148 (MQ=255)
                                                                                                                               tGGCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGAGGGGGTTTAGCGCTCTGCGGAGTGTCGGTAGGCTGGCGGCCCTGGAGCCCGCCACGCGCgggggc                    >  2:2810562/1‑146 (MQ=255)
                                                                                                                                ggCAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCg                   >  2:680578/1‑148 (MQ=255)
                                                                                                                                  cAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGca                 >  1:1625645/1‑148 (MQ=255)
                                                                                                                                  cAATATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGca                 >  1:333387/1‑148 (MQ=255)
                                                                                                                                   aaTATGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGcac                <  1:595685/148‑1 (MQ=255)
                                                                                                                                    atatGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGgaca               >  1:2609172/1‑148 (MQ=255)
                                                                                                                                     tatGCGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACAt              <  1:2392671/148‑1 (MQ=255)
                                                                                                                                         cGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGAt          <  1:2954529/148‑1 (MQ=255)
                                                                                                                                         cGTATAACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGAt          <  2:1344022/148‑1 (MQ=255)
                                                                                                                                           tataACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCt        >  2:1568583/1‑148 (MQ=255)
                                                                                                                                           tataACTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCt        <  1:498435/148‑1 (MQ=255)
                                                                                                                                              aaCTATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGt     >  1:1869371/1‑148 (MQ=255)
                                                                                                                                                 gATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGat  <  2:1199747/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
AATAGTTTAGCCTTTTTTAAGCTAAGTAAAGGGCTTTTTCTGCGACTTACGTTAAGAATTTGTAAATTCGCACCGCGTAATAAGTTGACAGTGATCACCCGGTTCGCGGTTATTTGATCAAGAAGAGTGGCAATATGCGTATAACGATTATTCTGGTCGCACCCGCCAGAGCAGAAAATATTGGGGCAGCGGCGCGGGCAATGAAAACGATGGGGTTTAGCGATCTGCGGATTGTCGATAGTCAGGCACACCTGGAGCCAGCCACCCGCTGGGTCGCACATGGATCTGGTGAT  >  CP009273/4630625‑4630917

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: