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 F309 I0 R1 1697 197.2 6627756 95.7% 6342762 147.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,723,809 G→T 83.0% E156* (GAA→TAA)  lhr → putative ATP‑dependent helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,723,8090GT83.0% 317.8 / 77.6 158E156* (GAA→TAA) lhrputative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base G (9/18);  new base T (66/65);  total (75/83)
Fisher's exact test for biased strand distribution p-value = 1.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 158 total aligned reads displayed.

CCCGTTGAAGGGTATTGCCGATGAACGGCGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCGAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCGATGCGCTG  >  CP009273/1723664‑1723952
                                                                                                                                                 |                                                                                                                                               
  cGTTGAAGGGTATTGCCGATGAACGGCGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATc                                                                                                                                             <  2:872861/148‑1 (MQ=255)
  cGTTGAAGGGTATTGCCGATGAACGGCGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATc                                                                                                                                             <  2:872833/148‑1 (MQ=255)
  cGTTGAAGGGTATTGCCGATGAACGGCGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATc                                                                                                                                             <  2:617176/148‑1 (MQ=255)
   gTTGAAGGGTATTGCCGATGAACGGCGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCa                                                                                                                                            <  2:1215771/148‑1 (MQ=255)
    ttGAAGGGTATTGCCGATGAACGGCGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCAc                                                                                                                                           >  2:3044102/1‑148 (MQ=255)
      gAAGGGTATTGCCGATGAACGGCGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCGAATCActc                                                                                                                                         <  1:2652831/148‑1 (MQ=255)
      gAAGGGTATTGCCGATGAACGGCGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCGAATCActc                                                                                                                                         <  1:1999581/148‑1 (MQ=255)
       aaGGGTATTGCCGATGAACGGCGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCActct                                                                                                                                        <  2:293919/148‑1 (MQ=255)
       aaGGGTATTGCCGATGAACGGCGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCGAATCActct                                                                                                                                        <  1:1149409/148‑1 (MQ=255)
         gggTATTGCCGATGAACGGCGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTAt                                                                                                                                      <  1:2379441/148‑1 (MQ=255)
         gggTATTGCCGATGAACGGCGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTAt                                                                                                                                      >  1:1649297/1‑148 (MQ=255)
         gggTATTGCCGATGAACGGCGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTAt                                                                                                                                      >  2:2091270/1‑148 (MQ=255)
         gggTATTGCCGATGAACGGCGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCGAATCACTCTAt                                                                                                                                      <  1:3070235/148‑1 (MQ=255)
          ggTATTGCCGATGAACGGCGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATc                                                                                                                                     <  2:482715/148‑1 (MQ=255)
           gTATTGCCGATGAACGGCGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCt                                                                                                                                    >  2:3008721/1‑148 (MQ=255)
               tGCCGATGAACGGCGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATg                                                                                                                                >  2:1796334/1‑148 (MQ=255)
                   gATGAACGGCGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGa                                                                                                                            <  2:1050390/148‑1 (MQ=255)
                    aTGAACGGCGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGAc                                                                                                                           >  1:2026418/1‑148 (MQ=255)
                    aTGAACGGCGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGAc                                                                                                                           >  1:2149922/1‑148 (MQ=255)
                      gAACGGCGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCGAATCACTCTATCTGATGCTGACCt                                                                                                                         <  2:493651/148‑1 (MQ=255)
                       aaCGGCGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTc                                                                                                                        <  1:597041/148‑1 (MQ=255)
                       aaCGGCGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTc                                                                                                                        <  1:597053/148‑1 (MQ=255)
                        aCGGCGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCGAATCACTCTATCTGATGCTGACCTcc                                                                                                                       <  1:1720630/148‑1 (MQ=255)
                                ggcgcggcgAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCcgcgcgcg                                                                                                               <  2:1004028/148‑1 (MQ=255)
                                  cgcggcgAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCcgcgcgcgcg                                                                                                             >  1:3030895/1‑148 (MQ=255)
                                   gcggcgAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGa                                                                                                            <  1:864318/148‑1 (MQ=255)
                                    cggcgAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGaa                                                                                                           <  1:3008721/148‑1 (MQ=255)
                                     ggcgAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGaaa                                                                                                          <  2:1327103/148‑1 (MQ=255)
                                       cgAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACg                                                                                                        <  2:2706452/148‑1 (MQ=255)
                                       cgAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACg                                                                                                        <  2:2430894/148‑1 (MQ=255)
                                       cgAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACg                                                                                                        <  2:2430873/148‑1 (MQ=255)
                                       cgAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACg                                                                                                        >  1:1267527/1‑148 (MQ=255)
                                       cgAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCGAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACg                                                                                                        >  2:2999138/1‑148 (MQ=255)
                                          aaCGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTa                                                                                                     >  2:2362992/1‑148 (MQ=255)
                                                 gTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACgcggcg                                                                                              >  1:1130040/1‑148 (MQ=255)
                                                  tCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATCTTCTGATCACCACACCCGAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTCCGCGGCGt                                                                                             >  2:2551229/1‑148 (MQ=255)
                                                     aTCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGa                                                                                          >  1:2145172/1‑148 (MQ=255)
                                                     aTCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGa                                                                                          >  2:1419815/1‑148 (MQ=255)
                                                      tCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGaa                                                                                         <  2:937520/148‑1 (MQ=255)
                                                       cTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGaaa                                                                                        <  1:454581/148‑1 (MQ=255)
                                                           gcgTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGt                                                                                    <  1:2942968/148‑1 (MQ=255)
                                                           gcgTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGt                                                                                    <  1:3085782/148‑1 (MQ=255)
                                                           gcgTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGt                                                                                    <  1:2901903/148‑1 (MQ=255)
                                                           gcgTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGt                                                                                    <  1:2942949/148‑1 (MQ=255)
                                                            cgTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTa                                                                                   <  1:3076261/148‑1 (MQ=255)
                                                              tAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAat                                                                                 >  2:2285604/1‑148 (MQ=255)
                                                              tAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGGAat                                                                                 >  2:3083642/1‑148 (MQ=255)
                                                               aGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAatt                                                                                >  2:882555/1‑148 (MQ=255)
                                                               aGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAatt                                                                                >  1:2466555/1‑148 (MQ=255)
                                                               aGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCGAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAatt                                                                                <  1:1672087/148‑1 (MQ=255)
                                                               aGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCGAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAatt                                                                                <  1:1672092/148‑1 (MQ=255)
                                                               aGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCGAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAatt                                                                                <  1:1672120/148‑1 (MQ=255)
                                                                   aTCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATtg                                                                            <  2:1431173/148‑1 (MQ=255)
                                                                   aTCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATtg                                                                            >  2:38808/1‑148 (MQ=255)
                                                                   aTCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATtg                                                                            <  1:1796334/148‑1 (MQ=255)
                                                                    tCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATtga                                                                           <  1:3262402/148‑1 (MQ=255)
                                                                      cGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATtgatg                                                                         <  2:2019731/148‑1 (MQ=255)
                                                                      cGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATtgatg                                                                         >  2:833271/1‑148 (MQ=255)
                                                                      cGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATtgatg                                                                         <  1:3056535/148‑1 (MQ=255)
                                                                          cTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGt                                                                     <  1:822274/148‑1 (MQ=255)
                                                                          cTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCGAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGt                                                                     <  2:1250887/148‑1 (MQ=255)
                                                                          cTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCGAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGt                                                                     <  2:1250109/148‑1 (MQ=255)
                                                                           tGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTc                                                                    >  2:1182068/1‑148 (MQ=255)
                                                                                aTACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCGAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCAcgg                                                               >  2:1573631/1‑147 (MQ=255)
                                                                                  aCGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCgg                                                             <  1:740922/148‑1 (MQ=255)
                                                                                       tGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCg                                                        >  1:874166/1‑148 (MQ=255)
                                                                                        gCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCgg                                                       >  1:1029964/1‑148 (MQ=255)
                                                                                         cacaGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCggg                                                      <  2:2885444/148‑1 (MQ=255)
                                                                                         cacaGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCggg                                                      <  1:1005960/148‑1 (MQ=255)
                                                                                            aGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAg                                                   >  2:1413873/1‑148 (MQ=255)
                                                                                              gAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTa                                                 >  2:2183612/1‑148 (MQ=255)
                                                                                              gAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTa                                                 >  1:1668540/1‑148 (MQ=255)
                                                                                                   cagcaAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACgt                                            <  2:3292555/148‑1 (MQ=255)
                                                                                                   cagcaAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACgt                                            <  1:771256/148‑1 (MQ=255)
                                                                                                     gcaAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACgtgg                                          >  2:425098/1‑148 (MQ=255)
                                                                                                         aCTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTgcg                                      <  2:918636/148‑1 (MQ=255)
                                                                                                         aCTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTgcg                                      >  1:2537159/1‑148 (MQ=255)
                                                                                                          cTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTgcgc                                     <  2:1431930/148‑1 (MQ=255)
                                                                                                          cTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTgcgc                                     >  1:125437/1‑148 (MQ=255)
                                                                                                                  tAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCg                             >  2:2662902/1‑148 (MQ=255)
                                                                                                                     tCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTa                          >  2:1519637/1‑148 (MQ=255)
                                                                                                                     tCCGCCGGATATTCTGATCACCACACCCGAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTa                          >  1:1752405/1‑148 (MQ=255)
                                                                                                                           ggATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTg                    >  2:220046/1‑148 (MQ=255)
                                                                                                                            gATATTCTGATCACCACACCCTACTCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTgg                   <  2:1251413/148‑1 (MQ=255)
                                                                                                                            gATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAGTTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTgg                   <  2:351128/148‑1 (MQ=255)
                                                                                                                              tatTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAg                 <  1:2878019/148‑1 (MQ=255)
                                                                                                                               atTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGc                >  2:918271/1‑148 (MQ=255)
                                                                                                                                  cTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGc             <  2:1649297/148‑1 (MQ=255)
                                                                                                                                  cTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGc             <  2:1987735/148‑1 (MQ=255)
                                                                                                                                     aTCACCACACCCTAATCCCTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCg          <  2:177453/148‑1 (MQ=255)
                                                                                                                                     aTCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTACACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCg          >  1:3188199/1‑148 (MQ=255)
                                                                                                                                      tCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCGa         <  1:562138/148‑1 (MQ=255)
                                                                                                                                      tCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCGa         >  2:1112187/1‑148 (MQ=255)
                                                                                                                                      tCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCGa         >  2:1173005/1‑148 (MQ=255)
                                                                                                                                         ccacACCCGAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCGATgc      >  2:489746/1‑148 (MQ=255)
                                                                                                                                          cacaCCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCGATgcg     <  1:3044122/148‑1 (MQ=255)
                                                                                                                                           acacCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTCAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCGATgcgc    >  1:323352/1‑148 (MQ=255)
                                                                                                                                            cacCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCGATGCgct   >  1:1954127/1‑148 (MQ=255)
                                                                                                                                            cacCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCGATGCgct   >  2:3091473/1‑148 (MQ=255)
                                                                                                                                             acCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCGACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCGATGCgctg  <  1:1258450/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
CCCGTTGAAGGGTATTGCCGATGAACGGCGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCGAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCGATGCGCTG  >  CP009273/1723664‑1723952

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
TCCCGTTGAAGGGTATTGCCGATGAACGGCGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCGAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCGATGCGCTGCTCCAC  >  CP009273/1723663‑1723958
                                                                                                                                                  |                                                                                                                                                     
TCCCGTTGAAGGGTATTGCCGATGAACGGCGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAAT                                                                                                                                                    <  SRR19973070.89514/150‑1 (MQ=60)
TCCCGTTGAAGGGTATTGCCGATGAACGGCGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAAT                                                                                                                                                    <  SRR19973070.89328/150‑1 (MQ=60)
  CCGTTGAAGGGTATTGCCGATGAACGGCGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCGAATCA                                                                                                                                                  <  SRR19973070.3089005/150‑1 (MQ=60)
       GAAGGGTATTGCCGATGAACGGCGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCGAATCACTCTA                                                                                                                                             <  SRR19973070.2013125/150‑1 (MQ=60)
       GAAGGGTATTGCCGATGAACGGCGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCGAATCACTCTA                                                                                                                                             <  SRR19973070.2670092/150‑1 (MQ=60)
        AAGGGTATTGCCGATGAACGGCGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTAT                                                                                                                                            >  SRR19973070.1660683/1‑150 (MQ=60)
        AAGGGTATTGCCGATGAACGGCGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCGAATCACTCTAT                                                                                                                                            <  SRR19973070.1157412/150‑1 (MQ=60)
          GGGTATTGCCGATGAACGGCGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCT                                                                                                                                          <  SRR19973070.2395164/150‑1 (MQ=60)
          GGGTATTGCCGATGAACGGCGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCGAATCACTCTATCT                                                                                                                                          <  SRR19973070.3089935/150‑1 (MQ=60)
          GGGTATTGCCGATGAACGGCGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCGAATCACTCTATCT                                                                                                                                          <  SRR19973070.3089950/150‑1 (MQ=60)
             TATTGCCGATGAACGGCGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGAT                                                                                                                                       <  SRR19973070.990501/150‑1 (MQ=60)
                   CGATGAACGGCGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGAC                                                                                                                                 >  SRR19973070.2040124/1‑150 (MQ=60)
                   CGATGAACGGCGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGAC                                                                                                                                 >  SRR19973070.2164159/1‑150 (MQ=60)
                   CGATGAACGGCGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGAC                                                                                                                                 >  SRR19973070.2164238/1‑150 (MQ=60)
                     ATGAACGGCGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCT                                                                                                                               >  SRR19973070.2457071/1‑150 (MQ=60)
                        AACGGCGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCC                                                                                                                            <  SRR19973070.601204/150‑1 (MQ=60)
                        AACGGCGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCC                                                                                                                            <  SRR19973070.601192/150‑1 (MQ=60)
                         ACGGCGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCGAATCACTCTATCTGATGCTGACCTCCCG                                                                                                                           <  SRR19973070.1732537/150‑1 (MQ=60)
                             CGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCG                                                                                                                       <  SRR19973070.578491/150‑1 (MQ=60)
                                 GGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCG                                                                                                                   >  SRR19973070.3050182/1‑150 (MQ=60)
                                    GCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAA                                                                                                                <  SRR19973070.870345/150‑1 (MQ=60)
                                     CGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAAC                                                                                                               <  SRR19973070.3027934/150‑1 (MQ=60)
                                      GGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACG                                                                                                              >  SRR19973070.1276398/1‑150 (MQ=60)
                                           AACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACG                                                                                                         <  SRR19973070.807935/150‑1 (MQ=60)
                                             CGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCGAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCG                                                                                                       >  SRR19973070.2911712/1‑150 (MQ=60)
                                                AAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCG                                                                                                    >  SRR19973070.1137838/1‑150 (MQ=60)
                                                AAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACATCCCGCGCGCGCGAAACGCTACGCGGCG                                                                                                    >  SRR19973070.1137831/1‑150 (MQ=60)
                                                    CAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGA                                                                                                >  SRR19973070.2159468/1‑150 (MQ=60)
                                                        CTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACG                                                                                            <  SRR19973070.457717/150‑1 (MQ=60)
                                                        CTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACG                                                                                            <  SRR19973070.457845/150‑1 (MQ=60)
                                                            GCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAA                                                                                        <  SRR19973070.3105629/150‑1 (MQ=60)
                                                            GCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAA                                                                                        <  SRR19973070.2920673/150‑1 (MQ=60)
                                                            GCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAA                                                                                        <  SRR19973070.2961913/150‑1 (MQ=60)
                                                            GCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAA                                                                                        <  SRR19973070.2961932/150‑1 (MQ=60)
                                                             CGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAAT                                                                                       <  SRR19973070.3096020/150‑1 (MQ=60)
                                                              GTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATT                                                                                      >  SRR19973070.2070925/1‑150 (MQ=60)
                                                              GTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATT                                                                                      >  SRR19973070.2482649/1‑150 (MQ=60)
                                                                AGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTAT                                                                                    >  SRR19973070.2932457/1‑150 (MQ=60)
                                                                AGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCGAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAGACGCTACGCGGCGTCGAAACGGTAATTAT                                                                                    <  SRR19973070.1683559/150‑1 (MQ=60)
                                                                AGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCGAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTAT                                                                                    <  SRR19973070.1683550/150‑1 (MQ=60)
                                                                AGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCGAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTAT                                                                                    <  SRR19973070.1683583/150‑1 (MQ=60)
                                                                AGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCGAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTAT                                                                                    <  SRR19973070.1683555/150‑1 (MQ=60)
                                                                    ATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGAT                                                                                <  SRR19973070.1808741/150‑1 (MQ=60)
                                                                     TCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATG                                                                               <  SRR19973070.3283399/150‑1 (MQ=60)
                                                                       CGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAA                                                                             <  SRR19973070.3076059/150‑1 (MQ=60)
                                                                           CTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCC                                                                         <  SRR19973070.827879/150‑1 (MQ=60)
                                                                               CGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGC                                                                     >  SRR19973070.1996315/1‑150 (MQ=60)
                                                                                   ACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTA                                                                 <  SRR19973070.745721/150‑1 (MQ=60)
                                                                                   ACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTA                                                                 <  SRR19973070.746158/150‑1 (MQ=60)
                                                                                      CCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCG                                                              >  SRR19973070.880307/1‑150 (MQ=60)
                                                                                       CTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGG                                                             >  SRR19973070.1037254/1‑150 (MQ=60)
                                                                                          CACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCA                                                          <  SRR19973070.1013129/150‑1 (MQ=60)
                                                                                             AGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTA                                                       >  SRR19973070.1679995/1‑150 (MQ=60)
                                                                                                    CAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGG                                                <  SRR19973070.776663/150‑1 (MQ=60)
                                                                                                        AAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCG                                            >  SRR19973070.2553745/1‑150 (MQ=60)
                                                                                                         AACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGC                                           >  SRR19973070.126539/1‑150 (MQ=60)
                                                                                                          ACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAAGTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCA                                          >  SRR19973070.2641368/1‑150 (MQ=60)
                                                                                                           CTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCAT                                         >  SRR19973070.3023409/1‑150 (MQ=60)
                                                                                                              ACCCGTAATCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTG                                      >  SRR19973070.1683123/1‑150 (MQ=60)
                                                                                                                    AATCCGCCGGATATTCTGATCACCACACCCGAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTA                                >  SRR19973070.1764546/1‑150 (MQ=60)
                                                                                                                    AATCCGCCGGATATTCTGATCACCACACCCGAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTA                                >  SRR19973070.1764505/1‑150 (MQ=60)
                                                                                                                    AATCCGCCGGATATTCTGATCACCACACCCGAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTA                                >  SRR19973070.1764567/1‑150 (MQ=60)
                                                                                                                      TCCGCCGGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAG                              <  SRR19973070.3108791/150‑1 (MQ=60)
                                                                                                                            GGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGA                        <  SRR19973070.1259023/150‑1 (MQ=60)
                                                                                                                            GGATATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGA                        >  SRR19973070.866652/1‑150 (MQ=60)
                                                                                                                               TATTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCG                     <  SRR19973070.2896659/150‑1 (MQ=60)
                                                                                                                                 TTCTGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGC                   >  SRR19973070.1536144/1‑150 (MQ=60)
                                                                                                                                    TGATCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTACACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCG                >  SRR19973070.3208826/1‑150 (MQ=60)
                                                                                                                                       TCACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCGATG             <  SRR19973070.565860/150‑1 (MQ=60)
                                                                                                                                         ACCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCGATGCG           >  SRR19973070.3267202/1‑150 (MQ=60)
                                                                                                                                          CCACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTCAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCGATGCGC          >  SRR19973070.325703/1‑150 (MQ=60)
                                                                                                                                           CACACCCTAATCACTCTATCTTATGCTGACCTCCCGCGCGCGCGAAACGCTCCGCGGCGTCTAAACGGTAATTATTGATGAAGCCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCGATGCGCT         <  SRR19973070.3063453/150‑1 (MQ=60)
                                                                                                                                           CACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCGATGCGCT         <  SRR19973070.3063473/150‑1 (MQ=60)
                                                                                                                                           CACACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCGATGCGCT         >  SRR19973070.1967272/1‑150 (MQ=60)
                                                                                                                                              ACCCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCGACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCGATGCGCTGCT      <  SRR19973070.1267271/150‑1 (MQ=60)
                                                                                                                                              ACCCGACTCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCGATGCGCTGCT      <  SRR19973070.1267403/150‑1 (MQ=60)
                                                                                                                                                CCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCGATGCGCTGCTCC    <  SRR19973070.817977/150‑1 (MQ=60)
                                                                                                                                                CCTAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCGATGCGCTGCTCC    <  SRR19973070.1994374/150‑1 (MQ=60)
                                                                                                                                                  TAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCGATGCGCTGCTCCAC  <  SRR19973070.1343373/150‑1 (MQ=60)
                                                                                                                                                  |                                                                                                                                                     
TCCCGTTGAAGGGTATTGCCGATGAACGGCGGCGGCGCGGCGAAACGGAAGTCAATCTTCGCGTAGGGATCCGTACTGGCGATACGCCTGCACAGGAACGCAGCAAACTCACCCGTAATCCGCCGGATATTCTGATCACCACACCCGAATCACTCTATCTGATGCTGACCTCCCGCGCGCGCGAAACGCTACGCGGCGTCGAAACGGTAATTATTGATGAAGTCCACGCGGTAGCGGGCAGTAAACGTGGTGCGCATCTGGCGTTAAGTCTGGAGCGGCTCGATGCGCTGCTCCAC  >  CP009273/1723663‑1723958

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