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 2,063,343 T→C A148A (GCA→GCG)  cobU ← cobinamide kinase/cobinamide phosphate guanylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,063,3430TC99.4% 493.7 / ‑4.9 155A148A (GCA→GCG) cobUcobinamide kinase/cobinamide phosphate guanylyltransferase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (84/70);  total (85/70)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.18e-01

Only 100 of 161 total aligned reads displayed.

TACTCCAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAA  >  NC_000913/2063253‑2063432
                                                                                          |                                                                                         
tACTCCAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCgcgc                                                                                         >  1:4010272/1‑93 (MQ=255)
  cgccAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCa                                                                                       <  1:386042/91‑1 (MQ=255)
  cTCCAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCa                                                                                       >  1:1039963/1‑93 (MQ=255)
  cTCCAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCa                                                                                       <  1:3983862/93‑1 (MQ=255)
   tCCAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAg                                                                                      >  2:192682/1‑93 (MQ=255)
   tCCAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAg                                                                                      >  1:245181/1‑93 (MQ=255)
    ccAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAga                                                                                     >  2:812638/1‑93 (MQ=255)
    ccAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAga                                                                                     >  2:2404602/1‑93 (MQ=255)
     cAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAgac                                                                                    <  2:3799308/93‑1 (MQ=255)
      aaTACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCAGGAAAATGTCGCGCCAgacg                                                                                   >  2:1169464/1‑93 (MQ=255)
       aTACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAgacga                                                                                  >  2:2828106/1‑93 (MQ=255)
       aTACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAgacga                                                                                  >  1:1620923/1‑93 (MQ=255)
       aTACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAgacga                                                                                  >  2:2135736/1‑93 (MQ=255)
         aCCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGAct                                                                                >  1:3445220/1‑93 (MQ=255)
          cccGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGActc                                                                               >  1:1681394/1‑93 (MQ=255)
          cccGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCACACGActc                                                                               >  1:2565478/1‑93 (MQ=255)
             gAACCCCCCAGCCATACTTCATTTGCCGCAGCGGCCGACTGCTGATTTGCCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTcc                                                                            <  2:1260843/93‑1 (MQ=255)
             gAAACCCCCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTcc                                                                            <  1:3876332/93‑1 (MQ=255)
             gAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACtctca                                                                            >  2:1788490/1‑92 (MQ=255)
             gAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTcc                                                                            >  2:3040528/1‑93 (MQ=255)
                accaccAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGc                                                                         >  2:4630950/1‑93 (MQ=255)
                   cccaGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGcaca                                                                      <  1:3100301/92‑1 (MQ=255)
                   cccaGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGcaca                                                                      <  2:1229704/92‑1 (MQ=255)
                   cccaGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGcaca                                                                      <  2:1229702/92‑1 (MQ=255)
                   accaGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGcaca                                                                      >  2:2471473/1‑93 (MQ=255)
                     caGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAAt                                                                    <  2:4453986/93‑1 (MQ=255)
                      aGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATc                                                                   <  2:2044425/93‑1 (MQ=255)
                      aGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATc                                                                   >  2:1779497/1‑93 (MQ=255)
                       gCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATcc                                                                  >  1:2610679/1‑93 (MQ=255)
                        ccATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATccc                                                                 >  1:2206810/1‑93 (MQ=255)
                        ccATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATccc                                                                 <  1:4576138/93‑1 (MQ=255)
                           tACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCAt                                                              <  2:3066031/93‑1 (MQ=255)
                            aCTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCAtc                                                             >  2:4019906/1‑92 (MQ=255)
                                  tttGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACt                                                       >  1:4409739/1‑93 (MQ=255)
                                   ttGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACtt                                                      <  1:2730754/93‑1 (MQ=255)
                                   ttGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACtt                                                      <  2:1500536/93‑1 (MQ=255)
                                   ttGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACtt                                                      >  2:3602751/1‑93 (MQ=255)
                                           gCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTc                                              <  1:3533109/93‑1 (MQ=255)
                                             ggCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCAc                                            <  1:3010892/93‑1 (MQ=255)
                                              gCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACt                                           >  1:1775847/1‑93 (MQ=255)
                                                cacctgctgATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTaa                                         <  2:4213634/90‑1 (MQ=255)
                                                cAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTaa                                         <  1:2159262/93‑1 (MQ=255)
                                                 aaCTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAAt                                        >  2:16202/1‑93 (MQ=255)
                                                   ctgctgATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATAc                                      >  1:2125759/1‑93 (MQ=255)
                                                   ctgctgATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATAc                                      <  2:4643175/93‑1 (MQ=255)
                                                    tgctgATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACa                                     >  1:763429/1‑93 (MQ=255)
                                                     gctgATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACaa                                    >  1:4704950/1‑93 (MQ=255)
                                                       tgATTTACCCGCCCGGCAATATCCCGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAAcc                                  <  2:2583710/93‑1 (MQ=255)
                                                       tgATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAAcc                                  <  1:4567926/93‑1 (MQ=255)
                                                       tgATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAAcc                                  >  2:3012717/1‑93 (MQ=255)
                                                        gATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCt                                 >  2:3524311/1‑93 (MQ=255)
                                                        gATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCt                                 >  2:2986742/1‑93 (MQ=255)
                                                         aTTTCCCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCtt                                <  2:3180356/93‑1 (MQ=255)
                                                          tttACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCttt                               >  1:2764866/1‑93 (MQ=255)
                                                           ttACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTg                              <  2:596030/93‑1 (MQ=255)
                                                              cccgcccgGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCgg                           <  1:1900202/93‑1 (MQ=255)
                                                              cccgcccgGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCgg                           <  2:1096927/93‑1 (MQ=255)
                                                              cccgcccgGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCgg                           <  2:1480110/93‑1 (MQ=255)
                                                              cccgcccgGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCgg                           <  1:3400807/93‑1 (MQ=255)
                                                              cccgcccgGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCgg                           <  1:3164283/93‑1 (MQ=255)
                                                               ccgcccgGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCggg                          <  2:1972551/93‑1 (MQ=255)
                                                               ccgcccgGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCggg                          <  2:4738482/93‑1 (MQ=255)
                                                                cgcccgGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCgggg                         <  2:1531628/93‑1 (MQ=255)
                                                                 gcccgGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTCGTCACTAATACCACCTTTGCGGGGc                        >  1:2302329/1‑93 (MQ=255)
                                                                 gcccgGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGc                        <  2:1282051/93‑1 (MQ=255)
                                                                   ccgGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCaa                      >  2:3383343/1‑93 (MQ=255)
                                                                   ccgGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCaa                      <  1:280502/93‑1 (MQ=255)
                                                                   ccgGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCaa                      <  2:3771810/93‑1 (MQ=255)
                                                                   ccgGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCaa                      >  2:1625370/1‑93 (MQ=255)
                                                                    cgGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAAc                     >  2:2135591/1‑93 (MQ=255)
                                                                    cgGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAAc                     <  2:248038/93‑1 (MQ=255)
                                                                    cgGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAAc                     <  1:711122/93‑1 (MQ=255)
                                                                    cgGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAAc                     <  2:3441621/93‑1 (MQ=255)
                                                                    cgGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAAc                     >  2:1288349/1‑93 (MQ=255)
                                                                     ggCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACg                    >  2:4527374/1‑93 (MQ=255)
                                                                      gCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGt                   <  2:420924/93‑1 (MQ=255)
                                                                       cAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACttt                  <  2:1503666/93‑1 (MQ=255)
                                                                       cAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGtt                  <  2:563466/93‑1 (MQ=255)
                                                                        aaTATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTg                 >  2:2303049/1‑93 (MQ=255)
                                                                        aaTATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTg                 >  2:1382297/1‑93 (MQ=255)
                                                                         atatCACGAAAATGTCGCGCCAGCCGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTgg                <  2:899754/93‑1 (MQ=255)
                                                                         atatCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTgg                <  2:2316081/93‑1 (MQ=255)
                                                                         atatCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTgg                >  1:376229/1‑93 (MQ=255)
                                                                          tatCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTggc               >  2:1707089/1‑93 (MQ=255)
                                                                          aatCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTggc               <  1:4079551/92‑1 (MQ=255)
                                                                           atCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTggca              >  2:2005214/1‑93 (MQ=255)
                                                                           atCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTggca              <  1:4338537/93‑1 (MQ=255)
                                                                           atCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTggca              <  1:2469048/93‑1 (MQ=255)
                                                                           atCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACATTggca              >  2:1765626/1‑93 (MQ=255)
                                                                               cGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTAACTAATACAACCTTTGCGGGGCAACGTTggcaggca          >  1:1872990/1‑93 (MQ=255)
                                                                                 aaaaTGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGgcagc        >  2:3337902/1‑93 (MQ=255)
                                                                                 aaaaTGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGgcagc        >  1:3534923/1‑93 (MQ=255)
                                                                                  aaaTGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGgcagca       >  1:505281/1‑93 (MQ=255)
                                                                                    aTGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCCCTAATACAACCTTTGCGGTGCAACGTTGGCAGGCAGCAAt     <  2:770811/93‑1 (MQ=255)
                                                                                    aTGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAAt     >  1:3119109/1‑93 (MQ=255)
                                                                                    aTGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCATTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAAt     >  1:1098692/1‑93 (MQ=255)
                                                                                     tGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATc    >  2:4593668/1‑93 (MQ=255)
                                                                                      gTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCa   <  2:1999887/93‑1 (MQ=255)
                                                                                       tCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCaa  <  1:2814400/93‑1 (MQ=255)
                                                                                       tCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCaa  <  2:2596175/93‑1 (MQ=255)
                                                                                          |                                                                                         
TACTCCAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAA  >  NC_000913/2063253‑2063432

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
TTTACTCCAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGAC  >  NC_000913/2063251‑2063436
                                                                                            |                                                                                             
TTTACTCCAATACCCGAACCCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCG                                                                                              <  SRR4302544.4537505/94‑1 (MQ=60)
 TTACTCCAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGC                                                                                             >  SRR4302544.4050167/1‑94 (MQ=60)
   ACTCCAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCA                                                                                           >  SRR4302544.2960693/1‑94 (MQ=60)
   ACTCCAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCA                                                                                           >  SRR4302544.1050288/1‑94 (MQ=60)
    CTCCAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAG                                                                                          >  SRR4302544.247744/1‑94 (MQ=60)
    CTCCAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAG                                                                                          >  SRR4302544.1188318/1‑94 (MQ=60)
    CTCCAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAG                                                                                          <  SRR4302544.4023516/94‑1 (MQ=60)
    CGCCAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAG                                                                                          <  SRR4302544.389967/94‑1 (MQ=60)
      CCAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGAC                                                                                        <  SRR4302544.2801353/94‑1 (MQ=60)
        AATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGA                                                                                      >  SRR4302544.1636888/1‑94 (MQ=60)
          TACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACT                                                                                    >  SRR4302544.3478138/1‑94 (MQ=60)
           ACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTC                                                                                   >  SRR4302544.1697961/1‑94 (MQ=60)
           ACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCACACGACTC                                                                                   >  SRR4302544.2590125/1‑94 (MQ=60)
              CGAACCCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCC                                                                                <  SRR4302544.1875579/94‑1 (MQ=60)
              CGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCC                                                                                >  SRR4302544.3748460/1‑94 (MQ=60)
              CGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCC                                                                                >  SRR4302544.2180948/1‑94 (MQ=60)
               GAAACCCCCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCG                                                                               <  SRR4302544.3914821/94‑1 (MQ=60)
                    CACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACA                                                                          >  SRR4302544.2049577/1‑94 (MQ=60)
                     CCCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAA                                                                         <  SRR4302544.3129814/94‑1 (MQ=60)
                       CAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATC                                                                       >  SRR4302544.4381063/1‑94 (MQ=60)
                        AGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCC                                                                      >  SRR4302544.2635754/1‑94 (MQ=60)
                         GCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCC                                                                     >  SRR4302544.2716452/1‑94 (MQ=60)
                         GCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCC                                                                     >  SRR4302544.2227588/1‑94 (MQ=60)
                          CCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCC                                                                    <  SRR4302544.4621738/94‑1 (MQ=60)
                                   ATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACT                                                           >  SRR4302544.4454268/1‑94 (MQ=60)
                                     TTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTC                                                         <  SRR4302544.2757073/94‑1 (MQ=60)
                                     TCGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTC                                                         >  SRR4302544.4051906/1‑94 (MQ=60)
                                             GCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCA                                                 <  SRR4302544.3566908/94‑1 (MQ=60)
                                               GGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACT                                               >  SRR4302544.1793230/1‑94 (MQ=60)
                                               GGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACT                                               <  SRR4302544.3039750/94‑1 (MQ=60)
                                                  CAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAAT                                            <  SRR4302544.2179778/94‑1 (MQ=60)
                                                    ACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATAC                                          >  SRR4302544.2146143/1‑94 (MQ=60)
                                                     CTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACA                                         >  SRR4302544.770870/1‑94 (MQ=60)
                                                      TGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAA                                        >  SRR4302544.4751444/1‑94 (MQ=60)
                                                       GCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTAACTAATACAAC                                       >  SRR4302544.3152019/1‑94 (MQ=60)
                                                         TGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCT                                     <  SRR4302544.4613463/94‑1 (MQ=60)
                                                          GATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTT                                    >  SRR4302544.3961852/1‑94 (MQ=60)
                                                           ATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTT                                   >  SRR4302544.2791457/1‑94 (MQ=60)
                                                            TTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTG                                  >  SRR4302544.930298/1‑94 (MQ=60)
                                                               ACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGG                               >  SRR4302544.4652615/1‑94 (MQ=60)
                                                                CCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGG                              <  SRR4302544.321225/94‑1 (MQ=60)
                                                                CCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGG                              <  SRR4302544.3433325/94‑1 (MQ=60)
                                                                CCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGG                              <  SRR4302544.1918561/94‑1 (MQ=60)
                                                                CCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGG                              <  SRR4302544.3194385/94‑1 (MQ=60)
                                                                 CCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGG                             <  SRR4302544.3049611/94‑1 (MQ=60)
                                                                 CCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGG                             <  SRR4302544.4485788/94‑1 (MQ=60)
                                                                  CGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTCGTCACTAATACCACCTTTGCGGGGC                            >  SRR4302544.2324072/1‑94 (MQ=60)
                                                                  CGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGC                            >  SRR4302544.1700643/1‑94 (MQ=60)
                                                                   GCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCA                           >  SRR4302544.349344/1‑94 (MQ=60)
                                                                     CCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAAC                         <  SRR4302544.283433/94‑1 (MQ=60)
                                                                      CGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACG                        <  SRR4302544.718171/94‑1 (MQ=60)
                                                                          AATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGG                    >  SRR4302544.380066/1‑94 (MQ=60)
                                                                          AATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGG                    <  SRR4302544.2478008/94‑1 (MQ=60)
                                                                            AATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCA                  <  SRR4302544.4120139/94‑1 (MQ=60)
                                                                             ATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAG                 <  SRR4302544.2492640/94‑1 (MQ=60)
                                                                             ATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAG                 <  SRR4302544.2277192/94‑1 (MQ=60)
                                                                             ATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAG                 <  SRR4302544.4382688/94‑1 (MQ=60)
                                                                               CACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGC               >  SRR4302544.3809247/1‑94 (MQ=60)
                                                                                ACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTAACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCA              >  SRR4302544.1891091/1‑94 (MQ=60)
                                                                                  GAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGC            >  SRR4302544.260628/1‑94 (MQ=60)
                                                                                  GAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGC            >  SRR4302544.3568736/1‑94 (MQ=60)
                                                                                   AAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCA           >  SRR4302544.510245/1‑94 (MQ=60)
                                                                                     AATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAAT         >  SRR4302544.3148772/1‑94 (MQ=60)
                                                                                     AATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCATTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAAT         >  SRR4302544.1109628/1‑94 (MQ=60)
                                                                                       TGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCA       <  SRR4302544.735621/94‑1 (MQ=60)
                                                                                         TCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAAC     <  SRR4302544.2841562/94‑1 (MQ=60)
                                                                                         TCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAAC     <  SRR4302544.503365/94‑1 (MQ=60)
                                                                                         TCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAAC     >  SRR4302544.4017944/1‑94 (MQ=60)
                                                                                          CGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACG    <  SRR4302544.213402/94‑1 (MQ=60)
                                                                                           GCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGA   >  SRR4302544.2842315/1‑94 (MQ=60)
                                                                                            CGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGAC  >  SRR4302544.1621363/1‑94 (MQ=60)
                                                                                            CGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGAC  >  SRR4302544.4207214/1‑94 (MQ=60)
                                                                                            CGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATAAAACCTTTGCGGGGCAACGTTGGCGGGAGGAAATCAACGAC  >  SRR4302544.1621246/1‑94 (MQ=60)
                                                                                            |                                                                                             
TTTACTCCAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAACGAC  >  NC_000913/2063251‑2063436

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