Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A1 F41 I0 R1 212 10.4 395026 92.3% 364608 148.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 8,903 T→G 82.5% V222V (GTT→GTG talB → transaldolase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092738,9030TG82.5% 146.7 / 16.6 56V222V (GTT→GTGtalBtransaldolase B
Reads supporting (aligned to +/- strand):  ref base T (7/3);  new base G (26/20);  total (33/23)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTGACCATCGCACCGGCACTGCTGAAAGAGCTGGCGGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAG  >  CP009273/8760‑9042
                                                                                                                                               |                                                                                                                                           
atttcGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATgg                                                                                                                                         <  2:38725/145‑1 (MQ=255)
  ttcGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGgcg                                                                                                                                      >  2:173326/2‑149 (MQ=255)
  ttcGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGgcg                                                                                                                                      >  1:156859/2‑149 (MQ=255)
    cGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAACGAGCACGGTTATGAAAACGTGGTTATGGgcgc                                                                                                                                     >  1:13499/1‑148 (MQ=255)
       cgttcgttGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCt                                                                                                                                <  2:145644/150‑1 (MQ=255)
       cgttcgttGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCt                                                                                                                                >  2:54252/1‑150 (MQ=255)
       cgttcgttGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGAGCAAGCt                                                                                                                                >  1:56432/1‑150 (MQ=255)
        gTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCtt                                                                                                                               <  1:192567/150‑1 (MQ=255)
         ttcgttgGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTc                                                                                                                              <  1:122670/147‑1 (MQ=255)
           cgttgGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTcc                                                                                                                             <  2:47141/148‑1 (MQ=255)
               ggCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAAc                                                                                                                        >  1:62989/1‑150 (MQ=255)
                gCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAAcg                                                                                                                       >  1:7906/1‑149 (MQ=255)
                 ccGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAAcgt                                                                                                                      >  1:157449/1‑148 (MQ=255)
                 ccGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAAcgt                                                                                                                      >  2:165241/1‑148 (MQ=255)
                 ccGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAAcgt                                                                                                                      >  1:48293/1‑148 (MQ=255)
                  cGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAAcgtt                                                                                                                     >  1:181283/1‑147 (MQ=255)
                     atcctcgACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGc                                                                                                                  >  2:139394/7‑150 (MQ=255)
                      tcctcgACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCg                                                                                                                 <  1:80813/145‑1 (MQ=255)
                                gTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGa                                                                                                       <  1:115272/150‑1 (MQ=255)
                                 tACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGaa                                                                                                      >  1:183717/1‑150 (MQ=255)
                                   cAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACt                                                                                                    >  1:3713/1‑150 (MQ=255)
                                     aaGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTgg                                                                                                  <  1:54252/150‑1 (MQ=255)
                                         caacaCCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAActggctgg                                                                                              <  2:52724/146‑1 (MQ=255)
                                            tACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTg                                                                                           >  2:132840/1‑150 (MQ=255)
                                              ccGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCGGGAACTGGCAGGCTGc                                                                                          <  2:13499/149‑1 (MQ=255)
                                               cGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGa                                                                                        <  2:65357/150‑1 (MQ=255)
                                               cGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTAGCGAAATTCTGGAACTGGCTGGCTGTGa                                                                                        >  2:27188/1‑150 (MQ=255)
                                                      aaagaGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTg                                                                                 <  2:26209/150‑1 (MQ=255)
                                                          agTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCa                                                                             <  1:165241/150‑1 (MQ=255)
                                                           gTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCAt                                                                            >  1:164040/1‑150 (MQ=255)
                                                            tACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATc                                                                           <  2:19165/150‑1 (MQ=255)
                                                            tACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATc                                                                           >  1:122567/1‑150 (MQ=255)
                                                            tACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATc                                                                           <  2:48293/150‑1 (MQ=255)
                                                            tACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATc                                                                           >  1:49438/1‑150 (MQ=255)
                                                            tACGCACCGACAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTGCAGACAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTGCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATc                                                                           <  2:149479/150‑1 (MQ=255)
                                                                    ggCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGc                                                                   >  1:160555/1‑150 (MQ=255)
                                                                     gCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGcg                                                                  <  2:88846/150‑2 (MQ=255)
                                                                      cAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGcgc                                                                 >  2:127249/1‑148 (MQ=255)
                                                                       agaagaTCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGct                                                                >  2:107126/1‑150 (MQ=255)
                                                                                gggCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTGACCATCGCACCGGCACTGCTGaaa                                                         >  2:160757/1‑148 (MQ=255)
                                                                                     tGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTGACCATCGCACCGGCACTGCTGAAAGAGCTgg                                                  >  2:145340/1‑150 (MQ=255)
                                                                                             gtcagcgAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGTCCGAAAGc                                          >  1:112853/7‑150 (MQ=255)
                                                                                                   gAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGTCCGAAAGCGAAggc                                    <  2:148319/150‑2 (MQ=255)
                                                                                                   gAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGTCCGAAAGCGAAggc                                    <  1:127249/150‑2 (MQ=255)
                                                                                                     taTCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGTCCGAAAGCGAAGgcgc                                  >  2:14118/2‑150 (MQ=255)
                                                                                                     aaTCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGTCCGAAAGCGAAGgcgc                                  >  1:82091/1‑150 (MQ=255)
                                                                                                     aaTCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGTCCGAAAGCGAAGgcgc                                  >  1:163373/1‑150 (MQ=255)
                                                                                                           cgaGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGTCCGAAAGCGAAGgcgccg                                <  1:173326/144‑3 (MQ=255)
                                                                                                           ccAGTACTACAAAGAGCACGGTTATGAAATCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTGACCATCGCACCGGCACTGCTGAAAGAGCCGGCGGAGAGCGAAGGGGCTATCg                             >  2:125336/1‑149 (MQ=255)
                                                                                                               tactacAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTGACCATCGCACCGGCACTGCTGAAAGAGCTGGCGGAGAGCGAAGGGGCTATCGAACg                         >  1:122902/1‑149 (MQ=255)
                                                                                                                  tacAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGTCCGAAAGCGAAGGCGCCGTTGAGCGTaag                     >  2:55266/1‑149 (MQ=255)
                                                                                                                    caagcagCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGTCCGAAAGCGAAGGCGCCGTTGAGCGTAAGCt                   >  1:76814/6‑150 (MQ=255)
                                                                                                                    cAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTGACCATCGCACCGGCACTGCTGAAAGAGCTGGCGGAGAGCGAAGGGGCTATCGAACGTAAACt                   <  1:132840/150‑1 (MQ=255)
                                                                                                                         agCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGTCCGAAAGCGAAGGCGCCGTTGAGCGTAAGCTggtat              <  1:169550/150‑5 (MQ=255)
                                                                                                                           cACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTGACCATCGCACCGGCACTGCTGAAAGAGCTGGCGGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTac            <  2:35894/150‑1 (MQ=255)
                                                                                                                                     aaaCCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTGACCATCGCACCGGCACTGCTGAAAGAGCTGGCGGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAg  >  2:72400/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                           
ATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTGACCATCGCACCGGCACTGCTGAAAGAGCTGGCGGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAG  >  CP009273/8760‑9042

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTGACCATCGCACCGGCACTGCTGAAAGAGCTGGCGGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGT  >  CP009273/8753‑9053
                                                                                                                                                      |                                                                                                                                                      
GTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGAACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTT                                                                                                                                                        >  SRR12117902.36543/1‑151 (MQ=60)
         TTCGCCGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGC                                                                                                                                                >  SRR12117902.159718/1‑150 (MQ=44)
           CGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAACGAGCACGGCTATGAAAACGTGGTTATGGGCGCA                                                                                                                                               >  SRR12117902.13718/1‑149 (MQ=60)
              CGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTT                                                                                                                                          <  SRR12117902.196060/151‑1 (MQ=47)
              CGTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGAGCAAGCTT                                                                                                                                          >  SRR12117902.57462/1‑151 (MQ=39)
               GTTCGTTGGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTC                                                                                                                                         <  SRR12117902.124908/151‑1 (MQ=48)
                      GGCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACG                                                                                                                                  >  SRR12117902.64128/1‑151 (MQ=53)
                       GCCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGT                                                                                                                                 >  SRR12117902.8054/1‑151 (MQ=51)
                        CCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACgtt                                                                                                                                >  SRR12117902.160320/1‑148 (MQ=50)
                        CCGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACgtt                                                                                                                                >  SRR12117902.49185/1‑148 (MQ=50)
                         CGTATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACgttg                                                                                                                               >  SRR12117902.184564/1‑147 (MQ=51)
                            ATCCTCGACTGGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCG                                                                                                                            <  SRR12117902.82279/151‑1 (MQ=50)
                                      GGTACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGA                                                                                                                  <  SRR12117902.117369/151‑1 (MQ=54)
                                        TACAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAAC                                                                                                                >  SRR12117902.187068/1‑151 (MQ=54)
                                          CAAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTG                                                                                                              >  SRR12117902.3764/1‑151 (MQ=54)
                                           AAAGCCAACACCGACAAGAAAGAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGG                                                                                                             <  SRR12117902.55250/151‑1 (MQ=56)
                                                                GAGTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCA                                                                                        <  SRR12117902.168263/151‑1 (MQ=54)
                                                                  GTACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATC                                                                                      >  SRR12117902.167038/1‑151 (MQ=54)
                                                                   TACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCG                                                                                     >  SRR12117902.50347/1‑151 (MQ=54)
                                                                   TACGCACCGGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCG                                                                                     >  SRR12117902.124804/1‑151 (MQ=54)
                                                                           GGCAGAAGATCCGGGCGTGGTTTCCGTCAGCGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCG                                                                             >  SRR12117902.163486/1‑151 (MQ=53)
                                                                                      cgggcgtggtttccgtcagcGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAA                                                                     >  SRR12117902.180746/21‑148 (MQ=56)
                                                                                              gtttccgtcagcGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACG                                                                                                                                  <  SRR12117902.97562/67‑1 (MQ=23)
                                                                                                   cgtcagcGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGtccgaaagc                                                     >  SRR12117902.179345/8‑142 (MQ=60)
                                                                                                    gtcagcGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGtccgaaagcg                                                    >  SRR12117902.114907/7‑141 (MQ=60)
                                                                                                         CGAAATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGtccgaaagcgaaggc                                               <  SRR12117902.129569/151‑16 (MQ=60)
                                                                                                            AATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGtccgaaagcgaaggcgcc                                            >  SRR12117902.166358/1‑133 (MQ=60)
                                                                                                            AATCTACGAGTATTACAAGCAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGtccgaaagcgaaggcgcc                                            >  SRR12117902.83582/1‑133 (MQ=60)
                                                                                                                 acgagtattacaagcAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGtccgaaagcgaaggcgccgt                                          <  SRR12117902.14825/133‑21 (MQ=60)
                                                                                                                 acgagtattacaagcAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGtccgaaagcgaaggcgccg                                           <  SRR12117902.176465/132‑20 (MQ=60)
                                                                                                                      TACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTGACCATCGCACCGGCACTGCTGAAAGAGCTGGCGGAGAGCGAAGGGGCTATCGAACGT                                   >  SRR12117902.125143/1‑150 (MQ=60)
                                                                                                                          ACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTGACCATCGCACCGGCACTGCTGAAAGAGCTGGCGGAGAGCGAAGGGGCTATCGAACGTAAACT                              <  SRR12117902.135267/151‑1 (MQ=60)
                                                                                                                           caagcAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGtccgaaagcgaaggcgccgtggagcgtaagctg                             >  SRR12117902.78205/6‑118 (MQ=60)
                                                                                                                               CAGCACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGtccgaaagcgaaggcgccgttgagcgtaagctggtat                         <  SRR12117902.172648/151‑38 (MQ=60)
                                                                                                                                   ACGGCTATGAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGtccgaaagcgaaggcgccgttgagcgtaagctggtatatac                     <  SRR12117902.109079/151‑42 (MQ=60)
                                                                                                                                           GAAACCGTCGTGATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGtccgaaagcgaaggcgccgttgagcgtaagctggtatataccggtga               >  SRR12117902.149108/1‑102 (MQ=60)
                                                                                                                                                      GATGGGCGCAAGCTTCCGTAACGTTGGCGAAATTCTGGAACTGGCTGGCTGTGACCGCCTGACCATCGCCCCGGCGCTGCTGAAAGAGCTGtccgaaagcgaaggcgccgttgagcgtaagctggtatataccggtgaagtgaaagcgcgt  <  SRR12117902.165009/151‑61 (MQ=60)
                                                                                                                                                      |                                                                                                                                                      
GTTCCTGATCTCGCCGTTTGTTGGCCGTATTCTTGACTGGTACAAAGCGAATACCGATAAGAAAGAGTACGCTCCGGCAGAAGATCCGGGCGTGGTTTCTGTATCTGAAATCTACCAGTACTACAAAGAGCACGGTTATGAAACCGTGGTTATGGGCGCAAGCTTCCGTAACATCGGCGAAATTCTGGAACTGGCAGGCTGCGACCGTCTGACCATCGCACCGGCACTGCTGAAAGAGCTGGCGGAGAGCGAAGGGGCTATCGAACGTAAACTGTCTTACACCGGCGAAGTGAAAGCGCGT  >  CP009273/8753‑9053

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: