Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,049,740 C→A 87.0% 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,7400CA87.0% 416.8 / 58.4 187L47L (CTC→CTAfau5‑formyltetrahydrofolate cyclo‑ligase family protein
Reads supporting (aligned to +/- strand):  ref base C (6/19);  new base A (62/100);  total (68/119)
Fisher's exact test for biased strand distribution p-value = 1.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 188 total aligned reads displayed.

CCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCA  >  CP009273/3049593‑3049883
                                                                                                                                                   |                                                                                                                                               
   cAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGac                                                                                                                                              >  2:2620810/1‑148 (MQ=255)
       aTGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACAccc                                                                                                                                          >  1:224260/1‑148 (MQ=255)
         gATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTTTTCCTCTCTTTTGAATGCGAACTAGACACCCAg                                                                                                                                        >  2:431903/1‑148 (MQ=255)
         gATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTTTTGGCGAACAAGAAACCCaa                                                                                                                                        >  1:238951/1‑147 (MQ=255)
         gATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAg                                                                                                                                        <  1:1956822/148‑1 (MQ=255)
         gATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGGATTGCTCTCTTTTGATGGCGAACTAGACACCCAg                                                                                                                                        >  2:1520265/1‑148 (MQ=255)
         gATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGAAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCGTTTTATGGCGTACTCGACACCCAg                                                                                                                                        >  2:939414/1‑148 (MQ=255)
          aTTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGTTGGCGAACTAGACACCCAGc                                                                                                                                       >  2:2709475/1‑148 (MQ=255)
            tCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGTTGGCGAACTAGACACCCAGCCa                                                                                                                                     >  2:2317326/1‑148 (MQ=255)
                cAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCa                                                                                                                                 <  2:1947396/148‑1 (MQ=255)
                cAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGTATTGCTCTCTTTTGGTGGCGAACTAGACACCCAGCCACTCa                                                                                                                                 >  2:2367708/1‑148 (MQ=255)
                 aaCGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCAt                                                                                                                                <  2:416143/148‑1 (MQ=255)
                     tcgtcgtGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTTATGGCGAACTAGACACCCAGCCACTCATAGaa                                                                                                                            >  2:2904429/1‑148 (MQ=255)
                      cgtcgtGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACCCCCAGCCACTCATAGaac                                                                                                                           >  2:1073048/1‑148 (MQ=255)
                        tcgtGCGTTAACGCCGGAACAACAGCAGAAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGaacaa                                                                                                                         <  2:2854899/148‑1 (MQ=255)
                        tcgtGCGTTAACGCCGGAACAACAGCAGAAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCAGAGaacaa                                                                                                                         >  2:2259609/1‑148 (MQ=255)
                            gCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGaacaaaact                                                                                                                     >  2:2375067/1‑144 (MQ=255)
                             cGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTCTTGCTCTCTTTTGTTGGCGAACTAGACACCCAGCCACTCATAGAACAACGCTg                                                                                                                    >  1:2851385/1‑148 (MQ=255)
                             cGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATGACGCAACTGTGCAGGCGTACAAGACACCCAGCCACCCCTACAAAAACACTg                                                                                                                    >  1:2851403/1‑148 (MQ=255)
                              gTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATCCGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTgg                                                                                                                   <  1:74408/148‑1 (MQ=255)
                              gTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTgg                                                                                                                   <  2:3011304/148‑1 (MQ=255)
                                tAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGgcg                                                                                                                 >  2:1028102/1‑148 (MQ=255)
                                   cGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGgcgcgc                                                                                                              >  2:2379339/1‑148 (MQ=255)
                                            aacaGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACCCCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAgcg                                                                                                     >  2:260721/1‑148 (MQ=255)
                                            aacaGCAGAAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAgcg                                                                                                     <  2:2978484/148‑1 (MQ=255)
                                              cagcagGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGcgcg                                                                                                   <  1:1498001/148‑1 (MQ=255)
                                              cagcagGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCACCCACTCATAGAACAACTCGGGCGCGCCGGTAAGcgcg                                                                                                   >  2:1340290/1‑148 (MQ=255)
                                               agcagGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGt                                                                                                  <  1:1964781/148‑1 (MQ=255)
                                               agcagGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGt                                                                                                  <  1:184674/148‑1 (MQ=255)
                                                gcagGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTTTTGCTGTCTTTTGATGGCGAACTAGAACCCCCCCCACACAAAGAACAACACAGGCGCGCCGGTAAGcgcgga                                                                                                 >  1:57729/1‑146 (MQ=255)
                                                 cagGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGTATTCCACTATTTTTATGGGGAACTATACACCCACCCACACAAAGAAAAAATCTGGCGCGCCGGGAAGCGCGgat                                                                                                >  1:2637260/1‑148 (MQ=255)
                                                 cagAAAATGGGTCAACACGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGtat                                                                                                <  2:1920021/148‑1 (MQ=255)
                                                   ggAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTatat                                                                                              <  2:433754/148‑1 (MQ=255)
                                                      aaTGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCtt                                                                                           <  2:349823/148‑1 (MQ=255)
                                                      aaTGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCGCATTTGATGGCGAACTAGACAACCAACCACTCATATAACAAATCAGGCGCGCCGGTAATCCCGtataaatt                                                                                           >  1:2220630/1‑144 (MQ=255)
                                                       aTGGTTCAACAAGCCTCTACCCGGATGATGAATTATCCCCCGGTGGTGATGGAACATACTGTCGATTTATTCCTCTCTTTTGTTGGCGAAATAGACACCCAGCAACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTc                                                                                          <  1:786776/148‑1 (MQ=255)
                                                       aTGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTc                                                                                          <  2:3173513/148‑1 (MQ=255)
                                                       aTGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTc                                                                                          <  1:1882286/148‑1 (MQ=255)
                                                          ggTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTGCTCTCTTTTTATGGCGAACTAGACACCCAGCCACACATAGAACAACACTGGCGCGCCGGTAAGCGCGTATATATTCCAg                                                                                       >  1:2434266/1‑148 (MQ=255)
                                                            tCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTTATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGtt                                                                                     <  1:3049245/148‑1 (MQ=255)
                                                            tCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGtt                                                                                     <  2:3034319/148‑1 (MQ=255)
                                                               acaaGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGTTGGGGAACTAGACACCCAGCCACTCATAGAAAAACTCTGGCGCGCCGGTAAGCGCGTATATATTCCAGTTTTg                                                                                  >  2:606848/1‑148 (MQ=255)
                                                               acaaGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGTTTTCCTCTCTTTTGATGGCGAACTAGACACCCACCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTg                                                                                  >  2:2934284/1‑148 (MQ=255)
                                                                caaGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGGATATCTTCCAGTTTTGc                                                                                 >  1:2186619/1‑148 (MQ=255)
                                                                caaGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGTATTGCTCTCTTTTGTTGGCGAACTAGACACCCAGCCACGCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGc                                                                                 >  1:2752124/1‑148 (MQ=255)
                                                                  aGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCAt                                                                               <  1:785410/148‑1 (MQ=255)
                                                                  aGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCAt                                                                               <  1:2251027/148‑1 (MQ=255)
                                                                  aGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGTTTTCCTCGCGTTTTTTTGCGAACTAGACACCCACCCCCACAAATAAAAAAAATGGCGCGCCGGGAAGCGCGGATATATTACAGtttttgat                                                                               >  1:952660/1‑144 (MQ=255)
                                                                  aGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGGATTCCTCTCTTTTGATTGCGAACTAGACACCCACCCACTCATAGAACAACACTGGCGCGCCGGGAAGCGCGGATATCTTCCAGTTTTGCAt                                                                               >  1:2199801/1‑148 (MQ=255)
                                                                   gCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATc                                                                              <  2:148422/148‑1 (MQ=255)
                                                                   gCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCGTTTTTTGGCGAACTAGACACCCAGCCACTCATAGAACAAAGCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATc                                                                              >  1:540854/1‑148 (MQ=255)
                                                                    ccGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATcc                                                                             <  1:1579182/148‑1 (MQ=255)
                                                                      gCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCt                                                                           <  1:1028102/148‑1 (MQ=255)
                                                                      gCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCt                                                                           <  1:2948631/148‑1 (MQ=255)
                                                                      gCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCt                                                                           <  2:1213481/148‑1 (MQ=255)
                                                                        tACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCttt                                                                         <  2:2207161/148‑1 (MQ=255)
                                                                        tACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCttt                                                                         <  2:816662/148‑1 (MQ=255)
                                                                              gatgatGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGGATTGCTCTCTTTTTATTGCGAACAAGAAACCCCCCCACACAAAGAACAACACTGGCGCGCCGGGAAGCGCCTATATCTTCCACTTTTTCATCCCATTatttcc                                                                   >  1:699700/1‑143 (MQ=255)
                                                                                tgatgaCTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCgg                                                                 <  2:2533581/148‑1 (MQ=255)
                                                                                  atgaCTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGTATTCCTCTCTTTTGATTGCGAACTAGACACCCAGCCACTCATAGAACAACTATGGCGCGCCGGTAAGCGCGTATATCTTACAGTTTTGCATCCCTTTAGTGCCGGtc                                                               >  2:1033577/1‑147 (MQ=255)
                                                                                   tgaCTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCTTAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTaa                                                              <  2:2011863/148‑1 (MQ=255)
                                                                                   tgaCTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTaa                                                              <  2:1405276/148‑1 (MQ=255)
                                                                                    gaCTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAAt                                                             <  1:3078283/148‑1 (MQ=255)
                                                                                         aTCCCCCGGTGGTGATGGCACATACGGTCGCTGTTTTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACACATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTGCAGTTTTTCATACCTTTAGTGCCGGTAATTTGCt                                                        >  2:744135/1‑148 (MQ=255)
                                                                                              ccGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTcc                                                   <  1:215492/148‑1 (MQ=255)
                                                                                              ccGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTcc                                                   <  1:215346/148‑1 (MQ=255)
                                                                                              ccGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGANCAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTcc                                                   <  1:1683090/148‑1 (MQ=255)
                                                                                              ccGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTcc                                                   <  2:3128205/148‑1 (MQ=255)
                                                                                                 gtggtgATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGa                                                <  2:1122734/148‑1 (MQ=255)
                                                                                                   ggtgATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAAt                                              <  1:2121249/148‑1 (MQ=255)
                                                                                                       aTGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTAcc                                          <  1:2298234/148‑1 (MQ=255)
                                                                                                        tGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCa                                         <  1:689377/148‑1 (MQ=255)
                                                                                                         ggCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCAt                                        <  1:2378385/148‑1 (MQ=255)
                                                                                                         ggCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCAt                                        <  1:902345/148‑1 (MQ=255)
                                                                                                             caTACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACAAGACACCCAGCCACACAAAGAACAAAACTGGCGCGCCGGTAAGCCCGTATATCTTACAGTTTTTCAACCCTTTAGTGCCGGTAATTTGATGTTCCTTAATTAACATCCGc                                    >  2:3188529/1‑148 (MQ=255)
                                                                                                             caTACGGTCGCGGGATTGCACTCTTTTGATTGCGAACTAGACACCCAGCCACGCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATACCGTTAGTGCCGGTAAATTTCTGTTCCTGAAATACCATCCGc                                    >  2:1368626/1‑148 (MQ=255)
                                                                                                                aCGGTCGCTGTATTCCTCTCTTTTGATGGCTAACTAGCCACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCaaa                                 <  1:2673087/148‑1 (MQ=255)
                                                                                                                 cGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAg                                <  1:2795481/148‑1 (MQ=255)
                                                                                                                   gTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCAGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTTCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCg                              >  2:1461470/1‑148 (MQ=255)
                                                                                                                      gCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAAc                           <  1:2557575/148‑1 (MQ=255)
                                                                                                                       cTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACt                          <  1:1051641/148‑1 (MQ=255)
                                                                                                                        tGTATTCCTCTCTTTTGAGGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTg                         <  2:2243054/148‑1 (MQ=255)
                                                                                                                          tATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGt                       >  1:400151/1‑148 (MQ=255)
                                                                                                                          tATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGt                       >  1:2834645/1‑148 (MQ=255)
                                                                                                                          tATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGt                       <  2:535473/148‑1 (MQ=255)
                                                                                                                              ccTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGtgatg                   <  1:2239425/148‑1 (MQ=255)
                                                                                                                               ctctctTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGtgatga                  <  2:884647/148‑1 (MQ=255)
                                                                                                                                 ctctTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAAc                <  2:3012124/148‑1 (MQ=255)
                                                                                                                                  tctTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACa               >  2:333160/1‑148 (MQ=255)
                                                                                                                                    ttttGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAgg             <  1:2259603/148‑1 (MQ=255)
                                                                                                                                     tttGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGt            <  1:373270/148‑1 (MQ=255)
                                                                                                                                     tttGATGGCGAACGCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGt            <  1:3062924/148‑1 (MQ=255)
                                                                                                                                          tGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAg       <  1:1334147/148‑1 (MQ=255)
                                                                                                                                            gCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGAt     <  1:1172797/148‑1 (MQ=255)
                                                                                                                                              gAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATcc   <  2:3068346/148‑1 (MQ=255)
                                                                                                                                               aaCTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCa  <  2:2509788/148‑1 (MQ=255)
                                                                                                                                               aaCTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCa  <  1:2620810/148‑1 (MQ=255)
                                                                                                                                               aaCTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCa  <  1:2347392/148‑1 (MQ=255)
                                                                                                                                               aaCTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCa  <  2:2509744/148‑1 (MQ=255)
                                                                                                                                               aaCTAGACACCCAGCCACTCATAGAACAACTCGGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCa  <  2:3258500/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
CCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCA  >  CP009273/3049593‑3049883

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