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
A6 F48 I1 R1 61 74.3 2613450 99.1% 2589928 139.8

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 447,009 A→T Y233N (TAT→AAT)  cyoE ← heme O synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913447,0090AT98.3% 203.6 / ‑3.7 59Y233N (TAT→AAT) cyoEheme O synthase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base T (33/25);  total (34/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

CACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGTTTGCCGCCTGGT  >  NC_000913/446875‑447143
                                                                                                                                         |                                                                                                                                      
caCCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATTTTtat                                                                                                                                   >  2:888228/1‑140 (MQ=255)
caCCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATTTTtat                                                                                                                                   >  1:834150/1‑140 (MQ=255)
 aCCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATTTTtata                                                                                                                                  <  1:391491/140‑1 (MQ=255)
     agagGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATTTTTATACCCa                                                                                                                              >  1:675492/1‑140 (MQ=255)
         ggCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATTTTTATACCCAGCGt                                                                                                                          >  2:736964/1‑140 (MQ=255)
          gCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATTTTTATACCCAGCGTa                                                                                                                         >  2:502381/1‑140 (MQ=255)
            aGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATTTTTATACCCAGCGTAAc                                                                                                                       >  2:266206/1‑140 (MQ=255)
                aTGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATTTTTATACCCAGCGTAACCGcc                                                                                                                   >  2:1190458/1‑140 (MQ=255)
                    ggataaaagagAAGCCAAACAGTTGGCGCCCCAAGTTTCTGTAATAAGAAACTTAAAACCCGCGCGGACCAATCCCAACCCCCAAGACCCTAACCGCCGCGGCGACCACCAGATTTTTATACCCAGCGTAACCGCCAaga                                                                                                               <  2:566873/133‑1 (MQ=255)
                    cGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATTTTTATACCCAGCGTAACCGCCAaga                                                                                                               >  2:573445/1‑140 (MQ=255)
                      agaatAGAGAACCAAAACAGCTTGGGCCCCAAAATTTGGTCATCGACAACTTTAAAACCGCGCAGACCCATCCCTACCCCCCAGACCATACCCGCCGCGGCGACCACCAGATTTTTATACCCAGCGTAACCGCCAagaga                                                                                                             <  2:772353/137‑1 (MQ=255)
                        gatAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATTTTTATACCCAGCGTAACCGCCAagagaga                                                                                                           >  2:133835/1‑140 (MQ=255)
                             agaAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATTTTTATACCCAGCGTAACCGCCAAGAGAGAGCATc                                                                                                      >  1:13499/1‑140 (MQ=255)
                             agaAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATTTTTATACCCAGCGTAACCGCCAAGAGAGAGCATc                                                                                                      >  2:820573/1‑140 (MQ=255)
                                aGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATTTTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGc                                                                                                   >  1:223945/1‑140 (MQ=255)
                                 gCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATTTTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCg                                                                                                  >  2:861942/1‑140 (MQ=255)
                                  ccGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATTTTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGt                                                                                                 >  1:1281725/1‑140 (MQ=255)
                                  ccGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGt                                                                                                 >  1:457558/1‑140 (MQ=255)
                                   cGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATTTTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTg                                                                                                >  1:990211/1‑140 (MQ=255)
                                   cGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATTTTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTg                                                                                                >  1:902765/1‑140 (MQ=255)
                                              ggcccccAGATTCGCGGATCAGAACTTTAAAACCCACGCAGCGCCAACCCAACCACCAGAACAATAACCGCCGCGGCGACCACCAGATTTTTATACCCCGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCaaa                                                                                     <  2:221154/136‑1 (MQ=255)
                                                cgcCAAAATTCTTCCATCAGCAACTTTAAACCCGCCAAGAGCCATACCAACCCACCAGACGCAAACCGCCGCGGCGACCACCAGATTTTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAgg                                                                                   <  2:5384/140‑1 (MQ=255)
                                                agcaaaGATTCTGTAACAAGAATTTTTAAAACCGCCAAAACCATCACCAACCCACCAGACGCTACCCGCCGCGGCGACCACCAGATTTTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAgg                                                                                   <  2:287438/135‑1 (MQ=255)
                                                agcCAAGATTCTGTCATCAGCAACTTAATACCCGCGAAGACCCATCCCAACCCACAAGACGCTACCCGCCGCGGCGACCACCAGATTTTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAgg                                                                                   <  2:29121/139‑1 (MQ=255)
                                                  ccaaGATTCTGTCATCAGCAACTTTATACCCGCGCAGAGCAATACCTACCCACCAGACGCTAACCGCCGCGGCGACCACCAGATTTTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCg                                                                                 <  2:675492/137‑1 (MQ=255)
                                                     aGATTCTGTAATCAGAAATTTTAAAACCGCGCAGAGCCATACCAACCCACCAGACGCAAACCGCCGCGGCGACCACCAGATTTTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCgatg                                                                              <  2:186617/140‑1 (MQ=255)
                                                             tcatcaGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCTGCGGCGACCACCAGATTTTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAg                                                                      >  1:38472/1‑140 (MQ=255)
                                                             gcatcaGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATTTTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAg                                                                      >  2:806908/2‑140 (MQ=255)
                                                               atcaGCAATTTTAAAACCCCCCAGAGCCATACCAAACCCCCAGACGCTAACCGCCGCGGCGACCACCAGATTTTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCg                                                                    <  2:1236652/140‑1 (MQ=255)
                                                               accaacAACTTTCAACCCGCGCAGAGACATCCCTAACCCCAAGCCCCCAACCGCCGCGGCGACCACCAGATATTTAACCCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCg                                                                    <  2:457558/135‑1 (MQ=255)
                                                                            taACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATTTTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGAttctt                                                       >  1:579222/1‑140 (MQ=255)
                                                                                   cgcAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATTTTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCAc                                                 >  2:229205/1‑140 (MQ=255)
                                                                                   cgcAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATTTTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCAc                                                 >  1:305587/1‑140 (MQ=255)
                                                                                        acccATACCAACCCACCAGACGATACCCGCCGCGGCGACCACAAGATTTTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGaaa                                            <  2:305587/138‑1 (MQ=255)
                                                                                             tACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATTTTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCt                                       >  2:507872/1‑140 (MQ=255)
                                                                                                  tAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATTTTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTa                                   <  1:229205/140‑1 (MQ=255)
                                                                                                    accaccAGACGCTAACCGCCGCGGCGACCACCAGATTTTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTAcc                                 >  2:806242/1‑140 (MQ=255)
                                                                                                    accacaAGCCGCAACCCGCCGCGGCGACCACCAGATTTTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTAcc                                 <  2:13499/140‑1 (MQ=255)
                                                                                                     cccccaGACGATAACCGCCGCGGCGACCCCCAGATTTTTATACAAAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCa                                <  2:820387/137‑1 (MQ=255)
                                                                                                     ccaccaGACGCTAACCGCCGCGGCGACCACCAGATTTTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCa                                >  2:1241523/1‑140 (MQ=255)
                                                                                                     ccaccaGACGCTAACCGCCGCGGCGACCACCAGATTTTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCa                                >  2:281007/1‑140 (MQ=255)
                                                                                                         caGACGCTAACCGCCGCGGCGACCACCAGATTTTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTgg                            >  1:995074/1‑140 (MQ=255)
                                                                                                         caGACGCTAACCGCCGCGGCGACCACCAGATTTTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTgg                            >  1:511008/1‑140 (MQ=255)
                                                                                                         caGACGCTAACCGCCGCGGCGACCACCAGATTTTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTgg                            >  1:620671/1‑140 (MQ=255)
                                                                                                          aGACGCTAACCGCCGCGGCGACCACCAGATTTTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGc                           >  2:257496/1‑140 (MQ=255)
                                                                                                           gACGCTAACCGCCGCGGCGACCACCAGATTTTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCa                          >  1:1296590/1‑140 (MQ=255)
                                                                                                           gACGCTAACCGCCGCGGCGACCACCAGATTTTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCa                          <  1:728563/140‑1 (MQ=255)
                                                                                                             aCGCTAACCGCCGCGGCGACCACCAGATTTTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCaa                         >  2:1001600/1‑140 (MQ=255)
                                                                                                             aCGCTAACCGCCGCGGCGACCACCAGATTTTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCaa                         <  1:133835/140‑1 (MQ=255)
                                                                                                                aaaaCCGCCGCGGCGACCACCAGATTTTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATAc                      <  2:355668/138‑1 (MQ=255)
                                                                                                                    ccgccgCGGCGACCACCAGATTTTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGa                  <  1:257496/140‑1 (MQ=255)
                                                                                                                      gccgcGGCGACCACCAGATTTTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAAt                >  1:992831/1‑140 (MQ=255)
                                                                                                                        cgcgGCGACCACCAGATTTTTATACCCAGCGAAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAAAGt              >  1:450729/1‑140 (MQ=255)
                                                                                                                           ggcgACCGCCAGATTTTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGTTTg           <  1:1138271/140‑1 (MQ=255)
                                                                                                                           ggcgACCACCAGATTTTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGTTTg           <  1:1001600/140‑1 (MQ=255)
                                                                                                                             cgACCACCAGATTTTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGTTTgcc         <  1:1190458/140‑1 (MQ=255)
                                                                                                                              gACCACCAGATTTTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGTTTgccg        <  1:575565/140‑1 (MQ=255)
                                                                                                                                    caGATTTTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGTTTGCCGCCTGGt  <  1:1116883/140‑1 (MQ=255)
                                                                                                                                    caGATTTTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGTTTGCCGCCTGGt  <  1:736964/140‑1 (MQ=255)
                                                                                                                                    caGATTTTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGTTTGCCGCCTGGt  <  1:861942/140‑1 (MQ=255)
                                                                                                                                    aaGATTTTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGTTTGCCGCCTGGt  <  2:620671/139‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                      
CACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAACCGCGCAGAGCCATACCTAACCACCAGACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATACCGGAATGTTTGCCGCCTGGT  >  NC_000913/446875‑447143

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

GATK/CNVnator alignment

N/A