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 4,203,198 Δ1 bp coding (312/1290 nt) purD ← phosphoribosylglycinamide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,203,1950T.94.6% 300.8 / 10.4 74coding (315/1290 nt)purDphosphoribosylglycinamide synthetase
Reads supporting (aligned to +/- strand):  ref base T (2/2);  new base . (28/42);  total (30/44)
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.97e-01

Only 100 of 103 total aligned reads displayed.

GCAGATACGCCAGCGCAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCTTTCACCA  >  NC_000913/4203104‑4203283
                                                                                            |                                                                                         
gCAGATACGCCAGCGCAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGGCAGGAAATCTTTTGTAAACGCtt                                                                                          >  1:4311194/1‑93 (MQ=255)
ccAGATACGCCAGCGCAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCtt                                                                                          <  2:3884558/92‑1 (MQ=255)
   gATACGCCAGCGCAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAGAAATTTGGTAAAAGCtttgg                                                                                       >  1:3743473/1‑91 (MQ=255)
   gATACGCCAGCGCAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGACAGGAAAGCTTTTGTTTAACGCtttg                                                                                        >  1:1897404/1‑92 (MQ=255)
     tACGCCAGCGCAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGGCAGGAAATCTTTTGTAAAAGCttttagc                                                                                     >  2:1812737/1‑90 (MQ=255)
      aCGCCAGCGCAGGTTCTACCTCGGTGAAGGTCTGGTATGCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGcc                                                                                   <  1:3373065/93‑1 (MQ=255)
           aGCGCAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCCGGGAATCTTTGGGAAACGCTTTTAGCCTTccc                                                                              >  1:844047/1‑92 (MQ=255)
           aGCGCAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCa                                                                              <  1:1302582/93‑1 (MQ=255)
           aGCGCAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAAACTTTGGTAAAACCTTTTAGCCTTCCa                                                                              >  1:1014776/1‑93 (MQ=255)
            gcgcAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGGCCGGAAATCTTTTGTAAACGCTTTTAGCCTTCCAg                                                                             >  1:74542/1‑93 (MQ=255)
            gcgcAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGGAATATTTGGTAAACGCTTTGAGGCTTcccg                                                                             >  1:4278475/1‑91 (MQ=255)
              gcAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGGCAGGAAATCTTTTGTAAACGCTTTGAGACTTTCAggt                                                                           >  2:250400/1‑91 (MQ=255)
               cAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGGAAACTTTGGTAAACGCTTTGAGCCCTCCCGttt                                                                          >  2:3531748/1‑92 (MQ=255)
                acgTTCTACCGCGGCGAAGTTCTGTTGTTCCGCCCTAGGGATCTTATGGCGGCCCAGTAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTgg                                                                         <  2:1090301/91‑1 (MQ=255)
                aGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCCGGAAATCTTTGGTAAACGCTTTGAGCCTTCCCGTTgg                                                                         >  2:2908173/1‑93 (MQ=255)
                aGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTgg                                                                         <  2:191762/93‑1 (MQ=255)
                  gTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGc                                                                       <  1:379072/93‑1 (MQ=255)
                      tACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTTGTAAACGCTTTGAGCCTTCCAGTTGGGGCGCa                                                                   >  2:3305666/1‑93 (MQ=255)
                        ccTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCAcc                                                                 <  2:1400619/93‑1 (MQ=255)
                        ccTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCAcc                                                                 <  1:3717824/93‑1 (MQ=255)
                        ccTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCAcc                                                                 <  1:3717823/93‑1 (MQ=255)
                         cTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTTGTAAACGCTTTTAGCCCTCCCGTTGGGCCGCAccc                                                                >  1:2457565/1‑92 (MQ=255)
                         cTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTTGAGCCTTCCAGTTGGGCCGCAcc                                                                 <  1:4024948/93‑1 (MQ=255)
                           cGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGc                                                              <  2:3834550/93‑1 (MQ=255)
                           cGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGc                                                              <  1:3237617/93‑1 (MQ=255)
                           cGGTGAAGTTCTGGTATTCCGCAGTAGGGATCTTATGGCGGGCCAGGGAATCTTTGGTAAAAGCTTTGAGCCTTCCAGTTGGGGCGCACCTGc                                                              >  2:2992005/1‑93 (MQ=255)
                            ggTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTTGTAAACGCTTTTAGCCTTCCAGTTGGGGCGCACCTGCg                                                             >  2:771734/1‑93 (MQ=255)
                            ggTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCg                                                             <  2:488185/93‑1 (MQ=255)
                             gTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCCGGAAATCTTTGGTAAACGCTTTGAGCCTTCCCGTTGGGCCGCACCTGCgg                                                            >  1:1938957/1‑93 (MQ=255)
                              tGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGt                                                           >  1:3315792/1‑93 (MQ=255)
                               gAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGGCGCACCTGCgggt                                                          >  1:4410036/1‑91 (MQ=255)
                                aaGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGGAATCTTTTGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCCGCGGGTg                                                         >  2:1217298/1‑93 (MQ=255)
                                  gTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCCGGGAATCTTTTGTAAACGCTTTTAGCCTTCCAGTTGGGGCGCACCTGCGGGTggg                                                       >  2:4713641/1‑93 (MQ=255)
                                   ttCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGGAATCTTTTGTAAAAGCTTTTAGCCTTCCCGTTGGGGCGCCCCCGCGGGTGGGc                                                      >  2:3796292/1‑93 (MQ=255)
                                   ttCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAAACTTTTGTAAACGCTTTTTGCCTTCCCGTTTGGCCGCACCTGCGGGTGGGc                                                      >  2:3558089/1‑93 (MQ=255)
                                    tCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTTGTAAACGCTTTGAGCCTTCCAGTTTGGCCGCACCTGCGGGTGGGcc                                                     >  2:3015185/1‑93 (MQ=255)
                                    tCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAAAGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGGTGGGcc                                                     >  1:3853607/1‑93 (MQ=255)
                                     cTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGGAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGGTGGGCCg                                                    >  1:555098/1‑93 (MQ=255)
                                         tATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGa                                                <  2:1485873/93‑1 (MQ=255)
                                              cgccgTAGGGATCTTATGGCGGGCCAGGAAATTTTTGGTAAACGCTTTTAGCCTTCCAGTTGGGGCGCACCTGCGGTTGGGCCGAAGAttttg                                          >  1:4164376/1‑92 (MQ=255)
                                               gccgTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCa                                         <  2:566001/93‑1 (MQ=255)
                                               gccgTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCa                                         <  2:4387895/93‑1 (MQ=255)
                                                ccgTAGGGATCTTATGGCGGGCCAGGAAATCTTTTGTAAACGCTTTGAGCCTTCCCGTTGGGGCGCACCTGCGGTTGGGCCGGAGATTTTCAg                                        >  2:2199057/1‑93 (MQ=255)
                                                  gTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGcc                                      <  1:2617035/93‑1 (MQ=255)
                                                  ctAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGcc                                      <  1:2943301/92‑1 (MQ=255)
                                                   tAGGGATCTTATGGCGGGCCCGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGGCGCACCTGCGGTTGGGGCGGAGATTTTTCAGcc                                      >  1:3776222/1‑93 (MQ=255)
                                                   tAGGGATCTTATGGCGGGCCAGGGAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGGCGAAGATTTTCAGccc                                     >  1:3263513/1‑93 (MQ=255)
                                                   tAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGccc                                     <  2:1964306/93‑1 (MQ=255)
                                                   tAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAAAGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGccc                                     >  1:963474/1‑93 (MQ=255)
                                                    aGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCcg                                    <  2:4023843/93‑1 (MQ=255)
                                                    aGGGATCTTATGGCGGGCCAGGAAATCTTTGGGAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTTGGCCGAAGATTTTCAGCCcg                                    >  2:587438/1‑93 (MQ=255)
                                                     gggATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCcgg                                   <  1:2159123/93‑1 (MQ=255)
                                                      ggATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCcggc                                  <  2:3947534/93‑1 (MQ=255)
                                                       gATCGTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCcggcg                                 <  1:792572/93‑1 (MQ=255)
                                                        aTCTTATGGCGGGCCAGGAAATCTTTTGTAAACGCTTTGAGCCTTCCAGTTGGGGCGCACCTGCGGGTGGGCCGAAGATTTTTAGCCcggcgg                                >  2:2991710/1‑93 (MQ=255)
                                                        aTCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCcggcgg                                <  2:2949552/93‑1 (MQ=255)
                                                         tCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCggcggc                               <  1:1211079/93‑1 (MQ=255)
                                                         tCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCggcggc                               <  2:421662/93‑1 (MQ=255)
                                                         tCTTATGGCGGGCCAGGAAATCTTTGGTAAAAGCTTTTAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCggcggc                               >  1:4416865/1‑93 (MQ=255)
                                                           ttATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGgcggcg                              <  2:4044362/93‑1 (MQ=255)
                                                           ttATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCCGCGgcgc                             >  2:399094/1‑93 (MQ=255)
                                                                 ggggCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAgg                       <  1:3716870/92‑1 (MQ=255)
                                                                 cGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAgg                       <  1:1010414/93‑1 (MQ=255)
                                                                  gggCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGt                      <  1:2489249/93‑1 (MQ=255)
                                                                  gggCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGt                      <  1:4091255/93‑1 (MQ=255)
                                                                    gCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTAt                    <  2:76824/93‑1 (MQ=255)
                                                                     ccAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATc                   <  2:798215/93‑1 (MQ=255)
                                                                      cAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCg                  >  2:3735711/1‑93 (MQ=255)
                                                                      cAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCg                  >  2:3067575/1‑93 (MQ=255)
                                                                      cAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCg                  >  2:3998893/1‑93 (MQ=255)
                                                                      cAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGGATCg                  >  2:3462530/1‑93 (MQ=255)
                                                                       aGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGa                 >  1:1350598/1‑93 (MQ=255)
                                                                       aGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGa                 >  1:797910/1‑93 (MQ=255)
                                                                        ggAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGAc                <  1:2981498/93‑1 (MQ=255)
                                                                        ggAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGAc                <  1:2772362/93‑1 (MQ=255)
                                                                        ggAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGAc                >  2:1124473/1‑93 (MQ=255)
                                                                        ggAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGAc                >  1:2376267/1‑93 (MQ=255)
                                                                         gAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGAcc               >  2:2810783/1‑93 (MQ=255)
                                                                          aaaTCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCTGCGGCGCCGAAGGTATCGACCa              >  2:4411245/1‑93 (MQ=255)
                                                                          aaaTCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCa              >  1:1219897/1‑93 (MQ=255)
                                                                          aaaTCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCa              >  2:3510126/1‑93 (MQ=255)
                                                                          aaaTCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCa              <  1:473406/93‑1 (MQ=255)
                                                                           aaTCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCAc             <  1:2455636/93‑1 (MQ=255)
                                                                            aTCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACg            >  1:3264799/1‑93 (MQ=255)
                                                                            aTCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACg            <  1:1825297/93‑1 (MQ=255)
                                                                             tCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGc           >  1:3850997/1‑93 (MQ=255)
                                                                              cTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGcc          >  1:2924713/1‑93 (MQ=255)
                                                                               gtttGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGTTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGcc          <  2:4580016/92‑1 (MQ=255)
                                                                                tttGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCt         >  1:1906376/1‑93 (MQ=255)
                                                                                tttGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCt         <  1:1642267/93‑1 (MQ=255)
                                                                                tttGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGACt         >  1:1359306/1‑93 (MQ=255)
                                                                                  tGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCttt       <  1:4297933/93‑1 (MQ=255)
                                                                                  tGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGCCGGCGCGGACGGTATCGACCACGCAttt       >  1:2284816/1‑93 (MQ=255)
                                                                                   ggTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCTTTc      <  2:1094343/93‑1 (MQ=255)
                                                                                   ggTAAACGCTTTGAGCCTTCCAGGTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCTTTc      <  1:4079797/93‑1 (MQ=255)
                                                                                     tAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCTTTCAc    <  2:1865011/93‑1 (MQ=255)
                                                                                      aaaCGCTTTGATCCTGCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCTTTCAcc   <  1:2241563/93‑1 (MQ=255)
                                                                                      aaaCGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCTTTCAcc   <  1:2241564/93‑1 (MQ=255)
                                                                                       aaCGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCTTTCACCa  <  2:1830996/93‑1 (MQ=255)
                                                                                       aaCGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGCGATCCACGACGCCTTTTAccg  >  2:3087934/1‑92 (MQ=255)
                                                                                            |                                                                                         
GCAGATACGCCAGCGCAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCTTTCACCA  >  NC_000913/4203104‑4203283

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 26 ≤ 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
AGATACGCCAGCGCAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCTTTCACCAGCGGCGCT  >  NC_000913/4203106‑4203291
                                                                                            |                                                                                              
AGATACGCCAGCGCAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAGAAATTTGGTAAAAGCTTTGG                                                                                               >  SRR4302544.3780667/1‑94 (MQ=60)
    ACGCCAGCGCAGGTTCTACCTCGGTGAAGGTCTGGTATGCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCT                                                                                          <  SRR4302544.3405270/94‑1 (MQ=60)
        CAGCGCAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCCGGGAATCTTTGGGAAACGCTTTTAGCCTTCCC                                                                                      >  SRR4302544.852380/1‑94 (MQ=60)
        CAGCGCAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAAACTTTGGTAAAACCTTTTAGCCTTCCA                                                                                      >  SRR4302544.1024897/1‑94 (MQ=60)
         AGCGCAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGGCCGGAAATCTTTTGTAAACGCTTTTAGCCTTCCAG                                                                                     >  SRR4302544.75257/1‑94 (MQ=60)
         AGCGCAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAG                                                                                     <  SRR4302544.1315390/94‑1 (MQ=60)
                GTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCC                                                                              <  SRR4302544.382933/94‑1 (MQ=60)
                      CCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTTGTAAACGCTTTTAGCCCTCCCGTTGGGCCGCACCC                                                                        >  SRR4302544.2481006/1‑94 (MQ=60)
                      CCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCT                                                                        <  SRR4302544.3754804/94‑1 (MQ=60)
                      CCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCT                                                                        <  SRR4302544.3754805/94‑1 (MQ=60)
                       CTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTTGAGCCTTCCAGTTGGGCCGCACCT                                                                        <  SRR4302544.4065038/94‑1 (MQ=60)
                         CGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCG                                                                     <  SRR4302544.3268604/94‑1 (MQ=60)
                          GGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCCGGAAATCTTTGGTAAACGCTTTGAGCCTTCCCGTTGGGCCGCACCTGCGG                                                                    >  SRR4302544.1957730/1‑94 (MQ=60)
                           GTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGT                                                                   >  SRR4302544.3347372/1‑94 (MQ=60)
                            TGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGGCGCACCTGCGGGT                                                                  >  SRR4302544.4454565/1‑94 (MQ=60)
                                 TTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAAAGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGGTGGGCC                                                             >  SRR4302544.3891857/1‑94 (MQ=60)
                                  TCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGGAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGGTGGGCCG                                                            >  SRR4302544.560530/1‑94 (MQ=60)
                                           CCGCCGTAGGGATCTTATGGCGGGCCAGGAAATTTTTGGTAAACGCTTTTAGCCTTCCAGTTGGGGCGCACCTGCGGTTGGGCCGAAGATTTTG                                                  >  SRR4302544.4205628/1‑94 (MQ=60)
                                                GTAGGGATCTTATGGCGGGCCCGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGGCGCACCTGCGGTTGGGGCGGAGATTTTTCAGCC                                              >  SRR4302544.3813645/1‑94 (MQ=60)
                                                GTAGGGATCTTATGGCGGGCCAGGGAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGGCGAAGATTTTCAGCCC                                             >  SRR4302544.3294679/1‑94 (MQ=60)
                                                GTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCC                                             <  SRR4302544.2642150/94‑1 (MQ=60)
                                                GTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAAAGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCC                                             >  SRR4302544.972935/1‑94 (MQ=60)
                                                CTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCC                                             <  SRR4302544.2971574/94‑1 (MQ=60)
                                                  AGGGATCTTATGGCGGGCCAGGAAATCTTTTGTAAACGCTTTGAGCCCTTCAGTTGGGGCGCACCCGCGGTTGGGGCGAAGATTTTCAGCCCGG                                           >  SRR4302544.4187310/1‑94 (MQ=60)
                                                  AGGGATCTTATGGCGGGCCAGGAAAACTTTAGTAAACGCTTTGAGCCTTCCAGTTGGGCCGGACCTGCGGTTGGGCCGGAGATTTTCAACCCGG                                           >  SRR4302544.2958663/1‑94 (MQ=60)
                                                   GGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGC                                          <  SRR4302544.2179637/94‑1 (MQ=60)
                                                     GATCGTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGG                                        <  SRR4302544.800244/94‑1 (MQ=60)
                                                      ATCTTATGGCGGGCCAGGAAATCTTTGGTAAAAGCTTTTAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGC                                       >  SRR4302544.4461448/1‑94 (MQ=60)
                                                       TCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCG                                      <  SRR4302544.1223071/94‑1 (MQ=60)
                                                               GGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGT                              <  SRR4302544.3753845/94‑1 (MQ=60)
                                                               CGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGT                              <  SRR4302544.1020480/94‑1 (MQ=60)
                                                                GGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTA                             <  SRR4302544.4131946/94‑1 (MQ=60)
                                                                GGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTA                             <  SRR4302544.2513092/94‑1 (MQ=60)
                                                                    CAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGA                         >  SRR4302544.805650/1‑94 (MQ=60)
                                                                    CAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGA                         >  SRR4302544.1363983/1‑94 (MQ=60)
                                                                     AGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGAC                        >  SRR4302544.2398766/1‑94 (MQ=60)
                                                                      GGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACC                       <  SRR4302544.3010092/94‑1 (MQ=60)
                                                                      GGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACC                       <  SRR4302544.2799031/94‑1 (MQ=60)
                                                                       GAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCA                      >  SRR4302544.1231964/1‑94 (MQ=60)
                                                                        AAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCAC                     <  SRR4302544.478033/94‑1 (MQ=60)
                                                                         AATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACG                    <  SRR4302544.2479063/94‑1 (MQ=60)
                                                                         AATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACG                    >  SRR4302544.3295972/1‑94 (MQ=60)
                                                                         AATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTGGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACG                    <  SRR4302544.2709366/94‑1 (MQ=60)
                                                                          ATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGC                   <  SRR4302544.1843075/94‑1 (MQ=60)
                                                                          ATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGC                   >  SRR4302544.3889223/1‑94 (MQ=60)
                                                                           TCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCC                  >  SRR4302544.2952826/1‑94 (MQ=60)
                                                                            CTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCT                 >  SRR4302544.1924806/1‑94 (MQ=60)
                                                                            CTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGACT                 >  SRR4302544.1372799/1‑94 (MQ=60)
                                                                             TTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCTT                <  SRR4302544.1658452/94‑1 (MQ=60)
                                                                              TTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGCCGGCGCGGACGGTATCGACCACGCATTT               >  SRR4302544.2306351/1‑94 (MQ=60)
                                                                               TGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCTTTC              <  SRR4302544.4341795/94‑1 (MQ=60)
                                                                                GGTAAACGCTTTGAGCCTTCCAGGTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCTTTCA             <  SRR4302544.4120386/94‑1 (MQ=60)
                                                                                   AAACGCTTTGATCCTGCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCTTTCACCA          <  SRR4302544.2262651/94‑1 (MQ=60)
                                                                                   AAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCTTTCACCA          <  SRR4302544.2262652/94‑1 (MQ=60)
                                                                                     aactcTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTGTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCcgccaccag         >  SRR4302544.4242432/6‑85 (MQ=60)
                                                                                      CGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCTTTCACCAGCG       >  SRR4302544.865149/1‑94 (MQ=60)
                                                                                      CGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCTTTCACCAGCG       <  SRR4302544.2848370/94‑1 (MQ=60)
                                                                                            TGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGGATCGACCACGCCTTTCACCAGCGGCGCT  >  SRR4302544.1988638/1‑94 (MQ=60)
                                                                                            |                                                                                              
AGATACGCCAGCGCAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCTTTCACCAGCGGCGCT  >  NC_000913/4203106‑4203291

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