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 I2 R1 1873 183.1 9376148 98.2% 9207377 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,011,430 A→G H263R (CAT→CGT)  metE → tetrahydropteroyltriglutamate methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,011,4300AG99.5% 654.3 / ‑5.1 209H263R (CAT→CGT) metEtetrahydropteroyltriglutamate methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (97/110);  total (98/111)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.60e-01

Only 100 of 216 total aligned reads displayed.

AGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGATCAATG  >  NC_000913/4011340‑4011519
                                                                                          |                                                                                         
aGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTg                                                                                         >  2:3089579/1‑93 (MQ=255)
 gggACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCgtgt                                                                                        >  1:3018301/1‑93 (MQ=255)
 gggACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCgtgt                                                                                        >  1:3163173/1‑93 (MQ=255)
 gggACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCgtgt                                                                                        >  2:3369584/1‑93 (MQ=255)
  ggACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGtt                                                                                       <  2:2953644/93‑1 (MQ=255)
    aCAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGCTCAGGGTCTGCGTGTTGa                                                                                     >  2:85969/1‑93 (MQ=255)
     cAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGAc                                                                                    >  1:3813105/1‑93 (MQ=255)
      aGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGAcc                                                                                   <  2:37474/93‑1 (MQ=255)
      aGGCGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGAcc                                                                                   >  1:4116648/1‑93 (MQ=255)
       ggTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCt                                                                                  <  1:4275330/93‑1 (MQ=255)
       ggTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCt                                                                                  >  1:3597120/1‑93 (MQ=255)
       ggTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCt                                                                                  >  2:58343/1‑93 (MQ=255)
         tGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCg                                                                                >  2:3438468/1‑93 (MQ=255)
          gAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGt                                                                               >  2:3271325/1‑93 (MQ=255)
          gAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTCCCTGTTCAGGGTCTGCGTGTTGACCTCGt                                                                               >  2:4602285/1‑93 (MQ=255)
           aaaCTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTa                                                                              >  1:2975352/1‑93 (MQ=255)
            aaCTGCTGCTGACCGCCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTAc                                                                             <  2:2787686/93‑1 (MQ=255)
            aaCTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTAc                                                                             <  1:939817/93‑1 (MQ=255)
            aaCTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTAc                                                                             >  1:1576181/1‑93 (MQ=255)
              ctgctgctgACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACAt                                                                           <  2:1315145/93‑1 (MQ=255)
              ctgctgctgACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACAt                                                                           <  1:3498714/93‑1 (MQ=255)
                gctgctgACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATgg                                                                         >  2:2259085/1‑93 (MQ=255)
                 ctgctgACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGTTCTGCGTGTTGACCTCGTACATGGt                                                                        <  1:4143703/93‑1 (MQ=255)
                   gctgACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTaa                                                                      >  2:1844678/1‑93 (MQ=255)
                    ctgACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTaaa                                                                     <  1:4355098/93‑1 (MQ=255)
                      gACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGa                                                                   <  2:372541/93‑1 (MQ=255)
                      gACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGa                                                                   <  1:2354056/93‑1 (MQ=255)
                      gACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGa                                                                   <  2:4331143/93‑1 (MQ=255)
                      gACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGa                                                                   <  2:4421812/93‑1 (MQ=255)
                       gccaccTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGAt                                                                  <  1:3242807/92‑1 (MQ=255)
                       accaccTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGAt                                                                  <  1:3833178/93‑1 (MQ=255)
                        ccaccTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATg                                                                 <  1:3440218/93‑1 (MQ=255)
                         caccTATTTTGAAGGCGTAACGCCAAGTCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGa                                                                <  1:2148120/93‑1 (MQ=255)
                         caccTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGa                                                                <  2:3298661/93‑1 (MQ=255)
                         caccTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGa                                                                <  1:1496651/93‑1 (MQ=255)
                          accTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGAc                                                               <  1:4637356/93‑1 (MQ=255)
                           ccTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACg                                                              >  2:3396529/1‑93 (MQ=255)
                            cTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGActt                                                             >  1:3888413/1‑91 (MQ=255)
                            cTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGt                                                             >  1:70091/1‑93 (MQ=255)
                            cTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGt                                                             >  2:352402/1‑93 (MQ=255)
                            cTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGt                                                             >  1:1388327/1‑93 (MQ=255)
                             tATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGtt                                                            >  2:2989078/1‑93 (MQ=255)
                              aTTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTg                                                           <  2:2052001/93‑1 (MQ=255)
                               ttttGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTTACCTCGTACATGGTAAAGATGACGTTGc                                                          <  2:1727263/93‑1 (MQ=255)
                               ttttGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGc                                                          >  2:3400102/1‑93 (MQ=255)
                               ttttGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGc                                                          <  2:506805/93‑1 (MQ=255)
                               ttttGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGc                                                          >  1:785923/1‑93 (MQ=255)
                                 ttGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTg                                                        <  1:2441697/93‑1 (MQ=255)
                                   gAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGaa                                                      <  2:1793258/93‑1 (MQ=255)
                                    aaGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAAc                                                     <  2:4324630/93‑1 (MQ=255)
                                     aGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACt                                                    <  1:1187589/93‑1 (MQ=255)
                                     aGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACt                                                    <  1:2067953/93‑1 (MQ=255)
                                       gCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGc                                                  >  2:101478/1‑93 (MQ=255)
                                         gTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGcac                                                >  2:1885887/1‑93 (MQ=255)
                                         gTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGcac                                                >  1:2827478/1‑93 (MQ=255)
                                           accGCCAAATCACGACGCGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCGCGTACATGGTAAAGATGACGTTGCTGAACTGCACaa                                              <  2:1089922/91‑1 (MQ=255)
                                            accacAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAg                                             <  2:1420999/89‑1 (MQ=255)
                                            aCGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAg                                             <  2:4575769/93‑1 (MQ=255)
                                              gCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAgcg                                           >  2:2691101/1‑93 (MQ=255)
                                               ccAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAgcgc                                          >  1:4093970/1‑93 (MQ=255)
                                                cAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCgcc                                         <  1:2397468/93‑1 (MQ=255)
                                                 aaaTCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCgcct                                        <  1:1737659/93‑1 (MQ=255)
                                                  aaTCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCgcctg                                       <  2:2449984/93‑1 (MQ=255)
                                                  aaTCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCgcctg                                       >  2:2859650/1‑93 (MQ=255)
                                                   aTCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCgcctgc                                      >  2:3230782/1‑93 (MQ=255)
                                                    tctcGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCgcctgcc                                     <  1:3402338/93‑1 (MQ=255)
                                                     ctcGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCgcctgcct                                    <  1:1714447/93‑1 (MQ=255)
                                                     cccGACACGATTACTGCGCGGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCgcctgcct                                    <  1:3343944/91‑1 (MQ=255)
                                                       cGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTc                                  <  2:4431199/93‑1 (MQ=255)
                                                       cGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTc                                  >  2:4134934/1‑93 (MQ=255)
                                                       cGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTc                                  >  1:2561619/1‑93 (MQ=255)
                                                       cGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTc                                  >  1:3917855/1‑93 (MQ=255)
                                                       cGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTc                                  <  2:4176807/93‑1 (MQ=255)
                                                       cGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTc                                  >  1:1281420/1‑93 (MQ=255)
                                                        gACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCt                                 >  2:1524440/1‑93 (MQ=255)
                                                        gACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCt                                 <  1:4573506/93‑1 (MQ=255)
                                                           acGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGAc                              >  2:2431817/1‑93 (MQ=255)
                                                            cGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACt                             <  1:262961/93‑1 (MQ=255)
                                                               ttACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACGTGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGt                          >  1:2521943/1‑93 (MQ=255)
                                                               ttACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGt                          <  1:1458531/93‑1 (MQ=255)
                                                                 aCTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTg                        >  2:4669671/1‑93 (MQ=255)
                                                                  cTGCGCTGCCTGTTCAGGGTCTGCGTGTTGCCCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGc                       <  1:1679227/93‑1 (MQ=255)
                                                                  cTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGc                       <  1:1824012/93‑1 (MQ=255)
                                                                  cTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGc                       <  2:4234763/93‑1 (MQ=255)
                                                                   tGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCt                      <  2:2353918/93‑1 (MQ=255)
                                                                    gcgcTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTg                     <  2:3186780/93‑1 (MQ=255)
                                                                     cgcTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGt                    <  2:1222610/93‑1 (MQ=255)
                                                                      gcTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTc                   >  1:4149053/1‑93 (MQ=255)
                                                                       cTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCt                  <  1:1042039/93‑1 (MQ=255)
                                                                        tGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTg                 >  1:551651/1‑93 (MQ=255)
                                                                        tGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTg                 >  1:3728312/1‑93 (MQ=255)
                                                                          ccTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCg               >  2:3908847/1‑93 (MQ=255)
                                                                            tGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCggg             >  2:4406525/1‑93 (MQ=255)
                                                                             gTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGt            <  1:1661842/93‑1 (MQ=255)
                                                                              ttCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTc           <  1:2707961/93‑1 (MQ=255)
                                                                               tCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCt          <  1:602639/93‑1 (MQ=255)
                                                                                  gggTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGAt       >  1:2426185/1‑93 (MQ=255)
                                                                                   ggTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGATc      >  1:1167474/1‑93 (MQ=255)
                                                                                    gTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGATCa     >  2:4214256/1‑93 (MQ=255)
                                                                                       tGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGATCAATg  <  2:2015337/93‑1 (MQ=255)
                                                                                          |                                                                                         
AGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGATCAATG  >  NC_000913/4011340‑4011519

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 30 ≤ 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
TCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGATCAATGGTCG  >  NC_000913/4011337‑4011523
                                                                                             |                                                                                             
TCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCG                                                                                               >  SRR4302535.2023820/1‑94 (MQ=60)
TCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCG                                                                                               <  SRR4302535.2340730/94‑1 (MQ=60)
  CAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTG                                                                                             <  SRR4302535.4115391/94‑1 (MQ=60)
  CAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTG                                                                                             <  SRR4302535.3853469/94‑1 (MQ=60)
   AGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGT                                                                                            >  SRR4302535.3039444/1‑94 (MQ=60)
   AGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGT                                                                                            >  SRR4302535.3185110/1‑94 (MQ=60)
    aggacaggttaaactgcgtgtgacccccttttttggaagcgttaccgCAAATCTCGACCCGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTT                                                                                           <  SRR4302535.350169/47‑1 (MQ=60)
     GGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTG                                                                                          <  SRR4302535.1746562/94‑1 (MQ=60)
       ACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGAC                                                                                        >  SRR4302535.3839419/1‑94 (MQ=60)
        CAGGCGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACC                                                                                       >  SRR4302535.4144824/1‑94 (MQ=60)
         AGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCT                                                                                      >  SRR4302535.3622081/1‑94 (MQ=60)
         AGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCT                                                                                      <  SRR4302535.1525845/94‑1 (MQ=60)
          GGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTC                                                                                     <  SRR4302535.4304595/94‑1 (MQ=60)
           GTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTCACGCCAAATCTCGACACGATTACTGAGCTGCCTGTTCCGTGTCTGCGTGTTGACCTCG                                                                                    >  SRR4302535.1292517/1‑94 (MQ=60)
             GAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTA                                                                                  >  SRR4302535.2996232/1‑94 (MQ=60)
              AAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTAC                                                                                 >  SRR4302535.1587551/1‑94 (MQ=60)
               AACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACA                                                                                <  SRR4302535.946473/94‑1 (MQ=60)
                 CTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATG                                                                              <  SRR4302535.3523031/94‑1 (MQ=60)
                 CTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATG                                                                              >  SRR4302535.1474307/1‑94 (MQ=60)
                 CTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATG                                                                              <  SRR4302535.336794/94‑1 (MQ=60)
                  TGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTAACCTCGTACATGG                                                                             >  SRR4302535.96932/1‑94 (MQ=60)
                    CTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGTTCTGCGTGTTGACCTCGTACATGGTA                                                                           <  SRR4302535.4172076/94‑1 (MQ=60)
                       CTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAG                                                                        >  SRR4302535.1465717/1‑94 (MQ=60)
                       CTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAG                                                                        <  SRR4302535.4384909/94‑1 (MQ=60)
                       CTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAG                                                                        <  SRR4302535.2214130/94‑1 (MQ=60)
                         GACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGAT                                                                      <  SRR4302535.944109/94‑1 (MQ=60)
                         GACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGAT                                                                      <  SRR4302535.3909205/94‑1 (MQ=60)
                         GACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGAT                                                                      <  SRR4302535.2370934/94‑1 (MQ=60)
                          GCCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATG                                                                     <  SRR4302535.3265332/94‑1 (MQ=60)
                          ACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATG                                                                     <  SRR4302535.3839450/94‑1 (MQ=60)
                          ACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATG                                                                     <  SRR4302535.3859613/94‑1 (MQ=60)
                          ACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATG                                                                     <  SRR4302535.4358653/94‑1 (MQ=60)
                           CCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGA                                                                    <  SRR4302535.3464142/94‑1 (MQ=60)
                           CCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGA                                                                    <  SRR4302535.4562810/94‑1 (MQ=60)
                            CACCTATTTTGAAGGCGTAACGCCAAGTCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGAC                                                                   <  SRR4302535.2163568/94‑1 (MQ=60)
                            CACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGAC                                                                   <  SRR4302535.1172147/94‑1 (MQ=60)
                            CACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGAC                                                                   <  SRR4302535.2509063/94‑1 (MQ=60)
                            CACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGAC                                                                   <  SRR4302535.1507414/94‑1 (MQ=60)
                             ACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACG                                                                  >  SRR4302535.2112029/1‑94 (MQ=60)
                             ACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACG                                                                  <  SRR4302535.4668810/94‑1 (MQ=60)
                             ACCTATGTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACG                                                                  <  SRR4302535.2819371/94‑1 (MQ=60)
                              CCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACTT                                                                 >  SRR4302535.3915081/1‑94 (MQ=60)
                              CCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGT                                                                 >  SRR4302535.1398237/1‑94 (MQ=60)
                              CCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGT                                                                 >  SRR4302535.70607/1‑94 (MQ=60)
                               CTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTT                                                                <  SRR4302535.3720708/94‑1 (MQ=60)
                                TATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTG                                                               <  SRR4302535.1678761/94‑1 (MQ=60)
                                 ATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGC                                                              >  SRR4302535.791432/1‑94 (MQ=60)
                                 ATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGC                                                              <  SRR4302535.3952468/94‑1 (MQ=60)
                                    TTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGA                                                           <  SRR4302535.2459287/94‑1 (MQ=60)
                                     TGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAA                                                          >  SRR4302535.2432233/1‑94 (MQ=60)
                                       AAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACT                                                        <  SRR4302535.689898/94‑1 (MQ=60)
                                       AAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACT                                                        <  SRR4302535.466815/94‑1 (MQ=60)
                                       AAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACT                                                        <  SRR4302535.4132953/94‑1 (MQ=60)
                                        AGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTG                                                       <  SRR4302535.2082858/94‑1 (MQ=60)
                                        AGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTG                                                       <  SRR4302535.1196082/94‑1 (MQ=60)
                                         GGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGC                                                      >  SRR4302535.348400/1‑94 (MQ=60)
                                          GCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGAGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCA                                                     <  SRR4302535.2524216/94‑1 (MQ=60)
                                           CGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCAC                                                    >  SRR4302535.2847543/1‑94 (MQ=60)
                                           CGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCAC                                                    <  SRR4302535.2007066/94‑1 (MQ=60)
                                            GTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACA                                                   <  SRR4302535.723888/94‑1 (MQ=60)
                                            GTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACA                                                   <  SRR4302535.1700529/94‑1 (MQ=60)
                                             TAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAA                                                  <  SRR4302535.4506062/94‑1 (MQ=60)
                                              AACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAG                                                 >  SRR4302535.3038759/1‑94 (MQ=60)
                                              AACGCCAAATCTCGACACGATTACCGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAG                                                 >  SRR4302535.831088/1‑94 (MQ=60)
                                               ACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGC                                                >  SRR4302535.1532100/1‑94 (MQ=60)
                                                CGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCG                                               <  SRR4302535.1831164/94‑1 (MQ=60)
                                                 GCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGC                                              >  SRR4302535.4121986/1‑94 (MQ=60)
                                                  CCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCC                                             >  SRR4302535.3453873/1‑94 (MQ=60)
                                                   CAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCT                                            <  SRR4302535.2414701/94‑1 (MQ=60)
                                                   CAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCT                                            >  SRR4302535.2610335/1‑94 (MQ=60)
                                                    AAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTG                                           <  SRR4302535.1750218/94‑1 (MQ=60)
                                                    AAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTG                                           >  SRR4302535.3181738/1‑94 (MQ=60)
                                                    AAATATCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTG                                           <  SRR4302535.2051977/94‑1 (MQ=60)
                                                       TCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCT                                        <  SRR4302535.3426015/94‑1 (MQ=60)
                                                        CTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTT                                       <  SRR4302535.1194351/94‑1 (MQ=60)
                                                        CTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTT                                       <  SRR4302535.1726845/94‑1 (MQ=60)
                                                        CCCGACACGATTACTGCGCGGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTT                                       <  SRR4302535.3367240/94‑1 (MQ=60)
                                                         TCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTC                                      >  SRR4302535.2580061/1‑94 (MQ=60)
                                                         TCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTC                                      >  SRR4302535.1290566/1‑94 (MQ=60)
                                                         TCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTC                                      <  SRR4302535.2453285/94‑1 (MQ=60)
                                                         TCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTC                                      >  SRR4302535.3944659/1‑94 (MQ=60)
                                                           GACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTG                                    <  SRR4302535.4604681/94‑1 (MQ=60)
                                                            ACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGA                                   >  SRR4302535.1602859/1‑94 (MQ=60)
                                                            ACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGA                                   >  SRR4302535.4687759/1‑94 (MQ=60)
                                                               CGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTG                                <  SRR4302535.264785/94‑1 (MQ=60)
                                                               CGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTG                                >  SRR4302535.3783354/1‑94 (MQ=60)
                                                                GATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGG                               >  SRR4302535.4164646/1‑94 (MQ=60)
                                                                GATTACTGCGCTGCCTGTTCAGGGCCGTCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGG                               <  SRR4302535.2681524/94‑1 (MQ=60)
                                                                 ATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACGTGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGT                              >  SRR4302535.2540116/1‑94 (MQ=60)
                                                                 ATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGT                              >  SRR4302535.3971990/1‑94 (MQ=60)
                                                                  TTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTT                             >  SRR4302535.567480/1‑94 (MQ=60)
                                                                  TTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTT                             <  SRR4302535.1468973/94‑1 (MQ=60)
                                                                     CTGCGCTGCCTGTTCAGGGTCTGCGTGTTGCCCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCT                          <  SRR4302535.1691388/94‑1 (MQ=60)
                                                                     CTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCT                          <  SRR4302535.3014696/94‑1 (MQ=60)
                                                                     CTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCT                          <  SRR4302535.1837182/94‑1 (MQ=60)
                                                                      TGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTG                         <  SRR4302535.4173829/94‑1 (MQ=60)
                                                                        CGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTC                       >  SRR4302535.4177461/1‑94 (MQ=60)
                                                                         GCTGCCTGTTCAGGGTCTGCTTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCT                      <  SRR4302535.322956/94‑1 (MQ=60)
                                                                          CTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTG                     >  SRR4302535.555514/1‑94 (MQ=60)
                                                                          CTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTG                     >  SRR4302535.3754136/1‑94 (MQ=60)
                                                                          CTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTG                     <  SRR4302535.1049426/94‑1 (MQ=60)
                                                                           TGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGC                    <  SRR4302535.3537486/94‑1 (MQ=60)
                                                                              CTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGG                 <  SRR4302535.2575189/94‑1 (MQ=60)
                                                                              CTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGGAAAGATGACGTTGCCGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGG                 >  SRR4302535.1216198/1‑94 (MQ=60)
                                                                                GTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTC               <  SRR4302535.1673850/94‑1 (MQ=60)
                                                                                 TTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCT              <  SRR4302535.2727311/94‑1 (MQ=60)
                                                                                  TCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTG             >  SRR4302535.3268638/1‑94 (MQ=60)
                                                                                  TCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTG             <  SRR4302535.606856/94‑1 (MQ=60)
                                                                                   CAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGA            >  SRR4302535.418585/1‑94 (MQ=60)
                                                                                    AGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGAT           >  SRR4302535.2443664/1‑94 (MQ=60)
                                                                                     GGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGATC          >  SRR4302535.1175808/1‑94 (MQ=60)
                                                                                       GTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGATCAA        >  SRR4302535.2411006/1‑94 (MQ=60)
                                                                                        TCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGATCAAT       <  SRR4302535.3929966/94‑1 (MQ=60)
                                                                                          TGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGATCAATGG     <  SRR4302535.2925758/94‑1 (MQ=60)
                                                                                          TGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGATCAATGG     >  SRR4302535.1847967/1‑94 (MQ=60)
                                                                                          TGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGATCAATGG     <  SRR4302535.354834/94‑1 (MQ=60)
                                                                                           GCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGATCAATGGT    <  SRR4302535.1712593/94‑1 (MQ=60)
                                                                                            CGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGATCAATGGTC   <  SRR4302535.3022804/94‑1 (MQ=60)
                                                                                            CGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGATCAATGGTC   <  SRR4302535.4619454/94‑1 (MQ=60)
                                                                                            CGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGATCAATGGTC   <  SRR4302535.873180/94‑1 (MQ=60)
                                                                                             GTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGATCAATGGTCG  >  SRR4302535.257165/1‑94 (MQ=60)
                                                                                             GTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGATCAATGGTCG  >  SRR4302535.2593161/1‑94 (MQ=60)
                                                                                             |                                                                                             
TCCAGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGATCAATGGTCG  >  NC_000913/4011337‑4011523

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