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: