Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A1 F1 I2 R1 1475 266.7 13469680 98.2% 13227225 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.6% 832.8 / ‑3.5 249intergenic (+145/+306)hemN/glnGO2‑independent coproporphyrinogen III oxidase/GlnG
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base . (122/126);  minor base T (0/1);  total (122/127)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 281 total aligned reads displayed.

GCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCCATCTGAGCATGATGCCTGATGCGGCGTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAG  >  NC_000913/4051051‑4051230
                                                                                            |                                                                                        
gCTGAGAACGAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGc                                                                                          <  1:6226796/93‑1 (MQ=255)
gCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGc                                                                                          <  2:1012273/93‑1 (MQ=255)
gCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGc                                                                                          <  1:5996119/93‑1 (MQ=255)
 cTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGca                                                                                         <  2:623287/93‑2 (MQ=255)
 cTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGGTTACGCCGCATCAGGca                                                                                         <  2:77570/93‑2 (MQ=255)
  tGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGcat                                                                                        <  2:4406103/93‑3 (MQ=255)
  tGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGcat                                                                                        >  1:5077156/1‑91 (MQ=255)
   gagaACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGcatc                                                                                       <  1:4361517/93‑4 (MQ=255)
   gagaACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGcatc                                                                                       <  1:5714675/93‑4 (MQ=255)
   gagaACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGcatc                                                                                       >  1:5268063/1‑90 (MQ=255)
      aaCAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGa                                                                                   >  1:1384029/1‑93 (MQ=255)
         aaGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCa                                                                                >  1:3285903/1‑93 (MQ=255)
         aaGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCa                                                                                >  2:4967535/1‑93 (MQ=255)
         aaGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCa                                                                                <  1:1952227/93‑1 (MQ=255)
          aGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCat                                                                               >  2:90256/1‑93 (MQ=255)
          aGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGACGAGGCGTTTACGACGCATCAGGCATCTGAGCat                                                                               >  2:3908446/1‑93 (MQ=37)
           gATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCatg                                                                              <  2:4959769/93‑1 (MQ=255)
           gATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCatg                                                                              >  1:1083733/1‑93 (MQ=255)
           gATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCatg                                                                              >  1:2698629/1‑93 (MQ=255)
            aTTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCatga                                                                             <  2:3520400/93‑1 (MQ=255)
            aTTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCatga                                                                             <  1:6355090/93‑1 (MQ=255)
             ttGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCatgat                                                                            >  1:4971106/1‑93 (MQ=255)
              tGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCatgatg                                                                           >  1:4078844/1‑93 (MQ=255)
                ccTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGcc                                                                         <  2:3421591/93‑1 (MQ=255)
                 cTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCt                                                                        >  1:2450545/1‑93 (MQ=255)
                 cTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCt                                                                        <  2:329196/93‑1 (MQ=255)
                 cTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCt                                                                        <  1:4855002/93‑1 (MQ=255)
                   gATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGa                                                                      <  1:1693955/93‑1 (MQ=255)
                   gATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGa                                                                      <  2:3743908/93‑1 (MQ=255)
                    aTGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCGTCAGGCATCTGAGCATGATGCCTGAt                                                                     <  2:1549958/93‑1 (MQ=255)
                    aTGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGAt                                                                     <  1:4745458/93‑1 (MQ=255)
                    aTGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGAt                                                                     <  1:5164768/93‑1 (MQ=255)
                    aTGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGACGCATCAGGCATCTGAGCATGATGCCTGAt                                                                     >  2:2090130/1‑93 (MQ=255)
                     tGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATg                                                                    <  2:5423806/93‑1 (MQ=255)
                     tGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGACGCATCAGGCATCTGAGCATGATGCCTGATg                                                                    >  2:143199/1‑93 (MQ=255)
                        gcTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATgcgg                                                                 <  2:3699623/93‑1 (MQ=255)
                          tGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATgcggcg                                                               >  1:2547055/1‑93 (MQ=255)
                          tGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATgcggcg                                                               >  2:3875163/1‑93 (MQ=255)
                                acaaAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTc                                                         >  2:5618688/1‑93 (MQ=21)
                                  aaaaCCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTcgc                                                       <  2:1981885/93‑1 (MQ=25)
                                   aaaCCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTcgcg                                                      >  2:2308186/1‑93 (MQ=25)
                                    aaCCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGACGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGt                                                     >  2:473452/1‑93 (MQ=21)
                                     aCCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTc                                                    <  2:2924870/93‑1 (MQ=25)
                                     aCCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTc                                                    <  2:2157429/93‑1 (MQ=25)
                                     aCCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTc                                                    <  1:3370670/93‑1 (MQ=25)
                                       cTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCtt                                                  <  1:5410478/93‑1 (MQ=21)
                                        tGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTa                                                 <  2:2247404/93‑1 (MQ=21)
                                         gTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTAt                                                <  2:4103077/93‑1 (MQ=21)
                                         gTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGATGCGCTGGTCGCGTCTTAt                                                >  1:4889004/1‑93 (MQ=21)
                                          ttGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATc                                               <  2:451296/93‑1 (MQ=21)
                                          ttGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATc                                               <  2:5782813/93‑1 (MQ=21)
                                           tGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCa                                              >  1:5878441/1‑93 (MQ=21)
                                           tGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCa                                              <  2:4773179/93‑1 (MQ=18)
                                            gAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAg                                             >  1:1572648/1‑93 (MQ=18)
                                            gAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAg                                             <  1:6551692/93‑1 (MQ=21)
                                              aaTTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGACGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGc                                           >  2:2299237/1‑93 (MQ=18)
                                               attattGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGcc                                          <  2:3612329/93‑1 (MQ=18)
                                                ttattGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCt                                         <  2:3145059/93‑1 (MQ=18)
                                                  attGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTAc                                       >  2:5964501/1‑93 (MQ=18)
                                                  attGCTTTTGTAGGCCGGATAAGGCCTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTAc                                       <  1:903962/93‑1 (MQ=17)
                                                   ttGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACg                                      >  2:188892/1‑93 (MQ=17)
                                                    tGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGa                                     <  2:2717518/93‑1 (MQ=17)
                                                     gCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGAc                                    <  1:3927514/93‑1 (MQ=17)
                                                      cTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACt                                   <  1:4484562/93‑1 (MQ=18)
                                                       ttttGTAGGCCGGATAAGGCGTTTACGGCGCATCAAGCATCTGAGCATGATGCCTGATTCGTCGCTGTTGGCGTCTTATCATGCCTACGACtt                                  >  2:26831/1‑93 (MQ=14)
                                                       ttttGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACtt                                  <  1:3121496/93‑1 (MQ=18)
                                                       ttttGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACtt                                  >  2:5340626/1‑93 (MQ=17)
                                                       ttttGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGGCGCGTCTTATCAGGCCTACGACtt                                  >  1:4395925/1‑93 (MQ=17)
                                                            tAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTcc                             >  2:5817656/1‑93 (MQ=18)
                                                            tAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTcc                             <  1:4641193/93‑1 (MQ=18)
                                                             aGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCg                            >  2:6713857/1‑93 (MQ=18)
                                                             aGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCg                            <  2:5040382/93‑1 (MQ=18)
                                                             aGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCg                            >  1:5612642/1‑93 (MQ=18)
                                                             aGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCg                            >  1:2592568/1‑93 (MQ=18)
                                                              ggCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGc                           >  2:2500277/1‑93 (MQ=18)
                                                              ggCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGc                           >  1:4438051/1‑93 (MQ=18)
                                                              ggCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGc                           >  2:1829810/1‑93 (MQ=18)
                                                              ggCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGc                           >  2:3817686/1‑93 (MQ=18)
                                                               gCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACtt                                  <  1:2237042/85‑1 (MQ=17)
                                                               gCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCa                          >  1:2422637/1‑93 (MQ=18)
                                                                 cGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATc                        >  1:426161/1‑93 (MQ=18)
                                                                   gATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCg                      >  2:1742815/1‑93 (MQ=18)
                                                                     tAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGc                    >  1:2918521/1‑93 (MQ=21)
                                                                       aGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCAt                  >  1:206781/1‑93 (MQ=21)
                                                                       aGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCAt                  <  1:1256246/93‑1 (MQ=21)
                                                                       aGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCAt                  >  1:3164872/1‑93 (MQ=21)
                                                                         gCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATaa                >  1:3667804/1‑93 (MQ=21)
                                                                            tttACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGAGCTGGTCGCGTCTTATCAGGACTACGACTTTCTCCGCATACGGCATAaaca             >  2:1065862/1‑93 (MQ=25)
                                                                             ttACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaacaa            <  1:1037087/93‑1 (MQ=21)
                                                                             ttACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaacaa            >  2:5094021/1‑93 (MQ=25)
                                                                             ttACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaacaa            <  1:2788122/93‑1 (MQ=21)
                                                                                cgccgcATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGa         <  1:18060/93‑1 (MQ=25)
                                                                                cgccgcATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGa         <  2:4918675/93‑1 (MQ=25)
                                                                                cgccgcATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGa         <  2:4649554/93‑1 (MQ=25)
                                                                                cgccgcATCAGTCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGCCTTTCTCCGCATCCGGCATAAACAACGa         <  2:6039688/93‑1 (MQ=21)
                                                                                 gccgcATCAGGCATCTGAGCATGATTCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTCCGACTTTCTCCGCATCCGGCATAAACAACGAg        <  2:6130599/93‑1 (MQ=25)
                                                                                     cATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTc    >  2:6544872/1‑93 (MQ=37)
                                                                                     cATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTc    <  1:828115/93‑1 (MQ=37)
                                                                                      aTCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCa   >  2:6331663/1‑93 (MQ=37)
                                                                                       tCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAg  >  2:2860072/1‑93 (MQ=37)
                                                                                            |                                                                                        
GCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCCATCTGAGCATGATGCCTGATGCGGCGTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAG  >  NC_000913/4051051‑4051230

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 28 ≤ 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
GCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCCATCTGAGCATGATGCCTGATGCGGCGTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAA  >  NC_000913/4051051‑4051235
                                                                                             |                                                                                            
GCTGAGAACGAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCA                                                                                              <  SRR4302536.6412324/94‑1 (MQ=60)
GCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCA                                                                                              <  SRR4302536.6174863/94‑1 (MQ=60)
  TGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATC                                                                                           >  SRR4302536.5425074/1‑94 (MQ=60)
   GAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCT                                                                                          <  SRR4302536.4491893/94‑1 (MQ=60)
   GAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCT                                                                                          <  SRR4302536.5884621/94‑1 (MQ=60)
   GAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCT                                                                                          <  SRR4302536.5190689/94‑1 (MQ=60)
   GAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCT                                                                                          <  SRR4302536.331349/94‑1 (MQ=60)
     GAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGA                                                                                        >  SRR4302536.1425251/1‑94 (MQ=60)
      AACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAG                                                                                       >  SRR4302536.3257443/1‑94 (MQ=60)
       ACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGC                                                                                      <  SRR4302536.6669655/94‑1 (MQ=60)
        CAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCA                                                                                     >  SRR4302536.2828003/1‑94 (MQ=60)
        CAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCA                                                                                     >  SRR4302536.3382557/1‑94 (MQ=60)
         AAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCAT                                                                                    >  SRR4302536.5091078/1‑94 (MQ=60)
         AAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCAT                                                                                    <  SRR4302536.4732693/94‑1 (MQ=60)
         AAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCAT                                                                                    <  SRR4302536.2010941/94‑1 (MQ=60)
          AGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATG                                                                                   >  SRR4302536.2779157/1‑94 (MQ=60)
          AGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATG                                                                                   >  SRR4302536.1116196/1‑94 (MQ=60)
            ATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGAT                                                                                 >  SRR4302536.5120427/1‑94 (MQ=60)
            ATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGAT                                                                                 <  SRR4302536.6544417/94‑1 (MQ=60)
             TTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATG                                                                                <  SRR4302536.2560395/94‑1 (MQ=60)
             TTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATG                                                                                >  SRR4302536.4199741/1‑94 (MQ=60)
             TTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATG                                                                                <  SRR4302536.4443632/94‑1 (MQ=60)
             TTGCCTGATGAGCTGGGCAACAAATCCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATG                                                                                <  SRR4302536.3382710/94‑1 (MQ=60)
               GCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCC                                                                              >  SRR4302536.5560440/1‑94 (MQ=60)
               GCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCC                                                                              >  SRR4302536.1489818/1‑94 (MQ=60)
                CCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCT                                                                             <  SRR4302536.2724217/94‑1 (MQ=60)
                CCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCT                                                                             >  SRR4302536.2524576/1‑94 (MQ=60)
                 CTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTG                                                                            <  SRR4302536.5001104/94‑1 (MQ=60)
                   GATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGAT                                                                          <  SRR4302536.6202505/94‑1 (MQ=60)
                   GATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGAT                                                                          <  SRR4302536.1742993/94‑1 (MQ=60)
                    ATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATG                                                                         <  SRR4302536.5318866/94‑1 (MQ=60)
                    ATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATG                                                                         <  SRR4302536.238066/94‑1 (MQ=60)
                    ATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATG                                                                         <  SRR4302536.4888273/94‑1 (MQ=60)
                    ATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATG                                                                         >  SRR4302536.4310159/1‑94 (MQ=60)
                    ATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATG                                                                         <  SRR4302536.2643993/94‑1 (MQ=60)
                    ATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATG                                                                         <  SRR4302536.123614/94‑1 (MQ=60)
                        GCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGC                                                                     <  SRR4302536.4126905/94‑1 (MQ=60)
                         CTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCG                                                                    <  SRR4302536.2318643/94‑1 (MQ=60)
                         CTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCG                                                                    <  SRR4302536.767251/94‑1 (MQ=60)
                         CTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCG                                                                    >  SRR4302536.2623770/1‑94 (MQ=60)
                         CTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCG                                                                    >  SRR4302536.2226656/1‑94 (MQ=60)
                          TGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGC                                                                  <  SRR4302536.4592399/94‑1 (MQ=60)
                              CAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGT                                                               <  SRR4302536.33847/94‑1 (MQ=60)
                              CAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGACGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGT                                                               >  SRR4302536.2222022/1‑94 (MQ=60)
                               AACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGctagtc                                                              >  SRR4302536.3103443/1‑88 (MQ=60)
                                ACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCG                                                             >  SRR4302536.2027906/1‑94 (MQ=60)
                                ACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCG                                                             <  SRR4302536.6436451/94‑1 (MQ=60)
                                 CAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGGCGCATCACGAATCTGAGCATGATGCCTGATGCGGggctgggcgg                                                            >  SRR4302536.936889/1‑84 (MQ=60)
                                  AAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCG                                                           >  SRR4302536.1413735/1‑94 (MQ=60)
                                   AAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGT                                                          >  SRR4302536.2961189/1‑94 (MQ=60)
                                    AACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTC                                                         >  SRR4302536.5105612/1‑94 (MQ=60)
                                    AACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTC                                                         >  SRR4302536.5841475/1‑94 (MQ=60)
                                     ACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCT                                                        <  SRR4302536.3470803/94‑1 (MQ=60)
                                     ACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCT                                                        <  SRR4302536.1078537/94‑1 (MQ=60)
                                     ACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCT                                                        >  SRR4302536.3150009/1‑94 (MQ=60)
                                     ACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCT                                                        >  SRR4302536.904124/1‑94 (MQ=60)
                                       CTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTA                                                      <  SRR4302536.5572147/94‑1 (MQ=60)
                                       CTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTA                                                      <  SRR4302536.5644530/94‑1 (MQ=60)
                                       CTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTA                                                      <  SRR4302536.2519778/94‑1 (MQ=60)
                                       CTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTA                                                      <  SRR4302536.6113298/94‑1 (MQ=60)
                                       CTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGAAGCATCAGGCATCTGAGCATGATGCCTGATGCGGAGCTGGTCGCGTCTTA                                                      >  SRR4302536.5381971/1‑94 (MQ=60)
                                        TGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTAT                                                     >  SRR4302536.480126/1‑94 (MQ=60)
                                        TGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTAT                                                     <  SRR4302536.6894516/94‑1 (MQ=60)
                                        TGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGATGCGCTGGTCGCGTCTTAT                                                     >  SRR4302536.5036057/1‑94 (MQ=60)
                                         GTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATC                                                    >  SRR4302536.263847/1‑94 (MQ=60)
                                          TTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCA                                                   >  SRR4302536.6053543/1‑94 (MQ=60)
                                           TGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAG                                                  >  SRR4302536.1618761/1‑94 (MQ=60)
                                            GAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGG                                                 <  SRR4302536.6746803/94‑1 (MQ=60)
                                             AAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGC                                                >  SRR4302536.5638399/1‑94 (MQ=60)
                                              AATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCC                                               >  SRR4302536.2834505/1‑94 (MQ=60)
                                               ATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCT                                              <  SRR4302536.3867501/94‑1 (MQ=60)
                                                 TATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTAC                                            >  SRR4302536.4967571/1‑94 (MQ=60)
                                                  ATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACG                                           <  SRR4302536.1595249/94‑1 (MQ=60)
                                                  ATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACG                                           >  SRR4302536.3631891/1‑94 (MQ=60)
                                                  ATTGCTTTTGTAGGCCGGATAAGGCCTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACG                                           <  SRR4302536.930599/94‑1 (MQ=60)
                                                   TTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGA                                          <  SRR4302536.1830678/94‑1 (MQ=60)
                                                   TTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGA                                          <  SRR4302536.3408313/94‑1 (MQ=60)
                                                     GCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATTAGGCCTACGACT                                        >  SRR4302536.4619093/1‑94 (MQ=60)
                                                     GCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACT                                        <  SRR4302536.4044219/94‑1 (MQ=60)
                                                      CTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTT                                       >  SRR4302536.5277161/1‑94 (MQ=60)
                                                      CTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTT                                       <  SRR4302536.4619255/94‑1 (MQ=60)
                                                      CTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGGCGCGTCTTATCAGGCCTACGACTT                                       >  SRR4302536.4527257/1‑94 (MQ=60)
                                                       ctaaggtGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTT                                      <  SRR4302536.2304716/87‑1 (MQ=60)
                                                       ctaaggtGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTT                                      <  SRR4302536.5149598/87‑1 (MQ=60)
                                                       TTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTT                                      <  SRR4302536.3212684/94‑1 (MQ=60)
                                                       TTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTT                                      <  SRR4302536.1411838/94‑1 (MQ=60)
                                                       TTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTT                                      <  SRR4302536.883640/94‑1 (MQ=60)
                                                         tggggggcgGGATAGGGCGTTTCCGCCGCATCAGCCATCTGAGCATGATGCTTGATGCGGCGCTGGTCGCGTTTTATCAGGCCTACGACTTTCT                                    <  SRR4302536.179440/85‑1 (MQ=60)
                                                           GTAGGCCGAATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCC                                  <  SRR4302536.1683711/94‑1 (MQ=60)
                                                            TAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCG                                 >  SRR4302536.446605/1‑94 (MQ=60)
                                                            TAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCG                                 <  SRR4302536.4780696/94‑1 (MQ=60)
                                                            TAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCG                                 >  SRR4302536.5779378/1‑94 (MQ=60)
                                                            TAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCG                                 >  SRR4302536.2240192/1‑94 (MQ=60)
                                                            TAGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCG                                 >  SRR4302536.2670404/1‑94 (MQ=60)
                                                            TAGGCCGGATAAGGCGTTTACGACGCATCAGGCATCTGAGAATGATGCCTGATGCGGAGCTGGTCCCGTCATATCAGGCCTACGACTTTttccc                                 >  SRR4302536.224641/1‑89 (MQ=41)
                                                             AGGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGC                                >  SRR4302536.4571245/1‑94 (MQ=60)
                                                              GGCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCA                               >  SRR4302536.2495889/1‑94 (MQ=60)
                                                               GCCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCAT                              <  SRR4302536.2094087/94‑1 (MQ=60)
                                                                CCGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATC                             >  SRR4302536.437382/1‑94 (MQ=60)
                                                                 CGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCC                            >  SRR4302536.4981458/1‑94 (MQ=60)
                                                                 CGGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCC                            >  SRR4302536.1969223/1‑94 (MQ=60)
                                                                  GGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCG                           <  SRR4302536.6310763/94‑1 (MQ=60)
                                                                  GGATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCG                           <  SRR4302536.4552761/94‑1 (MQ=60)
                                                                    ATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGC                         >  SRR4302536.1697445/1‑94 (MQ=60)
                                                                    ATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGC                         >  SRR4302536.529743/1‑94 (MQ=60)
                                                                    ATAAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGC                         >  SRR4302536.3004663/1‑94 (MQ=60)
                                                                      AAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCAT                       >  SRR4302536.5006599/1‑94 (MQ=60)
                                                                      AAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCAT                       >  SRR4302536.3257586/1‑94 (MQ=60)
                                                                      AAGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCAT                       >  SRR4302536.212236/1‑94 (MQ=60)
                                                                       AGGCTTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATA                      <  SRR4302536.5982715/94‑1 (MQ=60)
                                                                       AGGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATA                      <  SRR4302536.1293944/94‑1 (MQ=60)
                                                                        GGCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAA                     >  SRR4302536.3777456/1‑94 (MQ=60)
                                                                          CGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAAC                   <  SRR4302536.991211/94‑1 (MQ=60)
                                                                          CGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAAC                   >  SRR4302536.4653938/1‑94 (MQ=60)
                                                                          CGTTTACGACGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAAC                   >  SRR4302536.3314914/1‑94 (MQ=60)
                                                                           GTTTGCGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACA                  <  SRR4302536.3678211/94‑1 (MQ=60)
                                                                             TTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAAC                <  SRR4302536.1067868/94‑1 (MQ=60)
                                                                             TTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAAC                >  SRR4302536.273494/1‑94 (MQ=60)
                                                                             TTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAAC                <  SRR4302536.2870909/94‑1 (MQ=60)
                                                                              TGCGCCGCATCAGGCATCTGAGCATGATGCCTGAGGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACG               <  SRR4302536.1192583/94‑1 (MQ=60)
                                                                              TACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACG               <  SRR4302536.448352/94‑1 (MQ=60)
                                                                              TACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACG               >  SRR4302536.4054290/1‑94 (MQ=60)
                                                                              TACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCCGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACG               <  SRR4302536.5998984/94‑1 (MQ=60)
                                                                                CGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAG             >  SRR4302536.3095733/1‑94 (MQ=60)
                                                                                CGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAG             <  SRR4302536.18559/94‑1 (MQ=60)
                                                                                 GCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGC            >  SRR4302536.6664522/1‑94 (MQ=60)
                                                                                 GCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGC            <  SRR4302536.2617615/94‑1 (MQ=60)
                                                                                  CCGCATCAGGCATCTGAGCATGATGCCTGATGCGGGGCTGGTCGCGGCTTATCAGGCCTACGACTTTTTTCGCATCCGGCATAAACAACGAGCT           >  SRR4302536.380490/1‑94 (MQ=60)
                                                                                  CCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCT           >  SRR4302536.4523590/1‑94 (MQ=60)
                                                                                   CGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTT          >  SRR4302536.3247710/1‑94 (MQ=60)
                                                                                   CGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTT          <  SRR4302536.4875950/94‑1 (MQ=60)
                                                                                   CGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTT          <  SRR4302536.1867921/94‑1 (MQ=60)
                                                                                   CGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGCCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTT          <  SRR4302536.5272819/94‑1 (MQ=60)
                                                                                    GCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTC         <  SRR4302536.6172393/94‑1 (MQ=60)
                                                                                    GCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTC         >  SRR4302536.5838128/1‑94 (MQ=60)
                                                                                     CATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCA        <  SRR4302536.852492/94‑1 (MQ=60)
                                                                                      ATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAG       >  SRR4302536.4387671/1‑94 (MQ=60)
                                                                                        CAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCT     <  SRR4302536.6549248/94‑1 (MQ=60)
                                                                                        CAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCT     >  SRR4302536.737934/1‑94 (MQ=60)
                                                                                         AGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTA    >  SRR4302536.5681800/1‑94 (MQ=60)
                                                                                         AGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTA    <  SRR4302536.2217629/94‑1 (MQ=60)
                                                                                           GGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAA   <  SRR4302536.2265276/94‑1 (MQ=60)
                                                                                            GCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAA  >  SRR4302536.2335526/1‑94 (MQ=60)
                                                                                             |                                                                                            
GCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCAGGCCATCTGAGCATGATGCCTGATGCGGCGTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAA  >  NC_000913/4051051‑4051235

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