Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,538,785 C→T *239* (TAG→TAA)  nanC ← N‑acetylnuraminic acid outer membrane channel protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,538,7850CT100.0% 246.8 / NA 75*239* (TAG→TAA) nanCN‑acetylnuraminic acid outer membrane channel protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (34/41);  total (34/41)

TCAATTGCTCCGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTACTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACA  >  NC_000913/4538656‑4538920
                                                                                                                                 |                                                                                                                                       
tCAATTGCTCCGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTa                                                                                                                                <  1:1454934/139‑1 (MQ=255)
  aaTTGCTCCGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTaaa                                                                                                                              <  1:67250/139‑1 (MQ=255)
  aaTTGCTCCGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTaaa                                                                                                                              >  1:1430243/1‑139 (MQ=255)
  aaTTGCTCCGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTa                                                                                                                                >  2:331935/1‑137 (MQ=255)
  aaTTGCTCCGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTa                                                                                                                                <  1:331935/137‑1 (MQ=255)
     tGCTCCGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGa                                                                                                                           <  1:1182158/139‑1 (MQ=255)
       cTCCGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGaca                                                                                                                         >  1:1379450/1‑139 (MQ=255)
             tACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCaa                                                                                                                   <  1:1096521/139‑1 (MQ=255)
             tACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCaa                                                                                                                   >  2:1520003/1‑139 (MQ=255)
                 aCTTTTAAATGGCACAGGAGTTTCCGGTCATACAGCCGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTCTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCg                                                                                                               >  1:1457432/1‑139 (MQ=255)
                  cTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGAAACACCAATGcaa                                                                                                              >  1:85916/1‑137 (MQ=255)
                    tttAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTctatta                                                                                                                                       >  1:1464953/1‑112 (MQ=255)
                      tAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTa                                                                                                                                >  1:1872426/1‑117 (MQ=255)
                      tAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTa                                                                                                                                <  2:1872426/117‑1 (MQ=255)
                          tGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGtt                                                                                                      <  1:283066/139‑1 (MQ=255)
                              tcaGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACCCCAATGCGATAACTGTTTTcc                                                                                                  <  2:314371/138‑1 (MQ=255)
                                aGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGa                                                                                                <  2:15519/139‑1 (MQ=255)
                                               tACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGcc                                                                                 <  2:185451/139‑1 (MQ=255)
                                                 cAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGt                                                                               <  1:1711069/139‑1 (MQ=255)
                                                 cAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGt                                                                               >  2:1791155/1‑139 (MQ=255)
                                                    aCGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGt                                                                            >  1:631312/1‑139 (MQ=255)
                                                    aCGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTATCCTGTCTATTACAGTTTAAGTGa                                                                                                                           <  1:1082223/92‑1 (MQ=255)
                                                    aCGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTATCCTGTCTATTACAGTTTAAATGa                                                                                                                           >  2:1082223/1‑92 (MQ=255)
                                                     cgcgTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTAt                                                                                        >  2:779467/1‑126 (MQ=255)
                                                     cgcgTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTAt                                                                                        <  1:779467/126‑1 (MQ=255)
                                                     cgcgTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTa                                                                           >  1:412311/1‑139 (MQ=255)
                                                      gcgTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTa                                                                                                                                <  2:606167/85‑1 (MQ=255)
                                                      gcgTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTa                                                                                                                                >  1:606167/1‑85 (MQ=255)
                                                      gcgTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTaa                                                                          >  1:158333/1‑139 (MQ=255)
                                                              ggCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTg                                                                  <  2:412311/139‑1 (MQ=255)
                                                              ggCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTg                                                                  <  2:1375043/139‑1 (MQ=255)
                                                                        gCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAgg                                                        <  2:84318/139‑1 (MQ=255)
                                                                        gCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAgg                                                        <  1:96505/139‑1 (MQ=255)
                                                                            tcatGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGaca                                                                                                                         <  2:1796413/70‑1 (MQ=255)
                                                                            tcatGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGaca                                                                                                                         >  1:1796413/1‑70 (MQ=255)
                                                                             catGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTTTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCGGACGGTCAAGGTAGTc                                                   <  2:1457432/139‑1 (MQ=255)
                                                                               tGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCat                                                 <  2:1430243/139‑1 (MQ=255)
                                                                                 aTAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGt                                                             <  1:833059/125‑1 (MQ=255)
                                                                                 aTAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGt                                                             >  2:833059/1‑125 (MQ=255)
                                                                                   aGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTc                                             >  2:1074573/1‑139 (MQ=255)
                                                                                    gCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGaca                                                                                                                         >  2:341584/1‑62 (MQ=255)
                                                                                    gCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGaca                                                                                                                         <  1:341584/62‑1 (MQ=255)
                                                                                     cAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCta                                           >  1:642559/1‑139 (MQ=255)
                                                                                       aGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTctatta                                                                                                                                       <  1:401288/45‑1 (MQ=255)
                                                                                       aGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTctatta                                                                                                                                       <  2:760606/45‑1 (MQ=255)
                                                                                         cgcCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCtatgta                                       <  2:631312/139‑1 (MQ=255)
                                                                                         cgcCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCtatgta                                       <  1:1520003/139‑1 (MQ=255)
                                                                                         cgcCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCtatgta                                       <  2:1379450/139‑1 (MQ=255)
                                                                                           ccGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTc                                                   >  2:1613087/1‑125 (MQ=255)
                                                                                           ccGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTc                                                   <  1:1613087/125‑1 (MQ=255)
                                                                                              ttATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCaa                                                                                                                   <  1:859544/58‑1 (MQ=255)
                                                                                              ttATTGTTTTATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCaa                                                                                                                   >  2:859544/1‑58 (MQ=255)
                                                                                                aTTGTTTTATTCATTGTGACTGCCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTa                                <  1:1074573/139‑1 (MQ=255)
                                                                                                      ttATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTg                                                                                   <  2:1017420/82‑1 (MQ=255)
                                                                                                      ttATTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTg                                                                                   >  1:1017420/1‑82 (MQ=255)
                                                                                                        aTTCATTGTGACTGTCTCCTGTCTATTATAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCt                                                                   >  2:1788919/1‑96 (MQ=255)
                                                                                                        aTTCATTGTGACTGTCTCCTGTCTATTATAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCt                                                                   <  1:1788919/96‑1 (MQ=255)
                                                                                                        aTTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCtatgtat                                      <  2:1273918/125‑1 (MQ=255)
                                                                                                        aTTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCtatgtat                                      >  1:1273918/1‑125 (MQ=255)
                                                                                                        aTTCATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCggg                        <  2:642559/139‑1 (MQ=255)
                                                                                                           cATTGTGACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGt                     >  2:421147/1‑139 (MQ=255)
                                                                                                              tgtgACTGTCTCCTGTCTATTACAGTTTAAATGACACACCaa                                                                                                                   <  2:170523/42‑1 (MQ=255)
                                                                                                              tgtgACTGTCTCCTGTCTATTACAGTTTAAATGACACACCaa                                                                                                                   >  1:170523/1‑42 (MQ=255)
                                                                                                              tgtgACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCg                                                                                                 >  2:442743/1‑60 (MQ=255)
                                                                                                              tgtgACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCg                                                                                                 <  1:442743/60‑1 (MQ=255)
                                                                                                               gtgACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAAt                                                                                           <  2:540276/65‑1 (MQ=255)
                                                                                                               gtgACTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAAt                                                                                           >  1:540276/1‑65 (MQ=255)
                                                                                                                  aCTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCa                                                           >  2:81393/1‑94 (MQ=255)
                                                                                                                  aCTGTCTCCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCa                                                           <  1:81393/94‑1 (MQ=255)
                                                                                                                      tctcCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCtatgta                                       <  1:48590/110‑1 (MQ=255)
                                                                                                                      tctcCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCtatgta                                       >  2:48590/1‑110 (MQ=255)
                                                                                                                        tcCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTg                                                                                   >  1:796088/1‑64 (MQ=255)
                                                                                                                        tcCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTg                                                                                   <  2:796088/64‑1 (MQ=255)
                                                                                                                        tcCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCat                                                 <  2:45784/98‑1 (MQ=255)
                                                                                                                        tcCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCat                                                 >  1:45784/1‑98 (MQ=255)
                                                                                                                        tcCTGTCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGt        >  1:1872254/1‑139 (MQ=255)
                                                                                                                             tCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTc                                                   >  1:854570/1‑91 (MQ=255)
                                                                                                                             tCTATTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTc                                                   <  2:854570/91‑1 (MQ=255)
                                                                                                                              ctattaCAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACa  >  1:950715/1‑139 (MQ=255)
                                                                                                                                 |                                                                                                                                       
TCAATTGCTCCGGTACCACTTTTAAATGGCACAGGAGTTTCCGGTAATACAGACGCGTTTGCGGCAAATGAAGCCATCATGATAGCAAGCGCCGTTATTGTTTTATTCATTGTGACTGTCTCCTGTCTACTACAGTTTAAATGACACACCAATGCGATAACTGTTTTCCGATAAATTATCTCTGCCGTTGTAAACACCCTGACGGTCAAGGTAGTCATATTCTATGTATGGCGTAATATCGGGCGTCATATGGTATTGTAGAACA  >  NC_000913/4538656‑4538920

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: ‑