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: