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
A3 F55 I1 R1 61 67.4 2858014 98.3% 2809427 146.5

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947_CJ‑RC 5,677,568 G→C intergenic (+115/+25) 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,5680GC98.3% 199.2 / ‑6.7 58intergenic (+115/+25)rhlE‑II/yceIATP‑dependent RNA helicase/polyisoprene‑binding acidic stress response factor
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base C (28/29);  total (29/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGCGCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCA  >  NC_002947_CJ‑RC/5677427‑5677707
                                                                                                                                              |                                                                                                                                           
ccAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTGTTGTTCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGGcccctcg                                                                                                                                      >  3:168682/1‑145 (MQ=255)
ccAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCg                                                                                                                                       <  1:247861/149‑1 (MQ=255)
ccAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCg                                                                                                                                       >  8:262673/1‑149 (MQ=255)
ccAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCg                                                                                                                                       <  7:226638/149‑1 (MQ=255)
            aCCCGCAAGGGCTGATACAAGCACCGCGGCGGCGGCTTCGCAGCCCTATCACGACACACGGCCGCTCCTCCAGACAATTTCTTCATCTCAACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTgg                                                                                                                           <  4:140169/149‑1 (MQ=25)
                      gCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTAc                                                                                                                 >  8:362706/1‑149 (MQ=255)
                           tAGAAGTCCCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGTGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTg                                                                                                            <  2:295335/149‑1 (MQ=37)
                             gAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCACGGTTTTCAGCTGGCAGAAATTACTTCTGc                                                                                                           >  6:98765/1‑148 (MQ=255)
                              aaGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCtt                                                                                                         <  7:262673/149‑1 (MQ=255)
                                gaaCCGTGGGGGCTGCTTCGCAGCCCTAGCACGACACACGGCCGCTCCTGCATAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCTTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTc                                                                                                        <  5:186853/146‑1 (MQ=255)
                                       gggggCTGCTTCCCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTc                                                                                                <  4:22688/149‑1 (MQ=255)
                                         gggCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCGCCTCGGGTTTCCGCTGGCAGAAATTACTTCTGGGTCACACCTTCg                                                                                               >  5:165143/1‑148 (MQ=37)
                                          ggCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGa                                                                                              >  2:285123/1‑148 (MQ=255)
                                            ctgctTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGacac                                                                                           <  1:208574/149‑1 (MQ=255)
                                               ctTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACt                                                                                          <  8:306986/147‑1 (MQ=255)
                                                  cGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCGTCGACACTGata                                                                                      >  4:246492/1‑148 (MQ=255)
                                                   gCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATAtc                                                                                    <  4:50609/149‑1 (MQ=255)
                                                     atcccTATCACTACACAAGGCCGCTCCTGAAGAGAATTTCTTGATCTATACACGATCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATAtcc                                                                                   <  5:147495/146‑1 (MQ=38)
                                                      gCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATAtaca                                                                                  >  5:38399/1‑148 (MQ=255)
                                                       cccTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATAtccagg                                                                                >  8:169059/1‑149 (MQ=255)
                                                           aTCACGACACAAGGCCGCTCCTGTAGAGATCAAGAAATTCTCTGCAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATAtccaggtc                                                                              >  2:190740/1‑147 (MQ=14)
                                                            tCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATAtccaggtcca                                                                            >  5:173660/1‑148 (MQ=255)
                                                                     aaGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGGAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCg                                                                   >  3:33042/1‑148 (MQ=255)
                                                                       ggCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAg                                                                <  4:72045/149‑1 (MQ=255)
                                                                         ccGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTa                                                              <  1:96579/149‑1 (MQ=255)
                                                                               cTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCAccc                                                        >  3:142658/1‑149 (MQ=255)
                                                                                tGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCg                                                       >  7:11518/1‑149 (MQ=255)
                                                                                  cAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCggg                                                     >  4:262882/1‑149 (MQ=255)
                                                                                  cAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCgg                                                      >  5:108091/1‑148 (MQ=255)
                                                                                       aaTTTCTTTAACCCTACCGCCCCGGCCTTGGTTCGTGAAAGGGCCGCAAACCGGCCCCCGCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCttt                                                <  4:14368/149‑1 (MQ=255)
                                                                                           tCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATg                                            >  8:344991/1‑149 (MQ=255)
                                                                                            cTTGATCTCTACAGTAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGAt                                             <  5:206791/147‑1 (MQ=255)
                                                                                                    ctACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCg                                   >  3:95612/1‑149 (MQ=255)
                                                                                                      aCAGGCGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAg                                      <  5:292118/144‑1 (MQ=255)
                                                                                                       cAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCgttc                                <  8:153108/149‑1 (MQ=255)
                                                                                                        aGGAGCGGCCTTGTGTCGTGATAGGACCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAAATACTTCTGCTTCACACCTTCGACACTGATATACAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCgtt                                 >  3:29852/1‑147 (MQ=39)
                                                                                                         ggAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCgttcag                              <  8:233514/149‑1 (MQ=255)
                                                                                                            gCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCgttcaggt                            <  5:206087/148‑1 (MQ=255)
                                                                                                             cGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCgttcaggtt                           >  2:263803/1‑148 (MQ=255)
                                                                                                              ggCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCgttcaggttca                         <  1:285123/149‑1 (MQ=255)
                                                                                                              ggCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCgttcaggttc                          <  6:38399/148‑1 (MQ=255)
                                                                                                                 cTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCg                      >  8:342379/1‑149 (MQ=255)
                                                                                                                     tgtCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCgtgg                   <  1:263803/148‑1 (MQ=255)
                                                                                                                       tCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCgtggtgg                >  4:229331/1‑149 (MQ=255)
                                                                                                                        cGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCgtggtg                 <  5:244175/147‑1 (MQ=255)
                                                                                                                         gttaTAGGGCCGCAAAGCGGCCCCCTATTTTTTCATCTTTCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCCGGTTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCgtggtggt               <  2:70756/145‑1 (MQ=255)
                                                                                                                         gTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCgtggtggtg              >  4:201982/1‑149 (MQ=255)
                                                                                                                          tGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTTCTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGtggtggtgg             >  2:174230/1‑149 (MQ=255)
                                                                                                                          tGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCgtggtggtg              <  5:161292/148‑1 (MQ=255)
                                                                                                                          tGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGtggtggtgg             >  7:342190/1‑149 (MQ=255)
                                                                                                                            aTAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGc            >  6:329628/1‑148 (MQ=255)
                                                                                                                              aGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACCCTGATATCCAGGTCCAGGGTCTGCGAAGTAGCCCCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGtt         >  7:223339/1‑149 (MQ=255)
                                                                                                                               gggCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGtt         <  7:151310/148‑1 (MQ=255)
                                                                                                                                ggCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGa       <  7:228222/149‑1 (MQ=255)
                                                                                                                                 gCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGaa      <  3:12033/149‑1 (MQ=255)
                                                                                                                                  ccGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAg     >  3:225211/1‑149 (MQ=255)
                                                                                                                                  ccGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTT‑‑CACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAg     >  6:194672/1‑147 (MQ=255)
                                                                                                                                     cAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCa  <  4:138378/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                           
CCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGCGCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCA  >  NC_002947_CJ‑RC/5677427‑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
CCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCG‑‑CGCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGC  >  NC_002947_CJ‑RC/5677426‑5677705
                                                                                                                                                 |                                                                                                                                          
CCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTGTTGTTCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCG‑‑GGCCCCTCGG                                                                                                                                    >  NB501016:282:H2V7GAFX2:2:11310:18731:3132/1‑151 (MQ=60)
CCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑‑CCCCCTCGG                                                                                                                                    <  NB501016:282:H2V7GAFX2:4:21409:1694:2208/151‑1 (MQ=60)
CCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑‑CCCCCTCGG                                                                                                                                    <  NB501016:282:H2V7GAFX2:1:21201:20684:17703/151‑1 (MQ=60)
             ACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTAACACGACACAAGGCCGCTCCTGCAGAGAAATTATTGATCTCTAGAGGAGCGGCCTTGTGTAGTGATATGGACGCAAAGCGG‑‑CCCCCTCGGATTTCAGCTGGC                                                                                                                        >  NB501016:282:H2V7GAFX2:1:11102:3511:11053/1‑150 (MQ=60)
                              GAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGGTCTCTACAGGAGCGGCCTTGTGTCGTGCTAGGGCCGCAAAGCGG‑‑CCCCCTCGGGGTTCAGGTGGCAGAAATTACTTCTGCTTT                                                                                                      >  NB501016:282:H2V7GAFX2:4:21512:6349:5366/1‑151 (MQ=60)
                              GAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTC                                                                                                      <  NB501016:282:H2V7GAFX2:4:21504:18300:1237/151‑1 (MQ=60)
                               acgaaccgtGGGGGCTGCTTCGCAGCCCTAGCACGACACACGGCCGCTCCTGCATAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCTTGATAGGGCCGCAAAGCGG‑‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCA                                                                                                     <  NB501016:282:H2V7GAFX2:3:21401:11753:16911/142‑1 (MQ=60)
                                         GGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCGCCTCGGGTTTCCGCTGGCAGAAATTACTTCTGGGTCACACCTTCGAC                                                                                           >  NB501016:282:H2V7GAFX2:3:11609:10030:11670/1‑151 (MQ=60)
                                            GCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACT                                                                                        <  NB501016:282:H2V7GAFX2:1:21105:4383:15097/151‑1 (MQ=60)
                                                    GCATCCCTATCACTACACAAGGCCGCTCCTGAAGAGAATTTCTTGATCTATACACGATCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCA                                                                                <  NB501016:282:H2V7GAFX2:3:11606:6628:5332/151‑1 (MQ=60)
                                                      AGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGAAACTGATATACAGG                                                                              >  NB501016:282:H2V7GAFX2:3:11409:4339:2397/1‑151 (MQ=60)
                                                            ATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGG                                                                        >  NB501016:282:H2V7GAFX2:3:11611:9478:3953/1‑151 (MQ=60)
                                                                     CAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑‑CCCCCTCGGTTTTCAGCTGGGAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGA                                                                >  NB501016:282:H2V7GAFX2:2:11107:21486:15293/1‑150 (MQ=60)
                                                                     CAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑‑CCCCCGCGGTTTTCAGCTGGCAGAATTTACATCTGCTTCACACCTTCGACACTGATATCCAGTTCCAGGGTCTGCGA                                                                >  NB501016:282:H2V7GAFX2:2:21307:19828:12398/1‑150 (MQ=60)
                                                                         GCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAG                                                           <  NB501016:282:H2V7GAFX2:1:11207:20580:5410/151‑1 (MQ=60)
                                                                            gctcctgtagagatcaagaaatTCTCTGCAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCAC                                                        >  NB501016:282:H2V7GAFX2:4:21410:10201:8028/23‑151 (MQ=60)
                                                                               CCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCG                                                     >  NB501016:282:H2V7GAFX2:2:11305:9117:5824/1‑151 (MQ=60)
                                                                                CTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGG                                                    >  NB501016:282:H2V7GAFX2:4:11403:9143:6196/1‑151 (MQ=60)
                                                                                  GCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGC                                                  >  NB501016:282:H2V7GAFX2:3:11510:11431:18746/1‑151 (MQ=60)
                                                                                           TTCTTGATCTCTACAGTAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATG                                          <  NB501016:282:H2V7GAFX2:3:21406:26531:7743/150‑1 (MQ=60)
                                                                                                    TCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGT                                >  NB501016:282:H2V7GAFX2:2:11208:2327:4373/1‑151 (MQ=60)
                                                                                                     CTACAGGCGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGT                                   <  NB501016:282:H2V7GAFX2:3:21512:24913:6989/147‑1 (MQ=60)
                                                                                                        CAGGAGCGGCCTTGTGTCGTGATAGGACCGCAAAGCGG‑‑CCCCCTCGGTTTTCAGCTGGCAGAAAATACTTCTGCTTCACACCTTCGACACTGATATACAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTC                              >  NB501016:282:H2V7GAFX2:2:11107:20078:2809/1‑149 (MQ=60)
                                                                                                           GAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTT                         <  NB501016:282:H2V7GAFX2:3:21406:8570:4637/151‑1 (MQ=60)
                                                                                                              CGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAG                      <  NB501016:282:H2V7GAFX2:1:21209:10950:10154/151‑1 (MQ=60)
                                                                                                                     GTGTCGTGATAGGGCCGCAAAGCGG‑‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGT                <  NB501016:282:H2V7GAFX2:1:21205:5788:2121/150‑1 (MQ=60)
                                                                                                                       GTCGTGATAGGGCCGCAAAGCGG‑‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGG              <  NB501016:282:H2V7GAFX2:3:21502:5748:9286/150‑1 (MQ=60)
                                                                                                                         CGTGATAGGGCCGCAAAGCGG‑‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGG           <  NB501016:282:H2V7GAFX2:3:11608:6131:17020/151‑1 (MQ=60)
                                                                                                                          GTGATAGGGCCGCAAAGCGG‑‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGC          >  NB501016:282:H2V7GAFX2:4:21607:22801:16352/1‑151 (MQ=60)
                                                                                                                              TAGGGCCGCAAAGCGG‑‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACCCTGATATCCAGGTCCAGGGTCTGCGAAGTAGCCCCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTG      >  NB501016:282:H2V7GAFX2:4:21408:8621:9176/1‑151 (MQ=60)
                                                                                                                               AGGGCCGCAAAGCGG‑‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTG      <  NB501016:282:H2V7GAFX2:4:11606:20209:5525/150‑1 (MQ=60)
                                                                                                                                GGGCCGCAAAGCGG‑‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAA    <  NB501016:282:H2V7GAFX2:4:21409:23361:8083/151‑1 (MQ=60)
                                                                                                                                 GGCCGCAAAGCGG‑‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAG   <  NB501016:282:H2V7GAFX2:2:11103:9469:10711/151‑1 (MQ=60)
                                                                                                                                  GCCGCAAAGCGG‑‑CCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGC  >  NB501016:282:H2V7GAFX2:2:21109:13108:18011/1‑151 (MQ=60)
                                                                                                                                                 |                                                                                                                                          
CCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCG‑‑CGCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGC  >  NC_002947_CJ‑RC/5677426‑5677705

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: