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,7400CA98.2% 917.9 / 1.8 290L47L (CTC→CTAfau5‑formyltetrahydrofolate cyclo‑ligase family protein
Reads supporting (aligned to +/- strand):  ref base C (3/4);  new base A (129/153);  total (132/158)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 290 total aligned reads displayed.

CCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCATGA  >  CP009273/3049593‑3049886
                                                                                                                                                   |                                                                                                                                                  
 cGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGac                                                                                                                                                 <  2:251553/150‑1 (MQ=255)
   cAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGacac                                                                                                                                               >  1:965315/1‑150 (MQ=255)
          aTTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCa                                                                                                                                        >  1:4044459/1‑150 (MQ=255)
             cGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTc                                                                                                                                     <  2:1668694/150‑1 (MQ=255)
              ggCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCa                                                                                                                                    >  1:1240057/1‑150 (MQ=255)
                 aaCGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAg                                                                                                                                 <  2:1083915/150‑1 (MQ=255)
                   cgtcgtcgtGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGaa                                                                                                                               >  2:5126055/1‑150 (MQ=255)
                     tcgtcgtGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTTATGTCGAACTAGACACCCAGCCACACAGAGaaca                                                                                                                             >  1:461013/1‑150 (MQ=255)
                     tagtcgtGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGAGGTGATGGCACATACGGTCGCTGTATTCCACTCTTTTGATGGCGAACTATACACCCAGCCACTCATAGaaaa                                                                                                                             >  2:232862/3‑148 (MQ=255)
                       gtcgtGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGaacaac                                                                                                                           >  2:4429070/1‑150 (MQ=255)
                        tcgtGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAAct                                                                                                                          >  2:2886148/1‑150 (MQ=255)
                          gtGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAActct                                                                                                                        <  1:2080945/150‑1 (MQ=255)
                             cGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGgc                                                                                                                     <  2:2451891/150‑1 (MQ=255)
                               ttAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGgcgc                                                                                                                   >  1:324451/1‑150 (MQ=255)
                                tAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGgcgcg                                                                                                                  >  1:731469/1‑150 (MQ=255)
                                   cGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCg                                                                                                               <  2:2154505/150‑1 (MQ=255)
                                     ccGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGt                                                                                                             >  2:4043037/1‑150 (MQ=255)
                                     ccGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGTATTCCTCTCTTTTGATGGCGAACAAGACACCCAGCCACTCAGAGAACAACTCTGGCGCGCCGGt                                                                                                             >  2:546797/1‑150 (MQ=255)
                                      cGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTa                                                                                                            <  1:770896/150‑1 (MQ=255)
                                       ggAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTaa                                                                                                           <  1:5081860/150‑1 (MQ=255)
                                       ggAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTaa                                                                                                           <  2:4485054/150‑1 (MQ=255)
                                         aacaacAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAgc                                                                                                         <  1:3103955/150‑1 (MQ=255)
                                          acaacaGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAgcg                                                                                                        <  1:3818474/150‑1 (MQ=255)
                                           caacaGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAgcgc                                                                                                       <  2:1436455/150‑1 (MQ=255)
                                            aacaGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGcgcg                                                                                                      <  2:3949992/150‑1 (MQ=255)
                                            aacaGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGcgcg                                                                                                      <  2:2030228/150‑1 (MQ=255)
                                            aacaGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGcgcg                                                                                                      >  1:2107770/1‑150 (MQ=255)
                                            aacaGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGACGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGcgcg                                                                                                      <  2:259804/150‑1 (MQ=255)
                                              cagcagGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGta                                                                                                    <  2:2397731/150‑1 (MQ=255)
                                                gcagGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTGTTTTGATTGCGAACAAGAAACCCACCCCCAAAAATAACAACACTGGCGCGCCGGTAAGCGCGtaaa                                                                                                  >  1:920151/1‑148 (MQ=255)
                                                        tGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGACGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAg                                                                                          <  2:4218996/150‑1 (MQ=255)
                                                          ggTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGtt                                                                                        >  1:1579964/1‑150 (MQ=255)
                                                           gTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGttt                                                                                       <  2:953719/150‑1 (MQ=255)
                                                             caacaaGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTg                                                                                     <  2:303794/150‑1 (MQ=255)
                                                             caacaaGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGGATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTg                                                                                     >  2:3074888/1‑150 (MQ=255)
                                                              aacaaGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGc                                                                                    >  2:2845825/1‑150 (MQ=255)
                                                                caaGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATTGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCAt                                                                                  >  1:1019240/1‑150 (MQ=255)
                                                                caaGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCAt                                                                                  <  2:1532762/150‑1 (MQ=255)
                                                                   gCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATccc                                                                               >  1:3690521/1‑150 (MQ=255)
                                                                    ccGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATcccc                                                                              <  2:58126/150‑2 (MQ=255)
                                                                    ccGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCt                                                                              <  2:4743271/150‑1 (MQ=255)
                                                                     cGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGTTGGCGAACTAGACACCCAGCCACTCCTAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCtt                                                                             >  2:3097320/1‑150 (MQ=255)
                                                                     cGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCtt                                                                             <  2:3787619/150‑1 (MQ=255)
                                                                     cGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCtt                                                                             >  1:427176/1‑150 (MQ=255)
                                                                      gCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACCCCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCttt                                                                            >  2:3915418/1‑150 (MQ=255)
                                                                            cGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGcc                                                                      <  1:1890829/150‑1 (MQ=255)
                                                                              gatgatGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCgg                                                                    <  1:612421/150‑1 (MQ=255)
                                                                              gatgatGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGTATTCCTCTCGTTTTATTGCGAACAAGACACCCAGCCACACATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATATTCCAGTTTTGCATCCCGTTAGTGCCgg                                                                    >  1:1836945/1‑150 (MQ=255)
                                                                               atgatgACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGt                                                                   <  1:4056162/150‑1 (MQ=255)
                                                                                tgatgaCTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTa                                                                  <  1:4856801/150‑1 (MQ=255)
                                                                                 tatgaCTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCGTTTGATTGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTaa                                                                 >  1:4942112/2‑150 (MQ=255)
                                                                                 gatgaCTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTaa                                                                 <  2:1631691/150‑1 (MQ=255)
                                                                                   tgaCTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAAtt                                                               >  2:1141440/1‑150 (MQ=255)
                                                                                   tgaCTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAAtt                                                               <  2:4307778/150‑1 (MQ=255)
                                                                                    gaCTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAAttt                                                              <  1:4108623/150‑1 (MQ=255)
                                                                                     aCTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTg                                                             >  2:3843912/1‑150 (MQ=255)
                                                                                        tATCCCCCGGTGGTGATGGCACATACTGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACTCCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGAGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTg                                                          <  1:1920401/150‑1 (MQ=255)
                                                                                        tATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTg                                                          <  1:393894/150‑1 (MQ=255)
                                                                                          tCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTACTCTCTTTTGATTGCGTACTATACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTTCTGtt                                                        >  1:2444596/1‑150 (MQ=255)
                                                                                            ccccGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGTTGGCGAACTAGACACCCAGCCACACATAGAAAAACACTGGCGCGCCGGTAAGCGCGTATATATTCCAGTTTTTCATACCTTTATTGCCGGTAATTTGCTGTTcc                                                      >  1:3506397/1‑150 (MQ=255)
                                                                                               cGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGa                                                   <  2:4823421/150‑1 (MQ=255)
                                                                                               cGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGGTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCGGa                                                   >  2:126407/1‑150 (MQ=255)
                                                                                                ggtggtGATGGCACATACGGTCGCTGTATTACTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGaa                                                  >  2:609457/1‑150 (MQ=255)
                                                                                                     tgATGGCACATACGGTCGCGGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGGAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCGGAATTAcc                                             >  1:4147935/1‑150 (MQ=255)
                                                                                                      gATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACACTGGCGCGCCGGTAAGCGCGTATATCTTACAGTTTTGCATCCCTTTAGTGCCGGTAATTTTCTGTTCCTGAATTACCa                                            >  1:2853238/1‑150 (MQ=255)
                                                                                                        tGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACGCATAGAACAACACTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCGTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATc                                          >  1:646993/1‑150 (MQ=255)
                                                                                                         ggCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATcc                                         <  2:481206/150‑1 (MQ=255)
                                                                                                          gCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCg                                        <  1:1342329/150‑1 (MQ=255)
                                                                                                           cacaTACGGTCGCTGTATTCCTCTCTTTTTATGGCGAACAAGAAACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTACAGTTTTGCAAACCTTTAGTGCCGGTAATATGCTGTTACTGAATTAACAACCGc                                       >  2:4931523/1‑150 (MQ=255)
                                                                                                            acaTACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCa                                      >  2:1993968/1‑150 (MQ=255)
                                                                                                              aTACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCaaa                                    >  2:392781/1‑150 (MQ=255)
                                                                                                               tACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCATTAGTGCCGGTAATTTGCTGTTACTGAATTACCAAACGCAAAg                                   >  1:4241125/1‑150 (MQ=255)
                                                                                                                 cGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTACTTTGCTGTTCCTGAATTACCATCCGCAAAGCg                                 >  2:4104832/1‑150 (MQ=255)
                                                                                                                   gTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGaa                               <  1:2611795/150‑1 (MQ=255)
                                                                                                                    tCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAAc                              <  1:5073493/150‑1 (MQ=255)
                                                                                                                     cGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACt                             <  1:2935803/150‑1 (MQ=255)
                                                                                                                     cGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACt                             >  1:2795192/1‑150 (MQ=255)
                                                                                                                     cGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACt                             <  2:1821601/150‑1 (MQ=255)
                                                                                                                      gCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTg                            <  1:1705431/150‑1 (MQ=255)
                                                                                                                      gCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTg                            <  1:240166/150‑1 (MQ=255)
                                                                                                                       cTGTAATCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATATAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTTCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTgg                           >  2:801850/1‑150 (MQ=255)
                                                                                                                        tGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGt                          <  2:3450384/150‑1 (MQ=255)
                                                                                                                        tGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGt                          <  2:864191/150‑1 (MQ=255)
                                                                                                                         gTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGtg                         <  2:224203/150‑1 (MQ=255)
                                                                                                                         gTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGtg                         >  1:552329/1‑150 (MQ=255)
                                                                                                                           aTTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGtgat                       <  2:3575839/150‑1 (MQ=255)
                                                                                                                            ttCCTCTCTTTTGATGGCGAACTAGACACCCAGACACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGtgatg                      <  1:2532029/150‑1 (MQ=255)
                                                                                                                             tCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGtgatga                     >  2:137597/1‑150 (MQ=255)
                                                                                                                                   ctTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGt               <  1:4524038/150‑1 (MQ=255)
                                                                                                                                    ttttGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGtt              <  1:2957601/150‑1 (MQ=255)
                                                                                                                                    ttttGATGGCGAACTAGACACCCAGCCACTCAGAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGACTTACCATCCGCAAAGCGAACTGGTGATGAACAGGtt              <  1:2670342/150‑1 (MQ=255)
                                                                                                                                        gATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAg          <  2:338442/150‑1 (MQ=255)
                                                                                                                                           ggCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATATTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATc       >  1:2360153/1‑150 (MQ=255)
                                                                                                                                            gCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATcc      >  2:1117991/1‑150 (MQ=255)
                                                                                                                                              gAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTACTTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACGGGTGAGGACCAGGTTGAAGATCCAt    >  2:1870995/1‑150 (MQ=255)
                                                                                                                                              gAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCAt    <  1:3738727/150‑1 (MQ=255)
                                                                                                                                              gAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCAt    <  2:729218/150‑1 (MQ=255)
                                                                                                                                              gAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCAt    >  1:2229027/1‑150 (MQ=255)
                                                                                                                                              gAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCAt    <  1:181045/150‑1 (MQ=255)
                                                                                                                                              gAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCAGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCAt    >  1:1107941/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCATGA  >  CP009273/3049593‑3049886

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: