Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,542,093 G→A V563M (GTG→ATG)  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,0930GA98.0% 361.5 / ‑4.7 102V563M (GTG→ATG) feoBferrous iron transporter protein B and GTP‑binding protein; membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (49/51);  total (49/53)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 104 total aligned reads displayed.

GCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAG  >  NC_000913/3541998‑3542189
                                                                                               |                                                                                                
gCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCa                                                                                              <  2:618217/100‑1 (MQ=255)
gCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCa                                                                                              >  1:626939/1‑100 (MQ=255)
gCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGCCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCa                                                                                              <  1:1569991/100‑1 (MQ=255)
  ttCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGcatgcatg                                                                                            <  1:2496222/100‑1 (MQ=255)
   tctcGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGa                                                                                           <  2:869729/100‑1 (MQ=255)
    ctcGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGcatgcatgca                                                                                          >  1:1718963/1‑98 (MQ=255)
    ctcGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGaa                                                                                          <  2:457352/100‑1 (MQ=255)
      cGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGa                                                                                        <  2:1825337/100‑1 (MQ=255)
      cGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGa                                                                                        >  1:175607/1‑100 (MQ=255)
        ccgAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATa                                                                                      <  2:949599/98‑1 (MQ=255)
           aGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACt                                                                                   <  1:1972830/100‑1 (MQ=255)
           aGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACt                                                                                   <  1:1970712/100‑1 (MQ=255)
              gggAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTggc                                                                                >  2:819172/1‑100 (MQ=255)
              gggAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTggc                                                                                >  1:1259643/1‑100 (MQ=255)
               ggAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTggca                                                                               <  1:2370126/100‑1 (MQ=255)
                gAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTggcag                                                                              >  1:2571974/1‑100 (MQ=255)
                 aaaaTCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTggcagg                                                                             <  1:2073511/100‑1 (MQ=255)
                  aaaTCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTggcaggc                                                                            >  2:2378134/1‑100 (MQ=255)
                  aaaTCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTggcaggc                                                                            <  2:448362/100‑1 (MQ=255)
                    aTCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCaa                                                                          >  1:2207264/1‑100 (MQ=255)
                    aTCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCaa                                                                          >  1:1155467/1‑100 (MQ=255)
                       gtcgATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACgg                                                                       >  2:1318706/1‑100 (MQ=255)
                          gATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTg                                                                    <  2:308545/100‑1 (MQ=255)
                           aTAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTgg                                                                   <  1:2814409/100‑1 (MQ=255)
                            tAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGc                                                                  >  2:730914/1‑100 (MQ=255)
                            tAACAGCAACGACTCTGCGCTGGCGGGCGGCAGCCGGGTGATCACCCCGGTCTTCAAGTCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGc                                                                  <  1:1685459/100‑1 (MQ=255)
                               cATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTg                                                               <  2:2600351/100‑1 (MQ=255)
                               cATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTg                                                               <  2:2416848/100‑1 (MQ=255)
                               cATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTg                                                               >  1:2499796/1‑100 (MQ=255)
                               cATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTg                                                               <  2:58567/100‑1 (MQ=255)
                               cATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTg                                                               >  2:1937354/1‑100 (MQ=255)
                               cATCAACGACTCGACGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTg                                                               >  1:1749140/1‑100 (MQ=255)
                                 tCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGtt                                                             <  1:2405129/100‑1 (MQ=255)
                                  cAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGttt                                                            >  1:1218242/1‑100 (MQ=255)
                                   aaCGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTa                                                           >  2:1791578/1‑100 (MQ=255)
                                   aaCGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTa                                                           >  2:1595981/1‑100 (MQ=255)
                                   aaCGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTa                                                           <  1:28852/100‑1 (MQ=255)
                                     cGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACa                                                         >  2:423971/1‑100 (MQ=255)
                                     cGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACa                                                         <  1:1600712/100‑1 (MQ=255)
                                     cGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACa                                                         <  2:2065112/100‑1 (MQ=255)
                                     cGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACa                                                         <  1:1129986/100‑1 (MQ=255)
                                     cGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACa                                                         <  2:2084615/100‑1 (MQ=255)
                                      gACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAg                                                        >  2:1685167/1‑100 (MQ=255)
                                       aCTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAgg                                                       <  2:2383404/100‑1 (MQ=255)
                                       aCTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAgg                                                       <  2:2285119/100‑1 (MQ=255)
                                         tCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTg                                                     >  2:1819307/1‑100 (MQ=255)
                                         tCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTg                                                     <  2:152175/100‑1 (MQ=255)
                                         tCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGCCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTg                                                     >  1:1510748/1‑100 (MQ=255)
                                           ggCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGcc                                                   <  1:2256507/100‑1 (MQ=255)
                                              gcTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATg                                                <  1:730908/100‑1 (MQ=255)
                                               cTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATgg                                               >  1:1436427/1‑100 (MQ=255)
                                                 ggCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCg                                             >  1:2557177/1‑100 (MQ=255)
                                                 ggCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCg                                             <  1:2759679/100‑1 (MQ=255)
                                                  gCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGa                                            >  2:2251283/1‑100 (MQ=255)
                                                   cgtccgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGaa                                           >  1:2118922/1‑100 (MQ=255)
                                                    gtccgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGaaa                                          <  1:470267/100‑1 (MQ=255)
                                                    gtccgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGaaa                                          <  2:867149/100‑1 (MQ=255)
                                                    gtccgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGaa                                           <  2:2656224/99‑1 (MQ=255)
                                                       cgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAaagaa                                       >  2:2620362/1‑100 (MQ=255)
                                                       cgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAaagaa                                       >  2:1231595/1‑100 (MQ=255)
                                                       cgtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAaagaa                                       >  2:1438965/1‑100 (MQ=255)
                                                        gtcAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAATATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTTCCATTGCGAaagaa                                       >  2:127694/1‑99 (MQ=255)
                                                          cAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTg                                    <  2:1587684/100‑1 (MQ=255)
                                                            gCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTgg                                   >  1:61229/1‑99 (MQ=255)
                                                            gCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGt                                  <  1:586561/100‑1 (MQ=255)
                                                             ccGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTa                                 >  2:1881418/1‑100 (MQ=255)
                                                              cGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAg                                <  1:2747444/100‑1 (MQ=255)
                                                              cGGGTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCACGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGt                                  <  2:542099/98‑1 (MQ=255)
                                                                ggTGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTg                              <  2:2291658/100‑1 (MQ=255)
                                                                  tGATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTggg                            <  2:2702760/100‑1 (MQ=255)
                                                                   gATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGt                           >  1:1215793/1‑100 (MQ=255)
                                                                   gATCACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGt                           <  2:222200/100‑1 (MQ=255)
                                                                      cACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACg                        >  1:857399/1‑100 (MQ=255)
                                                                      cACCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACg                        >  1:1630200/1‑100 (MQ=255)
                                                                       aCCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGc                       >  1:960857/1‑100 (MQ=255)
                                                                       aCCCCGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGc                       <  1:714308/100‑1 (MQ=255)
                                                                        ccccGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCt                      <  2:403192/100‑1 (MQ=255)
                                                                        ccccGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCt                      >  1:71136/1‑100 (MQ=255)
                                                                        ccccGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCt                      >  2:1815930/1‑100 (MQ=255)
                                                                         cccGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTc                     <  1:1092338/100‑1 (MQ=255)
                                                                         cccGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTc                     <  2:175609/100‑1 (MQ=255)
                                                                          ccGGTTTTCAAGCCAATTGGCGCGCATTAAGATAACTGGCAGGCAACGGGTGGGCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTAc                         <  1:1404507/95‑1 (MQ=255)
                                                                          ccGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCa                    >  1:269705/1‑100 (MQ=255)
                                                                          ccGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCa                    >  1:2250741/1‑100 (MQ=255)
                                                                          ccGGTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCa                    >  2:962992/1‑100 (MQ=255)
                                                                            ggTCTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTTCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAac                  >  2:558570/1‑100 (MQ=255)
                                                                               cTTCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACAcc               >  2:1536914/1‑100 (MQ=255)
                                                                                ttCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGTAGTGGTAGTGGGTACGCTCAACAccc              <  2:1718976/100‑1 (MQ=255)
                                                                                ttCAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACAccc              <  2:1259651/100‑1 (MQ=255)
                                                                                 tCAAGCCAATTGGCACGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCct             >  1:700992/1‑100 (MQ=255)
                                                                                  cAAGCCAATTGGCGTGGATGAAGATAACTGGGAGGGAATGGTGGGCCTGTTTACAGGTGGCATGGCGAAAGAAGGGGGAGTGGGTACGCTCAACACCctc            <  1:2323159/100‑1 (MQ=255)
                                                                                  cAAGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCctc            >  1:804107/1‑100 (MQ=255)
                                                                                    aGCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTa          >  2:749247/1‑100 (MQ=255)
                                                                                     gCCAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTac         >  1:2806562/1‑100 (MQ=255)
                                                                                      ccAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTaca        >  1:851906/1‑100 (MQ=255)
                                                                                       cAATTGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTacac       <  2:269709/100‑1 (MQ=255)
                                                                                          ttGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGc    <  2:2562955/100‑1 (MQ=255)
                                                                                           tGGCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCa   <  2:1570056/100‑1 (MQ=255)
                                                                                            ggCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAg  >  1:748973/1‑100 (MQ=255)
                                                                                            ggCATGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAg  <  1:2253191/100‑1 (MQ=255)
                                                                                               |                                                                                                
GCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCAGCCGGGTGATCACCCCGGTCTTCAAGCCAATTGGCGTGCATGAAGATAACTGGCAGGCAACGGTTGGCCTGTTTACAGGTGCCATGGCGAAAGAAGTGGTAGTGGGTACGCTCAACACCCTCTACACCGCAG  >  NC_000913/3541998‑3542189

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