Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A1 F309 I0 R1 1697 197.2 6627756 95.7% 6342762 147.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 754,572 G→T 100% E137D (GAG→GAT sucA → 2‑oxoglutarate decarboxylase, thiamine triphosphate‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273754,5720GT99.4% 506.9 / ‑4.0 158E137D (GAG→GATsucA2‑oxoglutarate decarboxylase, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (82/75);  total (82/76)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01
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 158 total aligned reads displayed.

AGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAA  >  CP009273/754427‑754716
                                                                                                                                                 |                                                                                                                                                
   ttAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCtt                                                                                                                                             >  2:2122857/1‑148 (MQ=255)
   ttAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCtt                                                                                                                                             >  1:1197876/1‑148 (MQ=255)
     aaaGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCa                                                                                                                                           <  2:2676725/148‑1 (MQ=255)
     aaaGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCa                                                                                                                                           <  2:1393255/148‑1 (MQ=255)
      aaGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCaa                                                                                                                                          <  1:2441182/148‑1 (MQ=255)
       aGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAAc                                                                                                                                         >  1:1349963/1‑148 (MQ=255)
            tGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCgg                                                                                                                                    >  1:2100543/1‑148 (MQ=255)
            tGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCgg                                                                                                                                    >  1:2100876/1‑148 (MQ=255)
            tGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCgg                                                                                                                                    <  1:1938493/148‑1 (MQ=255)
                gCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCa                                                                                                                                <  2:1405018/148‑1 (MQ=255)
                gCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCa                                                                                                                                >  2:1952097/1‑148 (MQ=255)
                     ttAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGc                                                                                                                           >  2:245653/1‑148 (MQ=255)
                     ttAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGc                                                                                                                           <  2:1080454/148‑1 (MQ=255)
                        aCGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAg                                                                                                                        >  2:1099492/1‑148 (MQ=255)
                           cATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCgg                                                                                                                     >  2:2815478/1‑148 (MQ=255)
                              aCCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGcca                                                                                                                  >  2:1764781/1‑146 (MQ=255)
                              aCCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCaa                                                                                                                  >  1:60220/1‑148 (MQ=255)
                                 gCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGa                                                                                                               >  1:1283971/1‑148 (MQ=255)
                                   ttCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGaaa                                                                                                             <  1:1094899/148‑1 (MQ=255)
                                   ttCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGaaa                                                                                                             >  1:436856/1‑148 (MQ=255)
                                   ttCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGaaa                                                                                                             <  2:248523/148‑1 (MQ=255)
                                   ttCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGaaa                                                                                                             <  2:1918461/148‑1 (MQ=255)
                                     ccGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAAcc                                                                                                           <  2:174358/148‑1 (MQ=255)
                                         ggTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGa                                                                                                       <  1:2059521/148‑1 (MQ=255)
                                          gTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGaa                                                                                                      <  1:1068495/148‑1 (MQ=255)
                                            cACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAAc                                                                                                    >  1:371045/1‑148 (MQ=255)
                                              ccAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTc                                                                                                  >  1:2764555/1‑148 (MQ=255)
                                              ccAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTc                                                                                                  >  1:3213059/1‑148 (MQ=255)
                                                aGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCgg                                                                                                <  2:458382/148‑1 (MQ=255)
                                                  cATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCg                                                                                              <  1:284670/148‑1 (MQ=255)
                                                     gCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAgc                                                                                           >  2:2270640/1‑148 (MQ=255)
                                                     gCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAgc                                                                                           >  2:145309/1‑148 (MQ=255)
                                                     gCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAgc                                                                                           <  2:1068067/148‑1 (MQ=255)
                                                       gAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAgctg                                                                                         <  2:2674649/148‑1 (MQ=255)
                                                        aaTCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCCAAGTAACCATGAAACTCGGCGAgctgc                                                                                        >  1:2689005/1‑148 (MQ=255)
                                                           ctcGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTgg                                                                                     >  2:457742/1‑148 (MQ=255)
                                                           ctcGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTgg                                                                                     >  2:431021/1‑148 (MQ=255)
                                                              gATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAg                                                                                  <  2:3197164/148‑1 (MQ=255)
                                                               aTCCGCTGGGACTGTGTCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGc                                                                                 <  1:2087117/148‑1 (MQ=255)
                                                                tCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGcc                                                                                >  1:866295/1‑148 (MQ=255)
                                                                 ccGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGACACTCGGCGAGCTGCTGGAAGccc                                                                               >  2:3212343/1‑148 (MQ=255)
                                                                 ccGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGccc                                                                               <  1:2664452/148‑1 (MQ=255)
                                                                  cGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCt                                                                              >  2:68603/1‑148 (MQ=255)
                                                                       ggACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGc                                                                         >  2:509817/1‑148 (MQ=255)
                                                                           tgtgGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAAc                                                                     >  1:2850958/1‑148 (MQ=255)
                                                                            gtgGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAAcc                                                                    <  1:3276540/148‑1 (MQ=255)
                                                                               gcagcaAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTAc                                                                 >  1:1613760/1‑148 (MQ=255)
                                                                                  gcaAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGc                                                              <  2:2224751/148‑1 (MQ=255)
                                                                                      gATAAAGTGGCCGATCTGGTTCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGcc                                                          >  1:1373733/1‑148 (MQ=255)
                                                                                      gATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGcc                                                          >  2:741316/1‑148 (MQ=255)
                                                                                      gATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGcc                                                          <  2:2047873/148‑1 (MQ=255)
                                                                                      gATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGcc                                                          >  1:1771202/1‑148 (MQ=255)
                                                                                      gATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGcc                                                          >  1:1706348/1‑148 (MQ=255)
                                                                                      gATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGcc                                                          >  2:2369901/1‑148 (MQ=255)
                                                                                       aTAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGccc                                                         >  1:119720/1‑148 (MQ=255)
                                                                                        tAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCg                                                        <  1:3108776/148‑1 (MQ=255)
                                                                                        tAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCg                                                        <  1:3108778/148‑1 (MQ=255)
                                                                                          aaGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGAt                                                      >  2:737865/1‑148 (MQ=255)
                                                                                            gTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTg                                                    <  1:837656/148‑1 (MQ=255)
                                                                                            gTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTg                                                    >  2:1624899/1‑148 (MQ=255)
                                                                                            gTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTg                                                    >  2:568780/1‑148 (MQ=255)
                                                                                            gTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTg                                                    >  1:2767004/1‑148 (MQ=255)
                                                                                             tGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTgg                                                   >  1:2163460/1‑148 (MQ=255)
                                                                                               gCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTg                                                 <  2:920614/148‑1 (MQ=255)
                                                                                                ccGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGc                                                >  1:2637750/1‑148 (MQ=255)
                                                                                                  gATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCg                                              <  1:3186084/148‑1 (MQ=255)
                                                                                                   aTCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGa                                             >  2:2006858/1‑148 (MQ=255)
                                                                                                   aTCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGa                                             >  2:338715/1‑148 (MQ=255)
                                                                                                    tCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAg                                            >  2:2692731/1‑148 (MQ=255)
                                                                                                         aTCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTatat                                       >  2:1867605/1‑148 (MQ=255)
                                                                                                         aTCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTatat                                       <  2:736452/148‑1 (MQ=255)
                                                                                                         aTCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTatat                                       >  2:2400062/1‑148 (MQ=255)
                                                                                                          tCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATg                                      >  2:3268624/1‑148 (MQ=255)
                                                                                                          tCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGAGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATg                                      >  2:605387/1‑148 (MQ=255)
                                                                                                           ccGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGc                                     <  2:2596878/148‑1 (MQ=255)
                                                                                                           ccGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGc                                     <  2:2596885/148‑1 (MQ=255)
                                                                                                                 ttCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTa                               >  2:1293343/1‑148 (MQ=255)
                                                                                                                 ttCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTa                               <  2:1813196/148‑1 (MQ=255)
                                                                                                                 ttCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTa                               <  2:1813226/148‑1 (MQ=255)
                                                                                                                  tCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTAc                              <  1:2495920/148‑1 (MQ=255)
                                                                                                                  tCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTAc                              <  1:1154378/148‑1 (MQ=255)
                                                                                                                    cACGATCTGACCGAAGCAGACTTCCATGATACCTTCAACGTCGGTTCTTTTGACAGCGGCAAAGAAACCATGAAACTCGGCGAGCTTCTGGAAGACCTCAAGCAAACCTACTGCGGCCCGTTTGGTGCCGAGTATATGGACATTACCa                            <  2:2516814/148‑1 (MQ=255)
                                                                                                                    cACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCa                            <  2:2673322/148‑1 (MQ=255)
                                                                                                                    cACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCa                            >  1:1655725/1‑148 (MQ=255)
                                                                                                                    cACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCa                            >  1:3293892/1‑148 (MQ=255)
                                                                                                                        aTCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCTGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCTAGCTGCTGGAAGCCCTCACGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCAc                        >  2:3033541/1‑148 (MQ=255)
                                                                                                                             aCCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCgaag                   >  1:69351/1‑148 (MQ=255)
                                                                                                                              ccGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCgaaga                  >  2:1431864/1‑148 (MQ=255)
                                                                                                                                 aaGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGaaaa               <  1:3196302/148‑1 (MQ=255)
                                                                                                                                 aaGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGaaaa               >  2:1002110/1‑148 (MQ=255)
                                                                                                                                  aGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGaaaaa              <  1:114506/148‑1 (MQ=255)
                                                                                                                                    cAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACg            <  2:1197876/148‑1 (MQ=255)
                                                                                                                                    cAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACg            <  2:1483320/148‑1 (MQ=255)
                                                                                                                                      gACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCt          <  1:68603/148‑1 (MQ=255)
                                                                                                                                       aCTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTg         >  1:2127034/1‑148 (MQ=255)
                                                                                                                                       aCTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACAGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTg         >  2:1275017/1‑148 (MQ=255)
                                                                                                                                         ttCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGa       >  1:933138/1‑148 (MQ=255)
                                                                                                                                           ccAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATc     <  1:769982/148‑1 (MQ=255)
                                                                                                                                           ccAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATc     <  1:2562586/148‑1 (MQ=255)
                                                                                                                                              ggATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCaa  <  2:297055/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
AGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCAA  >  CP009273/754427‑754716

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
AAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCA  >  CP009273/754423‑754715
                                                                                                                                                     |                                                                                                                                               
AAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGAT                                                                                                                                                 <  SRR19973070.2318914/150‑1 (MQ=60)
   CAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACC                                                                                                                                              <  SRR19973070.1549028/150‑1 (MQ=60)
   CAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACC                                                                                                                                              <  SRR19973070.1549006/150‑1 (MQ=60)
     GGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTT                                                                                                                                            >  SRR19973070.1206293/1‑150 (MQ=60)
      GTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTC                                                                                                                                           <  SRR19973070.679294/150‑1 (MQ=60)
         AAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAAC                                                                                                                                        >  SRR19973070.3255238/1‑150 (MQ=60)
         AAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAAC                                                                                                                                        >  SRR19973070.1359213/1‑150 (MQ=60)
          AAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACG                                                                                                                                       <  SRR19973070.2457177/150‑1 (MQ=60)
              CCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGG                                                                                                                                   >  SRR19973070.2114639/1‑150 (MQ=60)
              CCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGG                                                                                                                                   >  SRR19973070.2114973/1‑150 (MQ=60)
                TGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTT                                                                                                                                 <  SRR19973070.1951524/150‑1 (MQ=60)
                   AGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCAT                                                                                                                              >  SRR19973070.1130161/1‑150 (MQ=60)
                                ATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAA                                                                                                                 >  SRR19973070.60906/1‑150 (MQ=60)
                                  ACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAG                                                                                                               >  SRR19973070.3030119/1‑150 (MQ=60)
                                   CCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGA                                                                                                              >  SRR19973070.1292932/1‑150 (MQ=60)
                                   CCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGA                                                                                                              >  SRR19973070.1292917/1‑150 (MQ=60)
                                     GCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAA                                                                                                            >  SRR19973070.440026/1‑150 (MQ=60)
                                       TTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACC                                                                                                          <  SRR19973070.2715407/150‑1 (MQ=60)
                                       TTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACC                                                                                                          <  SRR19973070.1102451/150‑1 (MQ=60)
                                             GGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAA                                                                                                    <  SRR19973070.2073426/150‑1 (MQ=60)
                                              GTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAAC                                                                                                   <  SRR19973070.1075970/150‑1 (MQ=60)
                                              GTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAAC                                                                                                   >  SRR19973070.373766/1‑150 (MQ=60)
                                                CACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTC                                                                                                 >  SRR19973070.2782265/1‑150 (MQ=60)
                                                CACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTC                                                                                                 >  SRR19973070.3233838/1‑150 (MQ=60)
                                                CACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTC                                                                                                 >  SRR19973070.1584351/1‑150 (MQ=60)
                                                   CAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGC                                                                                              >  SRR19973070.1793069/1‑150 (MQ=60)
                                                      CATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAG                                                                                           <  SRR19973070.286565/150‑1 (MQ=60)
                                                          CGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCCAAGTAACCATGAAACTCGGCGAGCTGC                                                                                       >  SRR19973070.2706428/1‑150 (MQ=60)
                                                                TCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAG                                                                                 <  SRR19973070.1312212/150‑1 (MQ=60)
                                                                  GATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCC                                                                               >  SRR19973070.204450/1‑150 (MQ=60)
                                                                  GATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCC                                                                               >  SRR19973070.872344/1‑150 (MQ=60)
                                                                   ATCCGCTGGGACTGTGTCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCC                                                                              <  SRR19973070.2101154/150‑1 (MQ=60)
                                                                     CCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTC                                                                            <  SRR19973070.2681766/150‑1 (MQ=60)
                                                                             ACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAAC                                                                    >  SRR19973070.2869421/1‑150 (MQ=60)
                                                                                GTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTA                                                                 <  SRR19973070.3297596/150‑1 (MQ=60)
                                                                                 TGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTAC                                                                >  SRR19973070.1624996/1‑150 (MQ=60)
                                                                                        AAGATAAAGTGGCCGATCTGGTTCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCC                                                         >  SRR19973070.1383069/1‑150 (MQ=60)
                                                                                        AAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCC                                                         >  SRR19973070.1783460/1‑150 (MQ=60)
                                                                                        AAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCC                                                         >  SRR19973070.1718124/1‑150 (MQ=60)
                                                                                         AGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCC                                                        >  SRR19973070.120795/1‑150 (MQ=60)
                                                                                          GATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACAGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCG                                                       <  SRR19973070.2605284/150‑1 (MQ=60)
                                                                                            TAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGAT                                                     <  SRR19973070.3128832/150‑1 (MQ=60)
                                                                                            TAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGAT                                                     <  SRR19973070.3128830/150‑1 (MQ=60)
                                                                                            TAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGAT                                                     <  SRR19973070.2149014/150‑1 (MQ=60)
                                                                                              AAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTG                                                   >  SRR19973070.2784723/1‑150 (MQ=60)
                                                                                               AGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGG                                                  <  SRR19973070.91780/150‑1 (MQ=60)
                                                                                               AGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGG                                                  >  SRR19973070.2177828/1‑150 (MQ=60)
                                                                                                GTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGT                                                 <  SRR19973070.843319/150‑1 (MQ=60)
                                                                                                GTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGT                                                 <  SRR19973070.844416/150‑1 (MQ=60)
                                                                                                  GGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGC                                               >  SRR19973070.2654942/1‑150 (MQ=60)
                                                                                                      GATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAG                                           <  SRR19973070.3206701/150‑1 (MQ=60)
                                                                                                      GATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAG                                           >  SRR19973070.2107113/1‑150 (MQ=60)
                                                                                                                      TCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCA                           <  SRR19973070.2512106/150‑1 (MQ=60)
                                                                                                                      TCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCA                           <  SRR19973070.1162433/150‑1 (MQ=60)
                                                                                                                      TCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCA                           >  SRR19973070.1667135/1‑150 (MQ=60)
                                                                                                                      TCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCA                           >  SRR19973070.3315014/1‑150 (MQ=60)
                                                                                                                      TCCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCA                           >  SRR19973070.3315017/1‑150 (MQ=60)
                                                                                                                       CCACGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAG                          <  SRR19973070.1716434/150‑1 (MQ=60)
                                                                                                                          CGATCTGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCAC                       <  SRR19973070.247342/150‑1 (MQ=60)
                                                                                                                               TGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAG                  <  SRR19973070.2686272/150‑1 (MQ=60)
                                                                                                                               TGACCGAAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAG                  >  SRR19973070.70116/1‑150 (MQ=60)
                                                                                                                                     AAGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAAC            <  SRR19973070.3216984/150‑1 (MQ=60)
                                                                                                                                      AGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACG           <  SRR19973070.116797/150‑1 (MQ=60)
                                                                                                                                      AGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACG           <  SRR19973070.2833682/150‑1 (MQ=60)
                                                                                                                                      AGCAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACG           <  SRR19973070.115557/150‑1 (MQ=60)
                                                                                                                                        CAGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCT         <  SRR19973070.2137065/150‑1 (MQ=60)
                                                                                                                                         AGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTG        >  SRR19973070.2141261/1‑150 (MQ=60)
                                                                                                                                         AGACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTG        <  SRR19973070.762963/150‑1 (MQ=60)
                                                                                                                                          GACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGG       <  SRR19973070.41138/150‑1 (MQ=60)
                                                                                                                                          GACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGG       <  SRR19973070.69364/150‑1 (MQ=60)
                                                                                                                                           ACTTCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGA      >  SRR19973070.939796/1‑150 (MQ=60)
                                                                                                                                              TCCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCC   <  SRR19973070.3212736/150‑1 (MQ=60)
                                                                                                                                               CCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCA  <  SRR19973070.775383/150‑1 (MQ=60)
                                                                                                                                               CCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCA  <  SRR19973070.3172553/150‑1 (MQ=60)
                                                                                                                                               CCAGGATACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCA  <  SRR19973070.2579332/150‑1 (MQ=60)
                                                                                                                                                     |                                                                                                                                               
AAGCAGGTTAAAGTCCTGCAGCTCATTAACGCATACCGCTTCCGTGGTCACCAGCATGCGAATCTCGATCCGCTGGGACTGTGGCAGCAAGATAAAGTGGCCGATCTGGATCCGTCTTTCCACGATCTGACCGAAGCAGACTTCCAGGAGACCTTCAACGTCGGTTCATTTGCCAGCGGCAAAGAAACCATGAAACTCGGCGAGCTGCTGGAAGCCCTCAAGCAAACCTACTGCGGCCCGATTGGTGCCGAGTATATGCACATTACCAGCACCGAAGAAAAACGCTGGATCCA  >  CP009273/754423‑754715

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: