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 2,665,772 G→T 100% E123* (GAA→TAA)  hcaD → phenylpropionate dioxygenase, ferredoxin reductase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,665,7720GT95.5% 430.5 / 12.8 155E123* (GAA→TAA) hcaDphenylpropionate dioxygenase, ferredoxin reductase subunit
Reads supporting (aligned to +/- strand):  ref base G (2/5);  new base T (70/78);  total (72/83)
Fisher's exact test for biased strand distribution p-value = 4.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.20e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 155 total aligned reads displayed.

TCTGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGA  >  CP009273/2665627‑2665916
                                                                                                                                                 |                                                                                                                                                
tCTGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATaa                                                                                                                                                <  2:2709580/148‑1 (MQ=255)
tCTGCAGTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATaa                                                                                                                                                <  2:2709577/148‑1 (MQ=255)
tCTGCAGTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATaa                                                                                                                                                <  2:2611380/148‑1 (MQ=255)
 cTGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCCGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAAc                                                                                                                                               >  2:2606518/1‑148 (MQ=255)
    cATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCGCGACCGCGGCCGGTGGGTGATGGACTGGGATAACgcg                                                                                                                                            >  1:2180037/1‑147 (MQ=255)
     aTTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACgctg                                                                                                                                           <  2:2858377/148‑1 (MQ=255)
      ttCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCGCGACCGCTGCCGTTGGTTGTTTCACGGGGATAAAgctgc                                                                                                                                          >  1:3190912/1‑148 (MQ=255)
       tCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCCTTGCTTGTTTCACAGGGATAAAgctgct                                                                                                                                         >  1:2037148/1‑148 (MQ=255)
        ccGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCtt                                                                                                                                        <  1:3306322/148‑1 (MQ=255)
         cGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCGGCCGTTGCGTGATGCACAGGGATAACGCGGCgtt                                                                                                                                       >  1:2882798/1‑148 (MQ=255)
           gtgtAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCCGCGCGACCGCAGCCGGTGCGTGTTGGACGGGGGTAACGCGGCGTTaa                                                                                                                                     >  1:1015945/1‑147 (MQ=255)
            tgtAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTAcc                                                                                                                                    <  2:689278/148‑1 (MQ=255)
             gtAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGGAGCTCGACCGCTGCCGTTGGGTGATGGCCTGGGGTAACGCTGCGTTAccc                                                                                                                                   >  1:2748811/1‑148 (MQ=255)
                aCCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTAACCTAc                                                                                                                                >  2:2645837/1‑148 (MQ=255)
                    tCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCa                                                                                                                            <  2:1009730/148‑1 (MQ=255)
                    tCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCa                                                                                                                            <  1:3142499/148‑1 (MQ=255)
                    tCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCa                                                                                                                            <  2:909528/148‑1 (MQ=255)
                     cAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCCGCGCGACCGCTGCCGTTTCTTGATTGACAGGGTTAACGCTTCTTTAACCAACGCCaa                                                                                                                           >  2:903922/1‑147 (MQ=255)
                         aCGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCg                                                                                                                       >  2:2031840/1‑148 (MQ=255)
                          cGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCgg                                                                                                                      >  2:1705182/1‑148 (MQ=255)
                          cGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCGTGATGCACTGGGATAACGCTGCTTTAACCTACGCCATGCCgg                                                                                                                      >  2:793250/1‑148 (MQ=255)
                          cGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCGCGACCGCGGCCGGTGCTTGATGCCCTGGGCTAACGCTGCGTTACCCTACGCCCTGCCgg                                                                                                                      >  2:3110485/1‑148 (MQ=255)
                             tGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGa                                                                                                                   >  2:2539586/1‑148 (MQ=255)
                                gTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATgc                                                                                                                <  1:1966320/148‑1 (MQ=255)
                                  cgcgACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATgccg                                                                                                              <  1:427400/148‑1 (MQ=255)
                                    cgACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATTCCGGTGATgccgcc                                                                                                            >  2:3239638/1‑148 (MQ=255)
                                       cacacGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGGTGGTTGATGGACTGGGGTAAAGCGGCTTTAACCAACCCCAAGCCGGTGATGCCGCCaca                                                                                                         >  1:1735560/1‑146 (MQ=255)
                                        acacGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGAc                                                                                                        <  1:2335977/148‑1 (MQ=255)
                                         cacGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCAACGCCATGCCGGTGATGCCGCCAGACt                                                                                                       >  2:182713/1‑148 (MQ=255)
                                          acGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTTCTTGATGCACTGGGATAACGCTGCTTTAACCTACGCCATGCCGGTGATGCCGCCAGAcag                                                                                                      >  2:2214817/1‑146 (MQ=255)
                                          acGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTg                                                                                                      >  2:3019240/1‑148 (MQ=255)
                                            gagaGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCg                                                                                                    >  2:2365347/1‑148 (MQ=255)
                                            gagaGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCg                                                                                                    <  1:927054/148‑1 (MQ=255)
                                                gTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGaa                                                                                                <  1:212341/148‑1 (MQ=255)
                                                gTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGaa                                                                                                <  2:1419247/148‑1 (MQ=255)
                                                     tgtTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCGCGACCCCGGCCGGTTGTTGTTGCCCAGGGGTAACCCTGCTTTAACCAACGCCAAGCCCGGGATGGCGCCCGACAGCGCGAAGTTCt                                                                                           >  1:2190449/1‑148 (MQ=255)
                                                      gtTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTg                                                                                          >  2:1450165/1‑148 (MQ=255)
                                                      gtTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTg                                                                                          >  2:3218905/1‑148 (MQ=255)
                                                      gtTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTg                                                                                          <  2:830892/148‑1 (MQ=255)
                                                      gtTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTg                                                                                          >  2:3218178/1‑148 (MQ=255)
                                                         aaCCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAg                                                                                       <  1:710147/148‑1 (MQ=255)
                                                               cGGCGACAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGaa                                                                                 <  1:1209462/148‑1 (MQ=255)
                                                               cGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGaa                                                                                 <  1:2564844/148‑1 (MQ=255)
                                                               cGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGaa                                                                                 <  2:3259859/148‑1 (MQ=255)
                                                                  cGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACgg                                                                              <  1:1532398/148‑1 (MQ=255)
                                                                     aaGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGGAGCGCGACCGCTGCCGTTTGGTGATGCACAGGGATAACGCGGCTTTAAACTACGCCAAGCCCGTGATGCCGCCAGACAGCGAGAAGTTCTGCAGCCCGCACGGtga                                                                           >  2:639321/1‑146 (MQ=255)
                                                                     aaGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCa                                                                           <  1:2642762/148‑1 (MQ=255)
                                                                     aaGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCa                                                                           <  2:3036243/148‑1 (MQ=255)
                                                                      aGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCGCGACCGCGGCCGGTGCTTGGTTCCCGGGGATAAAGCAGCTTTACCCTACCCCCAGCCCGGGATGCCGCCCGACGGCGAGAAGGTCTGCACCCCCAACCGTccg                                                                          >  2:968026/1‑146 (MQ=255)
                                                                        cTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCCGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCGACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTATGCAGCCCCAACGGTCAGTc                                                                        >  2:1884678/1‑148 (MQ=255)
                                                                         tGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCg                                                                       <  1:40667/148‑1 (MQ=255)
                                                                         tGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCg                                                                       <  2:2405848/148‑1 (MQ=255)
                                                                         tGGCACGGGGATCAGCTTTTGAGAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACGGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCg                                                                       <  2:2406559/148‑1 (MQ=255)
                                                                              cTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCAAGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGAtt                                                                  >  2:2043224/1‑148 (MQ=255)
                                                                               tGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTa                                                                 <  1:300125/148‑1 (MQ=255)
                                                                               tGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTa                                                                 <  1:1979516/148‑1 (MQ=255)
                                                                                gggATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCAGCCGTTTCTTGATTCACTGGGATAAAGCTGCTTTAACCTACGCCAAGCCCGTGATGCCGCCAGACTTCGAGAAGTTCTGCAGCCCCAACCGTCAGTCGGGATTAt                                                                >  2:3494740/1‑148 (MQ=255)
                                                                                  gATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCg                                                              <  1:2758597/148‑1 (MQ=255)
                                                                                     cAGCTTTTTATAGCAACCGGCGCGGCCGCTCGACCGCTGCCGGTGGTTGTTTGACGGGGATAACGCTGCGTTACCCTACGCCCAGCCGGTGATGCCGCCCGACGGCGCGAAGTTCTGCAGCCCGAACGGTCCGTCGTGGTTATAGGTg                                                           >  2:1258637/1‑148 (MQ=255)
                                                                                     cAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTg                                                           <  2:1805096/148‑1 (MQ=255)
                                                                                     cAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTg                                                           <  1:1390704/148‑1 (MQ=255)
                                                                                       gCTTTTTATAGCAACCGGCGCGGGAGCGCCACCGCGGCCGTTGCGTGATTGAAGGGGATAAAGCTGCGTTAACCTAAGCCAAGGCCGTGTTTCCGCCCCACAGCGAGAAGTTCTTCACCCCGAACGGTCACTCGGGCTTTTAGGTggc                                                         >  2:1270687/1‑146 (MQ=255)
                                                                                          ttttATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGa                                                      <  1:1457834/148‑1 (MQ=255)
                                                                                          ttttATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGa                                                      <  1:3303401/148‑1 (MQ=255)
                                                                                          ttttATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGa                                                      <  1:3171698/148‑1 (MQ=255)
                                                                                          ttttATAGCAACCGGCGCGGCAGCGCGACCGCTGCCGTTGCTTGATGCACAGGGCTAACGCTGCGTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGGGATTATCGGTGCCGGa                                                      >  2:497478/1‑148 (MQ=255)
                                                                                           tttATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGaa                                                     <  1:1842140/148‑1 (MQ=255)
                                                                                                       ggCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTgg                                         >  2:2667979/1‑148 (MQ=255)
                                                                                                         cgcgGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGaa                                       <  2:68508/148‑1 (MQ=255)
                                                                                                           cgGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACt                                     <  1:3034921/148‑1 (MQ=255)
                                                                                                             gCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTgg                                   <  2:3438801/148‑1 (MQ=255)
                                                                                                             gCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTgg                                   <  1:904443/148‑1 (MQ=255)
                                                                                                             gCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTgg                                   <  1:3218905/148‑1 (MQ=255)
                                                                                                             gCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTgg                                   <  1:3218178/148‑1 (MQ=255)
                                                                                                             gCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTgg                                   <  1:1450165/148‑1 (MQ=255)
                                                                                                              cAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGc                                  >  1:1098050/1‑148 (MQ=255)
                                                                                                              cAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGc                                  >  1:1388420/1‑148 (MQ=255)
                                                                                                                 cTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGc                               >  2:19565/1‑148 (MQ=255)
                                                                                                                 cTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGc                               >  1:3440000/1‑148 (MQ=255)
                                                                                                                  tCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGcc                              >  2:2415005/1‑148 (MQ=255)
                                                                                                                   cGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCa                             >  1:1364820/1‑148 (MQ=255)
                                                                                                                   cGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCa                             >  2:2284296/1‑148 (MQ=255)
                                                                                                                     aCCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAgc                           >  2:296845/1‑148 (MQ=255)
                                                                                                                          tGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAc                      >  2:2414932/1‑148 (MQ=255)
                                                                                                                            ccGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAcgc                    >  2:151282/1‑148 (MQ=255)
                                                                                                                              gTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAcgcag                  <  2:2357956/148‑1 (MQ=255)
                                                                                                                              gTTGCTTGATGCACTGGGATAACGCTGCTTTACCCGACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAcgcag                  <  2:37413/148‑1 (MQ=255)
                                                                                                                               ttGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAcgcagc                 >  2:1675084/1‑148 (MQ=255)
                                                                                                                               ttGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAcgcagc                 >  2:3082854/1‑148 (MQ=255)
                                                                                                                                tGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAcgcagcg                <  1:2942142/148‑1 (MQ=255)
                                                                                                                                tGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAcgcagcg                <  1:1375485/148‑1 (MQ=255)
                                                                                                                                 gCTTGATGCACTGGGATAACGCTGCTTTACCGTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAcgcagcgc               <  1:2812018/148‑1 (MQ=255)
                                                                                                                                 gCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAcgcagcgc               <  1:2811653/148‑1 (MQ=255)
                                                                                                                                    tGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGa            <  2:3298113/148‑1 (MQ=255)
                                                                                                                                    tGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGa            <  2:1159805/148‑1 (MQ=255)
                                                                                                                                    tGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGa            <  1:1140489/148‑1 (MQ=255)
                                                                                                                                      aTGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATg          >  2:180688/1‑148 (MQ=255)
                                                                                                                                      aTGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATg          >  2:3225391/1‑148 (MQ=255)
                                                                                                                                        gCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTa        <  1:430455/148‑1 (MQ=255)
                                                                                                                                             gggATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTg   <  2:2896581/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TCTGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGA  >  CP009273/2665627‑2665916

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
TTCATCTGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAG  >  CP009273/2665623‑2665919
                                                                                                                                                     |                                                                                                                                                   
TTCATCTGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAT                                                                                                                                                     <  SRR19973067.575285/150‑1 (MQ=60)
 TCATCTGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATA                                                                                                                                                    <  SRR19973067.1954343/150‑1 (MQ=60)
      TGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCT                                                                                                                                               <  SRR19973067.1536722/150‑1 (MQ=60)
         ATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCCTTGCTTGTTTCACAGGGATAAAGCTGCT                                                                                                                                            >  SRR19973067.2050785/1‑150 (MQ=60)
         ANTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCT                                                                                                                                            >  SRR19973067.3262881/1‑150 (MQ=60)
           TCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCGGCCGTTGCGTGATGCACAGGGATAACGCGGCGTT                                                                                                                                          >  SRR19973067.2901236/1‑150 (MQ=60)
            CCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTA                                                                                                                                         <  SRR19973067.3327182/150‑1 (MQ=60)
               GTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGGAGCTCGACCGCTGCCGTTGGGTGATGGCCTGGGGTAACGCTGCGTTACCC                                                                                                                                      >  SRR19973067.2766577/1‑150 (MQ=60)
                        TCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATG                                                                                                                             <  SRR19973067.3162537/150‑1 (MQ=60)
                        TCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATG                                                                                                                             <  SRR19973067.3162448/150‑1 (MQ=60)
                                    GTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCG                                                                                                                 <  SRR19973067.3074860/150‑1 (MQ=60)
                                    GTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCG                                                                                                                 <  SRR19973067.1979669/150‑1 (MQ=60)
                                      CGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCC                                                                                                               <  SRR19973067.430229/150‑1 (MQ=60)
                                      CGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCC                                                                                                               <  SRR19973067.430367/150‑1 (MQ=60)
                                         GACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGGTGGTTGATGGACTGGGGTAAAGCGGCTTTAACCAACCCCAAGCCGGTGATGCCGCCACA                                                                                                            >  SRR19973067.1747434/1‑150 (MQ=60)
                                            ACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTG                                                                                                         <  SRR19973067.2351059/150‑1 (MQ=60)
                                            ACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTG                                                                                                         <  SRR19973067.2351281/150‑1 (MQ=60)
                                              ACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTTATGCCGCCAGACTGCG                                                                                                       <  SRR19973067.2172992/150‑1 (MQ=60)
                                                GAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAG                                                                                                     <  SRR19973067.933499/150‑1 (MQ=60)
                                                GAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAG                                                                                                     <  SRR19973067.933312/150‑1 (MQ=60)
                                                    GTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGT                                                                                                 <  SRR19973067.213919/150‑1 (MQ=60)
                                                       AGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCTACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCT                                                                                              <  SRR19973067.168562/150‑1 (MQ=60)
                                                       AGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCCGCTCGACCGCGGCCGGTGCGTGTTGCACAGGGGAAACGCTGCTTTAgccaaaccccaacccggggtggcgccccacagcgccaacggca                                                                                              >  SRR19973067.2205513/1‑107 (MQ=60)
                                                           TTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAG                                                                                          <  SRR19973067.2556212/150‑1 (MQ=60)
                                                             AACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCC                                                                                        <  SRR19973067.715070/150‑1 (MQ=60)
                                                                 AACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAA                                                                                    <  SRR19973067.3443534/150‑1 (MQ=60)
                                                                   CGGCGACAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACG                                                                                  <  SRR19973067.1217421/150‑1 (MQ=60)
                                                                   CGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACG                                                                                  <  SRR19973067.2581596/150‑1 (MQ=60)
                                                                      CGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTC                                                                               <  SRR19973067.1542761/150‑1 (MQ=60)
                                                                         AAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGT                                                                            <  SRR19973067.2659831/150‑1 (MQ=60)
                                                                             TGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTG                                                                        <  SRR19973067.41235/150‑1 (MQ=60)
                                                                                   TGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATC                                                                  <  SRR19973067.1993719/150‑1 (MQ=60)
                                                                                   TGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATC                                                                  <  SRR19973067.1992921/150‑1 (MQ=60)
                                                                                   TGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATC                                                                  <  SRR19973067.302042/150‑1 (MQ=60)
                                                                                      GATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGT                                                               <  SRR19973067.2776424/150‑1 (MQ=60)
                                                                                       ATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTG                                                              <  SRR19973067.1266997/150‑1 (MQ=60)
                                                                                         CAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCC                                                            <  SRR19973067.1400157/150‑1 (MQ=60)
                                                                                              TTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAAC                                                       <  SRR19973067.1467589/150‑1 (MQ=60)
                                                                                              TTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAAC                                                       <  SRR19973067.3324228/150‑1 (MQ=60)
                                                                                              TTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAAC                                                       <  SRR19973067.3191780/150‑1 (MQ=60)
                                                                                              TTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAAC                                                       <  SRR19973067.2990619/150‑1 (MQ=60)
                                                                                               TTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACT                                                      <  SRR19973067.1854708/150‑1 (MQ=60)
                                                                                                             CGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACT                                        <  SRR19973067.2833736/150‑1 (MQ=60)
                                                                                                               CGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGG                                      <  SRR19973067.3054250/150‑1 (MQ=60)
                                                                                                                GGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGC                                     >  SRR19973067.1397860/1‑150 (MQ=60)
                                                                                                                GGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGC                                     >  SRR19973067.1397824/1‑150 (MQ=60)
                                                                                                                GGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGC                                     >  SRR19973067.1105527/1‑150 (MQ=60)
                                                                                                                 GCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCT                                    <  SRR19973067.909893/150‑1 (MQ=60)
                                                                                                                 GCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCT                                    <  SRR19973067.1459883/150‑1 (MQ=60)
                                                                                                                 GCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCT                                    <  SRR19973067.3238322/150‑1 (MQ=60)
                                                                                                                 GCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCT                                    <  SRR19973067.3238422/150‑1 (MQ=60)
                                                                                                                 GCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCT                                    <  SRR19973067.3239149/150‑1 (MQ=60)
                                                                                                                 GCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCT                                    <  SRR19973067.910415/150‑1 (MQ=60)
                                                                                                                   AGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGC                                  >  SRR19973067.3461898/1‑150 (MQ=60)
                                                                                                                     CTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCA                                >  SRR19973067.1374144/1‑150 (MQ=60)
                                                                                                                            GCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAC                         <  SRR19973067.1140468/150‑1 (MQ=60)
                                                                                                                                 CGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACGATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGC                    >  SRR19973067.1910423/1‑150 (MQ=60)
                                                                                                                                    TGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCA                 <  SRR19973067.2960789/150‑1 (MQ=60)
                                                                                                                                    TGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCA                 <  SRR19973067.1384872/150‑1 (MQ=60)
                                                                                                                                     GCTTGATGCACTGGGATAACGCTGCTTTACCGTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAG                <  SRR19973067.2830138/150‑1 (MQ=60)
                                                                                                                                     GCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAG                <  SRR19973067.2829770/150‑1 (MQ=60)
                                                                                                                                     GCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAG                <  SRR19973067.2829111/150‑1 (MQ=60)
                                                                                                                                        TGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATG             <  SRR19973067.1148153/150‑1 (MQ=60)
                                                                                                                                            GCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAG         <  SRR19973067.433289/150‑1 (MQ=60)
                                                                                                                                            GCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAG         <  SRR19973067.433449/150‑1 (MQ=60)
                                                                                                                                                 GGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGAC    <  SRR19973067.2662920/150‑1 (MQ=60)
                                                                                                                                                   GATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAG  >  SRR19973067.786722/1‑150 (MQ=60)
                                                                                                                                                   GATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAG  >  SRR19973067.786978/1‑150 (MQ=60)
                                                                                                                                                     |                                                                                                                                                   
TTCATCTGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAG  >  CP009273/2665623‑2665919

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