Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,184,471 A→T 26.5% I519I (ATT→ATA)  thiC ← phosphomethylpyrimidine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,184,4710AT26.5% 423.2 / 182.8 234I519I (ATT→ATA) thiCphosphomethylpyrimidine synthase
Reads supporting (aligned to +/- strand):  ref base A (73/99);  new base T (18/44);  total (91/143)
Fisher's exact test for biased strand distribution p-value = 6.96e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

Only 100 of 243 total aligned reads displayed.

TCATCGTGATAAGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGAAGTGGTCATAGCC  >  CP009273/4184325‑4184614
                                                                                                                                                  |                                                                                                                                               
 cATCGTGATAAGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTAt                                                                                                                                               <  2:1472214/148‑1 (MQ=255)
 cATCGTGATAAGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTAt                                                                                                                                               >  1:560690/1‑148 (MQ=255)
 cATCGTGATAAGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTAt                                                                                                                                               <  1:909651/148‑1 (MQ=255)
  aTCGTGATAAGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATc                                                                                                                                              <  2:156059/148‑1 (MQ=255)
   tCGTGATAAGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTAGCCAGGTCAGCGGCGTGGGCAGCAATCt                                                                                                                                             >  1:2904725/1‑148 (MQ=255)
    cGTGATAAGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCtt                                                                                                                                            <  2:564612/148‑1 (MQ=255)
    cGTGATAAGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCtt                                                                                                                                            >  1:1569837/1‑148 (MQ=255)
     gTGATAAGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTta                                                                                                                                           <  2:3216992/148‑1 (MQ=255)
     gTGATAAGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTta                                                                                                                                           <  2:2820055/148‑1 (MQ=255)
       gATAAGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTtata                                                                                                                                         >  2:708757/1‑148 (MQ=255)
        aTAAGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCACTCTTATAg                                                                                                                                        >  2:1057311/1‑148 (MQ=255)
        aTAAGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAg                                                                                                                                        <  1:2508539/148‑1 (MQ=255)
        aTAAGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAg                                                                                                                                        >  1:232836/1‑148 (MQ=255)
         tAAGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAgg                                                                                                                                       >  2:1094476/1‑148 (MQ=255)
            gcgcGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGa                                                                                                                                    <  1:2151338/148‑1 (MQ=255)
             cgcgGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGAt                                                                                                                                   >  2:239648/1‑148 (MQ=255)
             cgcgGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGAt                                                                                                                                   <  1:408756/148‑1 (MQ=255)
             cgcgGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGAt                                                                                                                                   <  2:1995237/148‑1 (MQ=255)
              gcgggcggTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATa                                                                                                                                  <  1:199261/148‑1 (MQ=255)
                gggcggTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAg                                                                                                                                <  2:86353/148‑1 (MQ=255)
                   cggTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGccc                                                                                                                             <  2:1466748/148‑1 (MQ=255)
                    ggTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGcccc                                                                                                                            <  1:2539582/148‑1 (MQ=255)
                     gTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCt                                                                                                                           <  1:1190934/148‑1 (MQ=255)
                      tAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTg                                                                                                                          <  2:3398108/148‑1 (MQ=255)
                           gggTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTaa                                                                                                                     <  1:1768828/148‑1 (MQ=255)
                            ggTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAAc                                                                                                                    >  2:2652466/1‑148 (MQ=255)
                               cGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACAtc                                                                                                                 <  1:727532/148‑1 (MQ=255)
                               cGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACAtc                                                                                                                 <  1:33501/148‑1 (MQ=255)
                                gAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACAtct                                                                                                                >  1:860184/1‑148 (MQ=255)
                                 aGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACAtctt                                                                                                               >  1:297165/1‑148 (MQ=255)
                                 aGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACAtctt                                                                                                               <  2:3357959/148‑1 (MQ=255)
                                   ggacAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACAtcttct                                                                                                             <  1:2400711/145‑1 (MQ=255)
                                   ggCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATGAGCCCCTGCTTAACAtcttct                                                                                                             >  1:1010837/1‑148 (MQ=255)
                                   ggCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACAtcttct                                                                                                             <  1:327135/148‑1 (MQ=255)
                                   ggCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACAtcttct                                                                                                             >  1:1594312/1‑148 (MQ=255)
                                      cAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTa                                                                                                          <  2:1540225/148‑1 (MQ=255)
                                       aGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTAt                                                                                                         <  2:297576/148‑1 (MQ=255)
                                            aaaCTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGTCCCTGCTTAACATCTTCTTTATTGGGc                                                                                                    <  1:912275/148‑1 (MQ=255)
                                                  gTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAAGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGAccc                                                                                              <  1:2603260/148‑1 (MQ=255)
                                                    cTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAg                                                                                            <  2:2722253/148‑1 (MQ=255)
                                                     ttCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGa                                                                                           <  2:2586003/148‑1 (MQ=255)
                                                       cccAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTTATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATg                                                                                         >  2:2555813/1‑148 (MQ=255)
                                                       cccAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATg                                                                                         >  1:1642405/1‑148 (MQ=255)
                                                       cccAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATg                                                                                         <  2:483911/148‑1 (MQ=255)
                                                        ccAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGc                                                                                        <  2:2793547/148‑1 (MQ=255)
                                                           gCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGctct                                                                                     <  2:493237/148‑1 (MQ=255)
                                                           gCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGctct                                                                                     >  2:1033148/1‑148 (MQ=255)
                                                             gAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCttt                                                                                   <  1:451151/148‑1 (MQ=255)
                                                              aaaTTCGACGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATATGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCtttt                                                                                  <  2:3030202/148‑1 (MQ=255)
                                                              aaaTTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATATGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCtttt                                                                                  <  2:3410033/148‑1 (MQ=255)
                                                              aaaTTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCtttt                                                                                  >  2:3038254/1‑148 (MQ=255)
                                                                  tCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCg                                                                              <  2:268125/148‑1 (MQ=255)
                                                                       gCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTAGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACg                                                                         <  2:2904725/148‑1 (MQ=255)
                                                                          ggCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTaa                                                                      >  2:3285166/1‑148 (MQ=255)
                                                                             tttCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACag                                                                   <  2:373084/148‑1 (MQ=255)
                                                                              ttCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACaga                                                                  >  1:41553/1‑148 (MQ=255)
                                                                               tCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACagag                                                                 <  1:2353028/148‑1 (MQ=255)
                                                                               tCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACagag                                                                 <  1:3039694/148‑1 (MQ=255)
                                                                                  aCATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCAt                                                              <  1:2570855/148‑1 (MQ=255)
                                                                                     tGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATcgc                                                           >  2:1641989/1‑148 (MQ=255)
                                                                                      ggCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATcgcg                                                          <  1:1355150/148‑1 (MQ=255)
                                                                                      ggCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATcgcg                                                          >  1:1285654/1‑148 (MQ=255)
                                                                                         gTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAg                                                       <  2:3463821/148‑1 (MQ=255)
                                                                                           tATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGcc                                                     >  2:541258/1‑148 (MQ=255)
                                                                                                cgAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAAcc                                                <  2:1928987/148‑1 (MQ=255)
                                                                                                cgAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAAcc                                                <  2:3232741/148‑1 (MQ=255)
                                                                                                   aTTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGc                                             <  2:2545744/148‑1 (MQ=255)
                                                                                                   aTTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGc                                             <  1:1378609/148‑1 (MQ=255)
                                                                                                   aTTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGc                                             >  1:1417162/1‑148 (MQ=255)
                                                                                                      tGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCaa                                          <  1:708757/148‑1 (MQ=255)
                                                                                                       gcgcgcCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAat                                         <  1:3404450/148‑1 (MQ=255)
                                                                                                         gcgcCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAatca                                       <  1:2964146/148‑1 (MQ=255)
                                                                                                          cgcCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAatcat                                      <  2:2262368/148‑1 (MQ=255)
                                                                                                          cgcCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAatcat                                      >  2:591466/1‑148 (MQ=255)
                                                                                                            cccGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATcg                                    <  2:1569837/148‑1 (MQ=255)
                                                                                                               ggATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATcgccg                                 >  2:787072/1‑148 (MQ=255)
                                                                                                                gATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATcgccgc                                >  1:217295/1‑148 (MQ=255)
                                                                                                                   gCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCAcc                             >  2:13895/1‑148 (MQ=255)
                                                                                                                    cccTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCa                            <  2:422521/148‑1 (MQ=255)
                                                                                                                      cTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAAt                          <  1:1546750/148‑1 (MQ=255)
                                                                                                                       tttCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATc                         >  2:2079909/1‑148 (MQ=255)
                                                                                                                        ttCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATcc                        <  1:645836/148‑1 (MQ=255)
                                                                                                                         tCGCCAGGTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATccc                       <  1:1099670/148‑1 (MQ=255)
                                                                                                                              aGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACg                  <  1:67006/148‑1 (MQ=255)
                                                                                                                               ggTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGt                 <  1:3146318/148‑1 (MQ=255)
                                                                                                                               ggTCAGCGGCGTGGGCAGCTATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGt                 >  1:2526825/1‑148 (MQ=255)
                                                                                                                                gTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTg                <  1:786707/148‑1 (MQ=255)
                                                                                                                                  cAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGaa              >  1:3467835/1‑148 (MQ=255)
                                                                                                                                  cAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGAGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGaa              <  2:3023483/148‑1 (MQ=255)
                                                                                                                                   aGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGGGAAg             >  2:343603/1‑148 (MQ=255)
                                                                                                                                     cggcgTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCCAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGAAGTg           >  2:2198315/1‑148 (MQ=255)
                                                                                                                                     cggcgTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGAAGTg           >  2:294740/1‑148 (MQ=255)
                                                                                                                                     cggcgTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGAAGTg           >  2:1239601/1‑148 (MQ=255)
                                                                                                                                          tGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGAAGTGGTCAt      <  1:833349/148‑1 (MQ=255)
                                                                                                                                            ggCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGAAGTGGTCATAg    >  2:1466090/1‑148 (MQ=255)
                                                                                                                                             gcagctATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGAAGTGGTCATAGc   <  2:1787487/148‑1 (MQ=255)
                                                                                                                                              cagcaATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGAAGTGGTCATAGcc  >  1:1970912/1‑148 (MQ=255)
                                                                                                                                              cagcaATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGAAGTGGTCATAGcc  <  1:2813563/148‑1 (MQ=255)
                                                                                                                                              cagcaATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGAAGTGGTCATAGcc  >  1:2618280/1‑148 (MQ=255)
                                                                                                                                              cagcaATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGAAGTGGTCATAGcc  <  1:1057311/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
TCATCGTGATAAGCGCGGGCGGTAAACGGGTCGAGGGCCAGATTAAACTGGTCTTCCCAGCGAAATTCGAAGCGGGCTTTCGACATGGCGTTATCGCGAATTTGCGCGCCCGGATGCCCTTTCGCCAGGTCAGCGGCGTGGGCAGCAATCTTATAGGTGATAAGCCCCTGCTTAACATCTTCTTTATTGGGCAGACCCAGATGCTCTTTTGGCGTTACGTAACAGAGCATCGCGCAGCCAAACCAGCCAATCATCGCCGCACCAATCCCCGACGTGAAGTGGTCATAGCC  >  CP009273/4184325‑4184614

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