Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A1 F1 I2 R1 1873 183.1 9376148 98.2% 9207377 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 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.92.9% 432.7 / 23.3 113coding (613/990 nt)araFL‑arabinose‑binding periplasmic protein
Reads supporting (aligned to +/- strand):  ref base C (2/6);  new base . (52/52);  total (54/59)
Fisher's exact test for biased strand distribution p-value = 2.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.77e-01

Only 100 of 118 total aligned reads displayed.

GCCCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTA  >  NC_000913/1983450‑1983627
                                                                                         |                                                                                        
gccCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGcccc                                                                                       <  1:2118751/93‑1 (MQ=255)
  ccAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGACGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGccccgg                                                                                     <  2:3586631/93‑3 (MQ=255)
   cAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCggg                                                                                    >  2:4416765/1‑93 (MQ=255)
   cAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCggg                                                                                    >  1:817278/1‑93 (MQ=255)
    aGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCggga                                                                                   <  1:1121165/93‑2 (MQ=255)
    aGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCggga                                                                                   <  2:757287/93‑2 (MQ=255)
    aGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCggga                                                                                   <  1:2451138/93‑2 (MQ=255)
    aGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCggga                                                                                   <  1:2745552/93‑2 (MQ=255)
    aGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCggga                                                                                   <  2:3564244/93‑2 (MQ=255)
    aGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCggga                                                                                   <  1:501266/93‑2 (MQ=255)
         ggTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCg                                                                             >  2:3696971/1‑93 (MQ=255)
          gTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGt                                                                            <  2:2395383/93‑1 (MQ=255)
            gCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTa                                                                          <  1:2326473/93‑1 (MQ=255)
            gCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTa                                                                          >  2:3515262/1‑93 (MQ=255)
             cTGTCGTTCATACCGACGATCAGCCCATGTTTAACTTCCGGATGTTGAACCAGCATTGCGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAg                                                                         >  2:721606/1‑93 (MQ=255)
             cTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAg                                                                         >  1:2961070/1‑93 (MQ=255)
               gTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGAt                                                                       >  1:1511968/1‑93 (MQ=255)
                tCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGAtt                                                                      <  1:48034/93‑1 (MQ=255)
                tCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGAtt                                                                      <  1:800817/93‑1 (MQ=255)
                    tCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTggt                                                                  >  1:1966886/1‑93 (MQ=255)
                    tCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTggt                                                                  >  1:4593101/1‑93 (MQ=255)
                    tCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTggt                                                                  <  1:1241962/93‑1 (MQ=255)
                     cATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTggta                                                                 >  2:1695959/1‑93 (MQ=255)
                      aTACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTggtag                                                                >  1:3983934/1‑93 (MQ=255)
                      aTACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTggtag                                                                >  1:3099703/1‑93 (MQ=255)
                       tACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTggtagg                                                               <  1:4243323/93‑1 (MQ=255)
                       tACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTggtagg                                                               <  1:4130585/93‑1 (MQ=255)
                        aCCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTggtaggt                                                              >  1:4602546/1‑93 (MQ=255)
                        aCCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTggtaggt                                                              >  2:2650775/1‑93 (MQ=255)
                        aCCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTggtaggt                                                              >  2:1497212/1‑93 (MQ=255)
                            acgaTCAGCCAATGTTTAACTTCCGTATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCt                                                          <  1:4425349/93‑1 (MQ=255)
                              gctCAGCCAATGTTTAACTTCCGGCTGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGa                                                        <  1:2100776/91‑1 (MQ=255)
                              gaTCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGa                                                        <  2:4659009/93‑1 (MQ=255)
                               aTCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGAt                                                       <  1:3127871/93‑1 (MQ=255)
                               aTCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGAt                                                       >  1:379282/1‑93 (MQ=255)
                                tCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATa                                                      <  2:538426/93‑1 (MQ=255)
                                tCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATa                                                      <  1:936991/93‑1 (MQ=255)
                                tCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATa                                                      <  2:2867761/93‑1 (MQ=255)
                                  aGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATaaa                                                    <  1:1945221/93‑1 (MQ=255)
                                  aGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATaaa                                                    >  1:1873172/1‑93 (MQ=255)
                                  aGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATaaa                                                    >  1:4379808/1‑93 (MQ=255)
                                  aGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATaaa                                                    >  2:2699198/1‑93 (MQ=255)
                                   gCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAAt                                                   >  1:308159/1‑93 (MQ=255)
                                    ccAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAAtt                                                  <  2:1890130/93‑1 (MQ=255)
                                    ccAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAAtt                                                  <  1:2061534/93‑1 (MQ=255)
                                     cAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAAttt                                                 <  1:98327/93‑1 (MQ=255)
                                      aaTGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTg                                                <  1:2590939/93‑1 (MQ=255)
                                      aaTGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGAATGTCGTTAGATTTGGTAGGTACCTGATAAATTTg                                                <  2:3050730/93‑1 (MQ=255)
                                       aTGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGt                                               >  2:884951/1‑93 (MQ=255)
                                        tGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGtt                                              <  1:2120261/93‑1 (MQ=255)
                                        tGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGtt                                              >  1:566526/1‑93 (MQ=255)
                                        tGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGtt                                              >  2:1018492/1‑93 (MQ=255)
                                         gTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGttt                                             <  1:2765232/93‑1 (MQ=255)
                                         gTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGttt                                             <  1:198131/93‑1 (MQ=255)
                                         gTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGttt                                             <  2:959168/93‑1 (MQ=255)
                                           ttAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGttttt                                           <  1:166478/93‑1 (MQ=255)
                                           ttAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGttttt                                           <  2:546627/93‑1 (MQ=255)
                                           ttAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGttttt                                           <  2:2683204/93‑1 (MQ=255)
                                            tAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTc                                          <  1:1895625/93‑1 (MQ=255)
                                            tAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTc                                          >  1:909371/1‑93 (MQ=255)
                                              aCTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCg                                        >  2:2767812/1‑93 (MQ=255)
                                              aCTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCg                                        <  2:4167521/93‑1 (MQ=255)
                                              aCTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTACg                                        >  1:2174592/1‑93 (MQ=255)
                                               cTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCgg                                       <  2:851360/93‑1 (MQ=255)
                                                ttCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCgcg                                      >  1:454700/1‑91 (MQ=255)
                                                 tCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGa                                     <  2:4155820/93‑1 (MQ=255)
                                                  ccGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGaa                                    <  2:4371008/93‑1 (MQ=255)
                                                    ggATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATc                                  >  2:543631/1‑93 (MQ=255)
                                                    ggATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATc                                  <  1:473280/93‑1 (MQ=255)
                                                     gATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAAttc                                 >  1:2165121/1‑91 (MQ=255)
                                                     gATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATcc                                 >  2:2413835/1‑93 (MQ=255)
                                                      aTGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCg                                <  2:2096664/93‑1 (MQ=255)
                                                       tGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCgg                               >  1:1954929/1‑93 (MQ=255)
                                                             aCCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCtt                         >  2:103599/1‑93 (MQ=255)
                                                             aCCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCtt                         >  1:259972/1‑93 (MQ=255)
                                                             aCCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCtt                         >  1:304190/1‑93 (MQ=255)
                                                              ccAGCATTGAGTTGGCAGCGTCAAATGCCCTGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCttt                        <  1:1106789/93‑1 (MQ=255)
                                                              ccAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCttt                        <  2:1230579/93‑1 (MQ=255)
                                                               cAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTc                       >  2:2510093/1‑93 (MQ=255)
                                                               cAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTc                       <  1:4099816/93‑1 (MQ=255)
                                                                aGCATTGAGTTGGCAGCGTCAAATGCGCCGGGATGTCGTTAGAGTTGGTAGGTACCTGATAAATTTGTTTTGCCGGGAATCCGGCCGCTTTCa                      <  2:1130494/93‑1 (MQ=255)
                                                                aGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCa                      <  1:4183760/93‑1 (MQ=255)
                                                                aGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCa                      <  1:3091173/93‑1 (MQ=255)
                                                                aGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCa                      >  2:3420148/1‑93 (MQ=255)
                                                                 ggaTTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGGACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAg                     >  2:2449421/3‑93 (MQ=255)
                                                                 gCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAg                     >  2:1953745/1‑93 (MQ=255)
                                                                 gCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAg                     <  1:2668252/93‑1 (MQ=255)
                                                                  cATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAgc                    >  1:1482699/1‑93 (MQ=255)
                                                                   aTTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAgcg                   >  1:471392/1‑93 (MQ=255)
                                                                    ttGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAgcga                  >  2:1187832/1‑92 (MQ=255)
                                                                        gTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATcc              <  2:2632082/93‑1 (MQ=255)
                                                                         ttGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCa             >  1:3267293/1‑93 (MQ=255)
                                                                         ttGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCa             >  1:289767/1‑93 (MQ=255)
                                                                          tGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCAt            >  2:4660643/1‑93 (MQ=255)
                                                                            gCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAg          >  2:3170864/1‑93 (MQ=255)
                                                                               gCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATc       >  2:1700001/1‑93 (MQ=255)
                                                                                cGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATcc      <  2:1926536/93‑1 (MQ=255)
                                                                                 gTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATccc     >  2:505083/1‑93 (MQ=255)
                                                                                    aaaTGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCgta  <  1:2109222/93‑1 (MQ=255)
                                                                                    aaaTGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCgta  >  2:111530/1‑93 (MQ=255)
                                                                                         |                                                                                        
GCCCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTA  >  NC_000913/1983450‑1983627

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GCCCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTAGTACGGCG  >  NC_000913/1983450‑1983635
                                                                                             |                                                                                            
GCCCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCG                                                                                              <  SRR4302535.2134034/94‑1 (MQ=60)
  CCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGG                                                                                            >  SRR4302535.823012/1‑94 (MQ=60)
    AGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGAT                                                                                         <  SRR4302535.1129160/94‑1 (MQ=60)
    AGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGAT                                                                                         <  SRR4302535.2765121/94‑1 (MQ=60)
    AGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGAT                                                                                         <  SRR4302535.2468804/94‑1 (MQ=60)
    AGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGAT                                                                                         <  SRR4302535.504776/94‑1 (MQ=60)
       ACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTC                                                                                      <  SRR4302535.2687318/94‑1 (MQ=60)
          GTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTT                                                                                   >  SRR4302535.1775199/1‑94 (MQ=60)
          GTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTT                                                                                   <  SRR4302535.595029/94‑1 (MQ=60)
            GCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAG                                                                                 >  SRR4302535.2981855/1‑94 (MQ=60)
            GCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAG                                                                                 <  SRR4302535.3755307/94‑1 (MQ=60)
            GCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAG                                                                                 <  SRR4302535.2343119/94‑1 (MQ=60)
              TGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGAT                                                                               >  SRR4302535.1522859/1‑94 (MQ=60)
                TCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTT                                                                             <  SRR4302535.48394/94‑1 (MQ=60)
                TCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTT                                                                             <  SRR4302535.806445/94‑1 (MQ=60)
                 CGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTG                                                                            >  SRR4302535.4248928/1‑94 (MQ=60)
                   TTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGT                                                                          >  SRR4302535.4624365/1‑94 (MQ=60)
                   TTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGT                                                                          >  SRR4302535.1981064/1‑94 (MQ=60)
                    TCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTA                                                                         <  SRR4302535.1250832/94‑1 (MQ=60)
                     CATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAG                                                                        >  SRR4302535.3121286/1‑94 (MQ=60)
                     CATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAG                                                                        >  SRR4302535.4011156/1‑94 (MQ=60)
                       TACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGT                                                                      <  SRR4302535.4272364/94‑1 (MQ=60)
                       TACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGT                                                                      >  SRR4302535.4633845/1‑94 (MQ=60)
                       TACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGT                                                                      <  SRR4302535.4158846/94‑1 (MQ=60)
                            ACGATCAGCCAATGTTTAACTTCCGTATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTG                                                                 <  SRR4302535.4455660/94‑1 (MQ=60)
                              GCTCAGCCAATGTTTAACTTCCGGCTGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGAT                                                               <  SRR4302535.2115924/94‑1 (MQ=60)
                              GATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGAT                                                               >  SRR4302535.381908/1‑94 (MQ=60)
                               ATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATA                                                              <  SRR4302535.3149654/94‑1 (MQ=60)
                                TCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAA                                                             <  SRR4302535.943627/94‑1 (MQ=60)
                                 CAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAA                                                            >  SRR4302535.1886672/1‑94 (MQ=60)
                                 CAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAA                                                            >  SRR4302535.4409774/1‑94 (MQ=60)
                                  AGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAAT                                                           <  SRR4302535.1959224/94‑1 (MQ=60)
                                  AGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAAT                                                           >  SRR4302535.310294/1‑94 (MQ=60)
                                    CCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTT                                                         <  SRR4302535.2076389/94‑1 (MQ=60)
                                     CAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTG                                                        <  SRR4302535.99021/94‑1 (MQ=60)
                                      AATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGT                                                       >  SRR4302535.2173006/1‑94 (MQ=60)
                                      AATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGT                                                       <  SRR4302535.2609541/94‑1 (MQ=60)
                                       ATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTT                                                      >  SRR4302535.570490/1‑94 (MQ=60)
                                       ATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATTCCTCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTT                                                      <  SRR4302535.4642173/94‑1 (MQ=60)
                                        TGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTT                                                     <  SRR4302535.2135553/94‑1 (MQ=60)
                                         GTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTT                                                    <  SRR4302535.2784927/94‑1 (MQ=60)
                                         GTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTT                                                    <  SRR4302535.199491/94‑1 (MQ=60)
                                           TTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTC                                                  >  SRR4302535.915830/1‑94 (MQ=60)
                                           TTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTC                                                  <  SRR4302535.167636/94‑1 (MQ=60)
                                            TAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCC                                                 <  SRR4302535.1909265/94‑1 (MQ=60)
                                             AACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTACG                                                >  SRR4302535.2190143/1‑94 (MQ=60)
                                               CTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGCG                                              >  SRR4302535.457885/1‑94 (MQ=60)
                                                    GGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATTC                                         >  SRR4302535.2180647/1‑94 (MQ=60)
                                                    GGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCC                                         <  SRR4302535.476592/94‑1 (MQ=60)
                                                      ATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGG                                       >  SRR4302535.1968999/1‑94 (MQ=60)
                                                            AACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTT                                 >  SRR4302535.261772/1‑94 (MQ=60)
                                                            AACCAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTT                                 >  SRR4302535.306302/1‑94 (MQ=60)
                                                              CCAGCATTGAGTTGGCAGCGTCAAATGCCCTGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTC                               <  SRR4302535.1114694/94‑1 (MQ=60)
                                                               CAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTaggggtggggggggtgtgtttggttttttggggggggggggggtgtgg                              >  SRR4302535.647637/1‑46 (MQ=60)
                                                               CAGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCA                              <  SRR4302535.4127886/94‑1 (MQ=60)
                                                                AGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAG                             <  SRR4302535.4212387/94‑1 (MQ=60)
                                                                AGCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAG                             <  SRR4302535.3112700/94‑1 (MQ=60)
                                                                 GCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGC                            <  SRR4302535.2687352/94‑1 (MQ=60)
                                                                 GCATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGC                            >  SRR4302535.1493302/1‑94 (MQ=60)
                                                                  CATTGAGTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCG                           >  SRR4302535.474691/1‑94 (MQ=60)
                                                                        GTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCA                     >  SRR4302535.3289989/1‑94 (MQ=60)
                                                                        GTTGGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCA                     >  SRR4302535.291795/1‑94 (MQ=60)
                                                                           GGCAGCGTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAG                  <  SRR4302535.203786/94‑1 (MQ=60)
                                                                                 GTCAAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCG            >  SRR4302535.306279/1‑94 (MQ=60)
                                                                                    AAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTAG         <  SRR4302535.2124433/94‑1 (MQ=60)
                                                                                    AAATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTAG         >  SRR4302535.197970/1‑94 (MQ=60)
                                                                                     AATGCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTAGT        <  SRR4302535.3542861/94‑1 (MQ=60)
                                                                                     AATGCCCCGGGATGTCGGTAGATTTGGTAGGTGCCTGATAAATTTGGTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTAGT        <  SRR4302535.3542860/94‑1 (MQ=60)
                                                                                         GCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTAGTACG     >  SRR4302535.4292235/1‑94 (MQ=60)
                                                                                         GCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTAGTACG     <  SRR4302535.2396738/94‑1 (MQ=60)
                                                                                          CCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTAGTACGG    <  SRR4302535.2832153/94‑1 (MQ=60)
                                                                                           CCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTAGTACGGC   >  SRR4302535.382281/1‑94 (MQ=60)
                                                                                            CCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTAGTACGGCG  <  SRR4302535.461207/94‑1 (MQ=60)
                                                                                            CCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTAGTACGGCG  >  SRR4302535.2766739/1‑94 (MQ=60)
                                                                                             |                                                                                            
GCCCAGCACGGTGCTGTCGTTCATACCGACGATCAGCCAATGTTTAACTTCCGGATGTTGAACCAGCATTGAGTTGGCAGCGTCAAATGCCCCCGGGATGTCGTTAGATTTGGTAGGTACCTGATAAATTTGTTTTTCCGGGAATCCGGCCGCTTTCAGCGCATCCATAGATCCCGTAGTACGGCG  >  NC_000913/1983450‑1983635

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