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,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.95.3% 449.9 / 12.6 107coding (315/1290 nt)purDphosphoribosylglycinamide synthetase
Reads supporting (aligned to +/- strand):  ref base T (5/0);  new base . (46/56);  total (51/56)
Fisher's exact test for biased strand distribution p-value = 2.21e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 163 total aligned reads displayed.

CGCAGATACGCCAGCGCAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCTTTCACCA  >  NC_000913/4203103‑4203283
                                                                                             |                                                                                        
cGCAGATACGCCAGCGCAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCt                                                                                          <  1:3566920/93‑1 (MQ=255)
    gATACGCCAGCGCAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGGAAACTTTTGTAAAAGCtttga                                                                                      >  2:1315299/1‑91 (MQ=255)
    gATACGCCAGCGCAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCtttga                                                                                      <  1:2921125/93‑3 (MQ=255)
    gATACGCCAGCGCAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCtttga                                                                                      <  1:4270271/93‑3 (MQ=255)
        cGCCAGCGCAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCt                                                                                 <  2:3336224/93‑1 (MQ=255)
           cAGCGCAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGGCCGGGAAACTTTTGTAAAAGCTTTGAGACTTcc                                                                              >  2:3475438/1‑93 (MQ=255)
            aGCGCAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTTGTAAAAGCTTTGTGCCTTccc                                                                             >  2:1889962/1‑92 (MQ=255)
            aGCGCAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCa                                                                             <  1:2764817/93‑1 (MQ=255)
             gcgcAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGGAATCTTTGGTAAACGCTTTGAGCCTTCCAg                                                                            >  2:2255016/1‑93 (MQ=255)
             gcgcAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAg                                                                            <  1:2228417/93‑1 (MQ=255)
             gcgcAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAAGGAATCTTTGGTAAACGCTTTGAGACTTcccg                                                                            >  2:2255015/1‑91 (MQ=255)
              cgcAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGGCAGGGAATATTTTGTTAACGCTTTGAGCCTTCCCGt                                                                           >  1:62651/1‑93 (MQ=255)
              cgcAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGGCAGGGAAACTTTTGTAAACGCTTTTAGGCTTCCCGt                                                                           >  1:2041475/1‑93 (MQ=255)
               gcAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGGCCGGGAAATTTTGGTAAAAGCTTTTGGGCTTCCCatt                                                                           >  1:3282856/1‑91 (MQ=255)
               gcAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGGAATCATTCGTAAACGCTTTGAGCCTTTCAGtt                                                                          >  2:1744483/1‑93 (MQ=255)
               gcAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGtt                                                                          <  2:1316998/93‑1 (MQ=255)
               gcAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTGCCAGtt                                                                          <  2:2552497/93‑1 (MQ=255)
                cAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTg                                                                         <  1:1696689/93‑1 (MQ=255)
                cAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTg                                                                         <  2:3801234/93‑1 (MQ=255)
                 aGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTgg                                                                        <  2:4338734/93‑1 (MQ=255)
                 aGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTgg                                                                        <  2:4469796/93‑1 (MQ=255)
                  ggTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGGCCGGGAATCTTTGGTAAACGCTTTGAGCCTTTCCGTTggg                                                                       >  2:1491673/1‑93 (MQ=255)
                  ggTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCCGGAAATCTTTGGTGAACGCTTTGAGACCTCCAGTTggg                                                                       >  2:3856801/1‑93 (MQ=255)
                  ggTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGGAATCTTTGGTAAAAGCTTTTAGCCTTCCCGTTggg                                                                       >  2:3424733/1‑93 (MQ=255)
                  ggTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTTGTAAACGCTTTGAGGCTTCCAGTTtgg                                                                       >  2:1885581/1‑93 (MQ=255)
                  ggTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAAAGCTTTTAGCCCTCCCGTTggg                                                                       >  2:3238697/1‑93 (MQ=255)
                   gTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGc                                                                      <  1:4641058/93‑1 (MQ=255)
                   gTTCTACCCCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGGGCCGTCCAGTTGGGc                                                                      <  2:1619981/93‑1 (MQ=255)
                    ttCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGcc                                                                     <  2:1425226/93‑1 (MQ=255)
                         ccTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGGCAGGGAAACTTTTGTAAACGCTTTGAGGCTTCCAGTTGGGGCGcccc                                                                >  2:3099674/1‑93 (MQ=255)
                         ccTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTTGTAAAAGCTTTTAGCCTTCCACTTGGGCCGCAcc                                                                >  1:1655517/1‑93 (MQ=255)
                         ccTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCAcc                                                                <  1:2361787/93‑1 (MQ=255)
                         ccTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCAcc                                                                <  2:806606/93‑1 (MQ=255)
                          cTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGGCAGGGAAACTTTTGTAAAAGCTTTTTGCCTTCCCGTTTGGGCGCAccc                                                               >  2:113496/1‑92 (MQ=255)
                           tCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTg                                                              <  2:1215811/93‑1 (MQ=255)
                           tCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAAAGCTTTGAGCCTTCCAGTTGGGGCGCACCTg                                                              >  1:4186709/1‑93 (MQ=255)
                            cGGTGAAGTTCTGGTATTCCGCCGTATGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGc                                                             <  1:753081/93‑1 (MQ=255)
                               tGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAAACTTTGGTAAAAGCTTTGAGCCTTCCCGTTGGGACGCACCTGCggg                                                          >  1:356665/1‑92 (MQ=255)
                                gAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGGCCGGGAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGtt                                                         >  1:2102194/1‑93 (MQ=255)
                                gAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGGCAGGGAATCTTTTGTAAAAGCTTTGAGGCTTCCCGTTTGGCCGCCCCCGCGGtt                                                         >  1:2424179/1‑93 (MQ=255)
                                gAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGGCAGGAAAACTTTGGGAAAAGCTTTTAGGCTTCCAGTTGGGCCGCCCCCGCgggt                                                         >  1:579467/1‑91 (MQ=255)
                                gAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGtt                                                         <  1:3192222/93‑1 (MQ=255)
                                   gTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGGAATCTTTTGTAAACGCTTTGAGCCTTCCAGTTGGGGCGCACCTGCGGGTggg                                                      >  2:3953737/1‑93 (MQ=255)
                                    ttCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTTGTAAACGCTTTGAGGCTTCCAGTTTGGGCGCACCTGCCGTTGGGc                                                     >  2:2050391/1‑93 (MQ=255)
                                    ttCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTTAGCCTTCCAGTTGGGCCGCACCCGCGGGTGGGc                                                     >  2:3496294/1‑93 (MQ=255)
                                      cTGTTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCg                                                   <  2:2880481/93‑1 (MQ=255)
                                        ggTATTCCGCCGTAGGGATCTTATGGCGGGCCCGGAAAACTTTGGTAAACGCTTTGAGCCTTCCCGTTGGGCCGCACCCGGGGTTGGGCCGaa                                                 >  2:4655256/1‑93 (MQ=255)
                                        ggTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGaa                                                 <  2:3972804/93‑1 (MQ=255)
                                         gTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAg                                                <  1:1128099/93‑1 (MQ=255)
                                          tATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTTGTAAACGCTTTTAGCCTTCCAGTTGGGGCGCACCTGCGGGTGGGGCGAAGa                                               >  2:68127/1‑93 (MQ=255)
                                          tATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGGCGCACCTGCGGTTGGGGCGgaga                                               >  1:1391146/1‑93 (MQ=255)
                                            ttCCGCCTTAGGGAGGTTAGGCCGGTCAAGGAATGCTTTGGTAAACGCTTTGAGCCTACAAGTTGGTCCGCACCTGCGGTTGGGCCGAAGAtt                                             <  1:2571128/93‑1 (MQ=255)
                                            ttCCGCCGTAGGGATCTTATGGCGGGCCCGGAAATCTTTGGTAAAAGCTTTGAGCCCTTCAGGTGGGGCGCACCTGCGGGTGGGGCGCAGAtt                                             >  2:1387299/1‑93 (MQ=255)
                                             tCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGAttt                                            <  1:2304265/93‑1 (MQ=255)
                                              ccgccgTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGAtttt                                           <  2:831273/93‑1 (MQ=255)
                                               cgccgTAGGGATCTTATGGCGGGCCAGGAAATCTTTTGTAAAACCTTTTAGACCTCCAGTTAGGCCGCACCTGCGGTTGGGCCGAATAttttt                                          >  2:731372/1‑92 (MQ=255)
                                                gccgTAGGGATCTTATGGCGGGCCAGGAAATCTTTTGTAAACGCTTTGAGCCTTCCAGTTGGGGCGCACCTGCGGGTGGGCCGGAGAttttta                                         >  2:1753762/1‑91 (MQ=255)
                                                  cgTAGGGATCTTATGGCGGGGCACGGAATCTTTGGTAAAAGCTTTTTGGCCTTCAGTTGGGGCGCACCCGCGGGTGGGGCGAAGAtttttcgc                                       >  2:42334/1‑89 (MQ=255)
                                                        gATCTTATGGCGGGCCAGGAAATCTTTTGTAAACGCTTTGAGCCTTCCCGTTGGGCCGCACCTGCGGTTGGGGCGAAGATTTTTAGCCcggcg                                 >  1:928909/1‑93 (MQ=255)
                                                         aTCTTATGGCGGGCCAGGAAATCTTTTGTAAACGCTTTGAGCCTTCCAGTTGGGGCGCACCTGCGGTTGGGCCCAAGATTTTTATCCcgggcg                                >  2:3551596/1‑90 (MQ=255)
                                                            ttATGGCGGGCCCGGAAATCTTTGGTAAAAGCTTTGAGCCTTCCAGTTGGGGCGCACCTGCGGGTGGGCCGAAGATTTTTAGCCCGGCGgcgc                             >  2:1990407/1‑93 (MQ=255)
                                                            ttATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGCAGATTTTCAGCCCGGCGgcgc                             >  1:2572290/1‑93 (MQ=255)
                                                             tATGGCGGGCCAGGAAATCTTTTGTAAACGCTTTTAGCCTTCCAGTTGGGCCGCACCTGCCGTTGGGGCGAAGATTTTCAGCCCGGCGGcgcg                            >  2:1694881/1‑93 (MQ=255)
                                                             tATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCCGCGGcgcg                            >  2:535017/1‑93 (MQ=255)
                                                               tGGCGGGCCAGGAAATCTTTGGTAAACGCTTTTAGCCTTCCAGTTGGGCCGCACCTGCGGGTTGGCCGAAGATTTTCAGCCCCGGGGCGCGGa                          >  2:2965079/1‑93 (MQ=255)
                                                               tGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGGCGCACCTGCGGTTGGGCCGAAGATTTTTAGCCCcgcggcgcggg                          >  2:24786/1‑92 (MQ=255)
                                                               tGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGa                          <  2:2191013/93‑1 (MQ=255)
                                                               tGGCGGGCCAGGAAATCTTTGGTAAAAGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGGTGGGCCGAAGGTTTTTCGCCCGGCGGCGCGGa                          >  2:1604737/1‑93 (MQ=255)
                                                                ggCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGaa                         >  2:1567124/1‑93 (MQ=255)
                                                                  cGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAgg                       <  1:4398837/93‑1 (MQ=255)
                                                                    gggCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGt                      >  1:2951755/1‑93 (MQ=255)
                                                                    gggCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGt                      <  1:2319618/93‑1 (MQ=255)
                                                                      gCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTAt                    <  2:596397/93‑1 (MQ=255)
                                                                        cAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCg                  >  2:4546038/1‑93 (MQ=255)
                                                                        cAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGGATCg                  >  2:4030699/1‑93 (MQ=255)
                                                                         aGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGa                 >  1:1778555/1‑93 (MQ=255)
                                                                         aGGAAATCTTTGGTAAACGTTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGa                 >  2:4547590/1‑93 (MQ=255)
                                                                          ggAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGAc                >  1:4482837/1‑93 (MQ=255)
                                                                          ggAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGAc                >  1:4433321/1‑93 (MQ=255)
                                                                          ggAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGAc                >  1:3797594/1‑93 (MQ=255)
                                                                          ggAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGAc                <  2:138458/93‑1 (MQ=255)
                                                                           gAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGAcc               >  1:3199746/1‑93 (MQ=255)
                                                                            aaaTCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCa              >  2:3620953/1‑93 (MQ=255)
                                                                            aaaTCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGGATCGAccc              >  2:959638/1‑92 (MQ=255)
                                                                             aaTCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCAc             >  1:3988462/1‑93 (MQ=255)
                                                                             aaTCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGGTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCAc             >  2:1494253/1‑93 (MQ=255)
                                                                              aTCTTTGGTAAACGCTTTGAGCCTTCCAGTTTGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACg            >  1:774395/1‑93 (MQ=255)
                                                                              aTCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACg            >  2:2765711/1‑93 (MQ=255)
                                                                              aTCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACg            <  1:4279863/93‑1 (MQ=255)
                                                                               tCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGc           >  1:1597040/1‑93 (MQ=255)
                                                                                cTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGcc          <  1:1401539/93‑1 (MQ=255)
                                                                                 tttGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCt         <  2:3739353/93‑1 (MQ=255)
                                                                                 tttGGTAAACGCTTTGAGCCTTCCAGTTGGGACGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCt         >  2:2064158/1‑93 (MQ=255)
                                                                                  ttGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCtt        <  2:2925177/93‑1 (MQ=255)
                                                                                    ggTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTGCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCTTTc      <  2:1152240/93‑1 (MQ=255)
                                                                                      tAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCTTTCAc    >  2:2750091/1‑93 (MQ=255)
                                                                                      tAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCTTTCAc    >  1:1808258/1‑93 (MQ=255)
                                                                                      tAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCTTTCAc    <  2:772721/93‑1 (MQ=255)
                                                                                      tAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCTTTCAc    >  1:1237278/1‑93 (MQ=255)
                                                                                        aaCGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCTTTCACCa  >  2:2381561/1‑93 (MQ=255)
                                                                                             |                                                                                        
CGCAGATACGCCAGCGCAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCTTTCACCA  >  NC_000913/4203103‑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
GATACGCCAGCGCAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCTTTCACCAGCGGCGCT  >  NC_000913/4203107‑4203291
                                                                                            |                                                                                             
GATACGCCAGCGCAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAG                                                                                            <  SRR4302535.4299506/94‑1 (MQ=60)
GATACGCCAGCGCAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAG                                                                                            <  SRR4302535.2941733/94‑1 (MQ=60)
       CAGCGCAGGTTCTACCTCGGTGAAGATCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCA                                                                                     <  SRR4302535.4559819/94‑1 (MQ=60)
        AGCGCAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAG                                                                                    <  SRR4302535.2784509/94‑1 (MQ=60)
         GCGCAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGGCAGGGAATATTTTGTTAACGCTTTGAGCCTTCCCGT                                                                                   >  SRR4302535.63119/1‑94 (MQ=60)
         GCGCAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGGCAGGGAAACTTTTGTAAACGCTTTTaggcttcccgt                                                                                    >  SRR4302535.2056180/1‑83 (MQ=60)
         GCGCAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGT                                                                                   <  SRR4302535.2244268/94‑1 (MQ=60)
          CGCAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGGCCGGGAAATTTTGGTAAAAGCTTTTGggcttcccatt                                                                                   >  SRR4302535.3305658/1‑83 (MQ=60)
          CGCAGGTTCTACCTCGGTGAAGTTCTGGTATGCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTT                                                                                  <  SRR4302535.737080/94‑1 (MQ=60)
           GCAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGGCAGGGAATCTTTTGTAAACGCTTTGAGGCTTCCCGTTG                                                                                 >  SRR4302535.57624/1‑94 (MQ=60)
            CAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGG                                                                                <  SRR4302535.1708953/94‑1 (MQ=60)
             AGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAAAGCTTTGAGCCCTCCAGTTGGG                                                                               >  SRR4302535.1521539/1‑94 (MQ=60)
              GGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCCGGGAATTTTTTGTAAACGCTTTTtgacctcccgttgggc                                                                               >  SRR4302535.1353686/1‑78 (MQ=60)
               GTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGGCAGGAAATCTTTTGTTAAAGCTTTGAGCCTTCCAGTTGGGGC                                                                             >  SRR4302535.2720892/1‑94 (MQ=60)
               GTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCC                                                                             <  SRR4302535.4672522/94‑1 (MQ=60)
                    ACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTTGTAAAAGCTTTTAGCCTTCCACTTGGGCCGCACC                                                                        >  SRR4302535.1667491/1‑94 (MQ=60)
                     CCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCT                                                                       <  SRR4302535.2378747/94‑1 (MQ=60)
                      CTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAAAGCTTTGAGCCTTCCAGTTGGGGCGCACCTG                                                                      >  SRR4302535.4215361/1‑94 (MQ=60)
                        CGGTGAAGTTCTGGTATTCCGCCGTATGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCG                                                                    <  SRR4302535.758382/94‑1 (MQ=60)
                          GTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAAACTTTGGTAAAAGCTTTGAGCCTTCCCGTTGGGACGCACCTGCGGG                                                                  >  SRR4302535.359130/1‑94 (MQ=60)
                           TGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGGCCGGGAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTT                                                                 >  SRR4302535.2117356/1‑94 (MQ=60)
                           TGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGGCAGGGAATCTTTTGTAAAAGCTTTGAGGCTTCCCGTTTGGCCGCCCCCGCGGTT                                                                 >  SRR4302535.2441647/1‑94 (MQ=41)
                           TGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGGCAGGAAAACTTTGGGAAAAGCTTTTAGGCTTCCAGTTGGGCCGCccccgcgggt                                                                 >  SRR4302535.583542/1‑84 (MQ=60)
                           TGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTT                                                                 <  SRR4302535.1727427/94‑1 (MQ=60)
                            GAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTG                                                                <  SRR4302535.3214306/94‑1 (MQ=60)
                                TTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTTGTAAACGCTTTGAGCCTTCCCGTTGGGCCGCACCTGCGGTTGGGGC                                                            >  SRR4302535.2391686/1‑94 (MQ=60)
                                     GTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGGCGCACCTGCGGTTGGGGCGGAGA                                                       >  SRR4302535.1401076/1‑94 (MQ=60)
                                     GTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGA                                                       <  SRR4302535.1136153/94‑1 (MQ=60)
                                       ATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTTGTAAACGCTTTGAGCCTTCCCGTTGGGCCGCACCTGCGGgtggggcggaggat                                                     >  SRR4302535.1039460/1‑80 (MQ=60)
                                       ATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAAAGCTTTGAGCCTTCCAGTTGGGCCGCCCCTGCGGTTGGGGCGAAGATT                                                     >  SRR4302535.2282978/1‑94 (MQ=60)
                                        TTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTTGTAAAAGCTTTGAGCCTTCCAGTTGGGGCGCACCTGCGGTTGGGGCGGAGATTT                                                    >  SRR4302535.1650622/1‑94 (MQ=60)
                                         ttccgccttagggaggttaggccggtcaaggaatgCTTTGGTAAACGCTTTGAGCCTACAAGTTGGTCCGCACCTGCGGTTGGGCCGAAGATTT                                                    <  SRR4302535.2589611/59‑1 (MQ=42)
                                         TCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTTGTAAACGCTTTTAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTT                                                   >  SRR4302535.2973211/1‑94 (MQ=60)
                                         TCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTT                                                   <  SRR4302535.2320744/94‑1 (MQ=60)
                                           CGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCA                                                 <  SRR4302535.91736/94‑1 (MQ=60)
                                             CCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTTGTAAACGCTTTTAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTTAGC                                               >  SRR4302535.1977619/1‑94 (MQ=60)
                                               GTAGGGATCTTATGGCGGGGCAGGCAAATTTTGGGAAACGCTTTGAGCCTTCCAGTTTGGGCGCACCTGCGGGTGGGGCGAAGATTTTGATCCC                                             >  SRR4302535.2950369/1‑94 (MQ=26)
                                                TAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCG                                            <  SRR4302535.2448403/94‑1 (MQ=60)
                                                 AGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGG                                           <  SRR4302535.4437382/94‑1 (MQ=60)
                                                   GGATCTTATGGCGGGCCAGGAAATCTTTTGTAAACGCTTTGAGCCTTCCCGTTGGGCCGCACCTGCGGTTGGGGCGAAGATTTTTAGCCCGGCG                                         >  SRR4302535.935483/1‑94 (MQ=60)
                                                    GATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGG                                        <  SRR4302535.2946552/94‑1 (MQ=60)
                                                    GATCTTATGGCGGGCCAGGAAAACTTTGGTAAACGCTTTTAGCCTTCCAGTTGGGCCGCACCTGCGGGTTGGCCGAAGATTTTCAGCCCGGCGG                                        >  SRR4302535.1278666/1‑94 (MQ=60)
                                                       CTTATGGCGGGCCAGGAAATCTTTTGGAAACGCTTTTAGCCTTCCCGGTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGC                                     >  SRR4302535.15592/1‑94 (MQ=60)
                                                       CTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGCAGATTTTCAGCCCGGCGGCGC                                     >  SRR4302535.2590778/1‑94 (MQ=60)
                                                              CGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGT                              >  SRR4302535.2972491/1‑94 (MQ=60)
                                                              CGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGT                              >  SRR4302535.1332406/1‑94 (MQ=60)
                                                              CGGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGT                              <  SRR4302535.4428944/94‑1 (MQ=60)
                                                                GGGCCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTA                             <  SRR4302535.2336223/94‑1 (MQ=60)
                                                                   CCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCG                          >  SRR4302535.4354203/1‑94 (MQ=60)
                                                                   CCAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATAG                          >  SRR4302535.3528048/1‑94 (MQ=60)
                                                                    CAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGA                         >  SRR4302535.1791433/1‑94 (MQ=60)
                                                                    CAGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGA                         >  SRR4302535.21165/1‑94 (MQ=60)
                                                                     AGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGAC                        >  SRR4302535.4513603/1‑94 (MQ=60)
                                                                     AGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGAC                        >  SRR4302535.4463699/1‑94 (MQ=60)
                                                                     AGGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGAC                        >  SRR4302535.3823772/1‑94 (MQ=60)
                                                                      GGAAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACC                       >  SRR4302535.3221900/1‑94 (MQ=60)
                                                                        AAATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCAC                     >  SRR4302535.4015738/1‑94 (MQ=60)
                                                                         AATCTTTGGTAAACGCTTTGAGCCTTCCAGTTTGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACG                    >  SRR4302535.779842/1‑94 (MQ=60)
                                                                          ATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGC                   >  SRR4302535.1608558/1‑94 (MQ=60)
                                                                          ATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGC                   <  SRR4302535.4309149/94‑1 (MQ=60)
                                                                          ATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGC                   >  SRR4302535.3794988/1‑94 (MQ=60)
                                                                          ATCTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGC                   >  SRR4302535.3237666/1‑94 (MQ=60)
                                                                            CTTTGGTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCT                 <  SRR4302535.1411542/94‑1 (MQ=60)
                                                                                 GTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCTTTCAC            >  SRR4302535.1646890/1‑94 (MQ=60)
                                                                                 GTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCTTTCAC            >  SRR4302535.1821346/1‑94 (MQ=60)
                                                                                 GTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCTTTCAC            >  SRR4302535.2332085/1‑94 (MQ=60)
                                                                                 GTAAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCTTTCAC            >  SRR4302535.1246114/1‑94 (MQ=60)
                                                                                   AAACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCCAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCTTTCACCA          <  SRR4302535.1393656/94‑1 (MQ=60)
                                                                                    AACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCTTTCACCAG         >  SRR4302535.4530632/1‑94 (MQ=60)
                                                                                    AACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCTTTCACCAG         >  SRR4302535.2019923/1‑94 (MQ=60)
                                                                                    AACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCTTTCACCAG         >  SRR4302535.1381658/1‑94 (MQ=60)
                                                                                     aacgctaTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCTTTCACCAG         <  SRR4302535.3056633/87‑1 (MQ=60)
                                                                                     ACGCTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATAGACCACGCCTTTCACCAGC        >  SRR4302535.574015/1‑94 (MQ=60)
                                                                                         CTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCTTTCACCAGCGGC     <  SRR4302535.1595529/94‑1 (MQ=60)
                                                                                            TGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCTTTCACCAGCGGCGCT  <  SRR4302535.1731437/94‑1 (MQ=60)
                                                                                            TGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCTTTCACCAGCGGCGCT  >  SRR4302535.2316260/1‑94 (MQ=60)
                                                                                            |                                                                                             
GATACGCCAGCGCAGGTTCTACCTCGGTGAAGTTCTGGTATTCCGCCGTAGGGATCTTATGGCGGGCCAGGAAATCTTTGGTAAACGCTTTTGAGCCTTCCAGTTGGGCCGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCTTTCACCAGCGGCGCT  >  NC_000913/4203107‑4203291

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