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 F1 I1 R1 1357 179.3 9500450 96.1% 9129932 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,953,902 T→C Q377R (CAG→CGG)  b1872 ← biotin sulfoxide reductase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,953,9020TC96.5% 347.0 / 1.5 118Q377R (CAG→CGG) b1872biotin sulfoxide reductase 2
Reads supporting (aligned to +/- strand):  ref base T (3/1);  new base C (54/58);  total (57/61)
Fisher's exact test for biased strand distribution p-value = 3.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01

Only 100 of 125 total aligned reads displayed.

TTTCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCAT  >  NC_000913/1953813‑1953990
                                                                                         |                                                                                        
tttCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGcc                                                                                       <  2:2245044/93‑1 (MQ=255)
tttCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGcc                                                                                       <  1:2138454/93‑1 (MQ=255)
tttCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGcc                                                                                       <  2:2245043/93‑1 (MQ=255)
tttCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGcc                                                                                       >  2:2685911/1‑93 (MQ=255)
tttCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGcc                                                                                       >  2:3164603/1‑93 (MQ=255)
tttCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGcc                                                                                       >  1:3226506/1‑93 (MQ=255)
    cGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTaa                                                                                   >  2:1688574/1‑93 (MQ=255)
     ggCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAAc                                                                                  >  1:4533467/1‑93 (MQ=255)
     ggCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAAc                                                                                  >  2:1572830/1‑93 (MQ=255)
       cagcaCGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACAt                                                                                <  2:4085077/93‑1 (MQ=255)
        gacaCGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAATCCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATc                                                                               <  2:3122183/91‑1 (MQ=255)
        agcaCGCCGCCAACACGCGTCGTTTTGCCACCGTTGGAGTAATGATAAGAGAACCCAAACCCGCCGCCCGGTGTACCAATCCGCCCTAACATc                                                                               <  2:4717585/93‑1 (MQ=255)
        agcaCGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATc                                                                               >  2:3522581/1‑93 (MQ=255)
        agcaCGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATc                                                                               <  2:355880/93‑1 (MQ=255)
        agcaCGCCGCCAACACGAGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATc                                                                               >  1:4296456/1‑93 (MQ=255)
         gcaCGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATcg                                                                              >  1:738633/1‑93 (MQ=255)
          caCGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATcgc                                                                             <  1:1331805/93‑1 (MQ=255)
           aCGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATcgcc                                                                            <  1:1783775/93‑1 (MQ=255)
            cgccgcCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATcgccg                                                                           >  1:4568294/1‑93 (MQ=255)
             gccgccAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATcgccgc                                                                          <  1:2841678/93‑1 (MQ=255)
               cgccAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCa                                                                        <  1:3574010/93‑1 (MQ=255)
                   aaCACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGt                                                                    <  2:69178/93‑1 (MQ=255)
                    acacGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTc                                                                   >  1:904858/1‑93 (MQ=255)
                     cacGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCa                                                                  <  2:3803073/93‑1 (MQ=255)
                      acGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCAc                                                                 >  1:1675169/1‑93 (MQ=255)
                      acGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCAc                                                                 <  1:340722/93‑1 (MQ=255)
                       cgcgTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCAcc                                                                <  1:1008799/93‑1 (MQ=255)
                        gcgTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCa                                                               >  2:2440505/1‑93 (MQ=255)
                           tCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCa                                                            >  1:1814200/1‑93 (MQ=255)
                           tCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCa                                                            <  2:2626385/93‑1 (MQ=255)
                              ggTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATcc                                                         <  1:1553002/93‑1 (MQ=255)
                                 tGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAgt                                                      <  1:4225181/93‑1 (MQ=255)
                                  gCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAgtg                                                     >  2:35996/1‑93 (MQ=255)
                                  gCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAgtg                                                     >  2:2900200/1‑93 (MQ=255)
                                   ccACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAgtgt                                                    <  2:1070704/93‑1 (MQ=255)
                                   ccACCGGTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAgtgt                                                    <  2:4534237/93‑1 (MQ=255)
                                    cACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGtt                                                   >  2:3708609/1‑93 (MQ=255)
                                      ccGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTc                                                 <  2:4599681/93‑1 (MQ=255)
                                       cGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCt                                                <  2:3904043/93‑1 (MQ=255)
                                        gTTGGAGTAATGATAAGAGAAACCGAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTg                                               >  2:4554996/1‑93 (MQ=255)
                                        gTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTg                                               <  1:3716307/93‑1 (MQ=255)
                                         ttGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGt                                              <  1:460033/93‑1 (MQ=255)
                                         ttGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGt                                              >  2:2885484/1‑93 (MQ=255)
                                          tGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGtt                                             <  2:847575/93‑1 (MQ=255)
                                           ggAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTc                                            <  1:2657613/93‑1 (MQ=255)
                                           ggAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTc                                            >  2:3177407/1‑93 (MQ=255)
                                            gAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGCGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCa                                           >  1:1891094/1‑93 (MQ=255)
                                              gTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCAcc                                         >  1:2781113/1‑93 (MQ=255)
                                              gTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCAcc                                         <  1:2233343/93‑1 (MQ=255)
                                               tAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCa                                        >  1:1156337/1‑93 (MQ=255)
                                                aaTGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCat                                       >  2:2644905/1‑93 (MQ=255)
                                                aaTGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCat                                       >  1:1756066/1‑93 (MQ=255)
                                                  tGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCatat                                     <  1:4326466/93‑1 (MQ=255)
                                                   gATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATAtt                                    >  1:1317585/1‑93 (MQ=255)
                                                    aTAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTg                                   <  2:10282/93‑1 (MQ=255)
                                                     aaaGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGc                                  <  1:3591295/92‑1 (MQ=255)
                                                      aaGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCt                                 <  2:127711/93‑1 (MQ=255)
                                                      aaGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCt                                 >  2:648307/1‑93 (MQ=255)
                                                      aaGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCt                                 <  1:2511149/93‑1 (MQ=255)
                                                      aaGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAAAATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCt                                 >  1:2430788/1‑93 (MQ=255)
                                                       agagAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTg                                <  2:1658021/93‑1 (MQ=255)
                                                       agagAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTg                                <  1:2398592/93‑1 (MQ=255)
                                                       agagAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTg                                >  2:2089836/1‑93 (MQ=255)
                                                         agaAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTCTCTGTTCACCATATTGCTGgc                              >  1:1752553/1‑93 (MQ=255)
                                                          gaAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGgcg                             >  1:2826875/1‑93 (MQ=255)
                                                            aaCCAAACCCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCt                           <  2:3351765/93‑1 (MQ=255)
                                                            aaCCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCt                           <  2:3595843/93‑1 (MQ=255)
                                                             aCCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCta                          >  1:1018494/1‑92 (MQ=255)
                                                             aCCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTg                          >  2:1566675/1‑93 (MQ=255)
                                                             aCCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTg                          >  1:1018493/1‑93 (MQ=255)
                                                               cAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGaa                        <  1:783479/93‑1 (MQ=255)
                                                                aaaaCCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAAt                       <  2:4380332/93‑1 (MQ=255)
                                                                 aaaCCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAAtt                      <  1:914275/93‑1 (MQ=255)
                                                                  aaCCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTc                     >  1:4309892/1‑93 (MQ=255)
                                                                   aCCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTcc                    >  2:3826033/1‑93 (MQ=255)
                                                                    ccgccgCCCGGTGTACCAATCGGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTccc                   <  1:1337102/93‑1 (MQ=255)
                                                                    ccgccgCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACTATATTGCTGGCGCTGAATTccc                   <  2:2990198/93‑1 (MQ=255)
                                                                    ccgccgCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGAATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTccc                   >  2:2567723/1‑93 (MQ=255)
                                                                      gccgccCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCa                 >  2:645055/1‑93 (MQ=255)
                                                                      gccgccCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCa                 <  1:2216718/93‑1 (MQ=255)
                                                                      cccgccCGGTGTACCAATCCGCCCTAGCATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCa                 <  1:1070634/92‑1 (MQ=255)
                                                                       ccgccCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAg                >  1:1625303/1‑93 (MQ=255)
                                                                       ccgccCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAg                <  1:3658545/93‑1 (MQ=255)
                                                                         gccCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGcc              <  1:1765622/93‑1 (MQ=255)
                                                                         gccCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGcc              <  2:475953/93‑1 (MQ=255)
                                                                         gccCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGcc              <  2:475951/93‑1 (MQ=255)
                                                                          cccGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCt             >  1:4071992/1‑93 (MQ=255)
                                                                           ccGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTg            >  1:2111999/1‑93 (MQ=255)
                                                                           ccGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTg            <  1:1497715/93‑1 (MQ=255)
                                                                            cGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGc           >  2:4006645/1‑93 (MQ=255)
                                                                             ggTGTCCCAATCCTCCCTAACATCGCCGCCAGCGTCACCAGCATCAAGTGTTTCTGTTCCCCATATTGCTGGCGCTGAATTCCCCAGCCTGcc          <  2:3709513/93‑1 (MQ=25)
                                                                                 tACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCa      <  1:3122384/93‑1 (MQ=255)
                                                                                   gcAATCCGCCCTAACATCCGCGCCAGCGTACCCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGc    <  1:1564988/92‑1 (MQ=25)
                                                                                   ccAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGc    >  1:3118841/1‑93 (MQ=255)
                                                                                   ccAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGc    >  2:1423065/1‑93 (MQ=255)
                                                                                   ccAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGc    <  1:1005745/93‑1 (MQ=255)
                                                                                    cAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCa   >  2:4051882/1‑93 (MQ=255)
                                                                                    cAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCa   <  1:2531912/93‑1 (MQ=255)
                                                                                     aaTCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCAt  >  2:3998792/1‑93 (MQ=255)
                                                                                     aaTCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCAt  <  2:2770844/93‑1 (MQ=255)
                                                                                         |                                                                                        
TTTCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCAT  >  NC_000913/1953813‑1953990

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
ATTTCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTG  >  NC_000913/1953812‑1953994
                                                                                          |                                                                                            
ATTTCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCC                                                                                           >  SRR4302544.3257356/1‑94 (MQ=60)
 TTTCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCC                                                                                          <  SRR4302544.2158880/94‑1 (MQ=60)
     CGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAAC                                                                                      >  SRR4302544.4578804/1‑94 (MQ=60)
        CAGCACGCCGCCAACACGAGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATC                                                                                   >  SRR4302544.4340306/1‑94 (MQ=60)
         AGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCG                                                                                  >  SRR4302544.745860/1‑94 (MQ=60)
           CACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCC                                                                                <  SRR4302544.1344912/94‑1 (MQ=60)
           CACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGACGCCCGGTGTACCAATCCGCCCTAACATCGCC                                                                                <  SRR4302544.2660235/94‑1 (MQ=60)
            ACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCG                                                                               <  SRR4302544.1801216/94‑1 (MQ=60)
            ACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCG                                                                               >  SRR4302544.4613834/1‑94 (MQ=60)
              GCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCC                                                                             <  SRR4302544.2869108/94‑1 (MQ=60)
                CGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAG                                                                           <  SRR4302544.3608286/94‑1 (MQ=60)
                CGCCAACACGCGTCGGGTGGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCAGGTGTACCAATCCGCCCTAACATCGCCGCCAG                                                                           <  SRR4302544.2502227/94‑1 (MQ=60)
                   CAACACGCGTCGGGTTGCCACCGTT‑GAGGAATGATAAGAGAAACCAAAACCGCCGCCCGGGGTACCAATCCGCCCTAACATCGCCGCCCGCGTC                                                                       >  SRR4302544.30728/1‑94 (MQ=60)
                    AACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTC                                                                       >  SRR4302544.913812/1‑94 (MQ=60)
                      CACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCAC                                                                     >  SRR4302544.1691669/1‑94 (MQ=60)
                       ACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACC                                                                    <  SRR4302544.344249/94‑1 (MQ=60)
                        CGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCA                                                                   <  SRR4302544.1018837/94‑1 (MQ=60)
                         GCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAG                                                                  <  SRR4302544.1781022/94‑1 (MQ=60)
                           GTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCA                                                                >  SRR4302544.1831910/1‑94 (MQ=60)
                             CGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGCTGTACCAATCCGCCCTAACATCGCCGACAGCGTCACCAGCATC                                                              >  SRR4302544.2388323/1‑94 (MQ=60)
                               GGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCA                                                            <  SRR4302544.1568459/94‑1 (MQ=60)
                                  TGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTG                                                         <  SRR4302544.4268135/94‑1 (MQ=60)
                                         GTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGT                                                  <  SRR4302544.3753274/94‑1 (MQ=60)
                                          TTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTT                                                 <  SRR4302544.464554/94‑1 (MQ=60)
                                           TGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCACCGTCACCAGCATCCAGTGTTTCTGTTC                                                >  SRR4302544.1608607/1‑94 (MQ=60)
                                            GGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCA                                               <  SRR4302544.2683103/94‑1 (MQ=60)
                                            GGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGCGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCA                                               >  SRR4302544.1909364/1‑94 (MQ=60)
                                              AGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACC                                             >  SRR4302544.2807882/1‑94 (MQ=60)
                                               GTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCA                                            <  SRR4302544.2254348/94‑1 (MQ=60)
                                               GTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCA                                            >  SRR4302544.1167845/1‑94 (MQ=60)
                                                TAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCAT                                           >  SRR4302544.1773255/1‑94 (MQ=60)
                                                 AATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATA                                          >  SRR4302544.4360617/1‑94 (MQ=60)
                                                   TGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATT                                        <  SRR4302544.4370525/94‑1 (MQ=60)
                                                   TGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATT                                        >  SRR4302544.1330540/1‑94 (MQ=60)
                                                      TAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAAAATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCT                                     >  SRR4302544.2453942/1‑94 (MQ=60)
                                                      AAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCT                                     <  SRR4302544.3625781/94‑1 (MQ=60)
                                                      AAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCT                                     <  SRR4302544.4002659/94‑1 (MQ=60)
                                                       AAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTG                                    <  SRR4302544.2535248/94‑1 (MQ=60)
                                                       AAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTG                                    >  SRR4302544.2115540/1‑94 (MQ=60)
                                                        AGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGG                                   <  SRR4302544.2421341/94‑1 (MQ=60)
                                                         GAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTCTCTGTTCACCATATTGCTGGC                                  >  SRR4302544.1769708/1‑94 (MQ=60)
                                                          AGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCG                                 >  SRR4302544.2854190/1‑94 (MQ=60)
                                                             AACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTG                              >  SRR4302544.1028662/1‑94 (MQ=60)
                                                             AACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTA                              >  SRR4302544.1028663/1‑94 (MQ=60)
                                                                CAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAAT                           <  SRR4302544.791065/94‑1 (MQ=60)
                                                                  AAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTC                         <  SRR4302544.923311/94‑1 (MQ=60)
                                                                  AAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTC                         <  SRR4302544.546651/94‑1 (MQ=60)
                                                                  AAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTC                         >  SRR4302544.4353822/1‑94 (MQ=60)
                                                                     CCGCCGCCCGGTGTACCAATCGGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCC                      <  SRR4302544.1350288/94‑1 (MQ=60)
                                                                     CCGCCGCCCGGTGTACCAATCGGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCC                      <  SRR4302544.1350291/94‑1 (MQ=60)
                                                                       GCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAG                    <  SRR4302544.3450273/94‑1 (MQ=60)
                                                                       GCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAG                    <  SRR4302544.2237568/94‑1 (MQ=60)
                                                                       GCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAG                    >  SRR4302544.1641317/1‑94 (MQ=60)
                                                                       CCCGCCCGGTGTACCAATCCGCCCTAGCATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAG                    <  SRR4302544.1081275/94‑1 (MQ=60)
                                                                        CCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGC                   <  SRR4302544.3695046/94‑1 (MQ=60)
                                                                          GCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCT                 <  SRR4302544.1782900/94‑1 (MQ=60)
                                                                          GCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCT                 >  SRR4302544.4112516/1‑94 (MQ=60)
                                                                           CCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTG                >  SRR4302544.2132304/1‑94 (MQ=60)
                                                                            CCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGC               <  SRR4302544.1512522/94‑1 (MQ=60)
                                                                             CGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCC              >  SRR4302544.106462/1‑94 (MQ=60)
                                                                                TGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATC           >  SRR4302544.3103074/1‑94 (MQ=60)
                                                                                  TACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAG         <  SRR4302544.3152066/94‑1 (MQ=60)
                                                                                   ACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGC        >  SRR4302544.3148502/1‑94 (MQ=60)
                                                                                    GCAATCCGCCCTAACATCCGCGCCAGCGTACCCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCA       <  SRR4302544.1580508/94‑1 (MQ=60)
                                                                                    CCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCA       <  SRR4302544.1015738/94‑1 (MQ=60)
                                                                                     CAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCAT      <  SRR4302544.2556239/94‑1 (MQ=60)
                                                                                         CCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTG  >  SRR4302544.3279925/1‑94 (MQ=60)
                                                                                          |                                                                                            
ATTTCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATTGTG  >  NC_000913/1953812‑1953994

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 2 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑