Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,511,146 C→T W104* (TGG→TGA fecE ← ferric citrate ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,511,1460CT100.0% 246.7 / NA 77W104* (TGG→TGAfecEferric citrate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (36/41);  total (36/41)

AAATGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCG  >  NC_000913/4511014‑4511272
                                                                                                                                    |                                                                                                                              
aaaTGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGa                                                                                                                          >  1:1274629/1‑139 (MQ=255)
  aTGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGAtt                                                                                                                        <  2:1879276/139‑1 (MQ=255)
     cgcgcTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCa                                                                                                                              >  1:1820961/1‑130 (MQ=255)
       cgctggcgctgACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGAc                                                                                                                   <  2:654222/139‑1 (MQ=255)
        gctggcgctgACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGAcc                                                                                                                  >  2:372230/1‑139 (MQ=255)
              gctgACCGCCGGAAAGCTCGGTTAACCGAGGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGa                                                                                                            <  2:1274629/139‑1 (MQ=255)
                      ccgGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAAcc                                                                                                        >  2:1505205/1‑135 (MQ=255)
                      ccgGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAAcc                                                                                                        <  1:1505205/135‑1 (MQ=255)
                      ccgGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCt                                                                                                    >  2:1858617/1‑139 (MQ=255)
                          aaaGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTg                                                                                                <  1:1673778/139‑1 (MQ=255)
                            aGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGa                                                                                              <  1:1511024/139‑1 (MQ=255)
                                 ggTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACtgtg                                                                                         <  1:1060542/139‑1 (MQ=255)
                                       ccGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATcccc                                                                                   >  1:1122707/1‑139 (MQ=255)
                                        cgacgaACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCct                                                                                  <  1:1726348/139‑1 (MQ=255)
                                            gaACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTc                                                                                                   >  1:1283673/1‑118 (MQ=255)
                                            gaACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTc                                                                                                   <  2:1283673/118‑1 (MQ=255)
                                             aaCGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGc                                                                             >  1:544732/1‑139 (MQ=255)
                                             aaCGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGc                                                                             >  2:1670580/1‑139 (MQ=255)
                                             aaCGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGc                                                                             >  1:1628083/1‑139 (MQ=255)
                                               cGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATc                                                                                      <  2:353216/128‑1 (MQ=255)
                                               cGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATc                                                                                      >  1:353216/1‑128 (MQ=255)
                                               cGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGAAATTAACTCGTTCATTGTCTTCAGCGGGGAGACGCCCCCAGAGTGACAGTCGGGGATCACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGt                                                                           >  2:1194778/1‑139 (MQ=255)
                                                   aaGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTaaa                                                                       <  2:550038/139‑1 (MQ=255)
                                                   aaGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTaaa                                                                       <  2:836630/139‑1 (MQ=255)
                                                     gatgatTGATCCGGGTCTGGTTCATGGCACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATgg                                                                    >  2:381466/1‑139 (MQ=255)
                                                               ccGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGgc                                                           >  2:681183/1‑139 (MQ=255)
                                                                 gggTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGAc                                                                                                                   <  2:1208638/81‑1 (MQ=255)
                                                                 gggTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGAc                                                                                                                   >  1:1208638/1‑81 (MQ=255)
                                                                 gggTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGgcag                                                         <  1:1670580/139‑1 (MQ=255)
                                                                  ggTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAAcc                                                                                                        <  2:900370/91‑1 (MQ=255)
                                                                  ggTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAAcc                                                                                                        >  1:900370/1‑91 (MQ=255)
                                                                   gTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAAcc                                                                                                        <  2:674622/90‑1 (MQ=255)
                                                                   gTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAAcc                                                                                                        >  1:674622/1‑90 (MQ=255)
                                                                   gTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCct                                                                                  >  2:1125352/1‑112 (MQ=255)
                                                                   gTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCct                                                                                  <  1:1125352/112‑1 (MQ=255)
                                                                       ggTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGa                                                   >  2:1144054/1‑139 (MQ=255)
                                                                        gTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAAcc                                                                                                        <  2:421994/85‑1 (MQ=255)
                                                                        gTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAAcc                                                                                                        >  1:421994/1‑85 (MQ=255)
                                                                        gTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCct                                                                                  <  1:425175/107‑1 (MQ=255)
                                                                        gTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCct                                                                                  >  2:425175/1‑107 (MQ=255)
                                                                             tGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAAcc                                                                                                        >  1:1263240/1‑80 (MQ=255)
                                                                             tGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAAcc                                                                                                        <  2:1263240/80‑1 (MQ=255)
                                                                               gCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGc                                           <  1:274786/139‑1 (MQ=255)
                                                                                             tGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCctc                                                                                 <  2:267691/87‑1 (MQ=255)
                                                                                             tGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCctc                                                                                 >  1:267691/1‑87 (MQ=255)
                                                                                              gCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGgcag                                                         >  2:1041998/1‑110 (MQ=255)
                                                                                              gCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGgc                                                           <  2:1565039/108‑1 (MQ=255)
                                                                                              gCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGgc                                                           >  1:1565039/1‑108 (MQ=255)
                                                                                              gCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGCTAAATGGTGCTGAGgcag                                                         <  1:1041998/110‑1 (MQ=255)
                                                                                               cATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGa                           <  1:1144054/139‑1 (MQ=255)
                                                                                                aTTGTCTTCAGCGGAGAGACGCCCCCCGAGTGACAGTCAGGGATCACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGAt                          <  1:1194778/139‑1 (MQ=255)
                                                                                                aTTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGAt                          <  1:130479/139‑1 (MQ=255)
                                                                                                  tGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATga                        >  2:1702877/1‑139 (MQ=255)
                                                                                                    tCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAAcc                                                                                                        <  1:1856070/57‑1 (MQ=255)
                                                                                                    tCTTCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAAcc                                                                                                        >  2:1856070/1‑57 (MQ=255)
                                                                                                      ttCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGc                    <  1:273883/139‑1 (MQ=255)
                                                                                                      ttCAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGc                    <  1:1014860/139‑1 (MQ=255)
                                                                                                        cAGGGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATga                        <  1:1291693/133‑1 (MQ=255)
                                                                                                        cAGCGGAGAGACGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATga                        >  2:1291693/1‑133 (MQ=255)
                                                                                                              agagaCGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCgg                                         <  1:32214/110‑1 (MQ=255)
                                                                                                              agagaCGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCgg                                         >  2:32214/1‑110 (MQ=255)
                                                                                                              agagaCGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGAt                          <  2:76686/125‑1 (MQ=255)
                                                                                                              agagaCGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGAt                          >  1:76686/1‑125 (MQ=255)
                                                                                                              agagaCGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTAt            <  2:1628083/139‑1 (MQ=255)
                                                                                                              agagaCGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTAt            <  1:681183/139‑1 (MQ=255)
                                                                                                                agaCGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTg                                                                                                <  2:465756/53‑1 (MQ=255)
                                                                                                                agaCGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTg                                                                                                >  1:465756/1‑53 (MQ=255)
                                                                                                                  aCGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCCGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGa        >  1:1914/1‑139 (MQ=255)
                                                                                                                   cGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATc                                                                                      <  2:1327137/60‑1 (MQ=255)
                                                                                                                   cGCCCCCAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATc                                                                                      >  1:1327137/1‑60 (MQ=255)
                                                                                                                     cccccAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTaaa                                                                       >  1:94182/1‑73 (MQ=255)
                                                                                                                     cccccAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTaaa                                                                       <  2:94182/73‑1 (MQ=255)
                                                                                                                     cccccAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGc                    >  2:1380089/1‑124 (MQ=255)
                                                                                                                     cccccAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGc                    <  1:1380089/124‑1 (MQ=255)
                                                                                                                        ccAGAGTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCATCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCg  <  1:1858617/139‑1 (MQ=255)
                                                                                                                            agTGACAGTCAGGGATTACGACCATACGAAACCAGGTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCg                                                    >  1:71753/1‑85 (MQ=255)
                                                                                                                            agTGACAGTCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCg                                                    <  2:71753/85‑1 (MQ=255)
                                                                                                                                    |                                                                                                                              
AAATGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCG  >  NC_000913/4511014‑4511272

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