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
A2 F24 I1 R1 60 68.7 2919236 98.6% 2878366 146.9

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.2% 185.7 / ‑3.6 54intergenic (+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 (23/30);  total (24/30)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGCGCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCA  >  NC_002947_CJ‑RC/5677418‑5677707
                                                                                                                                                    |                                                                                                                                              
ggTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAgcg                                                                                                                                               <  5:180705/149‑1 (MQ=255)
    aGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTCCCGCGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAgcggcc                                                                                                                                            <  7:340500/148‑4 (MQ=255)
      cgcCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAgcggccccc                                                                                                                                         >  1:226799/1‑143 (MQ=255)
      cgcCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCCATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGGTAGGGCCGCAAAgcggccccc                                                                                                                                         >  8:161672/1‑143 (MQ=255)
         ccAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGTCCGCTCCGGCAGAGAATATCATGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCCGCCCCCgcc                                                                                                                                     >  7:204920/1‑147 (MQ=37)
         ccAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCg                                                                                                                                       >  7:62377/1‑149 (MQ=255)
                  ctgACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGc                                                                                                                              >  8:106429/1‑149 (MQ=255)
                  ctgACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCGCGGGTTTAAGc                                                                                                                              >  3:116628/1‑149 (MQ=37)
                    gACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCTTCGGTTTTCAGCTg                                                                                                                            >  6:361436/1‑149 (MQ=255)
                       ccGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTACTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCACGGTTTTGAGCttg                                                                                                                           >  3:287630/1‑145 (MQ=37)
                         gCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAg                                                                                                                        >  4:299640/1‑148 (MQ=255)
                         gCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAg                                                                                                                        >  4:116178/1‑148 (MQ=255)
                               gCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTAc                                                                                                                 <  6:273764/149‑1 (MQ=255)
                                   aTAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTc                                                                                                              <  2:238666/148‑1 (MQ=255)
                                   aTAGAAGTACCGAGGGGGCGGCTTCGCAGCCCTATCACGACACACGGCCGCTCCTGCAGACAATTTCTTTATCTCAACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACtt                                                                                                               <  6:236230/147‑1 (MQ=37)
                                                     ctgctTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTTATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGacac                                                                                           <  5:153174/149‑1 (MQ=255)
                                                     ctgctTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGacac                                                                                           <  4:152506/149‑1 (MQ=255)
                                                     cggctCCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGacac                                                                                           <  4:31204/147‑1 (MQ=255)
                                                      tgctTCGCAGCCCTATCACGACACAAGGCCGCTACTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTATGCTTCACACCTTCGACACt                                                                                          >  7:20599/1‑149 (MQ=38)
                                                         ttCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGacac                                                                                           >  4:98720/1‑145 (MQ=255)
                                                         ttCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGacac                                                                                           <  3:98720/145‑1 (MQ=255)
                                                                  cTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATAtccaggtc                                                                              >  3:81947/1‑149 (MQ=255)
                                                                     tCACGAGACAAGGCCGCTCCTGCAGAGAATTTCTTTATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAATCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATAtccaggtcca                                                                            <  5:183652/148‑1 (MQ=38)
                                                                        cGACACAAGGCCGCTCCTGCAGAGAATTTCTTGCGCTCTACAGGCGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGGTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGGCCAGGGt                                                                        <  8:236187/149‑1 (MQ=38)
                                                                        cGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGt                                                                        >  2:151337/1‑149 (MQ=255)
                                                                         gACACACGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATAtccaggtccagg                                                                          <  2:340900/146‑1 (MQ=255)
                                                                         gACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGCCACTGATATCCAGGTCCAGGGTc                                                                       >  4:271545/1‑149 (MQ=255)
                                                                           cacaAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTg                                                                     <  8:282023/149‑1 (MQ=255)
                                                                            acaAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGc                                                                    <  1:151337/149‑1 (MQ=255)
                                                                                ggCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAg                                                                <  1:50114/149‑1 (MQ=255)
                                                                                  ccGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTa                                                              <  2:226799/149‑1 (MQ=255)
                                                                                        cTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCAccc                                                        >  7:17190/1‑149 (MQ=255)
                                                                                          gCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCgg                                                      >  6:319117/1‑149 (MQ=255)
                                                                                           cAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCggg                                                     >  7:9408/1‑149 (MQ=255)
                                                                                           cAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCggg                                                     >  6:144656/1‑149 (MQ=255)
                                                                                               gaATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCt                                                  <  4:290834/148‑1 (MQ=255)
                                                                                                         aTCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCg                                         >  2:194946/1‑147 (MQ=255)
                                                                                                          tctcTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAg                                      >  5:31782/1‑149 (MQ=255)
                                                                                                              tACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGaa                                       <  8:66298/144‑1 (MQ=255)
                                                                                                               aCAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATACAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCgtt                                 >  1:205668/1‑149 (MQ=255)
                                                                                                                cAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCgttc                                <  3:351329/149‑1 (MQ=255)
                                                                                                                cAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCgttc                                <  2:287095/149‑1 (MQ=255)
                                                                                                                   gAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCgttcagg                             <  7:106429/149‑1 (MQ=255)
                                                                                                                      cGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCgttcaggtt                           >  8:20500/1‑148 (MQ=255)
                                                                                                                       ggCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCgttcaggttca                         <  4:116628/149‑1 (MQ=255)
                                                                                                                        gCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCACGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCgttcaggttcag                        >  2:91399/1‑149 (MQ=255)
                                                                                                                          cTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCTTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCg                      <  5:361436/149‑1 (MQ=255)
                                                                                                                          cTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCGCGGTTTTAAGCTGGCAGAAATTAATTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCCCCCGGGCCTTTGATGCCGAAGGCgttcaggttaatcg                      >  1:325794/1‑144 (MQ=255)
                                                                                                                                tCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCgtggtg                 <  2:96994/148‑1 (MQ=255)
                                                                                                                                tCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCgtggtg                 <  1:358248/148‑1 (MQ=255)
                                                                                                                                   tGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGtggtggtgg             <  8:62377/149‑1 (MQ=255)
                                                                                                                                     aTAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGc            <  3:116178/148‑1 (MQ=255)
                                                                                                                                      tAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGGCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGt          <  3:178778/149‑1 (MQ=255)
                                                                                                                                        gggCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGTGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGtt         <  2:203422/148‑1 (MQ=255)
                                                                                                                                        gggCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCgttcag                              >  6:44828/1‑127 (MQ=255)
                                                                                                                                        gggCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTg        <  3:299640/149‑1 (MQ=255)
                                                                                                                                            cGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGc    <  8:204920/149‑1 (MQ=255)
                                                                                                                                              cAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCa  <  2:133010/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
GGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGCGCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCA  >  NC_002947_CJ‑RC/5677418‑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
CGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGCGCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGA  >  NC_002947_CJ‑RC/5677417‑5677702
                                                                                                                                                      |                                                                                                                                          
CGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG                                                                                                                                           <  NB501016:282:H2V7GAFX2:3:11611:14534:12799/151‑1 (MQ=60)
    GAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTCCCGCGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCC                                                                                                                                       <  NB501016:282:H2V7GAFX2:4:21605:8682:20183/150‑1 (MQ=60)
      GCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGgccccct                                                                                                                                     >  NB501016:282:H2V7GAFX2:1:21108:1978:10999/1‑144 (MQ=60)
         CCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCC‑‑TCGG                                                                                                                                 >  NB501016:282:H2V7GAFX2:4:11412:5610:10901/1‑151 (MQ=60)
                  GCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGgcccccgcgggtttaagct                                                                                                                         >  NB501016:282:H2V7GAFX2:2:11211:9313:18372/1‑132 (MQ=60)
                     cccccccaagggctgaaacacgaACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTCCAGATAATTTCTTGATCTCTACAGCAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCC‑‑TCGGTTTTCAGCTG                                                                                                                       <  NB501016:282:H2V7GAFX2:3:11407:9892:7416/126‑1 (MQ=60)
                       CCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTACTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGgcccccacggttttgagcttgc                                                                                                                      >  NB501016:282:H2V7GAFX2:2:21211:14164:13529/1‑127 (MQ=60)
                                    ATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCAGCAGAGAATTTCTTGATGTCTACAGGCGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCC‑‑TCGGTTTTCAGCTGGGAGAAATTTCTTCTGG                                                                                                      >  NB501016:282:H2V7GAFX2:3:11511:24976:8712/1‑151 (MQ=60)
                                        AAGTACCGCGGGGGCTGCTACCCACCCCTACCACGACACACGGCCGCACCTGCACATAATTTCTTGATCTCACCAGGCGCGGCCTTGTGTCGTTATAGGGCCGCAAAGCGGCCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCA                                                                                                  <  NB501016:282:H2V7GAFX2:2:21109:25164:6318/151‑1 (MQ=60)
                                                     GCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTTATAGGGCCGCAAAGCGGCCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACT                                                                                     <  NB501016:282:H2V7GAFX2:3:11606:8469:11497/151‑1 (MQ=60)
                                                      CTGCTTCGCAGCCCTATCACGACACAAGGCCGCTACTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTATGCTTCACACCTTCGACACTG                                                                                    >  NB501016:282:H2V7GAFX2:4:11404:21618:17170/1‑151 (MQ=60)
                                                         CTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGATAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACT                                                                                     <  NB501016:282:H2V7GAFX2:2:11208:10492:9474/147‑1 (MQ=60)
                                                                  CCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCC                                                                        >  NB501016:282:H2V7GAFX2:2:11205:22536:5171/1‑151 (MQ=60)
                                                                     ATCACGAGACAAGGCCGCTCCTGCAGAGAATTTCTTTATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAATCGGCCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAG                                                                      <  NB501016:282:H2V7GAFX2:3:11612:17335:4033/150‑1 (MQ=60)
                                                                       catgaaaaatggcgccaagtcctgagattttcttgaggtctaatggaggggcgtgttgtgGTGATAGGGCCGCAATGCGGCCCCC‑‑GCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGG                                                                    <  NB501016:282:H2V7GAFX2:1:21310:18728:16522/90‑1 (MQ=60)
                                                                            CACACGGCCGCTCCTGCAGATAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCG                                                              <  NB501016:282:H2V7GAFX2:1:11305:18930:10447/151‑1 (MQ=60)
                                                                                AGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGT                                                          <  NB501016:282:H2V7GAFX2:1:11110:22916:3324/151‑1 (MQ=60)
                                                                                        CCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCG                                                  >  NB501016:282:H2V7GAFX2:4:11404:17593:4620/1‑151 (MQ=60)
                                                                                           GCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGC                                               >  NB501016:282:H2V7GAFX2:4:11402:17266:15634/1‑151 (MQ=60)
                                                                                                          ATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGT                                >  NB501016:282:H2V7GAFX2:3:11407:11691:9937/1‑151 (MQ=60)
                                                                                                               TACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATACAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTC                           >  NB501016:282:H2V7GAFX2:1:21104:6736:8310/1‑151 (MQ=60)
                                                                                                                ACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCA                          <  NB501016:282:H2V7GAFX2:2:21312:21821:11884/151‑1 (MQ=60)
                                                                                                                   GGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGT                       <  NB501016:282:H2V7GAFX2:4:11508:22102:17015/151‑1 (MQ=60)
                                                                                                                          CCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCC‑‑GCGGTTTTAAGCTGGCAGAAATTAATTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCCCCCGGGCCTTTGATGCCGAAGGCGTTCAGGTTAATCGT                >  NB501016:282:H2V7GAFX2:1:21305:13791:1583/1‑151 (MQ=60)
                                                                                                                          CCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCTTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGT                <  NB501016:282:H2V7GAFX2:3:21611:12462:19733/151‑1 (MQ=60)
                                                                                                                             TGTGTCGTGTTAGGGCCGCAAAGCGGCCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGG              <  NB501016:282:H2V7GAFX2:4:11607:22771:13597/150‑1 (MQ=60)
                                                                                                                                GTCGTGATAGGGCCGCAAAGCGGCCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGT           <  NB501016:282:H2V7GAFX2:1:21311:14031:17697/150‑1 (MQ=60)
                                                                                                                                     GATAGGGCCGCAAAGCGGCCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCG      <  NB501016:282:H2V7GAFX2:2:11211:21616:16724/150‑1 (MQ=60)
                                                                                                                                      ATAGGGCCGCAAAGCGGCCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGGCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTT    <  NB501016:282:H2V7GAFX2:2:11311:2641:14328/151‑1 (MQ=60)
                                                                                                                                        AGGGCCGCAAAGCGGCCCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGA  <  NB501016:282:H2V7GAFX2:2:21302:4748:11801/151‑1 (MQ=60)
                                                                                                                                                      |                                                                                                                                          
CGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGCGCCC‑‑TCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGA  >  NC_002947_CJ‑RC/5677417‑5677702

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: