Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,542,503 G→A W699* (TGG→TGA feoB → ferrous iron transporter protein B and GTP‑binding protein; membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,542,5030GA97.5% 262.3 / ‑2.3 80W699* (TGG→TGAfeoBferrous iron transporter protein B and GTP‑binding protein; membrane protein
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (39/39);  total (40/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

CCTGATTTTCGTCCTGCTGTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGGGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGGTGTGCATTCTGGCGGT  >  NC_000913/3542407‑3542598
                                                                                                |                                                                                               
ccTGATTTTCGTCCTGCTGTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGTTTCTCCATCCTGTgaggg                                                                                              <  2:622521/100‑1 (MQ=255)
ccTGATTTTCGTCCTGCTGTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTgaggg                                                                                              <  2:1819357/100‑1 (MQ=255)
ccTGATTTTCGTCCTGCTGTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTgaggg                                                                                              >  2:346882/1‑100 (MQ=255)
ccTGATTTTCGTCCTGCTGTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTgaggg                                                                                              >  2:82648/1‑100 (MQ=255)
 cTGATTTTCGTCCTGCTGTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTgaggg                                                                                              >  1:2060084/1‑99 (MQ=255)
 cTGATTTTCGTCCTGCTGTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGc                                                                                             <  1:632748/100‑1 (MQ=255)
 cTGATTTTCGTCCTGCTGTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGc                                                                                             <  2:1650997/100‑1 (MQ=255)
 cTGATTTTCGTCCTGCTGTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGc                                                                                             >  2:431130/1‑100 (MQ=255)
 cTGATTTTCGTCCTGCTGTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGc                                                                                             >  1:1167327/1‑100 (MQ=255)
   gATTTTCGTCCTGCTGTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTg                                                                                           <  1:1179058/100‑1 (MQ=255)
   gATTTTCGTCCTGCTGTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTg                                                                                           >  1:1428539/1‑100 (MQ=255)
   gATTTGCGGCCTGCTGCATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTg                                                                                           <  1:1306045/100‑1 (MQ=255)
     ttttCGTCCTGCTGTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGaa                                                                                         <  2:966654/100‑1 (MQ=255)
     ttttCGTCCTGCTGTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGaa                                                                                         >  2:189544/1‑100 (MQ=255)
         cGTCCTGCTGTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATc                                                                                     <  2:460244/100‑1 (MQ=255)
           tCCTGCTGTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGc                                                                                   >  2:1239252/1‑100 (MQ=255)
            ccTGCTGTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCt                                                                                  >  2:1054867/1‑100 (MQ=255)
            ccTGCTGTATGTACCATGTATCTCAGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCt                                                                                  <  2:1693586/100‑1 (MQ=255)
               gctgTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGGGGGCTGAATATCGCTTAc                                                                               >  2:1519483/1‑100 (MQ=255)
               gctgTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTAc                                                                               <  2:1648876/100‑1 (MQ=255)
                    cttccgatcttATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACt                                                                          <  2:1936150/90‑1 (MQ=255)
                    atgtaCCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACt                                                                          <  1:505768/100‑1 (MQ=255)
                     tgtaCCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTg                                                                         >  1:1671024/1‑100 (MQ=255)
                      gtaCCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTgg                                                                        <  2:1516105/100‑1 (MQ=255)
                       taCCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGc                                                                       <  1:114905/100‑1 (MQ=255)
                       taCCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGc                                                                       >  2:1135760/1‑100 (MQ=255)
                       taCCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGc                                                                       >  1:2123966/1‑100 (MQ=255)
                       taCCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGc                                                                       >  2:997544/1‑100 (MQ=255)
                        aCCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGc                                                                       >  1:2142278/1‑99 (MQ=255)
                        aCCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCa                                                                      <  2:1490689/100‑1 (MQ=255)
                        aCCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCa                                                                      >  2:21369/1‑100 (MQ=255)
                        aCCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCa                                                                      <  1:1119830/100‑1 (MQ=255)
                        aCCATGTATCTCGGTGATGGGGGCTATCGCACGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCa                                                                      >  1:950684/1‑100 (MQ=255)
                         ccATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCaa                                                                     <  1:552366/100‑1 (MQ=255)
                         ccATGTATCTCGGTGATGGGGGCTATCGCCCGTGAACCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCaa                                                                     >  2:1459275/1‑100 (MQ=255)
                          cATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAAc                                                                    >  1:1877824/1‑100 (MQ=255)
                          cATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAAc                                                                    >  1:749367/1‑100 (MQ=255)
                              tATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACtagatcggaag                                                                >  1:1936140/1‑90 (MQ=255)
                              tATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTg                                                                >  1:492634/1‑100 (MQ=255)
                              tATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTg                                                                >  2:2044064/1‑100 (MQ=255)
                              tATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTg                                                                >  1:205255/1‑100 (MQ=255)
                                 ctcGGTGATGGGGGCTATTGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCTCTCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTc                                                             <  1:182679/100‑1 (MQ=255)
                                 ctcGGTGATGGGGGCTATCGGCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGt                                                               <  1:1541280/98‑1 (MQ=255)
                                 ctcGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTc                                                             <  2:1330959/100‑1 (MQ=255)
                                      tGATGGGTGCTAGCGCCCGGGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGGGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCa                                                        <  2:964542/100‑1 (MQ=255)
                                       gATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAg                                                       >  2:1388512/1‑100 (MQ=255)
                                           ggggCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCg                                                   >  1:979758/1‑100 (MQ=255)
                                                   cGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTAc                                           <  2:988490/100‑1 (MQ=255)
                                                    gCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACa                                          >  1:434030/1‑100 (MQ=255)
                                                      ccGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGt                                        <  1:123464/100‑1 (MQ=255)
                                                      ccGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGt                                        >  1:429732/1‑100 (MQ=255)
                                                       cGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTc                                       <  1:1648005/100‑1 (MQ=255)
                                                        gTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCa                                      <  2:1649337/100‑1 (MQ=255)
                                                         tGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAg                                     <  2:56446/100‑1 (MQ=255)
                                                          gAATCAAGCCGTGGCTGGATGGGCTTCTCCACCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGc                                    <  1:125647/100‑1 (MQ=255)
                                                           aaTCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCa                                   <  2:135486/100‑1 (MQ=255)
                                                           aaTCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCa                                   >  1:114864/1‑100 (MQ=255)
                                                            aTCAAGCCGTGGCTGTATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCa                                   <  1:1300363/99‑1 (MQ=255)
                                                               aaGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCa                               <  1:1074278/100‑1 (MQ=255)
                                                                aGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCaa                              >  1:1678296/1‑100 (MQ=255)
                                                                   cGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACtt                           <  2:1314681/100‑1 (MQ=255)
                                                                      ggctggctGGGCGTCTCCATGCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCATCATCCTACTTAc                         <  1:2079352/99‑1 (MQ=255)
                                                                      ggCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACa                        >  1:660220/1‑100 (MQ=255)
                                                                      ggCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACa                        >  1:895851/1‑100 (MQ=255)
                                                                        cTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGc                      >  1:2132208/1‑100 (MQ=255)
                                                                         tGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGcc                     >  2:2085095/1‑100 (MQ=255)
                                                                          ggATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCt                    >  2:1532832/1‑100 (MQ=255)
                                                                          ggATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCt                    >  2:484865/1‑100 (MQ=255)
                                                                              gggCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGgtg                <  1:1239250/100‑1 (MQ=255)
                                                                               ggCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGgtgt               >  1:1790991/1‑100 (MQ=255)
                                                                               ggCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGgtgt               >  2:93893/1‑100 (MQ=255)
                                                                                gCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGGtgtg              <  1:346881/100‑1 (MQ=255)
                                                                                gCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGGtgtg              <  2:1478155/100‑1 (MQ=255)
                                                                                    ctcCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGGTGTGCAtt          >  1:1155745/1‑100 (MQ=255)
                                                                                     tcCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGGTGTGCATTc         <  2:524921/100‑1 (MQ=255)
                                                                                     tcCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGGTGTGCATTc         >  1:2043509/1‑100 (MQ=255)
                                                                                     tcCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGGTGTGCATTc         <  2:1862281/100‑1 (MQ=255)
                                                                                      ccATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGGTGTGCATTCt        <  1:1809518/100‑1 (MQ=255)
                                                                                         tCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGGTGTGCATTCTGGc     >  1:70570/1‑100 (MQ=255)
                                                                                          ccTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGGTGTGCATTCTGGCg    >  2:1396563/1‑100 (MQ=255)
                                                                                          ccTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGGTGTGCATTCTGGCg    <  1:1068504/100‑1 (MQ=255)
                                                                                           cTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGGTGTGCATTCTGGCgg   >  1:1217774/1‑100 (MQ=255)
                                                                                           cTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGGTGTGCATTCTGGCgg   <  2:1167329/100‑1 (MQ=255)
                                                                                           cTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAACATCCAACTTACAGCCTGGTGTGCATTCTGGCgg   <  2:2123978/100‑1 (MQ=255)
                                                                                            tgtgAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGGTGTGCATTCTGGCGGt  <  2:677101/100‑1 (MQ=255)
                                                                                                |                                                                                               
CCTGATTTTCGTCCTGCTGTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGGGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGGTGTGCATTCTGGCGGT  >  NC_000913/3542407‑3542598

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