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 F1 I2 R1 593 264.2 13311096 97.8% 13018251 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,987,822 2 bp→CG coding (578‑579/963 nt) hemC ← porphobilinogen deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,987,8220GC100.0% 1012.3 / NA 275G193A (GGC→GCG) hemCporphobilinogen deaminase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (128/147);  total (128/147)
*NC_0009133,987,8230CG99.6% 998.5 / NA 272G193A (GGC→GCG) hemCporphobilinogen deaminase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (124/147);  minor base T (1/0);  total (125/147)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 284 total aligned reads displayed.

CAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAAGCCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTG  >  NC_000913/3987735‑3987909
                                                                                       ||                                                                                      
cagcTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGgc                                                                                    <  1:6446210/93‑1 (MQ=255)
 agcTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACgcggcg                                                                                   >  2:281763/1‑93 (MQ=255)
  gcTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGgcgc                                                                                  >  2:1487158/1‑93 (MQ=255)
     cgcgAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGaa                                                                               <  1:1844552/93‑1 (MQ=255)
     cgcgAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGaa                                                                               <  2:3366780/93‑1 (MQ=255)
      gcgAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAAt                                                                              <  2:4455281/93‑1 (MQ=255)
        gAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATAc                                                                            <  1:1660069/93‑1 (MQ=255)
        gAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATAc                                                                            >  1:3423537/1‑93 (MQ=255)
        gAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATAc                                                                            <  1:2610415/93‑1 (MQ=255)
         aGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACg                                                                           >  2:5751418/1‑93 (MQ=255)
         aGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACg                                                                           <  1:2346411/93‑1 (MQ=255)
           tGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTg                                                                         <  2:1312828/93‑1 (MQ=255)
             cGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGAc                                                                       >  2:561812/1‑93 (MQ=255)
              gTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACt                                                                      >  1:3583722/1‑93 (MQ=255)
                gAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTcc                                                                    >  1:533247/1‑93 (MQ=255)
                gAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTcc                                                                    >  2:6401895/1‑93 (MQ=255)
                 aaTCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATATGTGACTCCa                                                                   >  2:2139876/1‑93 (MQ=255)
                 aaTCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCa                                                                   <  2:3926975/93‑1 (MQ=255)
                 aaTCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCa                                                                   <  1:2096833/93‑1 (MQ=255)
                    catcaAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGAc                                                                <  1:3150023/93‑1 (MQ=255)
                      tcaAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCt                                                              >  1:1797651/1‑93 (MQ=255)
                         aGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTaaa                                                           >  2:5503741/1‑93 (MQ=255)
                         aGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTaaa                                                           >  2:6521042/1‑93 (MQ=255)
                         aGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTaaa                                                           <  1:5121114/93‑1 (MQ=255)
                          ggcggcATTCAATCCCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAAc                                                          <  1:1701225/93‑1 (MQ=255)
                          ggcggcATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAAc                                                          <  1:2207168/93‑1 (MQ=255)
                            cggcATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGt                                                        <  2:2013305/93‑1 (MQ=255)
                             ggcATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGtt                                                       >  1:6248671/1‑93 (MQ=255)
                               cATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGtttt                                                     >  2:5202778/1‑93 (MQ=255)
                               cATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGtttt                                                     <  2:6042631/93‑1 (MQ=255)
                                aTTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTa                                                    <  1:2823813/93‑1 (MQ=255)
                                 ttCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTat                                                   >  2:6264958/1‑92 (MQ=255)
                                 ttCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAg                                                   >  1:4676042/1‑93 (MQ=255)
                                 ttCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAg                                                   <  1:3626749/93‑1 (MQ=255)
                                 ttCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAg                                                   >  2:4085805/1‑93 (MQ=255)
                                  tCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGt                                                  >  1:700444/1‑93 (MQ=255)
                                  tCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGt                                                  <  2:3305219/93‑1 (MQ=255)
                                   cAATCCCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTc                                                 <  2:4549060/93‑1 (MQ=255)
                                   cAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTc                                                 <  2:1281667/93‑1 (MQ=255)
                                    aaTACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTcc                                                >  1:1105641/1‑93 (MQ=255)
                                        cccACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCt                                            <  2:374950/93‑1 (MQ=255)
                                         ccACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTa                                           >  1:2391496/1‑93 (MQ=255)
                                          cACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTAc                                          >  2:4242229/1‑93 (MQ=255)
                                          cACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTAc                                          <  1:6342779/93‑1 (MQ=255)
                                            ccGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACgg                                        <  2:2707366/93‑1 (MQ=255)
                                            ccGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACgg                                        >  1:1005532/1‑93 (MQ=255)
                                              gCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCa                                      >  1:4432546/1‑93 (MQ=255)
                                                 ccTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGa                                   >  2:5798115/1‑93 (MQ=255)
                                                  cTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGaa                                  <  2:2268313/93‑1 (MQ=255)
                                                   ttGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAat                                 >  2:5017588/1‑93 (MQ=255)
                                                   ttGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAat                                 >  1:4896122/1‑93 (MQ=255)
                                                    tGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAatg                                <  2:2394668/93‑1 (MQ=255)
                                                    tGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAatg                                <  1:527775/93‑1 (MQ=255)
                                                     gTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAatga                               >  1:4444753/1‑93 (MQ=255)
                                                     gTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAatga                               >  1:3724953/1‑93 (MQ=255)
                                                      tCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAatgat                              <  2:75466/93‑1 (MQ=255)
                                                      tCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAatgat                              <  2:3134591/93‑1 (MQ=255)
                                                      tCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAatgat                              >  2:650431/1‑93 (MQ=255)
                                                       ccTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAatgatg                             <  1:1940943/93‑1 (MQ=255)
                                                       ccTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAatgatg                             <  1:5623045/93‑1 (MQ=255)
                                                       ccTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAatgatg                             >  2:2313341/1‑93 (MQ=255)
                                                        cTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATgg                            >  1:3217876/1‑93 (MQ=255)
                                                         tACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGc                           >  2:2590379/1‑93 (MQ=255)
                                                         tACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGc                           >  1:4396317/1‑93 (MQ=255)
                                                         tACCGCCGGAAGAGAAATCTCGGGGGGCAACTCGGCGCGAATACGAGACTCCAGACCTAAACGTTTTAGCCCAGCTAGCGCAAGAATTATGGc                           >  1:3527109/1‑93 (MQ=255)
                                                           ccgccgGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCAt                         <  1:1236423/93‑1 (MQ=255)
                                                           ccgccgGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCAt                         <  2:4948670/93‑1 (MQ=255)
                                                            cgccgGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATc                        <  2:4809066/93‑1 (MQ=255)
                                                            cgccgGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATc                        <  2:1092249/93‑1 (MQ=255)
                                                             gccgGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCg                       >  1:5480176/1‑93 (MQ=255)
                                                             gccgGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCCACGGCAAGAATGATGGCATCg                       >  2:4178556/1‑93 (MQ=255)
                                                              ccgGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGt                      >  1:3914510/1‑93 (MQ=255)
                                                              ccgGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGt                      >  1:5170662/1‑93 (MQ=255)
                                                              ccgGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGt                      >  1:6142363/1‑93 (MQ=255)
                                                               cgGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATATGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTa                     >  1:6641389/1‑93 (MQ=255)
                                                               cgGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTa                     <  2:5951644/93‑1 (MQ=255)
                                                                 gAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTAtt                   <  1:5933610/93‑1 (MQ=255)
                                                                 gAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTAtt                   >  2:4691202/1‑93 (MQ=255)
                                                                  aaGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTc                  >  2:3677318/1‑93 (MQ=255)
                                                                  aaGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTc                  >  1:1485245/1‑93 (MQ=255)
                                                                   agagAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCg                 >  2:721057/1‑93 (MQ=255)
                                                                    gagaAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGc                <  1:3350434/93‑1 (MQ=255)
                                                                    gagaAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGc                >  1:884748/1‑93 (MQ=255)
                                                                       aaaTCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGt             >  1:6104981/1‑93 (MQ=255)
                                                                        aaTCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGtt            <  2:5643139/93‑1 (MQ=255)
                                                                          tctcGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTAt          <  2:6610356/93‑1 (MQ=255)
                                                                             cGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCa       <  1:3813581/93‑1 (MQ=255)
                                                                              gggTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAg      <  2:1022150/93‑1 (MQ=255)
                                                                              gggTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAg      >  2:169621/1‑93 (MQ=255)
                                                                               ggTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGt     <  1:235125/93‑1 (MQ=255)
                                                                               ggTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGt     >  2:2522313/1‑93 (MQ=255)
                                                                                gTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGtt    <  2:3516463/93‑1 (MQ=255)
                                                                                gTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGtt    >  1:1077634/1‑93 (MQ=255)
                                                                                gTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGtt    <  1:3858548/93‑1 (MQ=255)
                                                                                gTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGtt    <  1:4280011/93‑1 (MQ=255)
                                                                                gTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGtt    <  1:6227152/93‑1 (MQ=255)
                                                                                 tGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGttt   <  2:2942533/93‑1 (MQ=255)
                                                                                  ggCAACGCGGCGCGAATACGTGCCTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTg  <  1:2515298/93‑1 (MQ=255)
                                                                                  ggCAACGCGGCGCGAATACGTGCCTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTg  <  1:457218/93‑1 (MQ=255)
                                                                                  ggCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTg  >  1:6019344/1‑93 (MQ=255)
                                                                                       ||                                                                                      
CAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAAGCCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTG  >  NC_000913/3987735‑3987909

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAAGCCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCT  >  NC_000913/3987729‑3987911
                                                                                             |                                                                                         
GGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAAC                                                                                           >  SRR4302547.214196/1‑94 (MQ=60)
GGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAAC                                                                                           >  SRR4302547.2091553/1‑94 (MQ=60)
GGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAAA                                                                                           <  SRR4302547.3279782/94‑1 (MQ=60)
GGCAAGCAGCTCGCGAGGGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAAC                                                                                           <  SRR4302547.2975203/94‑1 (MQ=60)
     GCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGC                                                                                      <  SRR4302547.1962269/94‑1 (MQ=60)
      CAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCG                                                                                     <  SRR4302547.5472738/94‑1 (MQ=60)
      CAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCG                                                                                     <  SRR4302547.6580525/94‑1 (MQ=60)
       AGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGC                                                                                    <  SRR4302547.4381807/94‑1 (MQ=60)
       AGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGC                                                                                    >  SRR4302547.677849/1‑94 (MQ=60)
       AGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGC                                                                                    <  SRR4302547.2949230/94‑1 (MQ=60)
        GCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCG                                                                                   <  SRR4302547.1874339/94‑1 (MQ=60)
        GCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCG                                                                                   <  SRR4302547.3580903/94‑1 (MQ=60)
        GCTCGCGAGTGCGTGAAGCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCG                                                                                   <  SRR4302547.4379780/94‑1 (MQ=60)
         CTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTGCCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGA                                                                                  <  SRR4302547.6117542/94‑1 (MQ=60)
         CTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGA                                                                                  <  SRR4302547.460905/94‑1 (MQ=60)
          TCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAA                                                                                 >  SRR4302547.2599282/1‑94 (MQ=60)
           CGCGAGTGCGTGAATCATCAAGGCGGCATTCAATCCCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAAT                                                                                <  SRR4302547.1076822/94‑1 (MQ=60)
           CGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAAT                                                                                <  SRR4302547.1880993/94‑1 (MQ=60)
            ccgcGTGCGTGAATCATCAAGGCGGCATCCAATCCCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATA                                                                               <  SRR4302547.773942/90‑1 (MQ=60)
            GCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATA                                                                               >  SRR4302547.3750404/1‑94 (MQ=60)
             CGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATAC                                                                              <  SRR4302547.2398475/94‑1 (MQ=60)
             CGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATAC                                                                              >  SRR4302547.3904930/1‑94 (MQ=60)
             CGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATAC                                                                              >  SRR4302547.3493242/1‑94 (MQ=60)
             CGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATAC                                                                              >  SRR4302547.6347428/1‑94 (MQ=60)
              GAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACG                                                                             <  SRR4302547.1692806/94‑1 (MQ=60)
              GAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACG                                                                             <  SRR4302547.2663021/94‑1 (MQ=60)
              GAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACG                                                                             >  SRR4302547.2999073/1‑94 (MQ=60)
              GAGTGCGTGAAACATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACG                                                                             <  SRR4302547.3708680/94‑1 (MQ=60)
               AGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGT                                                                            <  SRR4302547.2393971/94‑1 (MQ=60)
               AGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGT                                                                            <  SRR4302547.6259083/94‑1 (MQ=60)
                GTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTG                                                                           <  SRR4302547.6765955/94‑1 (MQ=60)
                   CGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACT                                                                        >  SRR4302547.3131488/1‑94 (MQ=60)
                   CGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACT                                                                        <  SRR4302547.888135/94‑1 (MQ=60)
                   CGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACT                                                                        <  SRR4302547.6503523/94‑1 (MQ=60)
                   CGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACT                                                                        <  SRR4302547.3253194/94‑1 (MQ=60)
                   CGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACT                                                                        >  SRR4302547.1418137/1‑94 (MQ=60)
                   CGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACT                                                                        >  SRR4302547.3656842/1‑94 (MQ=60)
                    GTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTC                                                                       >  SRR4302547.6117352/1‑94 (MQ=60)
                    GTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTC                                                                       <  SRR4302547.3820264/94‑1 (MQ=60)
                     TGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCC                                                                      >  SRR4302547.542884/1‑94 (MQ=60)
                     TGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCC                                                                      <  SRR4302547.2181833/94‑1 (MQ=60)
                      GAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCA                                                                     >  SRR4302547.4324317/1‑94 (MQ=60)
                       AATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAG                                                                    <  SRR4302547.2138635/94‑1 (MQ=60)
                       AATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAG                                                                    >  SRR4302547.2317336/1‑94 (MQ=60)
                        ATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGA                                                                   <  SRR4302547.1139143/94‑1 (MQ=60)
                          CATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACC                                                                 >  SRR4302547.6566519/1‑94 (MQ=60)
                          CATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACC                                                                 <  SRR4302547.3212625/94‑1 (MQ=60)
                          CATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACC                                                                 <  SRR4302547.664532/94‑1 (MQ=60)
                          CATCAAGGCGGCATTCAATACCAACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACC                                                                 >  SRR4302547.1119927/1‑94 (MQ=60)
                           ATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCTGCGCGAATACGTGACTCCAGACCT                                                                <  SRR4302547.5451925/94‑1 (MQ=60)
                           ATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCT                                                                >  SRR4302547.1833211/1‑94 (MQ=60)
                           ATCAAGGCGGCATCCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCT                                                                >  SRR4302547.2035557/1‑94 (MQ=60)
                            TCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTA                                                               <  SRR4302547.4406470/94‑1 (MQ=60)
                              AAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAA                                                             <  SRR4302547.5998274/94‑1 (MQ=60)
                               AGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAAC                                                            <  SRR4302547.5226634/94‑1 (MQ=60)
                                GGCGGCATTCAATCCCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACG                                                           <  SRR4302547.1734653/94‑1 (MQ=60)
                                GGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACG                                                           <  SRR4302547.2251531/94‑1 (MQ=60)
                                 GCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGT                                                          >  SRR4302547.3409820/1‑94 (MQ=60)
                                  CGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTT                                                         >  SRR4302547.6378475/1‑94 (MQ=60)
                                  CGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTT                                                         <  SRR4302547.5664232/94‑1 (MQ=60)
                                   GGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTT                                                        <  SRR4302547.5401766/94‑1 (MQ=60)
                                   GGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTT                                                        <  SRR4302547.2742419/94‑1 (MQ=60)
                                    GCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCGCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTT                                                       <  SRR4302547.5328437/94‑1 (MQ=60)
                                     CATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATANGTGACTCCAGACCTAAACGTTTTA                                                      >  SRR4302547.2408284/1‑94 (MQ=60)
                                     CATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTA                                                      <  SRR4302547.5435028/94‑1 (MQ=60)
                                      ATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAG                                                     <  SRR4302547.416728/94‑1 (MQ=60)
                                      ATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAG                                                     >  SRR4302547.4772890/1‑94 (MQ=60)
                                      ATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAG                                                     <  SRR4302547.2880293/94‑1 (MQ=60)
                                       TTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGT                                                    <  SRR4302547.3700726/94‑1 (MQ=60)
                                       TTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGT                                                    >  SRR4302547.713387/1‑94 (MQ=60)
                                       TTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGT                                                    <  SRR4302547.6372490/94‑1 (MQ=60)
                                        TCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTC                                                   <  SRR4302547.5194355/94‑1 (MQ=60)
                                         CAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCC                                                  >  SRR4302547.1127891/1‑94 (MQ=60)
                                            TACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGC                                               >  SRR4302547.3460468/1‑94 (MQ=60)
                                              CCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTA                                             >  SRR4302547.2439943/1‑94 (MQ=60)
                                               CCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTAC                                            <  SRR4302547.4530191/94‑1 (MQ=60)
                                               CCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTAC                                            >  SRR4302547.4516862/1‑94 (MQ=60)
                                                CACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACG                                           >  SRR4302547.3673825/1‑94 (MQ=60)
                                                CACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACG                                           <  SRR4302547.3649564/94‑1 (MQ=60)
                                                CACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACG                                           <  SRR4302547.6474872/94‑1 (MQ=60)
                                                CACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACG                                           <  SRR4302547.6268172/94‑1 (MQ=60)
                                                 ACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGG                                          >  SRR4302547.1025152/1‑94 (MQ=60)
                                                  CCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGC                                         >  SRR4302547.2522112/1‑94 (MQ=60)
                                                   CGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCA                                        >  SRR4302547.4524174/1‑94 (MQ=60)
                                                    GCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAA                                       <  SRR4302547.2578560/94‑1 (MQ=60)
                                                     CACCTTGTGCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAG                                      >  SRR4302547.5483878/1‑94 (MQ=60)
                                                     CACCTTGTCCTGCCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAG                                      <  SRR4302547.3950630/94‑1 (MQ=60)
                                                        CTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGGGAATACGTGACTGCAGACCTAAACGTTTTAGTCCGGCTACGGCGAGAAT                                   >  SRR4302547.4498915/1‑94 (MQ=60)
                                                        CTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAAT                                   >  SRR4302547.4997519/1‑94 (MQ=60)
                                                        CTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAAT                                   >  SRR4302547.6412252/1‑94 (MQ=60)
                                                        CTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAAT                                   >  SRR4302547.3957579/1‑94 (MQ=60)
                                                         TTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATG                                  <  SRR4302547.5214739/94‑1 (MQ=60)
                                                         TTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATG                                  <  SRR4302547.5274794/94‑1 (MQ=60)
                                                         TTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATA                                  >  SRR4302547.6206242/1‑94 (MQ=60)
                                                          TGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGA                                 >  SRR4302547.3800995/1‑94 (MQ=60)
                                                          TGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGA                                 <  SRR4302547.537314/94‑1 (MQ=60)
                                                          TGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGA                                 >  SRR4302547.4536690/1‑94 (MQ=60)
                                                           GTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGAT                                <  SRR4302547.2089878/94‑1 (MQ=60)
                                                           GTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGAT                                >  SRR4302547.5268982/1‑94 (MQ=60)
                                                            TCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATG                               <  SRR4302547.6363768/94‑1 (MQ=60)
                                                            TCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATG                               >  SRR4302547.2138802/1‑94 (MQ=60)
                                                             CCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGG                              >  SRR4302547.2670621/1‑94 (MQ=60)
                                                             CCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGG                              <  SRR4302547.5739081/94‑1 (MQ=60)
                                                             CCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGG                              >  SRR4302547.4895936/1‑94 (MQ=60)
                                                             CCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGG                              <  SRR4302547.1979436/94‑1 (MQ=60)
                                                             CCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGG                              >  SRR4302547.3282064/1‑94 (MQ=60)
                                                             CCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGG                              <  SRR4302547.1819981/94‑1 (MQ=60)
                                                              CTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGC                             >  SRR4302547.4486825/1‑94 (MQ=60)
                                                              CTACCGCCGGAAGAGAAATCTCGGGGGGCAACTCGGCGCGAATACGAGACTCCAGACCTAAACGTTTTAGCCCAGCTAGCGCAAGAATTATGGC                             >  SRR4302547.3599126/1‑94 (MQ=55)
                                                                 CCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATC                          <  SRR4302547.3457000/94‑1 (MQ=60)
                                                                 CCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATC                          <  SRR4302547.1261373/94‑1 (MQ=60)
                                                                 CCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATC                          <  SRR4302547.1569021/94‑1 (MQ=60)
                                                                 CCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATC                          >  SRR4302547.3927025/1‑94 (MQ=60)
                                                                 CCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATC                          >  SRR4302547.11068/1‑94 (MQ=60)
                                                                  CGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCG                         <  SRR4302547.4761659/94‑1 (MQ=60)
                                                                  CGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCG                         >  SRR4302547.1711308/1‑94 (MQ=60)
                                                                  CGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCG                         >  SRR4302547.5593338/1‑94 (MQ=60)
                                                                   GCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGT                        >  SRR4302547.5277116/1‑94 (MQ=60)
                                                                   GCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGT                        >  SRR4302547.6269808/1‑94 (MQ=60)
                                                                   GCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGT                        >  SRR4302547.3994104/1‑94 (MQ=60)
                                                                   GCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGT                        >  SRR4302547.1900854/1‑94 (MQ=60)
                                                                    CCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATATGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTA                       >  SRR4302547.6779641/1‑94 (MQ=60)
                                                                    CCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTA                       <  SRR4302547.1434368/94‑1 (MQ=60)
                                                                     CGGAAGAGAAATCTCGGGTGGCACCGCGGCGCGAATACGTGACTCCAGACCTAACCGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTAT                      <  SRR4302547.1695049/94‑1 (MQ=60)
                                                                     CGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGAGTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAGGCATGATGGCATCGTAT                      >  SRR4302547.1085213/1‑94 (MQ=60)
                                                                     CGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTAT                      >  SRR4302547.1586257/1‑94 (MQ=60)
                                                                      GGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATT                     >  SRR4302547.3351967/1‑94 (MQ=60)
                                                                      GGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATT                     <  SRR4302547.4846948/94‑1 (MQ=60)
                                                                      GGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATT                     >  SRR4302547.3660729/1‑94 (MQ=60)
                                                                       GAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTC                    <  SRR4302547.907391/94‑1 (MQ=60)
                                                                       GAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTC                    <  SRR4302547.6056645/94‑1 (MQ=60)
                                                                       GAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTC                    >  SRR4302547.1514879/1‑94 (MQ=60)
                                                                        AAGAGAAATCTCGGGTTGCAACGCGGCGCGAATACGGGACTCCTGACCTAAACGTTTTAGTCCGGCTACGGCAAGTATGATGGCATCGCATTCG                   >  SRR4302547.786342/1‑94 (MQ=60)
                                                                        AAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCG                   >  SRR4302547.1876383/1‑94 (MQ=60)
                                                                         AGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGC                  >  SRR4302547.901749/1‑94 (MQ=60)
                                                                          GAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCC                 <  SRR4302547.3418369/94‑1 (MQ=60)
                                                                          GAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCC                 <  SRR4302547.6388589/94‑1 (MQ=60)
                                                                          GAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCC                 <  SRR4302547.3644709/94‑1 (MQ=60)
                                                                           AGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCG                <  SRR4302547.907469/94‑1 (MQ=60)
                                                                           AGAAATCTCGGGTGGCAACGCGGCGCGAATACATGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCG                >  SRR4302547.3503457/1‑94 (MQ=60)
                                                                            GAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGT               >  SRR4302547.6231618/1‑94 (MQ=60)
                                                                              AATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGACGTTA             >  SRR4302547.1929106/1‑94 (MQ=60)
                                                                               ATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTAT            <  SRR4302547.5706407/94‑1 (MQ=60)
                                                                                   CGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAG        <  SRR4302547.3891294/94‑1 (MQ=60)
                                                                                    GGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGT       >  SRR4302547.1217461/1‑94 (MQ=60)
                                                                                    GGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGT       >  SRR4302547.2477919/1‑94 (MQ=60)
                                                                                     GGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTT      <  SRR4302547.239542/94‑1 (MQ=60)
                                                                                     GGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTT      >  SRR4302547.1099139/1‑94 (MQ=60)
                                                                                     GGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTT      <  SRR4302547.4616191/94‑1 (MQ=60)
                                                                                     GGTGGCAACGCGCCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTT      <  SRR4302547.2391542/94‑1 (MQ=60)
                                                                                      GTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTT     <  SRR4302547.4367919/94‑1 (MQ=60)
                                                                                      GTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTT     <  SRR4302547.3937074/94‑1 (MQ=60)
                                                                                      GTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTT     <  SRR4302547.4755011/94‑1 (MQ=60)
                                                                                      GTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTT     <  SRR4302547.6356480/94‑1 (MQ=60)
                                                                                       TGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTG    >  SRR4302547.6144269/1‑94 (MQ=60)
                                                                                       TGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTG    <  SRR4302547.2404421/94‑1 (MQ=60)
                                                                                        GGCAACGCGGCGCGAATACGTGCCTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGC   <  SRR4302547.465677/94‑1 (MQ=60)
                                                                                        GGCAACGCGGCGCGAATACGTGCCTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGC   <  SRR4302547.2566152/94‑1 (MQ=60)
                                                                                         GCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCT  <  SRR4302547.5776818/94‑1 (MQ=60)
                                                                                             |                                                                                         
GGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAAGCCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCT  >  NC_000913/3987729‑3987911

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