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 F394 I0 R1 1342 215.0 7134094 96.2% 6862998 147.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,825,763 C→A 100% A159A (GCG→GCT astC ← succinylornithine transaminase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,825,7630CA100.0% 498.7 / ‑6.7 153A159A (GCG→GCTastCsuccinylornithine transaminase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (85/67);  total (86/67)
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 153 total aligned reads displayed.

TGTAAAAACGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGCGCC  >  CP009273/1825620‑1825907
                                                                                                                                               |                                                                                                                                                
tGTAAAAACGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCaa                                                                                                                                              >  1:1242896/1‑148 (MQ=255)
tGTAAAAACGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCaa                                                                                                                                              >  1:3250309/1‑148 (MQ=255)
 gTAAAAACGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCaaa                                                                                                                                             <  1:1011542/148‑1 (MQ=255)
   aaaaaCGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAAt                                                                                                                                           <  2:3559242/148‑1 (MQ=255)
   aaaaaCGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAAt                                                                                                                                           >  1:3100291/1‑148 (MQ=255)
   aaaaaCGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGGTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAAt                                                                                                                                           >  2:3382456/1‑148 (MQ=255)
    aaaaCGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATc                                                                                                                                          <  1:2870220/148‑1 (MQ=255)
       aCGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCCATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCtg                                                                                                                                       >  1:2104842/1‑148 (MQ=255)
       aCGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCCTCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCtg                                                                                                                                       >  1:1627252/1‑148 (MQ=255)
        cgcgTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCtgt                                                                                                                                      >  2:283036/1‑148 (MQ=255)
        cgcgTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCtgt                                                                                                                                      >  2:2045358/1‑148 (MQ=255)
         gcgTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCCGCGCGCTGGCCGAGTTAATATCGGTATATGGCGCATGACGCATCTCCGCCGGCAGGGGCGGGCAATCCtggg                                                                                                                                     >  1:2418178/1‑146 (MQ=255)
         gcgTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCtgtg                                                                                                                                     >  1:627682/1‑148 (MQ=255)
            ttGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAAt                                                                                                                                  <  1:3011365/148‑1 (MQ=255)
             tGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATa                                                                                                                                 <  2:1536134/148‑1 (MQ=255)
                tggctggcACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAGTATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAggc                                                                                                                              <  2:409373/148‑1 (MQ=255)
                tggctggcACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAggc                                                                                                                              <  1:2924483/148‑1 (MQ=255)
                tggctggcACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAggc                                                                                                                              <  1:2016809/148‑1 (MQ=255)
                tggctggcACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAggc                                                                                                                              <  1:418584/148‑1 (MQ=255)
                tggctggcACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAggc                                                                                                                              <  1:217431/148‑1 (MQ=255)
                tggctggcACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAggc                                                                                                                              <  1:3398612/148‑1 (MQ=255)
                  gctggcACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAggctg                                                                                                                            >  2:1054525/1‑148 (MQ=255)
                   ctggcACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAggctgg                                                                                                                           <  2:2566422/148‑1 (MQ=255)
                   ctggcACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAggctgg                                                                                                                           <  1:262654/148‑1 (MQ=255)
                      gcACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCCGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGgctggctg                                                                                                                        >  2:3114492/1‑148 (MQ=255)
                         ccacACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGccc                                                                                                                     >  2:2490707/1‑148 (MQ=255)
                              ccGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCg                                                                                                                >  1:2323479/1‑148 (MQ=255)
                               cGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGc                                                                                                               >  2:2535407/1‑148 (MQ=255)
                                gCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCa                                                                                                              >  1:637772/1‑148 (MQ=255)
                                gCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCa                                                                                                              >  1:1749532/1‑148 (MQ=255)
                                  cTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACt                                                                                                            >  1:3421260/1‑148 (MQ=255)
                                  cTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACt                                                                                                            >  1:2892689/1‑148 (MQ=255)
                                   ttCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTg                                                                                                           >  2:1291417/1‑148 (MQ=255)
                                      ccccTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACa                                                                                                        >  2:3412651/1‑148 (MQ=255)
                                       cccTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAg                                                                                                       >  2:462188/1‑148 (MQ=255)
                                       cccTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAg                                                                                                       <  1:2217888/148‑1 (MQ=255)
                                       cccTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGGTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAAACCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAg                                                                                                       >  2:1370282/1‑148 (MQ=255)
                                            gATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAAc                                                                                                  <  1:43201/148‑1 (MQ=255)
                                            gATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAAc                                                                                                  <  1:2713357/148‑1 (MQ=255)
                                             aTGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACa                                                                                                 <  2:3100291/148‑1 (MQ=255)
                                                ggTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCg                                                                                              <  2:43771/148‑1 (MQ=255)
                                                   tCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGc                                                                                           <  1:1931217/148‑1 (MQ=255)
                                                        aaTCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCa                                                                                      <  2:746861/148‑1 (MQ=255)
                                                           cACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGa                                                                                   >  2:1394656/1‑148 (MQ=255)
                                                            aCCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGaa                                                                                  >  1:1853511/1‑148 (MQ=255)
                                                             ccGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGaaa                                                                                 >  2:2426621/1‑148 (MQ=255)
                                                                cacaGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGc                                                                              >  1:1025537/1‑148 (MQ=255)
                                                                  caGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCAt                                                                            <  2:1219163/148‑1 (MQ=255)
                                                                  caGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCAt                                                                            <  2:2995798/148‑1 (MQ=255)
                                                                  caGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCAt                                                                            <  2:627682/148‑1 (MQ=255)
                                                                    ggTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGAAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCAttt                                                                          >  2:2081511/1‑148 (MQ=255)
                                                                     gTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCAtttt                                                                         >  1:1740239/1‑148 (MQ=255)
                                                                           tcgtcAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACg                                                                   <  2:386198/148‑1 (MQ=255)
                                                                            cgtcAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGc                                                                  <  1:1882025/148‑1 (MQ=255)
                                                                            cgtcAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGc                                                                  <  1:2384685/148‑1 (MQ=255)
                                                                              tcAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCa                                                                >  1:488391/1‑148 (MQ=255)
                                                                              tcAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCa                                                                <  2:2253109/148‑1 (MQ=255)
                                                                                aaTCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACTAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACg                                                              <  2:283287/148‑1 (MQ=255)
                                                                                aaTCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACg                                                              >  2:565359/1‑148 (MQ=255)
                                                                                   cAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATg                                                           <  1:3203445/148‑1 (MQ=255)
                                                                                    aGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGc                                                          <  2:2398558/148‑1 (MQ=255)
                                                                                      cgcgcTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCg                                                        <  2:1737451/148‑1 (MQ=255)
                                                                                          cTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGctct                                                    >  1:1246123/1‑148 (MQ=255)
                                                                                           tGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCtt                                                   <  2:1242896/148‑1 (MQ=255)
                                                                                           tGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCtt                                                   <  2:735595/148‑1 (MQ=255)
                                                                                           tGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCtt                                                   <  2:3020806/148‑1 (MQ=255)
                                                                                            ggCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTa                                                  >  2:3137603/1‑148 (MQ=255)
                                                                                              cAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATg                                                >  1:1786202/1‑148 (MQ=255)
                                                                                              cAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATg                                                >  1:969381/1‑148 (MQ=255)
                                                                                               agagTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATgg                                               <  2:3542824/148‑1 (MQ=255)
                                                                                                gagTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGc                                              >  1:698414/1‑148 (MQ=255)
                                                                                                gagTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGc                                              >  1:1681141/1‑148 (MQ=255)
                                                                                                 agTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGAt                                                            >  1:1317431/1‑133 (MQ=255)
                                                                                                        atCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAg                                      <  2:3080439/148‑1 (MQ=255)
                                                                                                          cGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCg                                    >  2:1048894/1‑148 (MQ=255)
                                                                                                             tataTGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTc                                 >  2:68195/1‑148 (MQ=255)
                                                                                                             tataTGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTc                                 >  1:1439484/1‑148 (MQ=255)
                                                                                                              atatGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCg                                <  2:2672613/148‑1 (MQ=255)
                                                                                                              atatGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCg                                <  2:1791122/148‑1 (MQ=255)
                                                                                                               tatGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGt                               >  2:855660/1‑148 (MQ=255)
                                                                                                               tatGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGt                               <  1:3171634/148‑1 (MQ=255)
                                                                                                                atGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTg                              <  2:2394513/148‑1 (MQ=255)
                                                                                                                atGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTg                              >  2:230098/1‑148 (MQ=255)
                                                                                                                   cagcaTGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGc                           >  1:1153408/1‑148 (MQ=255)
                                                                                                                     gcaTGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGa                         >  2:671307/1‑148 (MQ=255)
                                                                                                                     gcaTGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTCACAGTAAACAGCGTGCGACCATGAAACGCAGGTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGa                         >  2:2691232/1‑148 (MQ=255)
                                                                                                                             aaTATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTAcgc                 >  2:149937/1‑148 (MQ=255)
                                                                                                                             aaTATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTAcgc                 >  1:3560416/1‑148 (MQ=255)
                                                                                                                              atatCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTAcgcg                <  1:1292478/148‑1 (MQ=255)
                                                                                                                               tatCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACgcgc               >  1:2059690/1‑148 (MQ=255)
                                                                                                                               tatCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACgcgc               >  1:2244387/1‑148 (MQ=255)
                                                                                                                                 tCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCa             >  2:1574521/1‑148 (MQ=255)
                                                                                                                                 tCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCa             >  1:465133/1‑148 (MQ=255)
                                                                                                                                  ccgccgGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAg            <  2:1913574/148‑1 (MQ=255)
                                                                                                                                    gccgGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGtt          <  1:2490707/148‑1 (MQ=255)
                                                                                                                                     ccgGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGttt         <  1:1985521/148‑1 (MQ=255)
                                                                                                                                     ccgGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGttt         <  1:3497235/148‑1 (MQ=255)
                                                                                                                                       ggCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTa       >  1:2794387/1‑148 (MQ=255)
                                                                                                                                         cAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAgc     >  2:2161039/1‑148 (MQ=255)
                                                                                                                                            tGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGcgcc  >  2:3309096/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                                
TGTAAAAACGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGCGCC  >  CP009273/1825620‑1825907

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
GACCTTGTAAAAACGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGCGCCGCTTC  >  CP009273/1825615‑1825912
                                                                                                                                                    |                                                                                                                                                     
GACCTTGTAAAAACGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAG                                                                                                                                                      >  SRR19973066.1388077/1‑150 (MQ=60)
   CTTGTAAAAACGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAA                                                                                                                                                   >  SRR19973066.1251795/1‑150 (MQ=60)
   CTTGTAAAAACGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAA                                                                                                                                                   >  SRR19973066.3040957/1‑150 (MQ=60)
   CTTGTAAAAACGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAA                                                                                                                                                   >  SRR19973066.3271724/1‑150 (MQ=60)
      GTAAAAACGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAAT                                                                                                                                                <  SRR19973066.1018996/150‑1 (MQ=60)
      GTAAAAACGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAAT                                                                                                                                                >  SRR19973066.3120998/1‑150 (MQ=60)
         AAAACGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCT                                                                                                                                             <  SRR19973066.2889468/150‑1 (MQ=60)
         AAAACGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCT                                                                                                                                             <  SRR19973066.2889481/150‑1 (MQ=60)
          AAACGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCCATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTG                                                                                                                                            >  SRR19973066.2119347/1‑150 (MQ=60)
          AAACGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCCTCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTG                                                                                                                                            >  SRR19973066.1639054/1‑150 (MQ=60)
            ACGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCCGCGCGCTGGCCGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCCAAATCCTGTG                                                                                                                                          >  SRR19973066.632179/1‑150 (MQ=60)
                 TTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAG                                                                                                                                     <  SRR19973066.3031416/150‑1 (MQ=60)
                     TGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTG                                                                                                                                 <  SRR19973066.3421302/150‑1 (MQ=60)
                     TGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTG                                                                                                                                 <  SRR19973066.219215/150‑1 (MQ=60)
                     TGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTG                                                                                                                                 <  SRR19973066.421846/150‑1 (MQ=60)
                     TGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTG                                                                                                                                 <  SRR19973066.2943967/150‑1 (MQ=60)
                     TGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTG                                                                                                                                 <  SRR19973066.2030930/150‑1 (MQ=60)
                     TGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTG                                                                                                                                 <  SRR19973066.626079/150‑1 (MQ=60)
                        CTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCT                                                                                                                              <  SRR19973066.264554/150‑1 (MQ=60)
                        CTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCT                                                                                                                              <  SRR19973066.264625/150‑1 (MQ=60)
                                 CACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCG                                                                                                                     >  SRR19973066.2339347/1‑150 (MQ=60)
                                   CCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCA                                                                                                                   >  SRR19973066.642371/1‑150 (MQ=60)
                                   CCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCA                                                                                                                   >  SRR19973066.1762119/1‑150 (MQ=60)
                                     GCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACT                                                                                                                 >  SRR19973066.3444141/1‑150 (MQ=60)
                                     GCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACT                                                                                                                 >  SRR19973066.2912026/1‑150 (MQ=60)
                                     GCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACT                                                                                                                 >  SRR19973066.3444972/1‑150 (MQ=60)
                                            CCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTA                                                                                                          <  SRR19973066.2233236/150‑1 (MQ=60)
                                             CCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCCGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAAAATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAA                                                                                                         >  SRR19973066.2414738/1‑150 (MQ=60)
                                               TGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAAC                                                                                                       <  SRR19973066.2214777/150‑1 (MQ=60)
                                                 GATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAG                                                                                                     <  SRR19973066.43776/150‑1 (MQ=60)
                                                 GATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAG                                                                                                     <  SRR19973066.2731582/150‑1 (MQ=60)
                                                    GGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGT                                                                                                  <  SRR19973066.1634044/150‑1 (MQ=60)
                                                      GTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGC                                                                                                <  SRR19973066.1693860/150‑1 (MQ=60)
                                                        TCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGA                                                                                              <  SRR19973066.1944929/150‑1 (MQ=60)
                                                               TCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAA                                                                                       >  SRR19973066.1867288/1‑150 (MQ=60)
                                                               TCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAA                                                                                       >  SRR19973066.1866776/1‑150 (MQ=60)
                                                                   CGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGC                                                                                   >  SRR19973066.1033088/1‑150 (MQ=60)
                                                                     CACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCAT                                                                                 >  SRR19973066.2282749/1‑150 (MQ=60)
                                                                        AGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACAGACAGTAAACAGCGTGCGACCATGAAACGCATTTT                                                                              >  SRR19973066.1752790/1‑150 (MQ=60)
                                                                         GGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTT                                                                             >  SRR19973066.1137022/1‑150 (MQ=60)
                                                                                 CGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCA                                                                     <  SRR19973066.2400760/150‑1 (MQ=60)
                                                                                 CGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCA                                                                     >  SRR19973066.67714/1‑150 (MQ=60)
                                                                                 CGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCA                                                                     >  SRR19973066.492052/1‑150 (MQ=60)
                                                                                 CGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCA                                                                     <  SRR19973066.1895457/150‑1 (MQ=60)
                                                                                        CAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCC                                                              <  SRR19973066.3224704/150‑1 (MQ=60)
                                                                                             CGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCT                                                         >  SRR19973066.1254446/1‑150 (MQ=60)
                                                                                             CGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCT                                                         >  SRR19973066.1255060/1‑150 (MQ=60)
                                                                                                 GGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATG                                                     >  SRR19973066.976512/1‑150 (MQ=60)
                                                                                                 GGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATG                                                     <  SRR19973066.1787685/150‑1 (MQ=60)
                                                                                                 GGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATG                                                     >  SRR19973066.1798943/1‑150 (MQ=60)
                                                                                                   CAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGC                                                   >  SRR19973066.1693382/1‑150 (MQ=60)
                                                                                                   CAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGC                                                   >  SRR19973066.703566/1‑150 (MQ=60)
                                                                                                    AGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCagatcggaagagc                                                  >  SRR19973066.1326977/1‑137 (MQ=60)
                                                                                                                GTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTC                                      >  SRR19973066.1449733/1‑150 (MQ=60)
                                                                                                                    TATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGA                                  <  SRR19973066.3192752/150‑1 (MQ=60)
                                                                                                                      TGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGC                                >  SRR19973066.2952171/1‑150 (MQ=60)
                                                                                                                      TGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGC                                >  SRR19973066.1161620/1‑150 (MQ=60)
                                                                                                                          GCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAAT                            <  SRR19973066.2779499/150‑1 (MQ=60)
                                                                                                                           CATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATT                           >  SRR19973066.1984412/1‑150 (MQ=60)
                                                                                                                           CATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATT                           <  SRR19973066.465723/150‑1 (MQ=60)
                                                                                                                            ATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTT                          >  SRR19973066.202498/1‑150 (MQ=60)
                                                                                                                                CGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGC                      >  SRR19973066.3584027/1‑150 (MQ=60)
                                                                                                                                  AATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGC                    >  SRR19973066.2259866/1‑150 (MQ=60)
                                                                                                                                  AATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGC                    >  SRR19973066.2074119/1‑150 (MQ=60)
                                                                                                                                  AATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGC                    >  SRR19973066.2073976/1‑150 (MQ=60)
                                                                                                                                  AATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGC                    >  SRR19973066.383880/1‑150 (MQ=60)
                                                                                                                                   ATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCC                   <  SRR19973066.1301836/150‑1 (MQ=60)
                                                                                                                                    TATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCA                  >  SRR19973066.468681/1‑150 (MQ=60)
                                                                                                                                    TATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCA                  >  SRR19973066.1326970/1‑150 (MQ=60)
                                                                                                                                         GCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTT             <  SRR19973066.2507697/150‑1 (MQ=60)
                                                                                                                                          CCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTA            <  SRR19973066.1999520/150‑1 (MQ=60)
                                                                                                                                          CCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTA            <  SRR19973066.3520550/150‑1 (MQ=60)
                                                                                                                                          CCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTA            <  SRR19973066.3215477/150‑1 (MQ=60)
                                                                                                                                          CCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTA            >  SRR19973066.2813191/1‑150 (MQ=60)
                                                                                                                                           CGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAG           <  SRR19973066.3405020/150‑1 (MQ=60)
                                                                                                                                                    AGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGCGCCGCTTC  <  SRR19973066.1300766/150‑1 (MQ=60)
                                                                                                                                                    AGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGCGCCGCTTC  >  SRR19973066.1141403/1‑150 (MQ=60)
                                                                                                                                                    |                                                                                                                                                     
GACCTTGTAAAAACGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGCGCCGCTTC  >  CP009273/1825615‑1825912

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: