Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A1 F1 I1 R1 521 238.2 12299412 96.0% 11807435 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,051,170 +C intergenic (+172/+279) hemN → / ← glnG O2‑independent coproporphyrinogen III oxidase/GlnG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,051,1701.C99.5% 618.0 / ‑4.4 210intergenic (+172/+279)hemN/glnGO2‑independent coproporphyrinogen III oxidase/GlnG
Reads supporting (aligned to +/- strand):  ref base . (0/0);  major base C (105/104);  minor base T (1/0);  total (106/104)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.28e-01

Only 100 of 225 total aligned reads displayed.

TGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCCATCTGAGCATGATGCCTGATGCGGCGTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCAC  >  NC_000913/4051077‑4051255
                                                                                              |                                                                                     
tGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATgcggcg                                                                                        >  2:2006949/1‑93 (MQ=255)
 gggCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATgcggcgc                                                                                       >  2:4560067/1‑92 (MQ=255)
  ggCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATgcggcgcg                                                                                      >  1:882221/1‑91 (MQ=255)
    caacaaAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATgcggcgctgg                                                                                    <  1:4378246/93‑5 (MQ=255)
       caaAACCGGTTGAAATTATTGCTTTTGTAGACCGGATAAGGCGTTTACGCCGCATCAGGCATCTGCGCATGATGCCTGATGCGGCGCTGGTcg                                                                                 >  1:2072998/1‑93 (MQ=21)
         aaaCCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTcgcg                                                                               >  1:1121799/1‑93 (MQ=25)
           aCCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTc                                                                             >  1:1660173/1‑93 (MQ=21)
            ccTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCt                                                                            >  1:5418846/1‑93 (MQ=21)
             cTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTATTCCGCATCAGGCATCTTAGCATGATTTCTGATGCGGCGCCTGTCGCGTCtt                                                                           >  1:5452786/1‑93 (MQ=18)
                ttGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATAAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGGCGCGTCTTATc                                                                        >  2:5377606/1‑93 (MQ=18)
                 tGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCa                                                                       <  2:1163650/93‑1 (MQ=21)
                  gAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAg                                                                      >  2:1068044/1‑93 (MQ=18)
                    aaTTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGc                                                                    >  1:4230816/1‑93 (MQ=18)
                     attattGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGcc                                                                   >  2:318192/1‑93 (MQ=18)
                     attattGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGcc                                                                   <  1:3657102/93‑1 (MQ=18)
                       tattGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATAAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTa                                                                 <  2:4755639/93‑1 (MQ=17)
                         ttGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACg                                                               <  1:3131878/93‑1 (MQ=18)
                         ttGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACg                                                               <  1:454720/93‑1 (MQ=18)
                          tGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTCCGa                                                              <  1:5991476/93‑1 (MQ=17)
                          tGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGa                                                              <  1:2509987/93‑1 (MQ=18)
                          tGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGa                                                              <  2:741477/93‑1 (MQ=18)
                             ttttGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACtt                                                           <  2:5558758/93‑1 (MQ=18)
                             ttttGTAGGCCGGATAAGGCGTTTACGACGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACtt                                                           >  1:6096963/1‑93 (MQ=17)
                               ttGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCGGGCCTACGACTTtc                                                         <  1:2311048/93‑1 (MQ=18)
                                 gTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTtctc                                                       <  2:2463662/93‑1 (MQ=18)
                                 gTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTtctc                                                       <  1:2349533/93‑1 (MQ=18)
                                 gTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTtctc                                                       >  2:4546331/1‑93 (MQ=18)
                                 gTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTtctc                                                       <  1:4882898/93‑1 (MQ=18)
                                  tAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTcc                                                      <  2:2593774/93‑1 (MQ=18)
                                  tAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTcc                                                      <  2:4028703/93‑1 (MQ=18)
                                     gCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCa                                                   >  2:202775/1‑93 (MQ=18)
                                      ccGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCAt                                                  <  1:2730968/93‑1 (MQ=18)
                                        ggATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATcc                                                >  2:978032/1‑93 (MQ=18)
                                         gATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCg                                               >  2:2547964/1‑93 (MQ=18)
                                         gATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCg                                               >  1:2650023/1‑93 (MQ=18)
                                         gATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCg                                               >  2:1000482/1‑93 (MQ=18)
                                           tAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCTCATCCGGc                                             >  1:5184064/1‑93 (MQ=18)
                                           tAAGGCGTTTACGCCGCATCAGGCACCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGc                                             >  1:5862026/1‑93 (MQ=21)
                                             aGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCAt                                           <  2:5212749/93‑1 (MQ=21)
                                             aGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCAt                                           <  2:1685549/93‑1 (MQ=21)
                                              ggCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATa                                          >  2:2814573/1‑93 (MQ=21)
                                              ggCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATa                                          >  1:778604/1‑93 (MQ=21)
                                                cGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATaaa                                        <  1:1105391/93‑1 (MQ=21)
                                                 ttttACGCCGCATCAGGCATCTGAGCATGATGCTTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaac                                       <  2:1850820/92‑1 (MQ=21)
                                                 gTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaac                                       >  2:3486446/1‑93 (MQ=25)
                                                  tttACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaaca                                      >  1:4723611/1‑93 (MQ=25)
                                                  tttACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaaca                                      >  1:3383028/1‑93 (MQ=21)
                                                  tttACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaaca                                      >  1:4949473/1‑93 (MQ=25)
                                                  tttACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaaca                                      >  2:1277090/1‑93 (MQ=25)
                                                   ttACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaacaa                                     <  2:3695059/93‑1 (MQ=21)
                                                   ttACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaacaa                                     <  1:3310705/93‑1 (MQ=21)
                                                    tACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaacaac                                    >  1:3869664/1‑93 (MQ=25)
                                                    tACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaacaac                                    <  1:1342687/93‑1 (MQ=25)
                                                      cgccgcATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGa                                  <  2:2674632/93‑1 (MQ=25)
                                                       gccgcATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAg                                 >  2:5975257/1‑93 (MQ=25)
                                                        ccgcATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGc                                <  2:1987225/93‑1 (MQ=25)
                                                        ccgcATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGc                                >  1:3058832/1‑93 (MQ=37)
                                                          gcATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCtt                              >  1:1926554/1‑93 (MQ=37)
                                                           cATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTc                             <  2:5102086/93‑1 (MQ=37)
                                                           cATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTc                             <  2:233120/93‑1 (MQ=37)
                                                            aTCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGCCTTTCTCCGCATCCGGCATAAACAACGAGCTTCa                            <  2:3361270/93‑1 (MQ=25)
                                                             tCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAg                           <  2:236444/93‑1 (MQ=37)
                                                                aggcATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCt                         >  1:4972618/4‑93 (MQ=255)
                                                                 ggcATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTa                        >  1:919091/3‑93 (MQ=255)
                                                                  gcATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACACCGAGCTTCAGCTaa                       <  2:3852088/92‑1 (MQ=37)
                                                                    aTCTGAGCATGATTCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAAc                     <  2:1877100/93‑1 (MQ=255)
                                                                    aTCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAAc                     >  1:4011504/1‑93 (MQ=255)
                                                                    aTCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAAc                     >  1:3606235/1‑93 (MQ=255)
                                                                    aTCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAAc                     <  2:6147575/93‑1 (MQ=255)
                                                                     tCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACa                    >  1:2701982/1‑93 (MQ=255)
                                                                     tCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACa                    <  2:1362836/93‑1 (MQ=255)
                                                                     tCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACa                    <  1:4695420/93‑1 (MQ=255)
                                                                     tCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACa                    >  2:3940066/1‑93 (MQ=255)
                                                                     tCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACa                    >  1:649195/1‑93 (MQ=255)
                                                                     tCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACa                    >  2:34071/1‑93 (MQ=255)
                                                                      cTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAg                   <  2:291161/93‑1 (MQ=255)
                                                                       tGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGc                  >  1:5099297/1‑93 (MQ=255)
                                                                       tGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGc                  >  2:1973401/1‑93 (MQ=255)
                                                                         aGCATGATGCCTGATCCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGccc                >  1:225484/1‑93 (MQ=255)
                                                                          gCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCa               <  2:3603099/93‑1 (MQ=255)
                                                                           cATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCaa              >  1:1774692/1‑93 (MQ=255)
                                                                           cATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCaa              <  1:2874858/93‑1 (MQ=255)
                                                                           cATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCaa              >  2:4912453/1‑93 (MQ=255)
                                                                            atgatgCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAAt             <  1:2041858/93‑1 (MQ=255)
                                                                             tgatgCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATc            <  2:1710141/93‑1 (MQ=255)
                                                                               atgCCTGATGCGGCGTTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCAt          >  1:1824119/1‑93 (MQ=255)
                                                                               atgCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCAt          <  2:5640427/93‑1 (MQ=255)
                                                                                 gCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTg        <  1:4652271/93‑1 (MQ=255)
                                                                                 gCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTg        <  2:5932318/93‑1 (MQ=255)
                                                                                 gCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTg        <  2:2983551/93‑1 (MQ=255)
                                                                                 gCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCAGCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTg        <  1:970878/93‑1 (MQ=255)
                                                                                  ccTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGc       <  1:2058169/93‑1 (MQ=255)
                                                                                   cTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCg      >  2:5767294/1‑93 (MQ=255)
                                                                                   cTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCg      >  1:3698186/1‑93 (MQ=255)
                                                                                    tGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACTAGCTTCAGCTAAACAGCCCAATCATTGCgg     >  2:5083807/1‑93 (MQ=255)
                                                                                    tGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCgg     <  2:4195024/93‑1 (MQ=255)
                                                                                     gATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGc    <  2:280929/93‑1 (MQ=255)
                                                                                      aTGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGca   <  1:4343380/93‑1 (MQ=255)
                                                                                       tGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCACAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGcac  >  1:2138081/1‑93 (MQ=255)
                                                                                       ggcggcgcTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGcac  <  1:1218172/92‑1 (MQ=255)
                                                                                              |                                                                                     
TGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCCATCTGAGCATGATGCCTGATGCGGCGTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCAC  >  NC_000913/4051077‑4051255

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCCATCTGAGCATGATGCCTGATGCGGCG‑‑TGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACG  >  NC_000913/4051076‑4051263
                                                                                              |                                                                                               
CTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCG                                                                                                 <  SRR4302532.4204825/94‑1 (MQ=60)
  GGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCG‑‑CG                                                                                             >  SRR4302532.912382/1‑94 (MQ=60)
   GGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGctg                                                                                              <  SRR4302532.6277471/94‑4 (MQ=60)
     CAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGT                                                                                           >  SRR4302532.3643740/1‑94 (MQ=60)
     CAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGT                                                                                           <  SRR4302532.4526925/94‑1 (MQ=60)
      AACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTC                                                                                          >  SRR4302532.99082/1‑94 (MQ=60)
      AACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTC                                                                                          <  SRR4302532.5438056/94‑1 (MQ=60)
       ACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCG                                                                                         >  SRR4302532.2142088/1‑94 (MQ=60)
       ACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCG                                                                                         >  SRR4302532.5386803/1‑94 (MQ=60)
       ACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCG                                                                                         >  SRR4302532.2113127/1‑94 (MQ=60)
       ACAAAACCGGTTGAAATTATTGCTTTTGTAGACCGGATAAGGCGTTTACGCCGCATCAGGCATCTGCGCATGATGCCTGATGCGGCGCTGGTCG                                                                                         >  SRR4302532.2142089/1‑94 (MQ=60)
        CAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGC                                                                                        <  SRR4302532.3762150/94‑1 (MQ=60)
        CAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGC                                                                                        >  SRR4302532.2289966/1‑94 (MQ=60)
         AAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCG                                                                                       >  SRR4302532.1159257/1‑94 (MQ=60)
          AAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGT                                                                                      >  SRR4302532.2924219/1‑94 (MQ=60)
           AACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTC                                                                                     >  SRR4302532.1716432/1‑94 (MQ=60)
            ACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCT                                                                                    >  SRR4302532.5606660/1‑94 (MQ=60)
            ACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCT                                                                                    <  SRR4302532.6149424/94‑1 (MQ=60)
             CCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTATTCCGCATCAGGCATCTTAGCATGATTTCTGATGCGGCGCCTGTCGCGTCTT                                                                                   >  SRR4302532.5641788/1‑94 (MQ=37)
             CCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGGTTACGGCGGATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTT                                                                                   >  SRR4302532.5360582/1‑94 (MQ=60)
              CTGTTTAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACTCCGCATCAGGCATCTGAGCATGATGCCTGACGCGGCGCTGGTCGCGTCTTA                                                                                  <  SRR4302532.4365667/94‑1 (MQ=60)
               TGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGAATGATGCCTGATGCGGCGCTGGTCGCGTCTTAT                                                                                 <  SRR4302532.919752/94‑1 (MQ=60)
                 TTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCA                                                                               >  SRR4302532.4857582/1‑94 (MQ=60)
                  TGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAG                                                                              <  SRR4302532.4415761/94‑1 (MQ=60)
                  TGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAG                                                                              <  SRR4302532.5129028/94‑1 (MQ=60)
                    AAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGC                                                                            <  SRR4302532.1911414/94‑1 (MQ=60)
                    AAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGC                                                                            >  SRR4302532.4372839/1‑94 (MQ=60)
                    AAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGC                                                                            >  SRR4302532.3783412/1‑94 (MQ=60)
                     AATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCC                                                                           <  SRR4302532.258979/94‑1 (MQ=60)
                     AATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCC                                                                           >  SRR4302532.3197236/1‑94 (MQ=60)
                      ATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCT                                                                          <  SRR4302532.4767639/94‑1 (MQ=60)
                      ATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCT                                                                          <  SRR4302532.5099618/94‑1 (MQ=60)
                      ATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCT                                                                          <  SRR4302532.3779887/94‑1 (MQ=60)
                       TTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTC                                                                         >  SRR4302532.3315298/1‑94 (MQ=60)
                         cttcCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCAGGCATCTGAGCATGATGCCTGAGGCCGCGCTGGTCGCGTCTTATCAGGCCTACG                                                                       <  SRR4302532.4422190/90‑1 (MQ=49)
                          TTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGA                                                                      <  SRR4302532.470544/94‑1 (MQ=60)
                          TTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGA                                                                      <  SRR4302532.784452/94‑1 (MQ=60)
                          TTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGA                                                                      <  SRR4302532.3235844/94‑1 (MQ=60)
                           TGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTCCGAC                                                                     <  SRR4302532.6198479/94‑1 (MQ=60)
                           TGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGAC                                                                     <  SRR4302532.6361436/94‑1 (MQ=60)
                           TGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGAC                                                                     <  SRR4302532.2593728/94‑1 (MQ=60)
                           TGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGAC                                                                     >  SRR4302532.3621325/1‑94 (MQ=60)
                            GCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACT                                                                    <  SRR4302532.889552/94‑1 (MQ=60)
                             GTTTTATAGGCCGGATAAGGCGTTTCCGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTCCGACTT                                                                   <  SRR4302532.4053617/94‑1 (MQ=60)
                             CTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTT                                                                   >  SRR4302532.3982144/1‑94 (MQ=60)
                             CTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACAACTT                                                                   >  SRR4302532.5903372/1‑94 (MQ=60)
                             CTTTTGTAGGCCGGATAAGGCGTTTACGACGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTT                                                                   >  SRR4302532.6307408/1‑94 (MQ=60)
                              TTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTT                                                                  <  SRR4302532.3124108/94‑1 (MQ=60)
                              TTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTT                                                                  >  SRR4302532.299300/1‑94 (MQ=60)
                               TTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTC                                                                 <  SRR4302532.5004551/94‑1 (MQ=60)
                                TTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCGGGCCTACGACTTTCT                                                                <  SRR4302532.2387949/94‑1 (MQ=60)
                                TTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCT                                                                >  SRR4302532.4319498/1‑94 (MQ=60)
                                  GTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGTCCTACGACTTTCTCC                                                              <  SRR4302532.4472833/94‑1 (MQ=60)
                                  GTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCC                                                              <  SRR4302532.5050778/94‑1 (MQ=60)
                                  GTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCC                                                              <  SRR4302532.2427905/94‑1 (MQ=60)
                                    AGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGC                                                            >  SRR4302532.1005392/1‑94 (MQ=60)
                                     gcccggctaagcccTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATTAGACCTACGACTTTCTCCGCA                                                           <  SRR4302532.2707737/80‑1 (MQ=60)
                                     GGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCA                                                           >  SRR4302532.2381111/1‑94 (MQ=60)
                                       CCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATC                                                         >  SRR4302532.61678/1‑94 (MQ=60)
                                       CCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATC                                                         <  SRR4302532.2822785/94‑1 (MQ=60)
                                         GGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCG                                                       >  SRR4302532.2738924/1‑94 (MQ=60)
                                         GGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCG                                                       <  SRR4302532.4505950/94‑1 (MQ=60)
                                          AATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGG                                                      >  SRR4302532.1793156/1‑94 (MQ=60)
                                           ATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCTCATCCGGC                                                     >  SRR4302532.5362657/1‑94 (MQ=60)
                                           ATAAGGCGTTTACGCCGCATCAGGCACCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGC                                                     >  SRR4302532.6065267/1‑94 (MQ=60)
                                             AAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCAT                                                   <  SRR4302532.3486182/94‑1 (MQ=60)
                                             AAGGCGTTTACGCCGCATCAGGCATCTGAACATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCAT                                                   <  SRR4302532.4972392/94‑1 (MQ=60)
                                              AGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATA                                                  >  SRR4302532.805404/1‑94 (MQ=60)
                                                 CGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAAC                                               <  SRR4302532.1142344/94‑1 (MQ=60)
                                                  GTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACA                                              >  SRR4302532.3496170/1‑94 (MQ=60)
                                                  GTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACA                                              >  SRR4302532.5119317/1‑94 (MQ=60)
                                                  GTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACA                                              >  SRR4302532.4830585/1‑94 (MQ=60)
                                                  GTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACA                                              >  SRR4302532.4884781/1‑94 (MQ=60)
                                                  GTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACA                                              >  SRR4302532.2031896/1‑94 (MQ=60)
                                                    TTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAAC                                            <  SRR4302532.3421076/94‑1 (MQ=60)
                                                    TTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAAC                                            >  SRR4302532.4000247/1‑94 (MQ=60)
                                                     TACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACG                                           <  SRR4302532.5118013/94‑1 (MQ=60)
                                                     TACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACG                                           <  SRR4302532.1387577/94‑1 (MQ=60)
                                                     TACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACG                                           <  SRR4302532.2150119/94‑1 (MQ=60)
                                                      ACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGA                                          <  SRR4302532.1615901/94‑1 (MQ=60)
                                                       CGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAG                                         >  SRR4302532.2616211/1‑94 (MQ=60)
                                                       CGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAG                                         <  SRR4302532.3602556/94‑1 (MQ=60)
                                                        GCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGC                                        >  SRR4302532.3159831/1‑94 (MQ=60)
                                                         CCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCT                                       >  SRR4302532.1321485/1‑94 (MQ=60)
                                                          CGCATCAGGCATCTGAGCATGATGCCTGATGGGGGGCTGGTGGCGTCTTATCAGGCCTGCGACTTTCTCCGCATCCGGCAGAAACAATGAGCTT                                      >  SRR4302532.1990651/1‑94 (MQ=49)
                                                          CGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTT                                      >  SRR4302532.1990653/1‑94 (MQ=60)
                                                          CGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTT                                      >  SRR4302532.4066497/1‑94 (MQ=60)
                                                               CAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCCTAAACAACGAGCTTCAGCT                                 >  SRR4302532.1877879/1‑94 (MQ=60)
                                                               CAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCT                                 >  SRR4302532.5143085/1‑94 (MQ=60)
                                                                AGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTA                                >  SRR4302532.950545/1‑94 (MQ=60)
                                                                    GATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAAC                             <  SRR4302532.874590/94‑1 (MQ=60)
                                                                    CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAAC                             >  SRR4302532.3726774/1‑94 (MQ=60)
                                                                    CATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAAC                             >  SRR4302532.4146673/1‑94 (MQ=60)
                                                                     ATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACA                            >  SRR4302532.1430121/1‑94 (MQ=60)
                                                                     ATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACA                            >  SRR4302532.671720/1‑94 (MQ=60)
                                                                     ATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACA                            >  SRR4302532.2792763/1‑94 (MQ=60)
                                                                      TCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAG                           <  SRR4302532.4855571/94‑1 (MQ=60)
                                                                       CTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGC                          >  SRR4302532.5274618/1‑94 (MQ=60)
                                                                         GAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCC                        >  SRR4302532.233183/1‑94 (MQ=60)
                                                                         GAGCATGATGCCTGATCCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCC                        >  SRR4302532.233184/1‑94 (MQ=60)
                                                                          AGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCA                       >  SRR4302532.4946512/1‑94 (MQ=60)
                                                                          AGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCA                       <  SRR4302532.2809690/94‑1 (MQ=60)
                                                                           GCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAA                      >  SRR4302532.1834126/1‑94 (MQ=60)
                                                                           GCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAA                      >  SRR4302532.2739607/1‑94 (MQ=60)
                                                                            CATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAAT                     >  SRR4302532.532405/1‑94 (MQ=60)
                                                                            CATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAAT                     <  SRR4302532.2970981/94‑1 (MQ=60)
                                                                             ATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATC                    <  SRR4302532.2109859/94‑1 (MQ=60)
                                                                             ATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATC                    <  SRR4302532.2605234/94‑1 (MQ=60)
                                                                               GATGCCTGATGCGGCGTTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCAT                  >  SRR4302532.1884864/1‑94 (MQ=60)
                                                                               GATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCAT                  >  SRR4302532.5630129/1‑94 (MQ=60)
                                                                                 TGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTG                <  SRR4302532.1156579/94‑1 (MQ=60)
                                                                                 TGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTG                <  SRR4302532.4683859/94‑1 (MQ=60)
                                                                                  GCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGC               <  SRR4302532.4810884/94‑1 (MQ=60)
                                                                                  GCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCAGCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGC               <  SRR4302532.1003999/94‑1 (MQ=60)
                                                                                   CCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCG              <  SRR4302532.2126718/94‑1 (MQ=60)
                                                                                   CCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCG              >  SRR4302532.3822343/1‑94 (MQ=60)
                                                                                   CCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCG              <  SRR4302532.1993489/94‑1 (MQ=60)
                                                                                   CCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCG              >  SRR4302532.5263516/1‑94 (MQ=60)
                                                                                   CCTAATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCG              <  SRR4302532.1175732/94‑1 (MQ=60)
                                                                                      GATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCA           >  SRR4302532.4594746/1‑94 (MQ=60)
                                                                                      GATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCA           <  SRR4302532.2404787/94‑1 (MQ=60)
                                                                                       ATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCAC          <  SRR4302532.4490095/94‑1 (MQ=60)
                                                                                       ATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCACAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCAC          >  SRR4302532.2209484/1‑94 (MQ=60)
                                                                                          CGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACA       <  SRR4302532.5243499/94‑1 (MQ=60)
                                                                                          CGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACA       >  SRR4302532.2002920/1‑94 (MQ=60)
                                                                                           GGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAG      <  SRR4302532.37066/94‑1 (MQ=60)
                                                                                           GGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAG      >  SRR4302532.4064133/1‑94 (MQ=60)
                                                                                              C‑‑TGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACG  <  SRR4302532.1213642/94‑1 (MQ=60)
                                                                                              C‑‑TGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACG  <  SRR4302532.1292940/94‑1 (MQ=60)
                                                                                              C‑‑TGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACG  <  SRR4302532.410666/94‑1 (MQ=60)
                                                                                              C‑‑TGGTCGCGTCTTATCAGGCATACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACG  <  SRR4302532.3237354/94‑1 (MQ=60)
                                                                                              |                                                                                               
CTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCCATCTGAGCATGATGCCTGATGCGGCG‑‑TGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACG  >  NC_000913/4051076‑4051263

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