Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A1 F368 I0 R1 1472 214.4 7116750 95.9% 6824963 147.9

Breseq alignment

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

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

Only 100 of 170 total aligned reads displayed.

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

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCATGAGCC  >  CP009273/3049596‑3049889
                                                                                                                                                |                                                                                                                                                     
CAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTGCGCTCTTTTGATGGCGAACGAGACAC                                                                                                                                                  >  SRR19973067.2470010/1‑150 (MQ=60)
CAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACAC                                                                                                                                                  <  SRR19973067.562609/150‑1 (MQ=60)
CAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACAC                                                                                                                                                  <  SRR19973067.562747/150‑1 (MQ=60)
 AAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACC                                                                                                                                                 <  SRR19973067.692124/150‑1 (MQ=60)
 AAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGTATTCCTCTCGTTTGATGGCGAACTAGACACC                                                                                                                                                 >  SRR19973067.3254493/1‑150 (MQ=60)
   AATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCA                                                                                                                                               <  SRR19973067.486027/150‑1 (MQ=60)
   AATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCA                                                                                                                                               <  SRR19973067.491945/150‑1 (MQ=60)
   AATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCA                                                                                                                                               <  SRR19973067.2221238/150‑1 (MQ=60)
            GCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCAT                                                                                                                                      <  SRR19973067.1869104/150‑1 (MQ=60)
             CAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTGTTCCGCTCTTTTGTTGGCGACCTAGACACCCAGCCCCACATA                                                                                                                                     >  SRR19973067.402526/1‑150 (MQ=60)
               ACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGA                                                                                                                                   <  SRR19973067.2542690/150‑1 (MQ=60)
                 GTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAAC                                                                                                                                 <  SRR19973067.961166/150‑1 (MQ=60)
                    GTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTACTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAAC                                                                                                                              >  SRR19973067.527263/1‑150 (MQ=60)
                       GTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCT                                                                                                                           <  SRR19973067.2000232/150‑1 (MQ=60)
                       GTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCT                                                                                                                           <  SRR19973067.2366090/150‑1 (MQ=60)
                                   CGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTGCTCTCGTTTGTTTGCGAACTAGACACCCccccccacaaataaacaaaatggcgcgccgctc                                                                                                               >  SRR19973067.484334/1‑117 (MQ=60)
                                   CGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTA                                                                                                               <  SRR19973067.3102155/150‑1 (MQ=60)
                                    GGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGTTTTACTCTCTTTTTTTGGGGAACTAGAAACCCAGCCACACATAGAACAACACTGGCGCGCCGGTAA                                                                                                              >  SRR19973067.1730598/1‑150 (MQ=60)
                                     GAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAG                                                                                                             <  SRR19973067.2622117/150‑1 (MQ=60)
                                       ACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGTATTGCTCTCTTTTTTTTGCGAACTAGTCACCCCGCCACACAAAGAAAAACAATGGCGCGCCGGGAAGCG                                                                                                           >  SRR19973067.2735591/1‑150 (MQ=60)
                                        CAACAGCAGGAAATGGGTAAACAAGCCGCTACCCGGATTATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTTTTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGC                                                                                                          <  SRR19973067.328029/150‑1 (MQ=60)
                                         AACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTAATCCTCACTTTTGAAGGCCAAATAGAAACCCAACCACACAAATAACAACTCTGGCGCGCCCGTAAGCGCG                                                                                                         >  SRR19973067.3572196/1‑150 (MQ=60)
                                              CAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATAT                                                                                                    <  SRR19973067.2907525/150‑1 (MQ=60)
                                               AGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTACTCTCTTTTTATGGCGAACTAGAAACCCAGCCACACATAGAACAACACTGGCGCGCCGGTAAGCGCGGATATA                                                                                                   >  SRR19973067.1398488/1‑150 (MQ=60)
                                                     TGGGTCATCCAGCCGCTACCCGGATGATGAGTTATCCCCGGGTGGTGTTGGCACATAGGTTCGTTGTTTTCATCTCTTTTGTTGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAG                                                                                             <  SRR19973067.3361547/150‑1 (MQ=60)
                                                         TCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTT                                                                                         <  SRR19973067.1901524/150‑1 (MQ=60)
                                                         TCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGTATTGCGCTCTTTTGTTGGGGAACGAGAAACCCAGCCACACATAGAAAAAAACTGGCGCGCCGGTAAGCGCGGATATATTACAGTTTT                                                                                         >  SRR19973067.2817358/1‑150 (MQ=60)
                                                         TCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGTATTCCGCTCGTTTGATGGGGGACTAGACACCCCCCCACACAAAGAACAACACAGGCGCGCCGGTAAGCGCGGATATCTTCCAGTTTT                                                                                         >  SRR19973067.2817402/1‑150 (MQ=60)
                                                          CAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATTGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTG                                                                                        <  SRR19973067.3550945/150‑1 (MQ=60)
                                                          CAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGGTTTCGGCTCTTTTTTTTGGGAACGCGACACCCCGCCACACAAAGAACCACACGGGCGCGCCGGTAAGCGCGTATATCGTCCAGTTTTG                                                                                        >  SRR19973067.2458209/1‑150 (MQ=60)
                                                           AACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGC                                                                                       <  SRR19973067.3128673/150‑1 (MQ=60)
                                                            ACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCA                                                                                      <  SRR19973067.3151125/150‑1 (MQ=60)
                                                               AGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCC                                                                                   <  SRR19973067.445460/150‑1 (MQ=60)
                                                               AGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCC                                                                                   <  SRR19973067.1748977/150‑1 (MQ=60)
                                                               AGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGGATTCCTCTCTTTTGTTTGCGGACTAGACACCCACCCACGCAAAGAACAACAATGGCGCGCCGGGAAGCGCGGATATCTTGCAGTTTTGCATCC                                                                                   >  SRR19973067.381874/1‑150 (MQ=60)
                                                               AGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGTTTTGCGCGCTTTTGTTTGCGCACTAGACACCCACCCCCACCTAGAACAAAACAGGCGCGCCGGTAAGCGCGGATATCTTACAGTTTTTCATCC                                                                                   >  SRR19973067.157909/1‑150 (MQ=45)
                                                                  CGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTTTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTT                                                                                <  SRR19973067.1523502/150‑1 (MQ=60)
                                                                    CTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTGCTGTCTTTTTATGGGGAACTAGACACCCAcccccgcaaagaacaacacgggggcgcccgggagcgcggatatattgcagtttttcatcccgttt                                                                              >  SRR19973067.998006/1‑85 (MQ=60)
                                                                    CTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGTTTTGGTCTCTTTTTTTGGGGAACTAGAAACCCCGCCACACAAAGAAAAACACTGGCGCGCCGGGCAGCGCGGCTATCTTTCAGTTTTTGATCCCTTTA                                                                              >  SRR19973067.997970/1‑150 (MQ=36)
                                                                         CGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGTGTTGCTCTGTTTTTATGGGGAACGAGACACCCCCCCACACCTAATACAACACTGGCGCGCCGGGAAGCGCGGATATATTCCAGTTTTTCAAACCCTTAGTGCC                                                                         >  SRR19973067.2448195/1‑150 (MQ=29)
                                                                         CGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGTATTGCGCTCTTTTGATGGCGAACTAGACACCCACCCACACATAGAACAACACAGGCGCGCCGGTAAGCGCGGATATCTTGCAGTTTTGCATCCCTTTAGTGGC                                                                         >  SRR19973067.2619123/1‑150 (MQ=60)
                                                                          GGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCG                                                                        <  SRR19973067.402878/150‑1 (MQ=60)
                                                                          GGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTACTCTCTTTTGTTGGGGAACTAGACACCCACCCACACAAAGAACAACACTGGCGCGCCGGTAAGCGCGGATATATTCCAGTTTTGCATACCGTTAGTGGCG                                                                        >  SRR19973067.1778484/1‑150 (MQ=60)
                                                                          GGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCGGTATTGCTCGCTTTTGTTGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGGGCGCCGGGAAGCGCGGATTTCTTCCAGGTTTGCATCCCTTTTGTGCCG                                                                        >  SRR19973067.2811928/1‑150 (MQ=60)
                                                                               ATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAAT                                                                   <  SRR19973067.3011931/150‑1 (MQ=60)
                                                                               ATGAATTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGAGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAAT                                                                   <  SRR19973067.3013806/150‑1 (MQ=60)
                                                                                TGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATT                                                                  <  SRR19973067.1552352/150‑1 (MQ=60)
                                                                                 GACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTT                                                                 <  SRR19973067.728253/150‑1 (MQ=60)
                                                                                    TTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCT                                                              <  SRR19973067.2740970/150‑1 (MQ=60)
                                                                                               TGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATT                                                   <  SRR19973067.1077490/150‑1 (MQ=60)
                                                                                               TGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATT                                                   <  SRR19973067.1077425/150‑1 (MQ=60)
                                                                                                GGTGATGGCACATACGGTCGCGGTATTGCTGTCTTTTGTTGGCGAACAAGACACCCACCCCCACAGATAAAAACACTGGCGCGCCCGGGAGCGCGGATATCATCAAGTTTTGCATCCCTGTAGTGCCGGGAATTTGATGTTGCTGAAaaa                                                  >  SRR19973067.2810098/1‑147 (MQ=14)
                                                                                                GGTGATGGCACATACGGTCGCGGTATTGCTCTCTTTTTATTGCGAACTAGAAACCCAGCCACACATAGAACAAAACTGGCGCGCCGGTAAGCGCGGATATATTCCAGTTTTGCATCCCATTAGTTCCGGTGATTTGCTGGTGCTGAATTA                                                  >  SRR19973067.3207377/1‑150 (MQ=60)
                                                                                                  TGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACC                                                <  SRR19973067.1190031/150‑1 (MQ=60)
                                                                                                     TGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATC                                             <  SRR19973067.1158250/150‑1 (MQ=60)
                                                                                                     TGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATC                                             <  SRR19973067.1158208/150‑1 (MQ=60)
                                                                                                                      GTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTG                            <  SRR19973067.1618165/150‑1 (MQ=60)
                                                                                                                      GTATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTG                            <  SRR19973067.1618109/150‑1 (MQ=60)
                                                                                                                        ATTCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGAT                          >  SRR19973067.738053/1‑150 (MQ=60)
                                                                                                                          TCCTCTCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGA                        <  SRR19973067.868500/150‑1 (MQ=60)
                                                                                                                              ATCTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAG                    <  SRR19973067.657836/150‑1 (MQ=60)
                                                                                                                                CTTTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGT                  <  SRR19973067.2491402/150‑1 (MQ=60)
                                                                                                                                  TTTGATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTG                <  SRR19973067.282016/150‑1 (MQ=60)
                                                                                                                                     GATGGCGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAG             <  SRR19973067.1183325/150‑1 (MQ=60)
                                                                                                                                          CGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCA        >  SRR19973067.2937236/1‑150 (MQ=60)
                                                                                                                                          CGAACTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGCATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCA        <  SRR19973067.765946/150‑1 (MQ=60)
                                                                                                                                           GATCTAGACACCCAGCCACTCATAGAACAACTCTGGCGCGGCGGTAAGCGCGTATTTGTTCCTGTTTTGCATCCCTTTTGTGCCGGTTATTTGATGTTCGTGAATTACCATCCGCTAAGCGTTCTGGTGATGAACAGGTTGAAGATCCAT       <  SRR19973067.830696/150‑1 (MQ=60)
                                                                                                                                                AGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCATGAGCC  >  SRR19973067.1603320/1‑150 (MQ=60)
                                                                                                                                                AGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCATGAGCC  >  SRR19973067.950964/1‑150 (MQ=60)
                                                                                                                                                |                                                                                                                                                     
CAAAATGATTCGGCAACGTCGTCGTGCGTTAACGCCGGAACAACAGCAGGAAATGGGTCAACAAGCCGCTACCCGGATGATGACTTATCCCCCGGTGGTGATGGCACATACGGTCGCTGTATTCCTCTCTTTTGATGGCGAACTCGACACCCAGCCACTCATAGAACAACTCTGGCGCGCCGGTAAGCGCGTATATCTTCCAGTTTTGCATCCCTTTAGTGCCGGTAATTTGCTGTTCCTGAATTACCATCCGCAAAGCGAACTGGTGATGAACAGGTTGAAGATCCATGAGCC  >  CP009273/3049596‑3049889

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