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
A5 F25 I1 R1 61 74.4 3161000 98.6% 3116746 146.7

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947_CJ‑RC 5,677,566 +G intergenic (+113/+27) rhlE‑II → / ← yceI ATP‑dependent RNA helicase/polyisoprene‑binding acidic stress response factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947_CJ‑RC5,677,5651.G98.5% 238.0 / ‑3.2 68intergenic (+112/+28)rhlE‑II/yceIATP‑dependent RNA helicase/polyisoprene‑binding acidic stress response factor
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (41/26);  total (42/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

AGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGCGCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCA  >  NC_002947_CJ‑RC/5677422‑5677707
                                                                                                                                                |                                                                                                                                               
aGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAgcggcc                                                                                                                                             <  2:56736/148‑4 (MQ=255)
               tgACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCt                                                                                                                             >  4:159555/1‑149 (MQ=255)
               tgACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCGCGGTTTTCAGCt                                                                                                                             >  8:270009/1‑149 (MQ=255)
                 aCCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCAGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTCTGTCGTCATAGGGCCGCAAAGCGGCCCCCTCGGTTATCAGCTgg                                                                                                                           >  8:135129/1‑149 (MQ=37)
                             tGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACtt                                                                                                               >  8:305684/1‑149 (MQ=255)
                             tGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACtt                                                                                                               >  8:195089/1‑149 (MQ=255)
                               aTAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGGAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCGCGGGTTTCAGCTGGCAGAAATTACTTCt                                                                                                             >  8:147218/1‑149 (MQ=37)
                               aTAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACtt                                                                                                               <  8:187665/147‑1 (MQ=255)
                                tAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCGCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTg                                                                                                            <  8:353400/149‑1 (MQ=255)
                                        ccGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCacac                                                                                                    <  7:195089/149‑1 (MQ=255)
                                             ggggCTGCTTCGCAGCCCTATCACTACACAAGGCCGCTCCTGCACAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTc                                                                                                <  3:159555/148‑1 (MQ=38)
                                                gctgctTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACCGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGaca                                                                                            <  4:371573/149‑1 (MQ=255)
                                                 ctgctTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGacac                                                                                           <  8:147413/149‑1 (MQ=255)
                                                  tgctTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCTGAAATTACTTCTGCTTCACACCTTCGCCACt                                                                                          >  7:330209/1‑149 (MQ=37)
                                                    ctTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAATGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTg                                                                                         <  3:40090/148‑1 (MQ=255)
                                                     ttCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGGTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCATCGACACTGa                                                                                        >  4:269083/1‑148 (MQ=37)
                                                         ctgCCCTATCACGACACACGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGGGGCCTTGTGTCGTGATAGGGCGGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATAtcc                                                                                   <  5:35550/147‑1 (MQ=38)
                                                         cAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATAtcc                                                                                   <  7:393541/149‑1 (MQ=255)
                                                          aGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATAtc                                                                                    >  2:108945/1‑147 (MQ=255)
                                                           gCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGGCCCCGGGGGTTTCAGCTGGGAGAAATTACTTCTGCTTGACCCCTTCGACCCGGATTtacag                                                                                 >  3:128658/1‑149 (MQ=25)
                                                            cccTATCACGACACAAGGCCGCTCCTGTAGAGATCAAGAAATTCTCTGCAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATAtccagg                                                                                >  6:323569/1‑149 (MQ=14)
                                                            cccTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATAtccagg                                                                                >  6:348245/1‑149 (MQ=255)
                                                            cccTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATAtccagg                                                                                >  5:177136/1‑149 (MQ=255)
                                                              cTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCGCGGGTTTGAGCTGGCAGAAATTACTTCTGCTTCACACCGTCGACACTGATTtccaggtc                                                                              >  1:159140/1‑149 (MQ=37)
                                                                   aCGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATAtccaggtccagg                                                                          <  3:307499/148‑1 (MQ=255)
                                                                     gACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCGCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTc                                                                       <  4:157331/149‑1 (MQ=255)
                                                                        acaAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGc                                                                    <  8:73666/149‑1 (MQ=255)
                                                                          aaGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCg                                                                   >  7:11871/1‑148 (MQ=255)
                                                                           aGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGaa                                                                 <  1:90849/149‑1 (MQ=255)
                                                                             gCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAg                                                                >  2:128304/1‑148 (MQ=255)
                                                                             gCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAg                                                                >  8:38789/1‑148 (MQ=255)
                                                                              ccGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTa                                                              <  6:367954/149‑1 (MQ=255)
                                                                              ccGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTa                                                              <  8:116977/149‑1 (MQ=255)
                                                                               cGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAg                                                             >  4:118546/1‑149 (MQ=255)
                                                                               cGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAg                                                             <  7:135129/149‑1 (MQ=255)
                                                                               cGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAg                                                             >  6:249297/1‑149 (MQ=255)
                                                                                gCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGc                                                            <  4:99555/149‑1 (MQ=255)
                                                                                     tGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGc                                                                    <  2:1011/136‑1 (MQ=255)
                                                                                     tGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGc                                                                    >  1:1011/1‑136 (MQ=255)
                                                                                     tGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCg                                                       >  2:340102/1‑149 (MQ=255)
                                                                                      gCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCgg                                                      >  8:348215/1‑149 (MQ=255)
                                                                                         gagaATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGcc                                                   >  5:347248/1‑149 (MQ=255)
                                                                                             aTTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCttt                                                <  7:305684/148‑1 (MQ=255)
                                                                                             aTTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCttt                                                <  7:348215/148‑1 (MQ=255)
                                                                                             aTTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTg                                               <  7:304953/149‑1 (MQ=255)
                                                                                               ttCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGGTTTAAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATACAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGAt                                             >  1:383477/1‑149 (MQ=38)
                                                                                                tCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGATTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATg                                            >  1:185220/1‑149 (MQ=255)
                                                                                                    gATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCg                                         >  3:297175/1‑148 (MQ=255)
                                                                                                      tctcTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGaa                                       <  7:270009/148‑1 (MQ=255)
                                                                                                        tctACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTc                                    >  6:196519/1‑149 (MQ=255)
                                                                                                        tctACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTc                                    >  6:22889/1‑149 (MQ=255)
                                                                                                          tACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCgt                                  >  7:214102/1‑149 (MQ=255)
                                                                                                          tACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACAGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCgt                                  >  2:361468/1‑149 (MQ=255)
                                                                                                            cAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCgtt                                 >  8:50585/1‑148 (MQ=255)
                                                                                                             aGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCCGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCACGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCgttca                               >  2:237413/1‑149 (MQ=39)
                                                                                                               gAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCgttcagg                             >  8:75604/1‑149 (MQ=255)
                                                                                                                       ttGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCgt                     >  7:56383/1‑149 (MQ=255)
                                                                                                                       ttGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGAGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCgt                     >  8:255022/1‑149 (MQ=255)
                                                                                                                              gTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCgtggtggtg              <  5:323569/149‑1 (MQ=255)
                                                                                                                               tGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGtggtggtgg             >  1:368142/1‑149 (MQ=255)
                                                                                                                               tGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCATCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGtggtggtgg             >  1:195062/1‑149 (MQ=255)
                                                                                                                                 aTAGGGCCGCAAAGCGGCCCCCTCGGGTTTCAGCTGGGAGAAATTACTTCTGCTTCACACCTTCGACACTGATATGCAGGGCCGGGGTCTGCGAAGTAGCACCGGGGCGGGTGGTGCCGGGGTCGTGCGGGTGCGGGGGGGTGGGGGCg           >  4:322792/1‑149 (MQ=255)
                                                                                                                                   aGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGtt         >  3:216199/1‑149 (MQ=255)
                                                                                                                                      gCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTg        >  5:337084/1‑147 (MQ=255)
                                                                                                                                      gCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGaa      <  4:359729/149‑1 (MQ=255)
                                                                                                                                      gCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGaa      <  6:337084/149‑1 (MQ=255)
                                                                                                                                        cGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGc    >  1:56332/1‑149 (MQ=255)
                                                                                                                                          cAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCa  <  6:177136/149‑1 (MQ=255)
                                                                                                                                          cAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTACACACCTTCGACACTCATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGCGCCTTTGATCCCGAAGTCGTACAGCTTCAGCGTGGTGGTGGCGTTGAAGcc   >  3:70380/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                               
AGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGCGCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCA  >  NC_002947_CJ‑RC/5677422‑5677707

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGCGCCCTC‑‑GGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCG  >  NC_002947_CJ‑RC/5677417‑5677716
                                                                                                                                                      |                                                                                                                                                        
CGGTGCGCCCCCAGCGCTGCTCACCCGCAAGGCCTGATAGACGAACCGAGGGGGCTGCTTCCCAGCCCTACCACGACACAAGCCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG                                                                                                                                                         <  NB501016:282:H2V7GAFX2:4:21601:12967:6263/151‑1 (MQ=60)
 GGTGAGCGCCCCGCGCGGCTGACCCCCACGGGCTTATAGACGTCCCGCGGGGGCGGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGgc                                                                                                                                                        <  NB501016:282:H2V7GAFX2:2:21212:2937:18447/151‑3 (MQ=60)
        GCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGCGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGgcccccgcg                                                                                                                                                 >  NB501016:282:H2V7GAFX2:2:11108:24889:12179/1‑142 (MQ=60)
                                            ACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTC‑‑GGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACC                                                                                                            <  NB501016:282:H2V7GAFX2:4:11611:14463:10711/151‑1 (MQ=60)
                                             CCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCGTGTGGCGTGGTAGGGCCGCACAGCGGCCCCCTC‑‑GGTGTTCAGCTGGgtgaaattagtggtgagtgtcacct                                                                                                           >  NB501016:282:H2V7GAFX2:1:21102:11959:7072/1‑126 (MQ=60)
                                                 GGGGGCTGATTCGCAGCCCTATCACTACACAAGGCCGCTCCTGCACAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTC‑‑GGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCG                                                                                                        <  NB501016:282:H2V7GAFX2:2:11306:13012:5434/150‑1 (MQ=60)
                                                      CTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTC‑‑GGTTTTCAGCTGGCTGAAATTACTTCTGCTTCACACCTTCGCCACTG                                                                                                  >  NB501016:282:H2V7GAFX2:4:21510:16124:15536/1‑151 (MQ=60)
                                                        GCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAATGCGGCCCCCTC‑‑GGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGA                                                                                                 <  NB501016:282:H2V7GAFX2:2:11108:1758:11174/150‑1 (MQ=60)
                                                             GCTGCCCTATCACGACACACGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGGGGCCTTGTGTCGTGATAGGGCGGCAAAGCGGCCCCCTC‑‑GGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCA                                                                                           <  NB501016:282:H2V7GAFX2:3:11407:10733:10636/151‑1 (MQ=60)
                                                             GCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTC‑‑GGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCA                                                                                           <  NB501016:282:H2V7GAFX2:4:21610:6309:4953/151‑1 (MQ=60)
                                                               AGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGGCCCCGGGGGTTTCAGCTGGGAGAAATTACTTCTGCTTGACCCCTTCGACccggatttacagg                                                                                         >  NB501016:282:H2V7GAFX2:2:11212:6905:15872/1‑138 (MQ=60)
                                                                GCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTC‑‑GGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGT                                                                                        >  NB501016:282:H2V7GAFX2:3:11608:6071:10225/1‑151 (MQ=60)
                                                                  CCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCGC‑‑GGGTTTGAGCTGGCAGAAATTACTTCTGCTTCACACCGTCGACACTGATTTCCAGGTCC                                                                                      >  NB501016:282:H2V7GAFX2:1:11305:14753:1083/1‑151 (MQ=60)
                                                                       CACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTC‑‑GGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGG                                                                                  <  NB501016:282:H2V7GAFX2:2:21211:22542:11094/150‑1 (MQ=60)
                                                                              CAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTC‑‑GGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGA                                                                           >  NB501016:282:H2V7GAFX2:4:11403:21698:1736/1‑150 (MQ=60)
                                                                               AAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTC‑‑GGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAG                                                                         <  NB501016:282:H2V7GAFX2:1:11204:24207:12407/151‑1 (MQ=60)
                                                                                   CCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTC‑‑GGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGC                                                                     <  NB501016:282:H2V7GAFX2:4:11512:13274:10164/151‑1 (MQ=60)
                                                                                         CTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTC‑‑GGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCG                                                                            >  NB501016:282:H2V7GAFX2:1:11101:10642:4510/1‑138 (MQ=60)
                                                                                             AGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTC‑‑GGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCT                                                           >  NB501016:282:H2V7GAFX2:3:21603:22971:14833/1‑151 (MQ=60)
                                                                                                 AATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTC‑‑GGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGA                                                       <  NB501016:282:H2V7GAFX2:4:21506:14584:6749/151‑1 (MQ=60)
                                                                                                 AATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTC‑‑GGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTG                                                        <  NB501016:282:H2V7GAFX2:4:21601:16623:18488/150‑1 (MQ=60)
                                                                                                 AATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTC‑‑GGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTG                                                        <  NB501016:282:H2V7GAFX2:4:21506:4980:9181/150‑1 (MQ=60)
                                                                                                   TTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTC‑‑GGGTTTAAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATACAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATG                                                     >  NB501016:282:H2V7GAFX2:1:21311:24890:5181/1‑151 (MQ=60)
                                                                                                    TTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTC‑‑GATTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGC                                                    >  NB501016:282:H2V7GAFX2:1:11309:12196:19973/1‑151 (MQ=60)
                                                                                                        TGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTC‑‑GGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGA                                                 >  NB501016:282:H2V7GAFX2:2:21209:15569:12493/1‑150 (MQ=60)
                                                                                                          ATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTC‑‑GGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAG                                               <  NB501016:282:H2V7GAFX2:4:21412:23622:4437/150‑1 (MQ=60)
                                                                                                              CTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTC‑‑GGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTT                                          >  NB501016:282:H2V7GAFX2:4:21402:17583:18018/1‑151 (MQ=60)
                                                                                                                           CTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTC‑‑GGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTG                             >  NB501016:282:H2V7GAFX2:4:11410:25865:11468/1‑151 (MQ=60)
                                                                                                                                  CGTGATAGGGCCGCAAAGCGGCCCCCTC‑‑GGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGG                      <  NB501016:282:H2V7GAFX2:3:21511:12521:13245/151‑1 (MQ=60)
                                                                                                                                   GTGATAGGGCCGCAAAGCGGCCCCCTC‑‑GGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGC                     >  NB501016:282:H2V7GAFX2:1:21308:12705:6399/1‑151 (MQ=60)
                                                                                                                                   GTGATAGGGCCGCAAAGCGGCCCCCTC‑‑GGTTTTCAGCTGGCAGAAATTACTTCTGCATCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGC                     >  NB501016:282:H2V7GAFX2:1:11311:24852:16644/1‑151 (MQ=60)
                                                                                                                                       TAGGGCCGCAAAGCGGCCCCCTC‑‑GGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTT                 >  NB501016:282:H2V7GAFX2:2:21104:2434:7282/1‑151 (MQ=60)
                                                                                                                                          GGCCGCAAAGCGGCCCCCTC‑‑GGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGA                >  NB501016:282:H2V7GAFX2:3:21602:12310:1296/1‑149 (MQ=60)
                                                                                                                                            CCGCAAAGCGGCCCCCTC‑‑GGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCC            >  NB501016:282:H2V7GAFX2:1:11110:17041:10828/1‑151 (MQ=60)
                                                                                                                                              GCAAAGCGGCCCCCTC‑‑GGTTTTCAGCTGGCAGAAATTACTTCTGCTACACACCTTCGACACTCATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGCGCCTTTGATCCCGAAGTCGTACAGCTTCAGCGTGGTGGTGGCGTTGAAGCCA           >  NB501016:282:H2V7GAFX2:2:11202:1342:5194/1‑150 (MQ=60)
                                                                                                                                                      GCCCCCTC‑‑GGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCG  <  NB501016:282:H2V7GAFX2:1:11211:14197:7069/151‑1 (MQ=60)
                                                                                                                                                      |                                                                                                                                                        
CGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGCGCCCTC‑‑GGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCG  >  NC_002947_CJ‑RC/5677417‑5677716

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: