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 1,983,543 (C)5→4 coding (609/990 nt) araF ← L‑arabinose‑binding periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,983,5390C.97.2% 951.7 / 13.4 218coding (613/990 nt)araFL‑arabinose‑binding periplasmic protein
Reads supporting (aligned to +/- strand):  ref base C (1/5);  new base . (110/102);  total (111/107)
Fisher's exact test for biased strand distribution p-value = 1.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.82e-01

Only 100 of 232 total aligned reads displayed.

CCGCCCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTA  >  NC_000913/1983448‑1983627
                                                                                           |                                                                                        
 cgccCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGccc                                                                                        <  1:2031299/93‑1 (MQ=255)
 cgccCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGccc                                                                                        <  1:2533588/93‑1 (MQ=255)
  gccCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGcccc                                                                                       >  2:2029910/1‑93 (MQ=255)
   cccAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGccccg                                                                                      <  2:62029/93‑2 (MQ=255)
    ccAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGccccgg                                                                                     >  2:3794291/1‑91 (MQ=255)
      aGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCggga                                                                                   <  2:702669/93‑2 (MQ=255)
      aGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCggga                                                                                   <  1:2336463/93‑2 (MQ=255)
      aGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCggga                                                                                   >  2:2178910/1‑92 (MQ=255)
       gCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGAt                                                                                 <  1:5400789/93‑1 (MQ=255)
       gCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGAt                                                                                 <  2:2578724/93‑1 (MQ=255)
       gCACGGGGCTGTCTTACATACCGACGATCAGCGAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGAt                                                                                 <  1:5267445/93‑1 (MQ=255)
          cGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTc                                                                              >  2:733304/1‑93 (MQ=255)
          cGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTc                                                                              <  2:3980240/93‑1 (MQ=255)
           ggTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCg                                                                             <  2:1389989/93‑1 (MQ=255)
            gTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGt                                                                            <  2:2689980/93‑1 (MQ=255)
              gCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTa                                                                          <  1:1140004/93‑1 (MQ=255)
               cTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAg                                                                         >  1:2154452/1‑93 (MQ=255)
                tGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGa                                                                        >  2:1225763/1‑93 (MQ=255)
                 gTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGAt                                                                       <  1:6435194/93‑1 (MQ=255)
                 gTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGAt                                                                       >  1:19353/1‑93 (MQ=255)
                  tCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGAtt                                                                      >  2:2219388/1‑93 (MQ=255)
                    gTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTg                                                                    >  2:6335101/1‑93 (MQ=255)
                     ttCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTgg                                                                   <  2:3330472/93‑1 (MQ=255)
                     ttCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTgg                                                                   >  1:4530919/1‑93 (MQ=255)
                      tCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTggt                                                                  <  1:3581299/93‑1 (MQ=255)
                       cATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTggta                                                                 <  2:4357982/93‑1 (MQ=255)
                       cATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTggta                                                                 >  1:4923984/1‑93 (MQ=255)
                       cATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTggta                                                                 <  2:388314/93‑1 (MQ=255)
                        aTACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTggtag                                                                <  2:5760570/93‑1 (MQ=255)
                        aTACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTggtag                                                                <  1:1557664/93‑1 (MQ=255)
                         tACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTggtagg                                                               >  2:5078131/1‑93 (MQ=255)
                          aCCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTggtaggt                                                              <  2:5665792/93‑1 (MQ=255)
                          aCCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGACCCGGGATGTCGTTAGATTTggtaggt                                                              >  1:5367640/1‑93 (MQ=255)
                           ccGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTggtaggta                                                             >  1:1466663/1‑93 (MQ=255)
                             gacgaTCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTAcc                                                           <  2:1291917/93‑1 (MQ=255)
                              acgaTCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCt                                                          >  2:4711414/1‑93 (MQ=255)
                              acgaTCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCt                                                          >  1:3079744/1‑93 (MQ=255)
                               cgaTCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTg                                                         >  1:3698178/1‑93 (MQ=255)
                                gaTCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGa                                                        <  2:3177618/93‑1 (MQ=255)
                                  tCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATa                                                      >  1:5345270/1‑93 (MQ=255)
                                      ccAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAAtt                                                  <  1:3030683/93‑1 (MQ=255)
                                       cAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAAttt                                                 >  2:2977867/1‑93 (MQ=255)
                                        aaTGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTg                                                >  2:265391/1‑93 (MQ=255)
                                         aTGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGt                                               <  2:772374/93‑1 (MQ=255)
                                         aTGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGt                                               <  2:3640744/93‑1 (MQ=255)
                                         aTGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGt                                               <  1:2654387/93‑1 (MQ=255)
                                         aTGTGTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGt                                               <  1:2438213/93‑1 (MQ=255)
                                           gTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGttt                                             >  2:6627627/1‑93 (MQ=255)
                                            tttAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGtttt                                            <  2:4249911/93‑1 (MQ=255)
                                            tttAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGtttt                                            <  2:1927154/93‑1 (MQ=255)
                                             ttAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGttttt                                           >  2:5053894/1‑93 (MQ=255)
                                             ttAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGttttt                                           >  2:3240871/1‑93 (MQ=255)
                                               aaCTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTcc                                         >  1:333979/1‑93 (MQ=255)
                                               aaCTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTcc                                         <  1:3954584/93‑1 (MQ=255)
                                                aCTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCg                                        >  2:1350159/1‑93 (MQ=255)
                                                aCTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCg                                        <  1:5898867/93‑1 (MQ=255)
                                                 cTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCgg                                       <  1:6472280/93‑1 (MQ=255)
                                                   tCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGa                                     <  2:4923835/93‑1 (MQ=255)
                                                   tCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGa                                     <  1:3559695/93‑1 (MQ=255)
                                                    ccGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGaa                                    <  2:16150/93‑1 (MQ=255)
                                                    ccGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGaa                                    <  2:2946427/93‑1 (MQ=255)
                                                      ggATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATc                                  <  2:5245598/93‑1 (MQ=255)
                                                        aTGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTTCGGGAATCCg                                >  1:5626552/1‑93 (MQ=255)
                                                        aTGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCg                                >  2:549265/1‑93 (MQ=255)
                                                         tGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCgg                               <  2:590047/93‑1 (MQ=255)
                                                         tGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCgg                               <  2:4536068/93‑1 (MQ=255)
                                                          gTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGc                              <  2:3128538/93‑1 (MQ=255)
                                                          gTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGc                              <  2:1144090/93‑1 (MQ=255)
                                                           ttGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGcc                             >  2:790029/1‑93 (MQ=255)
                                                           ttGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGcc                             >  1:5697497/1‑93 (MQ=255)
                                                             gAACCAGCATTGATTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGc                           <  2:6709214/93‑1 (MQ=255)
                                                             gAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGc                           >  1:6581008/1‑93 (MQ=255)
                                                              gaCCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCt                          <  2:2288831/92‑1 (MQ=255)
                                                              aaCCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCt                          >  2:3424641/1‑93 (MQ=255)
                                                              aaCCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCt                          <  1:3414164/93‑1 (MQ=255)
                                                               aCCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCtt                         <  1:6516414/93‑1 (MQ=255)
                                                               aCCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCtt                         >  1:3776714/1‑93 (MQ=255)
                                                                  aGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCa                      >  1:1792232/1‑93 (MQ=255)
                                                                   gCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAg                     >  2:4680075/1‑93 (MQ=255)
                                                                   gCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAg                     >  1:5806026/1‑93 (MQ=255)
                                                                   gCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAg                     >  2:6164224/1‑93 (MQ=255)
                                                                   gCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCAGCTTTCAg                     >  1:1644359/1‑93 (MQ=255)
                                                                    cATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAgc                    <  1:1962227/93‑1 (MQ=255)
                                                                     aTTGAGTTGGCAGCGTCAAATGCCCCGGGGTGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTTCGGGAATCCGGCCGCTTTCAgcg                   >  2:6428348/1‑93 (MQ=255)
                                                                      ttGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAgcgc                  <  1:2301521/93‑1 (MQ=255)
                                                                      ttGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAgcgc                  >  2:1839528/1‑93 (MQ=255)
                                                                         acttGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCAGGAATCCGGCCGCTTTCAGCGCATc               >  2:6484162/3‑93 (MQ=255)
                                                                           ttGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCa             <  1:692618/93‑1 (MQ=255)
                                                                            tGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCAt            <  2:5611750/93‑1 (MQ=255)
                                                                             ggCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATa           <  1:99077/93‑1 (MQ=255)
                                                                             ggCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATa           >  1:2918867/1‑93 (MQ=255)
                                                                             ggCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATa           <  2:5425136/93‑1 (MQ=255)
                                                                             ggCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATa           >  1:5055999/1‑93 (MQ=255)
                                                                               cAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGa         <  1:2162373/93‑1 (MQ=255)
                                                                               cAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGa         <  2:4444626/93‑1 (MQ=255)
                                                                                 gCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTCATAACTTTGTTTTTCCGGGAATCGGGCCGCTTTCAGCGCATCCATAGATc       >  2:3560407/1‑93 (MQ=255)
                                                                                   gTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATccc     >  2:3504796/1‑93 (MQ=255)
                                                                                   gTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATccc     <  1:609170/93‑1 (MQ=255)
                                                                                    tCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCg    >  1:94055/1‑93 (MQ=255)
                                                                                      aaaTGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCACgta  >  1:5865003/1‑93 (MQ=255)
                                                                                           |                                                                                        
CCGCCCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTA  >  NC_000913/1983448‑1983627

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
CCCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTAGTACGGCG  >  NC_000913/1983451‑1983635
                                                                                            |                                                                                            
CCCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGG                                                                                             >  SRR4302536.3671677/1‑94 (MQ=60)
  CAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGA                                                                                          <  SRR4302536.1012438/94‑1 (MQ=60)
   AGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGAT                                                                                         <  SRR4302536.2407350/94‑1 (MQ=60)
   AGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGAT                                                                                         <  SRR4302536.2302782/94‑1 (MQ=60)
   AGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGAT                                                                                         >  SRR4302536.350124/1‑94 (MQ=60)
    GCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATG                                                                                        <  SRR4302536.5562114/94‑1 (MQ=60)
    GCACGGGGCTGTCTTACATACCGACGATCAGCGAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATG                                                                                        <  SRR4302536.5424433/94‑1 (MQ=60)
      CCGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTC                                                                                      <  SRR4302536.5520245/94‑1 (MQ=60)
      ACGGTGCTGTCGTTCGTACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTC                                                                                      <  SRR4302536.5657205/94‑1 (MQ=60)
      ACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGTCCCGGGATGTC                                                                                      <  SRR4302536.4853141/94‑1 (MQ=60)
      ACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTC                                                                                      <  SRR4302536.6712872/94‑1 (MQ=60)
        GGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGT                                                                                    >  SRR4302536.5117823/1‑94 (MQ=60)
        GGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGT                                                                                    >  SRR4302536.3465103/1‑94 (MQ=60)
         GTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTT                                                                                   <  SRR4302536.2666951/94‑1 (MQ=60)
         GTGCTGGCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTT                                                                                   <  SRR4302536.4552408/94‑1 (MQ=60)
           GCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAG                                                                                 <  SRR4302536.1174283/94‑1 (MQ=60)
           GCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAG                                                                                 >  SRR4302536.2219190/1‑94 (MQ=60)
             TGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGAT                                                                               <  SRR4302536.1555164/94‑1 (MQ=60)
             TGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGAT                                                                               <  SRR4302536.2226039/94‑1 (MQ=60)
             TGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGAT                                                                               >  SRR4302536.19896/1‑94 (MQ=60)
             TGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGGTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGAT                                                                               >  SRR4302536.1806674/1‑94 (MQ=60)
              GTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATT                                                                              <  SRR4302536.6627522/94‑1 (MQ=60)
                CGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTG                                                                            <  SRR4302536.4602583/94‑1 (MQ=60)
                 GTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGG                                                                           >  SRR4302536.4666738/1‑94 (MQ=60)
                   TCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTA                                                                         >  SRR4302536.5932890/1‑94 (MQ=60)
                   TCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTA                                                                         >  SRR4302536.5071943/1‑94 (MQ=60)
                   TCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTA                                                                         <  SRR4302536.3688384/94‑1 (MQ=60)
                    CATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAG                                                                        <  SRR4302536.2392218/94‑1 (MQ=60)
                    CATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAG                                                                        >  SRR4302536.2045932/1‑94 (MQ=60)
                    CATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAG                                                                        <  SRR4302536.605313/94‑1 (MQ=60)
                     ATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGG                                                                       >  SRR4302536.1701141/1‑94 (MQ=60)
                     ATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGG                                                                       <  SRR4302536.1603413/94‑1 (MQ=60)
                      TACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGACCCGGGATGTCGTTAGATTTGGTAGGT                                                                      >  SRR4302536.5527707/1‑94 (MQ=60)
                       ACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTA                                                                     >  SRR4302536.1510052/1‑94 (MQ=60)
                         CGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACC                                                                   <  SRR4302536.629526/94‑1 (MQ=60)
                          GACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCT                                                                  >  SRR4302536.3169889/1‑94 (MQ=60)
                          GACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCT                                                                  >  SRR4302536.3406245/1‑94 (MQ=60)
                           ACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCGGCGTCAAATGCCCCGGGATGTCGTGAGATTTGGTAGGTACCTG                                                                 <  SRR4302536.2008432/94‑1 (MQ=60)
                           ACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTG                                                                 >  SRR4302536.3808760/1‑94 (MQ=60)
                             GATCAGCCAATGTTTAACTTCCGGATGTTGAACCCGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGAT                                                               >  SRR4302536.5994949/1‑94 (MQ=60)
                              ATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATA                                                              >  SRR4302536.5504897/1‑94 (MQ=60)
                                  GCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATT                                                          >  SRR4302536.5497279/1‑94 (MQ=60)
                                   CCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTTATAAATTT                                                         >  SRR4302536.1987518/1‑94 (MQ=60)
                                   CCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTT                                                         <  SRR4302536.3119738/94‑1 (MQ=60)
                                   CCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTT                                                         >  SRR4302536.2580295/1‑94 (MQ=60)
                                    CAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTG                                                        >  SRR4302536.3799572/1‑94 (MQ=60)
                                     AATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGT                                                       >  SRR4302536.5577455/1‑94 (MQ=60)
                                     AATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGT                                                       <  SRR4302536.2068099/94‑1 (MQ=60)
                                      ATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTT                                                      >  SRR4302536.5797977/1‑94 (MQ=60)
                                      ATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTT                                                      <  SRR4302536.2733860/94‑1 (MQ=60)
                                      ATGTGTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTT                                                      <  SRR4302536.2511916/94‑1 (MQ=60)
                                        GTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTT                                                    <  SRR4302536.2976243/94‑1 (MQ=60)
                                         TTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTT                                                   >  SRR4302536.3609891/1‑94 (MQ=60)
                                          TTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTC                                                  >  SRR4302536.2573782/1‑94 (MQ=60)
                                           TAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCC                                                 >  SRR4302536.342842/1‑94 (MQ=60)
                                            AACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCG                                                <  SRR4302536.4072014/94‑1 (MQ=60)
                                             ACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGG                                               <  SRR4302536.6074506/94‑1 (MQ=60)
                                              CTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGG                                              <  SRR4302536.6665436/94‑1 (MQ=60)
                                              CTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGG                                              >  SRR4302536.2474487/1‑94 (MQ=60)
                                              CTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGG                                              <  SRR4302536.3122177/94‑1 (MQ=60)
                                              CTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGG                                              >  SRR4302536.3885574/1‑94 (MQ=60)
                                                TCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAA                                            <  SRR4302536.3666128/94‑1 (MQ=60)
                                                   TGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCC                                         <  SRR4302536.6653622/94‑1 (MQ=60)
                                                   GGATGTTGAACCGGCATTGAGTTGGCAGCGTCAAATGCGCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCC                                         >  SRR4302536.6709150/1‑94 (MQ=60)
                                                   GGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCC                                         >  SRR4302536.6214335/1‑94 (MQ=60)
                                                    GATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTTCGGGAATCCG                                        >  SRR4302536.5793686/1‑94 (MQ=60)
                                                      TGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGC                                      >  SRR4302536.1316613/1‑94 (MQ=60)
                                                       GTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCC                                     >  SRR4302536.5866908/1‑94 (MQ=60)
                                                       GTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCC                                     <  SRR4302536.1646208/94‑1 (MQ=60)
                                                         TGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGC                                   >  SRR4302536.6776830/1‑94 (MQ=60)
                                                           GACCAACATTGAGTTGGCAGCCTCAAATGCCCCGGGATGTGGTTAGGTTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTT                                 <  SRR4302536.5540407/94‑1 (MQ=60)
                                                           AACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTT                                 >  SRR4302536.3889439/1‑94 (MQ=60)
                                                           AACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTT                                 <  SRR4302536.3515750/94‑1 (MQ=60)
                                                            ACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTT                                <  SRR4302536.6710571/94‑1 (MQ=60)
                                                             CCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTC                               <  SRR4302536.6757274/94‑1 (MQ=60)
                                                              CAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCA                              >  SRR4302536.6008215/1‑94 (MQ=60)
                                                              CAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCA                              >  SRR4302536.1846358/1‑94 (MQ=60)
                                                               AGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAG                             >  SRR4302536.3104165/1‑94 (MQ=60)
                                                               AGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAG                             >  SRR4302536.5978943/1‑94 (MQ=60)
                                                               AGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCAGCTTTCAG                             >  SRR4302536.1692188/1‑94 (MQ=60)
                                                                 CATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCG                           <  SRR4302536.2021226/94‑1 (MQ=60)
                                                                   TTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCA                         <  SRR4302536.2371347/94‑1 (MQ=60)
                                                                    TGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCAT                        >  SRR4302536.2332289/1‑94 (MQ=60)
                                                                     GAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATC                       <  SRR4302536.5234358/94‑1 (MQ=60)
                                                                        TTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCAT                    <  SRR4302536.712856/94‑1 (MQ=60)
                                                                         TGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATA                   >  SRR4302536.5207454/1‑94 (MQ=60)
                                                                         TGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATA                   >  SRR4302536.3005013/1‑94 (MQ=60)
                                                                          GGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAG                  <  SRR4302536.6335945/94‑1 (MQ=60)
                                                                          GGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAG                  <  SRR4302536.101713/94‑1 (MQ=60)
                                                                            CAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGAT                <  SRR4302536.2227415/94‑1 (MQ=60)
                                                                            CAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGAT                <  SRR4302536.3744277/94‑1 (MQ=60)
                                                                            CAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGAT                >  SRR4302536.934108/1‑94 (MQ=60)
                                                                             AGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATC               >  SRR4302536.686615/1‑94 (MQ=60)
                                                                               CGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCC             <  SRR4302536.2145555/94‑1 (MQ=60)
                                                                               CGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCC             >  SRR4302536.3965461/1‑94 (MQ=60)
                                                                                GTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCG            >  SRR4302536.96556/1‑94 (MQ=60)
                                                                                GTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCG            <  SRR4302536.627373/94‑1 (MQ=60)
                                                                                  CAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTA          <  SRR4302536.1009857/94‑1 (MQ=60)
                                                                                  CAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCACGTA          >  SRR4302536.6039763/1‑94 (MQ=60)
                                                                                    AATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTAGT        >  SRR4302536.3940489/1‑94 (MQ=60)
                                                                                     ATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTAGTA       <  SRR4302536.6864998/94‑1 (MQ=60)
                                                                                        gcgCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTGCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTAGTACG     <  SRR4302536.4330689/91‑1 (MQ=60)
                                                                                        GCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTAGTACG     <  SRR4302536.6614183/94‑1 (MQ=60)
                                                                                        GCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTAGTACG     <  SRR4302536.698395/94‑1 (MQ=60)
                                                                                        GCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTAGTACG     <  SRR4302536.835403/94‑1 (MQ=60)
                                                                                          CCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTAGTACGGC   >  SRR4302536.4088295/1‑94 (MQ=60)
                                                                                           CCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTAGTACGGCG  <  SRR4302536.5699457/94‑1 (MQ=60)
                                                                                           CCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTAGTACGGCG  <  SRR4302536.4459444/94‑1 (MQ=60)
                                                                                            |                                                                                            
CCCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTAGTACGGCG  >  NC_000913/1983451‑1983635

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