Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,931,725 A→G N108D (AAT→GAT)  rbsA → ATP‑binding component of D‑ribose high‑affinity transport system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,931,7250AG82.9% 696.2 / 146.1 322N108D (AAT→GAT) rbsAATP‑binding component of D‑ribose high‑affinity transport system
Reads supporting (aligned to +/- strand):  ref base A (27/28);  new base G (135/132);  total (162/160)
Fisher's exact test for biased strand distribution p-value = 8.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.89e-01

Only 100 of 325 total aligned reads displayed.

GAAGCCGGGATTGGGATTATCCATCAGGAACTGAACCTGATCCCGCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTAATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCGGATAAATTGCTGGCTAAACTTAACCTGCGCTTTAAAAGCGACAAGCT  >  NC_000913/3931635‑3931816
                                                                                          |                                                                                           
gAAGCCGGGATTGGGATTATCCATCAGGAACTGAACCTGATCCCGCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGAt                                                                                           <  1:4983444/93‑1 (MQ=255)
gAAGCCGGGATTGGGATTATCCATCAGGAACTGAACCTGATCCCGCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGAt                                                                                           <  2:7470226/93‑1 (MQ=255)
gAAGCCGGGATTGGGATTATCCATCAGGAACTGAACCTGATCCCGCAGTTGACCATTGCCAAAAACATTTTCCTCGGTCGTGAGTTTGTTGAt                                                                                           <  1:1063853/93‑1 (MQ=255)
 aaGCCGGGATTGGGATTATCCATCAGGAACTGAACCTGATCCCGCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATc                                                                                          <  1:2877439/93‑1 (MQ=255)
 aaGCCGGGATTGGGATTATCCATCAGGAACTGAACCTGATCCCGCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATc                                                                                          <  1:6989676/93‑1 (MQ=255)
  aGCCGGGATTGGGATTATCCATCAGGAACTGAACCTGATCCCGCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCg                                                                                         >  1:3901012/1‑93 (MQ=255)
   gCCGGGATTGGGATTATCCATCAGGAACTGAACCTGATCCCGCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGc                                                                                        >  1:309529/1‑93 (MQ=255)
     cGGGATTGGGATTATCCATCAGGAACTGAACCTGATCCCGCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCtt                                                                                      <  2:1859197/93‑1 (MQ=255)
      gggattgggattATCCATCAGGAACTGAACCTGATCCCGCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCttt                                                                                     <  2:516280/93‑1 (MQ=255)
      gggattgggattATCCATCAGGAACTGAACCTGATCCCGCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCttt                                                                                     <  1:5301535/93‑1 (MQ=255)
      gggattgggattATCCATCAGGAACTGAACCTGATCCCGCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCttt                                                                                     <  2:2362163/93‑1 (MQ=255)
       ggattgggattATCCATCAGGAACTGAACCTGATCCCGCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTg                                                                                    <  1:4531689/93‑1 (MQ=255)
        gattgggattATCCATCAGGAACTGAACCTGATCCCGCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTgg                                                                                   >  1:5765894/1‑93 (MQ=255)
        gattgggattATCCATCAGGAACTGAACCTGATCCCGCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTgg                                                                                   <  2:3340733/93‑1 (MQ=255)
        gattgggattATCCATCAGGAACTGAACCTGATCCCGCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTgg                                                                                   <  2:5312606/93‑1 (MQ=255)
        gattgggattATCCATCAGGAACTGAACCTGATCCCGCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTAATCGCTTTgg                                                                                   <  1:6740358/93‑1 (MQ=255)
         attgggattATCCATCAGGAACTGAACCTGATCCCGCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTGGc                                                                                  <  1:3170825/93‑1 (MQ=255)
          ttgggattATCCATCAGGAACTGAACCTGATCCCGCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCa                                                                                 >  2:3080774/1‑93 (MQ=255)
            gggattATCCATCAGGAACTGAACCTGATCCCGCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCaaa                                                                               <  1:3267743/93‑1 (MQ=255)
             ggattATCCATCAGGAACTGAACCTGATCCCGCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCaaaa                                                                              >  1:7471187/1‑93 (MQ=255)
             ggattATCCATCAGGAACTGAACCTGATCCCGCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCaaaa                                                                              <  2:337397/93‑1 (MQ=255)
             ggattATCCATCAGGAACTGAACCTGATCCCGCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCaaaa                                                                              <  2:3539348/93‑1 (MQ=255)
              gattATCCATCAGGAACTGAACCTGATCCCGCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAAt                                                                             <  2:1318676/93‑1 (MQ=255)
              gattATCCATCAGGAACTGAACCTGATCCCGCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAAt                                                                             >  1:4676374/1‑93 (MQ=255)
                ttATCCATCAGGAACTGAACCTGATCCCGCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTg                                                                           >  1:2692245/1‑93 (MQ=255)
                 tATCCATCAGGAACTGAACCTGATCCCGCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGa                                                                          <  1:3857805/93‑1 (MQ=255)
                 tATCCATCAGGAACTGAACCTGATCCCGCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGa                                                                          <  2:4791947/93‑1 (MQ=255)
                  aTCCATCAGGAACTGAACCTGATCCCGCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTAATCGCTTTGGCAAAATTGAc                                                                         <  1:6478629/93‑1 (MQ=255)
                   tCCATCAGGAACTGAACCTGATCCCGCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACt                                                                        >  2:810332/1‑93 (MQ=255)
                    ccATCAGGAACTGAACCTGATCCCGCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTg                                                                       >  2:7439541/1‑93 (MQ=255)
                     cATCAGGAACTGAACCTGATCCCGCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTgg                                                                      <  2:6297283/93‑1 (MQ=255)
                     cATCAGGAACTGAACCTGATCCCGCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTgg                                                                      >  1:6843107/1‑93 (MQ=255)
                       tCAGGAACTGAACCTGATCCCGCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGaa                                                                    >  2:256973/1‑93 (MQ=255)
                         aGGAACTGAACCTGATCCCGCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGaaaa                                                                  >  1:5629474/1‑93 (MQ=255)
                           gAACTGAACCTGATCCCGCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGaaaann                                                                >  1:3390788/1‑91 (MQ=255)
                            aaCTGAACCTGATCCCGCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCa                                                               >  1:5943776/1‑93 (MQ=255)
                             aCTGAACCTGATCCCGCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCAt                                                              >  1:4622702/1‑93 (MQ=255)
                             aCTGAACCTGATCCCGCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCAt                                                              >  1:887749/1‑93 (MQ=255)
                                gAACCTGATCCCGCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCATGTa                                                           <  1:4843661/93‑1 (MQ=255)
                                gAACCTGATCCCGCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCATGTa                                                           <  2:6199047/93‑1 (MQ=255)
                                gAACCTGATCCCGCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCATGTa                                                           <  2:6943188/93‑1 (MQ=255)
                                gAACCTGATCCCGCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTAATCGCTTTGGCAAAATTGACTGGAAAACCATGTa                                                           <  1:3774022/93‑1 (MQ=255)
                                 aaCCTGATCCCGCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCATGTAt                                                          >  2:2267916/1‑93 (MQ=255)
                                 aaCCTGATCCCGCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCATGTAt                                                          <  2:3651520/93‑1 (MQ=255)
                                    cTGATCCCGCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGcc                                                       >  2:5434496/1‑93 (MQ=255)
                                      gATCCCGCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGa                                                     <  2:5093037/93‑1 (MQ=255)
                                       aTCCCGCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGaa                                                    <  2:5677269/93‑1 (MQ=255)
                                       aTCCCGCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTAATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGaa                                                    <  1:4260430/93‑1 (MQ=255)
                                        tCCCGCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAg                                                   <  1:2163911/93‑1 (MQ=255)
                                          ccGCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCg                                                 >  2:6120015/1‑93 (MQ=255)
                                          ccGCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTAATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCg                                                 >  1:5516639/1‑93 (MQ=255)
                                            gCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCGGa                                               <  2:4459690/93‑1 (MQ=255)
                                            gCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTAATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCGGa                                               <  2:382024/93‑1 (MQ=255)
                                              aGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCGGATa                                             >  2:4343154/1‑93 (MQ=255)
                                               gTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCGGATaa                                            >  1:2541210/1‑93 (MQ=255)
                                               gTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCGGATaa                                            >  2:2631078/1‑93 (MQ=255)
                                                ttGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCGGATaaa                                           <  1:2308943/93‑1 (MQ=255)
                                                 tGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCGGATAAAt                                          >  1:7365532/1‑93 (MQ=255)
                                                  gACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCGGATAAAtt                                         >  2:7377557/1‑93 (MQ=255)
                                                  gACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCGGATAAAtt                                         <  1:1527598/93‑1 (MQ=255)
                                                  gACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCGGATAAAtt                                         >  1:6142879/1‑93 (MQ=255)
                                                  gACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTAATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCGGATAAAtt                                         <  2:6757988/93‑1 (MQ=255)
                                                   aCCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCGGATAAATTg                                        >  2:1999170/1‑93 (MQ=255)
                                                    ccATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCGGATAAATCGc                                       <  2:318087/93‑1 (MQ=255)
                                                    ccATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTAATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCGGATAAATTGc                                       >  1:515918/1‑93 (MQ=255)
                                                     cATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCGGATAAATTGCt                                      <  2:6432065/93‑1 (MQ=255)
                                                     cATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCGGATAAATTGCt                                      >  2:5078165/1‑93 (MQ=255)
                                                      aTTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCGGATAAATTGCTg                                     <  1:2404447/93‑1 (MQ=255)
                                                      aTTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCGGATAAATTGCTg                                     >  2:4742539/1‑93 (MQ=255)
                                                      aTTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCGGATAAATCGCTg                                     <  2:551542/93‑1 (MQ=255)
                                                       ttGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCGGATAAATTGCTgg                                    >  1:3713265/1‑93 (MQ=255)
                                                            gAAAACATTTTCCTCGGTCGTGAGTTTGTTAATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCGGATAAATTGCTGGCTaaa                               >  2:6002438/1‑93 (MQ=255)
                                                             aaaaCATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCGGATAAATTGCTGGCTAAAc                              >  2:4173098/1‑93 (MQ=255)
                                                              aaaCATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCGGATAAATCGCTGGCTAAACt                             <  1:2771806/93‑1 (MQ=255)
                                                                 cATTTTCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCGGATAAATCGCTGGCTAAACTTaa                          <  2:6553385/93‑1 (MQ=255)
                                                                   ttttCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCGGATAAATTGCTGGCTAAACTTAAcc                        <  1:6905763/93‑1 (MQ=255)
                                                                   ttttCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCGGATAAATTGCTGGCTAAACTTAAcc                        <  1:4160540/93‑1 (MQ=255)
                                                                   ttttCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCGGATAAATCGCTGGCTAAACTTAAcc                        >  1:92907/1‑93 (MQ=255)
                                                                    tttCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCGGATAAATTGCTGGCTAAACTTAACCt                       <  1:5378643/93‑1 (MQ=255)
                                                                     ttCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCGGATAAATTGCTGGCTAAACTTAACCTg                      <  1:352946/93‑1 (MQ=255)
                                                                     ttCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCGGATAAATTGCTGGCTAAACTTAACCTg                      >  1:2962759/1‑93 (MQ=255)
                                                                     ttCCTCGGTCGTGAGTTTGTTAATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCGGATAAATTGCTGGCTAAACTTAACCTg                      >  2:1081817/1‑93 (MQ=255)
                                                                      tCCTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCGGATAAATTGCTGGCTAAACTTAACCTgc                     <  1:1814817/93‑1 (MQ=255)
                                                                        cTCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCGGATAAATCGCTGGCTAAACTTAACCTgcgc                   >  1:1321545/1‑93 (MQ=255)
                                                                         tCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCGGATAAATTGCTGGCTAAACTTAACCTGCGCt                  >  1:1940463/1‑93 (MQ=255)
                                                                         tCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCGGATAAATTGCTGGCTAAACTTAACCTGCGCt                  <  2:3894100/93‑1 (MQ=255)
                                                                         tCGGTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCGGATAAATTGCTGGCTAAACTTAACCTGCGCt                  >  2:24369/1‑93 (MQ=255)
                                                                           ggTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCGGATAAATTGCTGGCTAAACTTAACCTGCGCttt                >  1:1227954/1‑93 (MQ=255)
                                                                           ggTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCGGATAAATTGCTGGCTAAACTTAACCTGCGCttt                >  1:6661733/1‑93 (MQ=255)
                                                                            gTCGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCGGATAAATTGCTGGCTAAACTTAACCTGCGCTTTa               >  2:1544744/1‑93 (MQ=255)
                                                                              cGTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCGGATAAATTGCTGGCTAAACTTAACCTGCGCTTTaaa             <  1:6019997/93‑1 (MQ=255)
                                                                               gTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCGGATAAATTGCTGGCTAAACTTAACCTGCGCTTTaaaa            <  2:4977840/93‑1 (MQ=255)
                                                                               gTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCGGATAAATTGCTGGCTAAACTTAACCTGCGCTTTaaaa            <  1:6431820/93‑1 (MQ=255)
                                                                               gTGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCGGATAAATTGCTGGCTAAACTTAACCTGCGCTTTaaaa            <  1:4826930/93‑1 (MQ=255)
                                                                                tGAGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCGGATAAATTGCTGGCTAAACTTAACCTGCGCTTTAAAAg           <  2:4084414/93‑1 (MQ=255)
                                                                                  aGTTTGTTGATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCGGATAAATTGCTGGCTAAACTTAACCTGCGCTTTAAAAGCg         >  1:6286415/1‑93 (MQ=255)
                                                                                   gTTTGTTAATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCGGATAAATTGCTGGCTAAACTTAACCTGCGCTTTAAAAGCGa        >  2:483316/1‑93 (MQ=255)
                                                                                     ttgttgATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCGGATAAATTGCTGGCTAAACTTAACCTGCGCTTTAAAAGCGACa      <  2:4688978/93‑1 (MQ=255)
                                                                                       gTTGATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCGGATAAATTGCTGGCTAAACTTAACCTGCGCTTTAAAAGCGACAAg    <  1:7228647/93‑1 (MQ=255)
                                                                                         tAATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCGGATAAATTGCTGGCTAAACTTAACCTGCGCTTTAAAAGCGACAAGCt  <  2:290767/93‑1 (MQ=255)
                                                                                          |                                                                                           
GAAGCCGGGATTGGGATTATCCATCAGGAACTGAACCTGATCCCGCAGTTGACCATTGCCGAAAACATTTTCCTCGGTCGTGAGTTTGTTAATCGCTTTGGCAAAATTGACTGGAAAACCATGTATGCCGAAGCGGATAAATTGCTGGCTAAACTTAACCTGCGCTTTAAAAGCGACAAGCT  >  NC_000913/3931635‑3931816

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