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 F434 I0 R1 1130 215.9 7031662 97.1% 6827743 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,7720GT97.9% 375.8 / ‑1.0 127E123* (GAA→TAA) hcaDphenylpropionate dioxygenase, ferredoxin reductase subunit
Reads supporting (aligned to +/- strand):  ref base G (3/1);  new base T (45/77);  total (49/78)
Fisher's exact test for biased strand distribution p-value = 1.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 128 total aligned reads displayed.

TCTGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACA  >  CP009273/2665627‑2665918
                                                                                                                                                 |                                                                                                                                                   
tCTGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATaa                                                                                                                                                   <  1:1101183/148‑1 (MQ=255)
tCTGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATaa                                                                                                                                                   <  1:1101465/148‑1 (MQ=255)
 cTGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAAc                                                                                                                                                  >  2:2775716/1‑148 (MQ=255)
  tGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACg                                                                                                                                                 <  2:1742/148‑1 (MQ=255)
    cATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCGCGACCGCGGCCGGTGGTGGTTGGACGGGGGTAACgct                                                                                                                                               >  1:3009759/1‑148 (MQ=255)
    cATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCGCGACCCCAGCCGGTGGGTGTTTCACTGGGTTAACgcg                                                                                                                                               >  1:3009848/1‑147 (MQ=255)
      ttCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACgctgc                                                                                                                                             <  2:1673080/148‑1 (MQ=255)
       tCCGGTGTAACGATCAACCCGCTGGGTCGCGACACACGAGAGTGAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCGGCCGTTGCTTGATGCACTGGGATAACgctgc‑t                                                                                                                                           <  2:806372/148‑1 (MQ=255)
       tCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACgctgc‑t                                                                                                                                           <  2:2566902/148‑1 (MQ=255)
        ccGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGcttc‑tt                                                                                                                                          <  1:1868734/148‑1 (MQ=255)
          ggTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGC‑TTTa                                                                                                                                        <  2:3026622/148‑1 (MQ=255)
          ggTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGC‑TTTa                                                                                                                                        <  1:2823649/148‑1 (MQ=255)
               aaCCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCCCTGCCGTTGCTTGATTGAAGGGGTTAACGCGGC‑TTTaaacaa                                                                                                                                   >  2:1310741/1‑143 (MQ=255)
               aaCCATCAAAACGCTGGCTCGCGACACACGAGAGTTATTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCACCTCCACCCCAGCCGTTTCATGATGGACGGGGTTACCACTGC‑ATTAcccaa                                                                                                                                   >  2:1310303/1‑146 (MQ=255)
                  cATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAAAACGCTGC‑GTTACCCAACGc                                                                                                                                >  2:431941/1‑148 (MQ=255)
                   aTCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGcc                                                                                                                               <  1:1230962/148‑1 (MQ=255)
                   aTCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGcc                                                                                                                               <  1:1171079/148‑1 (MQ=255)
                      aaaaCGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATg                                                                                                                            <  1:892245/148‑1 (MQ=255)
                            ttGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGC‑TTTACCCTACGCCATGCCGGTg                                                                                                                      <  1:901854/147‑1 (MQ=255)
                            cTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCCCTGCCGTTGCTTGATGCACGGGGATAACGCTGC‑TTTAACCAAAGCCAAGCCGGTg                                                                                                                      >  2:1235959/1‑148 (MQ=255)
                              gggTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGAt                                                                                                                    <  2:848543/148‑1 (MQ=255)
                               ggTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGGCGCGCGACCGCTGCCCGTTCTTGATTGACGGGGGGCACGCTGC‑GTTTAACAACGCCCAGCCGGGGAta                                                                                                                   >  2:901854/1‑147 (MQ=255)
                                       cacacGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCCGCACCCCCCCGGCCGGTGCGTGATGCACTGGGCTAACCCTGC‑TTTAACCAACCCCCAGCCCGGGATGCCGCCAGa                                                                                                           >  1:2168460/1‑148 (MQ=255)
                                       cacacGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGa                                                                                                           <  1:3489016/148‑1 (MQ=255)
                                        acacGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCCGCGCGACCGCGGCCGGTGGTTGGTGCACGGGGATCACGCTGG‑GTTACCCAACGCCATGCCGGGGGTggcgcccgcc                                                                                                          >  1:2559219/1‑144 (MQ=255)
                                          acGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGGGCGACCCCGGCCGGTGGTTGTTGGACGGGGGGAACGCAGC‑GTTAACCTAAGCCATGCCGGTGATTCCGCCAGAcag                                                                                                        >  1:1732896/1‑146 (MQ=255)
                                          acGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTg                                                                                                        >  2:2585576/1‑148 (MQ=255)
                                           cGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCGCGACCGCGGCCGTTGGTTGTTGCACTGGGATTAAGCGGC‑TTTAACCAACGCCCAGCCGGTGATTCCGCCACACAGc                                                                                                       >  1:3460414/1‑148 (MQ=255)
                                            gagaGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCg                                                                                                      <  1:2579739/148‑1 (MQ=255)
                                                 ttAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAg                                                                                                 >  2:1904370/1‑148 (MQ=255)
                                                  tAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGt                                                                                                <  2:533857/148‑1 (MQ=255)
                                                    gtgtTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGGCGCGCGACCGCGGCCGGTGGTTGTTGGAAAGGGGGAAAGCTGC‑GGTTCCCCACCCCCAGCCGGGGGTGCCGCCCGACTGCGCGGAGTTc                                                                                              >  1:993252/1‑148 (MQ=255)
                                                    gtgtTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTc                                                                                              <  1:991438/148‑1 (MQ=255)
                                                       ttAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGc                                                                                           <  1:2111384/148‑1 (MQ=255)
                                                        tAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCa                                                                                          <  2:1530008/148‑1 (MQ=255)
                                                        tAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCa                                                                                          <  2:3271537/148‑1 (MQ=255)
                                                        tAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCCCGGCCGTTTCGTGTTGCACTGGGATAAAGCTGC‑TTTAACCAACCCCATGCCGGTGATGCCGCCAGACAGCGAGAAGTTCTTCa                                                                                          >  2:1663137/1‑148 (MQ=255)
                                                               cGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGaa                                                                                   <  2:2011626/148‑1 (MQ=255)
                                                               cGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCGTGATGCCCGGGGTTAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGATAAGTTCTGCAGCCCGaa                                                                                   >  2:2131309/1‑148 (MQ=255)
                                                               cGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCGTGATGCACGGGGATAACGCTGC‑TTTACCCAACGCCAAGCCGGTGATGCCGCCCGACTGCGAGCAGTTATGCACCCCGaa                                                                                   >  2:2798732/1‑148 (MQ=255)
                                                               cGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCGGCCGTTGCTTGTTGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGaa                                                                                   >  2:1561909/1‑148 (MQ=255)
                                                                 gCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACg                                                                                 <  2:795721/148‑1 (MQ=255)
                                                                 gCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACg                                                                                 <  2:349359/148‑1 (MQ=255)
                                                                         tGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCg                                                                         <  2:1441345/148‑1 (MQ=255)
                                                                            cACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGa                                                                      <  2:211515/148‑1 (MQ=255)
                                                                            cACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGa                                                                      <  2:3282877/148‑1 (MQ=255)
                                                                              cTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGAtt                                                                    <  2:2385430/148‑1 (MQ=255)
                                                                              cTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGAtt                                                                    <  2:1280794/148‑1 (MQ=255)
                                                                              cTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGAtt                                                                    <  2:2386033/148‑1 (MQ=255)
                                                                               tGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTa                                                                   <  1:2847491/148‑1 (MQ=255)
                                                                               tGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTa                                                                   <  1:2779306/148‑1 (MQ=255)
                                                                               tGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTa                                                                   <  1:2698207/148‑1 (MQ=255)
                                                                               tGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTa                                                                   <  1:3300213/148‑1 (MQ=255)
                                                                                   aTCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCgg                                                               <  2:1049342/148‑1 (MQ=255)
                                                                                    tCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGt                                                              <  1:1235959/148‑1 (MQ=255)
                                                                                     cAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGGTGCGTGATGCACTGGGCTAACGCTTC‑TTTACCCTACGCCAAGCCGGGGATGGCCCCCGACTTCGAGAAGTTGTGCAGCCCGAAAGGGCCGGCGTGCTTTTCGGTg                                                             >  2:2472285/1‑148 (MQ=255)
                                                                                         tttttATAGCAACCGGCGCGGGAGCGCGACCGCGGCCGGTGCGTGTTGCACAGGGCTAAAGCTGCGTTTACCC‑ACGCCCTGCCGGGGATGCCGCCCGACAGCGAGAAGTTCTTCCCCCCGCCCGGGGAGTCGTGATTATCGGTGCCgg                                                         >  2:909658/1‑148 (MQ=255)
                                                                                         tttttATAGCAACCGGCGCGGCCGCTCGACCGCTGCCGGTGCGTTATTCACTGGGTTAAAGCGGC‑GTTACACGAAGCCAAGCCGGTGATTCCGCCCCACAGCGAGAAGTTATGCAGCCCCAACCGGCAGTCGGGATTATAGGTGGCgg                                                         >  2:909642/1‑148 (MQ=255)
                                                                                          ttttATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGa                                                        <  2:3052942/148‑1 (MQ=255)
                                                                                           tttATAGCAACCGGCGCGGCATCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGC‑TTGACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGaa                                                       <  2:3053583/148‑1 (MQ=255)
                                                                                            ttATAGCAACCGGCGCGGCCGCTCGACCGCTGCCGGTGGTTGATGCACGGGGCTAACGCGGC‑TTTTACCTACGCCACGCCGGGGGTGCCGCCCGACTGCGGGAAGTTCTGGAGCCCGAACGGTCAGGCGTGGTTATCGGTGGCGGAAc                                                      >  2:1562978/1‑148 (MQ=255)
                                                                                              ataGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTa                                                    <  2:1267767/148‑1 (MQ=255)
                                                                                               taGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTAt                                                   <  2:2356219/148‑1 (MQ=255)
                                                                                               taGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTAt                                                   <  2:592917/148‑1 (MQ=255)
                                                                                                aGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTAtt                                                  <  1:2166698/148‑1 (MQ=255)
                                                                                                aGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTAtt                                                  <  2:2033008/148‑1 (MQ=255)
                                                                                                aGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTAtt                                                  <  2:3320254/148‑1 (MQ=255)
                                                                                                aGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTAtt                                                  <  1:398074/148‑1 (MQ=255)
                                                                                                      cGGCGCGGCAGCTCGACCGCTGCCGATGCTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTg                                            <  1:3477668/148‑1 (MQ=255)
                                                                                                        gcgcGGGCGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCGGC‑GTTACCCAACGCCCGGCCGGTGATGCCGCCAGACTGCGAGACGTTCTGCAGCCCGAACCGGCCGTCGGGCTTATCGGTGGCGGAACAAATGGTCTggc                                          >  2:1536235/1‑147 (MQ=255)
                                                                                                         cgcgGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGaa                                         <  2:2814802/148‑1 (MQ=255)
                                                                                                             gCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTgg                                     <  1:2131460/148‑1 (MQ=255)
                                                                                                                gCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTg                                  >  2:2069265/1‑148 (MQ=255)
                                                                                                                 cTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTggcggc                                 <  2:873559/148‑1 (MQ=255)
                                                                                                                 cTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGc                                 >  2:1070231/1‑148 (MQ=255)
                                                                                                                   cGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCa                               >  2:2586159/1‑148 (MQ=255)
                                                                                                                   cGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCCGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCa                               >  1:2208378/1‑148 (MQ=255)
                                                                                                                      ccGCTGCCGTTGCTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAgcg                            >  2:1105345/1‑148 (MQ=255)
                                                                                                                       cGCTGCCGTTGCTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAgcgc                           <  1:1620279/148‑1 (MQ=255)
                                                                                                                       cGCTGCCGTTGCTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAgcgc                           <  1:2262861/148‑1 (MQ=255)
                                                                                                                            ccGTTGCTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAcgc                      >  2:2022358/1‑148 (MQ=255)
                                                                                                                             cGTTGCTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAcgca                     <  1:3178609/148‑1 (MQ=255)
                                                                                                                              gTTGCTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAcgcag                    >  1:3281608/1‑148 (MQ=255)
                                                                                                                              gTTGCTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAcgcag                    <  2:565375/148‑1 (MQ=255)
                                                                                                                              gTTGCTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAcgcag                    >  1:1541203/1‑148 (MQ=255)
                                                                                                                                  cTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAcgcagcgca                >  2:281092/1‑148 (MQ=255)
                                                                                                                                  cTTGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAcgcagcgca                >  2:1841820/1‑148 (MQ=255)
                                                                                                                                   ttGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAcgcagcgcag               <  1:808662/148‑1 (MQ=255)
                                                                                                                                   ttGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAcgcagcgcag               <  1:2943399/148‑1 (MQ=255)
                                                                                                                                   ttGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCAcgcagcgcag               <  1:2943694/148‑1 (MQ=255)
                                                                                                                                    tGATGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGa              <  1:2320221/148‑1 (MQ=255)
                                                                                                                                      aTGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATg            <  1:2585576/148‑1 (MQ=255)
                                                                                                                                       tGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGt           >  2:2985638/1‑148 (MQ=255)
                                                                                                                                       tGCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGt           <  1:1487388/148‑1 (MQ=255)
                                                                                                                                        gCACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTa          <  2:51721/148‑1 (MQ=255)
                                                                                                                                         cACTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTaa         <  2:1500335/148‑1 (MQ=255)
                                                                                                                                           cTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAgg       >  1:1006900/1‑148 (MQ=255)
                                                                                                                                           cTGGGATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAgg       <  2:1544419/148‑1 (MQ=255)
                                                                                                                                              ggATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGa    >  2:703110/1‑148 (MQ=255)
                                                                                                                                              ggATAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGa    >  2:242853/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TCTGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGC‑TTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACA  >  CP009273/2665627‑2665918

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
TCTGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGT  >  CP009273/2665627‑2665920
                                                                                                                                                 |                                                                                                                                                    
TCTGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACG                                                                                                                                                  <  SRR19973065.1108418/150‑1 (MQ=60)
TCTGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACG                                                                                                                                                  <  SRR19973065.1108702/150‑1 (MQ=60)
        CCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTTCTTTA                                                                                                                                          <  SRR19973065.1881033/150‑1 (MQ=60)
          GGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACC                                                                                                                                        <  SRR19973065.2841085/150‑1 (MQ=60)
                   ATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCAT                                                                                                                               <  SRR19973065.1178542/150‑1 (MQ=60)
                   ATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCAT                                                                                                                               <  SRR19973065.1239015/150‑1 (MQ=60)
                      AAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCC                                                                                                                            <  SRR19973065.897908/150‑1 (MQ=60)
                            TTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGAT                                                                                                                      <  SRR19973065.907648/150‑1 (MQ=60)
                                      ACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCCGCGCGACCGCGGCCGGTGGTTGGTGCACGGGGATCACGCTGGGTTACCCAACGCCATGCCGGgggtggcgcccgcc                                                                                                            >  SRR19973065.2575225/1‑136 (MQ=60)
                                       CACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACT                                                                                                           <  SRR19973065.3510234/150‑1 (MQ=60)
                                        ACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGGGCGACCCCGGCCGGTGGTTGTTGGACGGGGGGAACGCAGCGTTAACCTAAGCCATGCCGGTGATTCCGCCAGACAG                                                                                                          >  SRR19973065.1744364/1‑150 (MQ=60)
                                            GAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAG                                                                                                      <  SRR19973065.2595823/150‑1 (MQ=60)
                                                    GTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTG                                                                                              <  SRR19973065.998025/150‑1 (MQ=60)
                                                    GTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTG                                                                                              <  SRR19973065.998024/150‑1 (MQ=60)
                                                     TGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGC                                                                                             <  SRR19973065.1991490/150‑1 (MQ=60)
                                                       TTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAG                                                                                           <  SRR19973065.2124969/150‑1 (MQ=60)
                                                                    AAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAG                                                                              <  SRR19973065.2197535/150‑1 (MQ=60)
                                                                           GCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGAT                                                                       <  SRR19973065.2937712/150‑1 (MQ=60)
                                                                             ACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTA                                                                     <  SRR19973065.3121932/150‑1 (MQ=60)
                                                                               TGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATC                                                                   <  SRR19973065.2865012/150‑1 (MQ=60)
                                                                               TGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATC                                                                   <  SRR19973065.2714965/150‑1 (MQ=60)
                                                                               TGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATC                                                                   <  SRR19973065.2796392/150‑1 (MQ=60)
                                                                               TGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATC                                                                   <  SRR19973065.2796554/150‑1 (MQ=60)
                                                                               TGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATC                                                                   <  SRR19973065.3320314/150‑1 (MQ=60)
                                                                                GGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCG                                                                  <  SRR19973065.1238526/150‑1 (MQ=60)
                                                                                    TCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGC                                                              <  SRR19973065.1244063/150‑1 (MQ=60)
                                                                                                AGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGG                                                  <  SRR19973065.2180619/150‑1 (MQ=60)
                                                                                                AGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGG                                                  <  SRR19973065.400849/150‑1 (MQ=60)
                                                                                                      CGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGA                                            <  SRR19973065.3498827/150‑1 (MQ=60)
                                                                                                      CGGCGCGGCAGCTCGACCGCTGCCGATGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGA                                            <  SRR19973065.3498828/150‑1 (MQ=60)
                                                                                                             GCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCT                                     <  SRR19973065.2145175/150‑1 (MQ=60)
                                                                                                               AGCTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGC                                   >  SRR19973065.518153/1‑150 (MQ=60)
                                                                                                                 CTCGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCCGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCA                                 >  SRR19973065.2222500/1‑150 (MQ=60)
                                                                                                                   CGACCGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGC                               >  SRR19973065.2480183/1‑150 (MQ=60)
                                                                                                                       CGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCA                           <  SRR19973065.1631143/150‑1 (MQ=60)
                                                                                                                       CGCTGCCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCA                           <  SRR19973065.2277166/150‑1 (MQ=60)
                                                                                                                            CCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAG                      >  SRR19973065.3301558/1‑150 (MQ=60)
                                                                                                                            CCGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAG                      >  SRR19973065.1551392/1‑150 (MQ=60)
                                                                                                                             CGTTGCTTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGC                     <  SRR19973065.3197818/150‑1 (MQ=60)
                                                                                                                                   TTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGAT               <  SRR19973065.2961575/150‑1 (MQ=60)
                                                                                                                                   TTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGAT               <  SRR19973065.814003/150‑1 (MQ=60)
                                                                                                                                   TTGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGAT               <  SRR19973065.2961277/150‑1 (MQ=60)
                                                                                                                                    TGATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATG              <  SRR19973065.2334761/150‑1 (MQ=60)
                                                                                                                                      ATGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTA            <  SRR19973065.2601681/150‑1 (MQ=60)
                                                                                                                                       TGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAA           <  SRR19973065.1580629/150‑1 (MQ=60)
                                                                                                                                       TGCACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAA           <  SRR19973065.1496900/150‑1 (MQ=60)
                                                                                                                                         CACTGGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGG         >  SRR19973065.1013584/1‑150 (MQ=60)
                                                                                                                                             GGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGAC     >  SRR19973065.2287187/1‑150 (MQ=60)
                                                                                                                                             GGGATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCAAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGAC     >  SRR19973065.2287204/1‑150 (MQ=60)
                                                                                                                                               GATAACGCTGCTTTACCCTACGCCATTCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAG   >  SRR19973065.898081/1‑150 (MQ=60)
                                                                                                                                               GATAACGCTGCTTTACCCTACGCCATTCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAG   >  SRR19973065.898213/1‑150 (MQ=60)
                                                                                                                                               GATAACGCTGCTTTACCCTACGCCATTCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAG   >  SRR19973065.2844072/1‑150 (MQ=60)
                                                                                                                                                ATAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGT  <  SRR19973065.1665058/150‑1 (MQ=60)
                                                                                                                                                 |                                                                                                                                                    
TCTGCATTCCGGTGTAACCATCAAAACGCTGGGTCGCGACACACGAGAGTTAGTGTTAACCAACGGCGAAAGCTGGCACTGGGATCAGCTTTTTATAGCAACCGGCGCGGCAGCTCGACCGCTGCCGTTGCTTGATGCACTGGGAGAACGCTGCTTTACCCTACGCCATGCCGGTGATGCCGCCAGACTGCGAGAAGTTCTGCAGCCCGAACGGTCAGTCGTGATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGT  >  CP009273/2665627‑2665920

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: ‑