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,144 Δ1 bp intergenic (+146/+305) 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,1430C.99.5% 731.9 / ‑3.5 217intergenic (+145/+306)hemN/glnGO2‑independent coproporphyrinogen III oxidase/GlnG
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base . (106/110);  minor base G (1/0);  total (107/110)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 239 total aligned reads displayed.

GCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCCATCTGAGCATGATGCCTGATGCGGCGTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAG  >  NC_000913/4051051‑4051230
                                                                                            |                                                                                        
gCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGc                                                                                          >  2:2676742/1‑93 (MQ=255)
gCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGc                                                                                          <  2:269573/93‑1 (MQ=255)
 cTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGca                                                                                         >  1:1029349/1‑92 (MQ=255)
 cTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGca                                                                                         <  1:5750598/93‑2 (MQ=255)
 cTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGACGCATCAGGca                                                                                         >  1:2534181/1‑92 (MQ=255)
     gaACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTg                                                                                    <  1:4079956/93‑1 (MQ=255)
     gaACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTg                                                                                    <  1:1383514/93‑1 (MQ=255)
     gaACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTg                                                                                    <  2:4453611/93‑1 (MQ=255)
     gaACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTg                                                                                    <  2:6109883/93‑1 (MQ=255)
     gaACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTg                                                                                    >  1:2810498/1‑93 (MQ=255)
      aaCAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGa                                                                                   >  2:1713680/1‑93 (MQ=255)
       aCAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAg                                                                                  >  2:4693577/1‑93 (MQ=255)
        cAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGc                                                                                 <  1:5048327/93‑1 (MQ=255)
        cAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGc                                                                                 <  2:5179165/93‑1 (MQ=255)
           gATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCatg                                                                              <  1:2330625/93‑1 (MQ=255)
           gATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCatg                                                                              >  1:3726594/1‑93 (MQ=255)
            aTTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCatga                                                                             <  1:4910410/93‑1 (MQ=255)
            aTTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCatga                                                                             <  1:1255930/93‑1 (MQ=255)
             ttGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCatgat                                                                            <  1:5997059/93‑1 (MQ=255)
               gCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGc                                                                          >  2:5991160/1‑93 (MQ=255)
               gCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGc                                                                          >  2:5381737/1‑93 (MQ=255)
                ccTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGcc                                                                         <  2:5079334/93‑1 (MQ=255)
                ccTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGcc                                                                         >  2:4929624/1‑93 (MQ=255)
                ccTGAGGCGCCGGGCAACAAAACCTGTTGAAACTATTGCTTTTGTAGGCCGGAGAAGGCGTGTACGCCGCATCAGGCATCTGAGCAGGAGGcc                                                                         >  2:5489982/1‑93 (MQ=21)
                 cTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCt                                                                        >  1:1087628/1‑93 (MQ=255)
                 cTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGGTTACGCCGCATCAGGCATCTGATCATGATGCCt                                                                        >  2:312928/1‑93 (MQ=37)
                   gATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGa                                                                      >  2:4857391/1‑93 (MQ=255)
                   gATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGa                                                                      >  1:2865409/1‑93 (MQ=255)
                   gATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGa                                                                      <  1:2611414/93‑1 (MQ=255)
                   gATGAGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGa                                                                      >  1:4036467/1‑93 (MQ=255)
                    aTGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGAt                                                                     >  2:5343380/1‑93 (MQ=255)
                    aTGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGAt                                                                     <  2:2768395/93‑1 (MQ=255)
                      gcgcTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATgc                                                                   <  2:1165776/93‑1 (MQ=255)
                       cgcTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATgcg                                                                  <  2:227182/93‑1 (MQ=255)
                        gcTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATgcgg                                                                 >  2:3540850/1‑93 (MQ=255)
                        gcTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATgcgg                                                                 >  2:1462215/1‑93 (MQ=255)
                        gcTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATgcgg                                                                 <  2:3854159/93‑1 (MQ=255)
                         cTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATgcggc                                                                <  2:860672/93‑1 (MQ=255)
                          tGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATgcggcg                                                               >  2:2006949/1‑93 (MQ=255)
                            ggCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATgcggcgct                                                             <  2:4348855/93‑3 (MQ=255)
                                acaaAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTc                                                         >  1:95743/1‑93 (MQ=25)
                                 caaAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTcg                                                        <  1:3640011/93‑1 (MQ=25)
                                 caaAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTcg                                                        >  1:2072997/1‑93 (MQ=21)
                                  aaaaCCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTcgc                                                       <  2:1066579/93‑1 (MQ=25)
                                   aaaCCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTcgcg                                                      >  1:1121799/1‑93 (MQ=25)
                                     aCCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTc                                                    >  1:1660173/1‑93 (MQ=21)
                                      ccTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCt                                                   >  1:5418846/1‑93 (MQ=21)
                                        tGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGAATGATGCCTGATGCGGCGCTGGTCGCGTCTTa                                                 <  1:889368/93‑1 (MQ=21)
                                           tGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCa                                              <  1:4271921/93‑1 (MQ=18)
                                           tGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCa                                              <  2:5939640/93‑1 (MQ=18)
                                             aaaTTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAgg                                            <  1:1850013/93‑1 (MQ=18)
                                               attattGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGcc                                          <  1:4930316/93‑1 (MQ=18)
                                               attattGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGcc                                          >  2:318192/1‑93 (MQ=18)
                                               attattGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGcc                                          >  1:3094645/1‑93 (MQ=18)
                                                 tattGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCtc                                        >  1:3208752/1‑92 (MQ=18)
                                                 tattGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATAAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTa                                        <  2:4755639/93‑1 (MQ=17)
                                                   ttGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACg                                      <  1:758414/93‑1 (MQ=17)
                                                   ttGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACg                                      <  1:454720/93‑1 (MQ=18)
                                                    tGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGa                                     <  1:2509987/93‑1 (MQ=18)
                                                     gCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGAc                                    <  1:860012/93‑1 (MQ=18)
                                                       ttttGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACtt                                  >  1:3852088/1‑93 (MQ=17)
                                                       ttttGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACtt                                  >  2:2938317/1‑93 (MQ=18)
                                                       ttttGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACtt                                  <  2:5558758/93‑1 (MQ=18)
                                                       ttttGTAGGCCGGATAAGGCGTTTACGACGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACtt                                  >  1:6096963/1‑93 (MQ=17)
                                                        tttGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACttt                                 <  2:1075122/93‑1 (MQ=17)
                                                           gTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTtctc                              <  2:2463662/93‑1 (MQ=18)
                                                           gTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTtctc                              >  2:4546331/1‑93 (MQ=18)
                                                            tAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTcc                             <  2:4028703/93‑1 (MQ=18)
                                                               gCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCa                          >  1:2304490/1‑93 (MQ=18)
                                                               gCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCa                          >  2:202775/1‑93 (MQ=18)
                                                                 cGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATc                        >  1:59531/1‑93 (MQ=18)
                                                                  ggATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATcc                       >  2:978032/1‑93 (MQ=18)
                                                                  ggATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATcc                       <  1:4358315/93‑1 (MQ=18)
                                                                  ggATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATcc                       <  2:3703274/93‑1 (MQ=18)
                                                                   gATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCg                      >  2:2547964/1‑93 (MQ=18)
                                                                   gATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCg                      <  2:1661634/93‑1 (MQ=18)
                                                                   gATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCg                      >  2:1000482/1‑93 (MQ=18)
                                                                    aTAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCgg                     >  2:1857551/1‑93 (MQ=18)
                                                                     tAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCTCATCCGGc                    >  1:5184064/1‑93 (MQ=18)
                                                                     tAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGc                    >  2:5872823/1‑93 (MQ=21)
                                                                     tAAGGCGTTTACGCCGCATCAGGCACCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGc                    >  1:5862026/1‑93 (MQ=21)
                                                                       aGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCAt                  <  2:1685549/93‑1 (MQ=21)
                                                                        ggCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATa                 >  2:2814573/1‑93 (MQ=21)
                                                                           gTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaac              >  2:3486446/1‑93 (MQ=25)
                                                                            tttACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaaca             >  2:1277090/1‑93 (MQ=25)
                                                                            tttACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaaca             >  1:4723611/1‑93 (MQ=25)
                                                                            tttACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaaca             >  1:4671316/1‑93 (MQ=25)
                                                                            tttACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaaca             >  1:3383028/1‑93 (MQ=21)
                                                                            tttACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaaca             >  1:1966575/1‑93 (MQ=25)
                                                                              tACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaacaac           <  1:2080754/93‑1 (MQ=25)
                                                                              tACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaacaac           <  1:1342687/93‑1 (MQ=25)
                                                                                cgccgcATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGa         <  1:3485777/93‑1 (MQ=25)
                                                                                 gccgcATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAg        >  2:5975257/1‑93 (MQ=25)
                                                                                  ccgcATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGc       >  1:3058832/1‑93 (MQ=37)
                                                                                   cgcATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCt      >  2:2573437/1‑93 (MQ=37)
                                                                                    gcATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCtt     >  1:3933943/1‑93 (MQ=37)
                                                                                    gcATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCtt     >  2:1940487/1‑93 (MQ=37)
                                                                                     cATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTc    <  2:233120/93‑1 (MQ=37)
                                                                                      aTCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGCCTTTCTCCGCATCCGGCATAAACAACGAGCTTCa   <  2:3361270/93‑1 (MQ=25)
                                                                                       tCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAg  <  2:236444/93‑1 (MQ=37)
                                                                                            |                                                                                        
GCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCCATCTGAGCATGATGCCTGATGCGGCGTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAG  >  NC_000913/4051051‑4051230

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCCATCTGAGCATGATGCCTGATGCGGCG‑‑TGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAAC  >  NC_000913/4051051‑4051236
                                                                                             |                                                                                              
GCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCA                                                                                                >  SRR4302532.1064140/1‑94 (MQ=60)
GCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGACGCATCAGGCA                                                                                                >  SRR4302532.2618725/1‑94 (MQ=60)
    AGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTG                                                                                           <  SRR4302532.728226/94‑1 (MQ=60)
    AGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTG                                                                                           >  SRR4302532.1272992/1‑94 (MQ=60)
    AGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTG                                                                                           >  SRR4302532.2904600/1‑94 (MQ=60)
     GAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGA                                                                                          <  SRR4302532.4216954/94‑1 (MQ=60)
     GAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGA                                                                                          <  SRR4302532.1429872/94‑1 (MQ=60)
     GAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGA                                                                                          <  SRR4302532.1429873/94‑1 (MQ=60)
      AACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAG                                                                                         <  SRR4302532.2618827/94‑1 (MQ=60)
      AACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAG                                                                                         <  SRR4302532.2023226/94‑1 (MQ=60)
        CAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCA                                                                                       <  SRR4302532.5221540/94‑1 (MQ=60)
         AAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCAT                                                                                      >  SRR4302532.1111399/1‑94 (MQ=60)
          AGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATG                                                                                     >  SRR4302532.3851765/1‑94 (MQ=60)
           GATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGA                                                                                    <  SRR4302532.2408254/94‑1 (MQ=60)
            ATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGAT                                                                                   <  SRR4302532.1297938/94‑1 (MQ=60)
            ATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGAT                                                                                   <  SRR4302532.318287/94‑1 (MQ=60)
            ATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGAT                                                                                   <  SRR4302532.5079134/94‑1 (MQ=60)
            ATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGAT                                                                                   <  SRR4302532.6196946/94‑1 (MQ=60)
             TTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATG                                                                                  <  SRR4302532.144093/94‑1 (MQ=60)
             TTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATG                                                                                  <  SRR4302532.6204213/94‑1 (MQ=60)
               TCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCC                                                                                <  SRR4302532.3978157/94‑1 (MQ=60)
               GCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCC                                                                                >  SRR4302532.4495488/1‑94 (MQ=60)
                CCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCT                                                                               >  SRR4302532.1124045/1‑94 (MQ=60)
                 CTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTG                                                                              >  SRR4302532.4073237/1‑94 (MQ=60)
                  TGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGA                                                                             >  SRR4302532.2961247/1‑94 (MQ=60)
                  TGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGA                                                                             >  SRR4302532.986849/1‑94 (MQ=60)
                  TGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGAAGGATAAGGCGTTTACGCCGCATCAGACATCTGAGAATGATGTCTGA                                                                             >  SRR4302532.6011105/1‑94 (MQ=60)
                  TGATGAGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGA                                                                             >  SRR4302532.4172440/1‑94 (MQ=60)
                   GATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGAT                                                                            <  SRR4302532.2698775/94‑1 (MQ=60)
                   GATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGAT                                                                            <  SRR4302532.2363633/94‑1 (MQ=60)
                   GATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGAT                                                                            <  SRR4302532.963296/94‑1 (MQ=60)
                   GATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCGGGCATCTGAGCATGATGCCTGAT                                                                            >  SRR4302532.3599643/1‑94 (MQ=60)
                    ATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATG                                                                           <  SRR4302532.5080592/94‑1 (MQ=60)
                    ATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATG                                                                           <  SRR4302532.1287875/94‑1 (MQ=60)
                       CGCTGGCCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCTTCAGGCATCTGAGCATGATGCCTGATGCGG                                                                        <  SRR4302532.2762515/94‑1 (MQ=60)
                         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.2113127/1‑94 (MQ=60)
                                ACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCG                                                              >  SRR4302532.5386803/1‑94 (MQ=60)
                                ACAAAACCGGTTGAAATTATTGCTTTTGTAGACCGGATAAGGCGTTTACGCCGCATCAGGCATCTGCGCATGATGCCTGATGCGGCGCTGGTCG                                                              >  SRR4302532.2142089/1‑94 (MQ=60)
                                 CAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGC                                                             >  SRR4302532.2289966/1‑94 (MQ=60)
                                 CAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGC                                                             <  SRR4302532.3762150/94‑1 (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.5129028/94‑1 (MQ=60)
                                           TGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAG                                                   <  SRR4302532.4415761/94‑1 (MQ=60)
                                             AAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGC                                                 >  SRR4302532.3783412/1‑94 (MQ=60)
                                             AAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGC                                                 <  SRR4302532.1911414/94‑1 (MQ=60)
                                             AAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGC                                                 >  SRR4302532.4372839/1‑94 (MQ=60)
                                              AATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCC                                                <  SRR4302532.258979/94‑1 (MQ=60)
                                              AATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCC                                                >  SRR4302532.3197236/1‑94 (MQ=60)
                                               ATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCT                                               <  SRR4302532.3779887/94‑1 (MQ=60)
                                               ATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCT                                               <  SRR4302532.5099618/94‑1 (MQ=60)
                                               ATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCT                                               <  SRR4302532.4767639/94‑1 (MQ=60)
                                                TTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTC                                              >  SRR4302532.3315298/1‑94 (MQ=60)
                                                  cttcCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCAGGCATCTGAGCATGATGCCTGAGGCCGCGCTGGTCGCGTCTTATCAGGCCTACG                                            <  SRR4302532.4422190/90‑1 (MQ=49)
                                                   TTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGA                                           <  SRR4302532.3235844/94‑1 (MQ=60)
                                                   TTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGA                                           <  SRR4302532.470544/94‑1 (MQ=60)
                                                   TTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGA                                           <  SRR4302532.784452/94‑1 (MQ=60)
                                                    TGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTCCGAC                                          <  SRR4302532.6198479/94‑1 (MQ=60)
                                                    TGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGAC                                          <  SRR4302532.2593728/94‑1 (MQ=60)
                                                    TGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGAC                                          <  SRR4302532.6361436/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.2427905/94‑1 (MQ=60)
                                                           GTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCC                                   <  SRR4302532.5050778/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.2031896/1‑94 (MQ=60)
                                                                           GTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACA                   >  SRR4302532.4830585/1‑94 (MQ=60)
                                                                           GTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACA                   >  SRR4302532.4884781/1‑94 (MQ=60)
                                                                           GTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACA                   >  SRR4302532.5119317/1‑94 (MQ=60)
                                                                             TTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAAC                 <  SRR4302532.3421076/94‑1 (MQ=60)
                                                                             TTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAAC                 >  SRR4302532.4000247/1‑94 (MQ=60)
                                                                              TACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACG                <  SRR4302532.2150119/94‑1 (MQ=60)
                                                                              TACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACG                <  SRR4302532.1387577/94‑1 (MQ=60)
                                                                              TACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACG                <  SRR4302532.5118013/94‑1 (MQ=60)
                                                                               ACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGA               <  SRR4302532.1615901/94‑1 (MQ=60)
                                                                                CGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAG              <  SRR4302532.3602556/94‑1 (MQ=60)
                                                                                CGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAG              >  SRR4302532.2616211/1‑94 (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)
                                                                                             |                                                                                              
GCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCCATCTGAGCATGATGCCTGATGCGGCG‑‑TGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAAC  >  NC_000913/4051051‑4051236

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