Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,049,716 A→T 6.4% V39V (GTA→GTT 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,7160AT 6.4% 1052.8 / 6.4 337V39V (GTA→GTTfau5‑formyltetrahydrofolate cyclo‑ligase family protein
Reads supporting (aligned to +/- strand):  ref base A (105/194);  new base T (17/2);  total (141/196)
Fisher's exact test for biased strand distribution p-value = 3.12e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 342 total aligned reads displayed.

CCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGAT  >  CP009273/3049567‑3049865
                                                                                                                                                     |                                                                                                                                                     
      aCATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCctc                                                                                                                                                 <  1:1241747/150‑1 (MQ=255)
       cATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGGATTCctgg                                                                                                                                                >  1:2653284/1‑148 (MQ=255)
       cATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGAGGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCctct                                                                                                                                                <  2:4730410/150‑1 (MQ=255)
          tATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCtt                                                                                                                                             <  2:31515/150‑1 (MQ=255)
           aTCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCttt                                                                                                                                            <  1:204524/150‑1 (MQ=255)
           aTCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCttt                                                                                                                                            <  1:2825456/150‑1 (MQ=255)
           aTCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGTGTTccgcgcgtt                                                                                                                                            >  2:888116/1‑143 (MQ=255)
           aTCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGTGGTCCTCTCgtt                                                                                                                                            >  1:3156242/1‑150 (MQ=255)
            tCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCtttt                                                                                                                                           <  1:5196351/150‑1 (MQ=255)
             cccGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGGATTCCTCTCTTTTg                                                                                                                                          >  2:2983992/1‑150 (MQ=255)
              ccGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGa                                                                                                                                         <  1:4140420/150‑1 (MQ=255)
               cGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGTATTACGCTCTGTTTAt                                                                                                                                        >  1:3625474/1‑150 (MQ=255)
                gACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTTTTACTCTCTTTTGATg                                                                                                                                       >  2:1799675/1‑150 (MQ=255)
                gACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATg                                                                                                                                       <  1:4581407/150‑1 (MQ=255)
                gACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATg                                                                                                                                       <  2:163442/150‑1 (MQ=255)
                gACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATg                                                                                                                                       <  1:729108/150‑1 (MQ=255)
                   aaGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGGATTCCTCTCTTTTGATGGCg                                                                                                                                    >  2:1916384/1‑150 (MQ=255)
                   aaGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGGGTTGCTCTCATTTGATTGCg                                                                                                                                    >  2:4025611/1‑150 (MQ=255)
                      aaaTCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGGATTCCTCTAGTTTGTTTGCGAAc                                                                                                                                 >  1:3302639/1‑150 (MQ=255)
                             cAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGGATTCCTCTCTTTTGATGGCGAACTCGacac                                                                                                                          >  2:2048275/1‑150 (MQ=255)
                                aaTGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCa                                                                                                                       <  2:3397145/150‑1 (MQ=255)
                                aaTGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCa                                                                                                                       <  2:1244528/150‑1 (MQ=255)
                                aaTGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCa                                                                                                                       <  2:4588573/150‑1 (MQ=255)
                                aaTGAGTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCa                                                                                                                       <  2:4316403/150‑1 (MQ=255)
                                 aTGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTTATGGCGTACTCGACACCCAg                                                                                                                      >  1:4289940/1‑150 (MQ=255)
                                  tGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGc                                                                                                                     >  2:2013812/1‑150 (MQ=255)
                                  tGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGGGTTGGGCTCTTTTGTTGGCGCACTCGCCCCCCAGc                                                                                                                     >  1:1783573/1‑150 (MQ=255)
                                   gATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGTATTAAGCTCTTTTTTTTGGGAACACGACACCCCGcc                                                                                                                    >  1:4032961/1‑150 (MQ=255)
                                       cGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGTATTACTATCATTTTTTTGGGAACACGACacccacccccac                                                                                                                >  1:1410545/1‑146 (MQ=255)
                                        ggCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCa                                                                                                               <  1:4612904/150‑1 (MQ=255)
                                          cAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGGGTTGCTCGCTTTTGTTGGCGAACTCGACACCCAGCCACACATa                                                                                                             >  1:618262/1‑150 (MQ=255)
                                           aaCGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGGATTGCTCTCTTTTTTTTGCGTACTCGACCCCCACCCACACAAAg                                                                                                            >  2:1316142/1‑150 (MQ=255)
                                            aCGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGa                                                                                                           <  2:1840174/150‑1 (MQ=255)
                                            aCGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGa                                                                                                           <  1:3523382/150‑1 (MQ=255)
                                               tcgtcgtGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGaaca                                                                                                        <  1:3430454/150‑1 (MQ=255)
                                                  tcgtGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAAct                                                                                                     <  1:1698535/150‑1 (MQ=255)
                                                   cgtGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGGGGTCCTCTCGTTTGATTGCGGACTCGACACCCACCCCCGCATCTAACAActc                                                                                                    >  1:1893428/1‑150 (MQ=255)
                                                    gtGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCACTCTTTTTATGGCGAACTCGACACCCACCCACACAAATaacaacaaa                                                                                                   >  1:3769720/1‑147 (MQ=255)
                                                      gCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTgg                                                                                                 >  2:5033606/1‑150 (MQ=255)
                                                       cGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGgc                                                                                                <  2:3101598/150‑1 (MQ=255)
                                                             cGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCg                                                                                          <  2:3035073/150‑1 (MQ=255)
                                                               ccGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGt                                                                                        <  2:2601565/150‑1 (MQ=255)
                                                                cGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTTTATTCCTCTCTTTTTATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTa                                                                                       <  1:22557/150‑1 (MQ=255)
                                                                cGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTa                                                                                       <  1:2365505/150‑1 (MQ=255)
                                                                cGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGTGTTTCTCACTGTTGTTGGGGCACACGACACCCACCCACACAAAGAAAAACACAGGGGCGCCCGTa                                                                                       >  1:4399069/1‑150 (MQ=255)
                                                                cGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGGTTTGCTGTCGTTTTTTGGGGAACTCGACAACCAGCCACACAAATAAAAAATAAGGCGCGCCCGTa                                                                                       >  1:32073/1‑150 (MQ=255)
                                                                 ggAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTaa                                                                                      <  2:2202310/150‑1 (MQ=255)
                                                                  gAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGTTGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAg                                                                                     <  1:4445179/150‑1 (MQ=255)
                                                                  gAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGGCTTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGGAAg                                                                                     >  2:2526385/1‑150 (MQ=255)
                                                                   aacaacAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAAAAAGACACCCACCCACTCAAATAAAAACTCTGGCGCGCCGGTAAgc                                                                                    >  2:3137286/1‑150 (MQ=255)
                                                                        cagcagGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGACGCGCCGGTAAGCGCGta                                                                               <  2:2819282/150‑1 (MQ=255)
                                                                           cagGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGCACTCGACAACCAGCCCCTCATAGAAAAACACTGGCGCGCCGGTAAGCGCGtataa                                                                            >  2:648127/1‑149 (MQ=255)
                                                                            agGAAATGGGTCAACAAGCCGCTACCCGGAGGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTGTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATc                                                                           <  2:3670657/150‑1 (MQ=255)
                                                                             ggATATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCt                                                                          <  1:4922303/150‑1 (MQ=255)
                                                                             ggAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCt                                                                          <  2:4473452/150‑1 (MQ=255)
                                                                             ggAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCt                                                                          <  1:807711/150‑1 (MQ=255)
                                                                             ggAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCt                                                                          <  2:4089588/150‑1 (MQ=255)
                                                                              gAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATTCGGTCGCTGTATTCCTCTCTTTTGATGGCTAACTCTACACCCAGCCACTCATAGAACAACTCTGGCGCTCCGGTAAGCTCGTATATCtt                                                                         <  2:4176758/150‑1 (MQ=255)
                                                                              gAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCtt                                                                         <  1:2512410/150‑1 (MQ=255)
                                                                              gAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCtt                                                                         <  2:107133/150‑1 (MQ=255)
                                                                                aaTGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGTATTCCTCTCGTTTGATGGCGAACTCGACACCCAGCCACACATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTcc                                                                       >  2:2856305/1‑150 (MQ=255)
                                                                                aaTGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGGATTCCGCGCTTTTGATGGCGAACACGACACCCAGCCACACCAAGAACAACACTGGCGCGCCGGTAAGCGCGTATATCTTcc                                                                       >  2:232981/1‑150 (MQ=255)
                                                                                 aTGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCa                                                                      <  2:4141227/150‑1 (MQ=255)
                                                                                  tGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAg                                                                     >  2:5216909/1‑150 (MQ=255)
                                                                                    ggTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGGGTTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCCCTCCTAGAACAACGCTGGCGCGCCGGTAAGCGCGTAAATCTTCCAGtt                                                                   >  1:4770563/1‑150 (MQ=255)
                                                                                       caacaaGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATTGCGAACTCGACACCCAGCCACTCATAGTACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTg                                                                >  2:668020/1‑150 (MQ=255)
                                                                                        aacaaGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCGCTTTTGATGGCGAACTCGACACCCAGCCACGCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTGCAGGTTTGc                                                               >  2:1520288/1‑150 (MQ=255)
                                                                                           aaGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATc                                                            <  1:848700/150‑1 (MQ=255)
                                                                                           aaGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACACATAGAACAACTCTGGCGCGCCGGGAAGCGCGGATATCTTCCAGTTTTGCATc                                                            >  2:3283736/1‑150 (MQ=255)
                                                                                              ccGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCt                                                         <  1:1011615/150‑1 (MQ=255)
                                                                                               cGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTGCTCTATTTTGATGGCGAACTCGACACCCACCCACACATAGAACAACACTGGCGCGCCGGTGATCGCGTATATAATACAGTTTTTCATaccgt                                                        >  1:4591250/1‑148 (MQ=255)
                                                                                                gCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCttt                                                       <  1:2643819/150‑1 (MQ=255)
                                                                                                gCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCttt                                                       <  1:3937563/150‑1 (MQ=255)
                                                                                                      cGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGcc                                                 >  2:1709088/1‑150 (MQ=255)
                                                                                                       ggATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCg                                                <  1:2167028/150‑1 (MQ=255)
                                                                                                         atgatgACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGt                                              <  1:1185691/150‑1 (MQ=255)
                                                                                                         atgatgACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGt                                              <  2:3613750/150‑1 (MQ=255)
                                                                                                           gatgaCTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTaa                                            <  1:5039252/150‑1 (MQ=255)
                                                                                                           gatgaCTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTaa                                            <  2:3578731/150‑1 (MQ=255)
                                                                                                           gatgaCTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGTTTTCCTCTCTTTTGATGGCGGACGCGACACCCACCCACACATATAAAAACACTGGCGCGCCGGTAAGCGCGGATATCTTACACTTTTGCATACCGGTAGTGGCGGGaa                                            >  1:1622297/1‑150 (MQ=255)
                                                                                                                 ttATCCCCCGGTGGTGATGGCACATACGGTCGCTGGGTTCCTCTCTTTTGATGGGGAACTCGACACCCAGCCACTCATAGAACCACGCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGGAATTTGCt                                      >  2:2427625/1‑150 (MQ=255)
                                                                                                                  tATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTg                                     <  2:3961947/150‑1 (MQ=255)
                                                                                                                  tATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTg                                     <  2:213460/150‑1 (MQ=255)
                                                                                                                        ccGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTg                               <  2:5134210/150‑1 (MQ=255)
                                                                                                                        ccGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTg                               <  1:2877414/150‑1 (MQ=255)
                                                                                                                          ggtggtGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACACATAGAACAACTCGGGCGCGCCGGGAAGCGCGGATATCTTCCAGGTTTGCATCCCATTAGTGCCGGTAATTTTCTGTTCCTGaa                             >  2:2927786/1‑150 (MQ=255)
                                                                                                                           gtggtgATGGCACATACGGTCGCGGGTTTCCTCTCATTTTTTGGGGCACTCGACACCCACCCACACAAAGAAAAAAACGGGCGCGCCGGTAAGCGCGTATATATTCCACGTTTGGATCCCGTTAGTTCCGGGACTTTGCTGTTGCTGAAt                            >  2:3511585/1‑150 (MQ=255)
                                                                                                                              gtgATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTAc                         <  1:561026/150‑1 (MQ=255)
                                                                                                                              gtgATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTAc                         <  1:1480327/150‑1 (MQ=255)
                                                                                                                                    gCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATcca                   <  2:71518/150‑2 (MQ=255)
                                                                                                                                          aCGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGc             <  2:4254761/150‑1 (MQ=255)
                                                                                                                                          aCGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGc             <  2:7986/150‑1 (MQ=255)
                                                                                                                                           cGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCg            <  2:4115677/150‑1 (MQ=255)
                                                                                                                                               cGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACt        >  1:4210259/1‑150 (MQ=255)
                                                                                                                                                gCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTg       >  2:3091681/1‑150 (MQ=255)
                                                                                                                                                  tGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGt     >  2:110757/1‑150 (MQ=255)
                                                                                                                                                  tGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGt     <  2:2653284/150‑1 (MQ=255)
                                                                                                                                                   gTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGtg    <  2:3790811/150‑1 (MQ=255)
                                                                                                                                                    tATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTCATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGtga   >  1:3868847/1‑150 (MQ=255)
                                                                                                                                                     aTTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGtgat  >  1:4848173/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCACTGACATTATCCCGACAAGAAATCCGCAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGAT  >  CP009273/3049567‑3049865

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: