Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,049,740 C→A 100% L47L (CTC→CTA fau → 5‑formyltetrahydrofolate cyclo‑ligase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,049,7400CA96.5% 476.6 / 6.9 170L47L (CTC→CTAfau5‑formyltetrahydrofolate cyclo‑ligase family protein
Reads supporting (aligned to +/- strand):  ref base C (4/3);  new base A (77/86);  total (81/89)
Fisher's exact test for biased strand distribution p-value = 7.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 170 total aligned reads displayed.

CAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCA  >  CP009273/3049596‑3049883
                                                                                                                                                 |                                                                                                                                               
cAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGac                                                                                                                                              <  1:558951/148‑1 (MQ=255)
 aaaaTGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGaca                                                                                                                                             <  1:687394/148‑1 (MQ=255)
  aaaTGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTGCGCTCTTTTGATGGCGAACGAGacac                                                                                                                                            >  1:2453881/1‑148 (MQ=255)
  aaaTGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGacac                                                                                                                                            <  2:2670484/148‑1 (MQ=255)
   aaTGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACAcc                                                                                                                                           <  1:488589/148‑1 (MQ=255)
   aaTGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACAcc                                                                                                                                           <  1:482706/148‑1 (MQ=255)
   aaTGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACAcc                                                                                                                                           <  1:2206523/148‑1 (MQ=255)
   aaTGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGTATTCCTCTCGTTTGATGGCGAACTAGACAcc                                                                                                                                           >  1:3234202/1‑148 (MQ=255)
            gCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTc                                                                                                                                  <  1:1856453/148‑1 (MQ=255)
             cAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGGATTGCGCTCTTTTGTTGGCGACCTAGACAACCACCcacaca                                                                                                                                 >  1:2395873/1‑148 (MQ=255)
              aaCGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCAt                                                                                                                                <  2:218765/148‑1 (MQ=255)
              aaCGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCAt                                                                                                                                <  2:1277587/148‑1 (MQ=255)
               aCGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTGTTCCGCTCTTTTGTTGGCGACCTAGACACCCAGCCCCACATa                                                                                                                               >  1:399737/1‑148 (MQ=255)
               aCGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTAGTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATa                                                                                                                               >  2:2768226/1‑148 (MQ=255)
               aCGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGGATTTCTCTCTTTTTTTTGCGAACTAGACACCCCCCcacacaaa                                                                                                                               >  1:399695/1‑146 (MQ=255)
                cgtcgtcgtGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAg                                                                                                                              >  2:3033016/1‑148 (MQ=255)
                 gtcgtcgtGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGa                                                                                                                             <  1:954629/148‑1 (MQ=255)
                 gtcgtcgtGCGTTAACGCCGGAACAACAGCAGAAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGa                                                                                                                             <  2:2713725/148‑1 (MQ=255)
                   cgtcgtGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTGCTCTCTTTTGATGGCGAACTAGACACCCAGCCACACATAGaac                                                                                                                           >  2:1699084/1‑148 (MQ=255)
                       gtGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAAct                                                                                                                       <  1:1986798/148‑1 (MQ=255)
                       gtGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAAct                                                                                                                       <  2:1004758/148‑1 (MQ=255)
                        tGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAActc                                                                                                                      >  2:1175546/1‑148 (MQ=255)
                          cGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTg                                                                                                                    >  2:2949948/1‑148 (MQ=255)
                          cGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGGATTCCGCTCTTTTGTTGGCGAACTAGACACCCAGCCACTCAGAGAACAACTCTg                                                                                                                    >  2:1178646/1‑148 (MQ=255)
                             tAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGGAGTGCTCGCTTTTGTTGGGGTACTAGACACCCCCCCACACAAAGAACAAAACTGgcg                                                                                                                 >  1:3004308/1‑148 (MQ=255)
                               aCGCCGGAACAACAGCAGAAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGgcgcg                                                                                                               >  2:2588289/1‑148 (MQ=255)
                                cGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGgcgcgc                                                                                                              <  2:1484238/148‑1 (MQ=255)
                                cGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGgcgcgc                                                                                                              <  2:3250655/148‑1 (MQ=255)
                                   cGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCgg                                                                                                           <  1:3082541/148‑1 (MQ=255)
                                     gAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTa                                                                                                         <  1:2605186/148‑1 (MQ=255)
                                        caacaGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAgc                                                                                                      >  2:3358858/1‑148 (MQ=255)
                                        caacaGCAGGAAATGGGTAAACAAGCCGCTACCCGGATTATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTTTTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAgc                                                                                                      <  1:325806/148‑1 (MQ=255)
                                         aacaGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGTATTGCTCTCTTTTTTTTGCGAACTAGTCACCCCGCCACACAAAGAAAAACAATGGCGCGCCGGGAAgcg                                                                                                     >  1:2718072/1‑148 (MQ=255)
                                           cagcagGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTAATCCTCACTTTTGAAGGCCAAATAGAAACCCAACCACACAAATAACAACTCTGGCGCGCCCGTAAGcgcg                                                                                                   >  1:3549806/1‑148 (MQ=255)
                                              cagGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGtat                                                                                                <  1:2889059/148‑1 (MQ=255)
                                                ggAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTatat                                                                                              >  2:2320984/1‑148 (MQ=255)
                                                 gAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTACTCTCTTTTTATGGCGAACTAGAAACCCAGCCACACATAGAACAACACTGGCGCGCCGGTAAGCGCGGatata                                                                                             >  1:1389046/1‑147 (MQ=255)
                                                    aTGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTc                                                                                          <  2:479303/148‑1 (MQ=255)
                                                    aTGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTc                                                                                          <  2:481148/148‑1 (MQ=255)
                                                     tGGGTCATCCAGCCGCTACCCGGATGATGAGTTATCCCCGGGTGGTGTTGGCACATAGGTTCGTTGTTTTCATCTCTTTTGTTGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTcc                                                                                         <  1:3340394/148‑1 (MQ=255)
                                                     tGGGTCAACAATCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGTATTCCTCTAGTGTGAATTGGCAAAACACACCCAGCACCACATATAAAAACTCATGGGGGCCCGTAAGCCCGTATATATTca                                                                                         >  1:2570457/1‑147 (MQ=255)
                                                        gTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTTTTCCTATCTTTTGTTTGCGAACTAGACACCCCGCCACTCATAGAACAACACTGGCGCGCCGGTCAGCGCGGATAAATTCCAGt                                                                                      >  2:312771/1‑148 (MQ=255)
                                                        gTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACAATACACCCAGCCACTCATAGAACAACACTGGCGCGCCGGGAAGCGCGGATATCTTACAgg                                                                                      >  2:312709/1‑147 (MQ=255)
                                                          caacaaGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTGCACTCTTTTGATTGCGAACTAGACACCCAGCCCCTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTAAATCGTCCAGttt                                                                                    >  2:755670/1‑148 (MQ=255)
                                                           aacaaGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGTATTCCGCTCTTTTGATGGCGGACTAGACACCCAGCCCCGCATCGAACAACTCGGGCGCGCCGGTAAGCGCGTATATCTTCCAGgttt                                                                                   >  2:1667245/1‑148 (MQ=255)
                                                           aacaaGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGTATTCCGCTCGTTTGATGGGGGACTAGACACCCCCCCACACAAAGAACAACACAGGCGCGCCGGTAAGCGCGGATATCTTCCAGtttt                                                                                   >  1:2799344/1‑148 (MQ=255)
                                                            acaaGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTg                                                                                  <  1:3131230/148‑1 (MQ=255)
                                                               aGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCAt                                                                               <  1:1737096/148‑1 (MQ=255)
                                                               aGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCAt                                                                               <  1:442388/148‑1 (MQ=255)
                                                                gCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATc                                                                              <  2:2595484/148‑1 (MQ=255)
                                                                gCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATc                                                                              <  2:3251969/148‑1 (MQ=255)
                                                                      aCCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTa                                                                        <  2:886853/148‑1 (MQ=255)
                                                                            atgatgACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGTATTGCTCGCTTTTGTTGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGGGCGCCGGGAAGCGCGGATTTCTTCCAGGTTTGCATCCCTTTTGTGCCg                                                                  >  1:2793900/1‑148 (MQ=255)
                                                                             tgatgaCTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACGCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGGTTTGCATCCCTTTAGTGCCgg                                                                 >  2:1677694/1‑148 (MQ=255)
                                                                              gatgaCTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTTTTACTCTCTTTTTATTGCGAAAAAGAAACCCCCCCCCTCATATAACAACACTTGCGCGCCGGGAAGCGCGGATATATTACAGTTTTGCATCACGTTAGTGGCGGt                                                                >  2:1522031/1‑148 (MQ=255)
                                                                              gatgaCTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGTTTTGCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAAAAACTCTGGCGCGCCGGTAAGCGCGGATATCTTACAGTTTTGCATCCCTTTAGTGGCGGt                                                                >  2:1522147/1‑148 (MQ=255)
                                                                               atgaATTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGAGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTa                                                               <  1:2994822/148‑1 (MQ=255)
                                                                                tgaCTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTaa                                                              <  1:1541893/148‑1 (MQ=255)
                                                                                 gaCTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAAt                                                             <  2:2970414/148‑1 (MQ=255)
                                                                                    ttATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTg                                                          <  2:1941171/148‑1 (MQ=255)
                                                                                      aTCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCt                                                        >  2:2940644/1‑148 (MQ=255)
                                                                                       tCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTg                                                       <  2:2532367/148‑1 (MQ=255)
                                                                                         ccccGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCGTTAGTGCCGGTAATTTGCTGtt                                                     >  2:547989/1‑148 (MQ=255)
                                                                                         ccccGGTGGTGATGGCACATACGGTCGCTGTAGTCCTCTCTTTTGATGGCGAACGAGACACCCAGCCACTCATAGAACAACACTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGtt                                                     <  2:2363455/148‑1 (MQ=255)
                                                                                             ggtggtGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTg                                                 <  2:1212077/148‑1 (MQ=255)
                                                                                               tggtgATGGCACATACGGTCGCTGTATTCCTCTCTTTTGTTGGCGAACTAGACACCCAGCCACTCCTAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTTCAACCCTTTAGTGCCGGTACTTTGCTGTTCCTGaa                                               >  2:801511/1‑148 (MQ=255)
                                                                                               tggtgATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGaa                                               <  1:1070154/148‑1 (MQ=255)
                                                                                               tggtgATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGaa                                               <  1:1070089/148‑1 (MQ=255)
                                                                                                ggtgATGGCACATACGGTCGCTGTTTTGCTCTCTATTGTTGGCGAACTAGACACCCACCCCCACAAACAACAACAGTGGCGCGCCGGGAAGCGCGGATATCTTACAGTTTTGCATCCCTTTAGTGCCGGGAATTTGCTGTTACTGAAt                                              >  2:3327441/1‑148 (MQ=255)
                                                                                                  tgATGGCACATACGGTCGCTGTATTCCTGTCGTTTTATGGCGAACAAGACACCCAGCCACGCATAGAACAACTCAGGCGCGCCGGTAAGCGCGTATATATTACAGTTTTTCATCCCGTTAGTGCCGGGAATTTGCTGTTACGGAATTa                                            >  2:3388461/1‑148 (MQ=255)
                                                                                                  tgATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTa                                            <  1:1182231/148‑1 (MQ=255)
                                                                                                    aTGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTAcc                                          <  2:156596/148‑1 (MQ=255)
                                                                                                     tGGCACATACGGTCGCTGTATTGCTCTGTTTTGTTGGCGAACTAGACACCCACCCACTCATATAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCATTAGTGCCGGTCATTTGCTGTTCCTGAATTAccc                                         >  2:2520894/1‑147 (MQ=255)
                                                                                                     tGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCa                                         <  1:1150552/148‑1 (MQ=255)
                                                                                                      ggCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCAt                                        <  2:3519959/148‑1 (MQ=255)
                                                                                                      ggCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCAt                                        <  2:181634/148‑1 (MQ=255)
                                                                                                       gCACATACGGTCGCTGTATTGCTCTCTTTTTATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGGATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTTATTTGCTGTTCCTGAATTAACATc                                       >  2:3090982/1‑148 (MQ=255)
                                                                                                              cGGTCGCTGTTTTACTCTCTTTTTATGGCGAACAAGAAACCCAGCCACACAAATAAAAAAAATGGCGCGCCGGTAAGCGCGTATATCTTACAGTTTTTCATCCCTTTATTTGCGGTTATTTGATGTTACTGAATTAAAAAACGCAAAg                                >  2:2142749/1‑148 (MQ=255)
                                                                                                              cGGTCGCTGTATTCCTCTCTTTTGATTGCGAACAATACAACCAGCCACACATAAAACAAAACAGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCAACCCTTTAGTTCCGGTAATTTTCTGTTCCTGAAATACCAAACCCAAAg                                >  2:2142785/1‑148 (MQ=255)
                                                                                                              cGGTCGCTGTATTCCTCTCTTTTGATGGCGAACAATACACCCAGCCACTCATAGAACAAAACTGGCGCGCCGGTAAGCGCGTATATATTCCAGTTTTTCATCCCTTTAGTGCCGGTAATTTTCTGTTACTGAATTAACATCCGCAAAg                                >  2:214780/1‑148 (MQ=255)
                                                                                                                    cTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACt                          >  2:1461736/1‑148 (MQ=255)
                                                                                                                    cTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACt                          <  2:1271846/148‑1 (MQ=255)
                                                                                                                    cTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACt                          >  2:1044238/1‑148 (MQ=255)
                                                                                                                      gTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTgg                        <  2:653630/148‑1 (MQ=255)
                                                                                                                      gTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTgg                        <  1:1607076/148‑1 (MQ=255)
                                                                                                                      gTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTgg                        <  1:1607131/148‑1 (MQ=255)
                                                                                                                       tATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGt                       <  2:742850/148‑1 (MQ=255)
                                                                                                                          tCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGtgat                    <  1:862713/148‑1 (MQ=255)
                                                                                                                          tCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGtgat                    >  1:732962/1‑148 (MQ=255)
                                                                                                                            ctctctTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGtgatga                  <  2:3234202/148‑1 (MQ=255)
                                                                                                                            ctctctTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGtgatga                  >  2:1993243/1‑148 (MQ=255)
                                                                                                                             tctctTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGaa                 <  2:3410938/148‑1 (MQ=255)
                                                                                                                             tctctTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGaa                 <  2:1756586/148‑1 (MQ=255)
                                                                                                                              atctTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAAc                <  1:653399/147‑1 (MQ=255)
                                                                                                                               tctTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACa               <  2:3484304/148‑1 (MQ=255)
                                                                                                                                ctTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAg              <  1:2475130/148‑1 (MQ=255)
                                                                                                                                           cGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATc    >  2:2834893/1‑148 (MQ=255)
                                                                                                                                           cGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATc    >  2:2540983/1‑148 (MQ=255)
                                                                                                                                           cGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGATTGAAGATc    <  2:3282400/148‑1 (MQ=255)
                                                                                                                                             aaCTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCa  >  1:2918670/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
CAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCA  >  CP009273/3049596‑3049883

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