Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,541,962 T→G 100% I519S (ATC→AGC)  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,541,9620TG99.3% 410.2 / ‑5.7 117I519S (ATC→AGC) feoBferrous iron transporter protein B and GTP‑binding protein; membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (51/65);  total (51/66)
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.53e-01

Only 100 of 124 total aligned reads displayed.

CATGGAGTTGCCGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCA  >  NC_000913/3541864‑3542057
                                                                                                  |                                                                                               
cATGGAGCTGCCGGTCTACCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGatt                                                                                                >  1:357744/1‑99 (MQ=255)
 aTGGAGTTGCCGGTCTATCATGTACCACACGTTAAAAGCCTTATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGagca                                                                                               <  1:1160141/100‑1 (MQ=255)
 aTGGAGTTGCCGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGagca                                                                                               <  2:172463/100‑1 (MQ=255)
    gAGTTGCCGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGagcatca                                                                                            >  2:325925/1‑100 (MQ=255)
    gAGTTGCCGGTCTATCATGTACCACACGCTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGagcatc                                                                                             <  2:480450/99‑1 (MQ=255)
        tGCCGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGt                                                                                        >  1:191329/1‑100 (MQ=255)
         gCCGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTc                                                                                       <  2:49796/100‑1 (MQ=255)
          ccGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCa                                                                                      >  1:2379819/1‑100 (MQ=255)
          ccGGTCTATCATGTAACACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCa                                                                                      <  1:686971/100‑1 (MQ=255)
           cGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAg                                                                                     >  2:1249800/1‑100 (MQ=255)
           cGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAg                                                                                     <  2:1103869/100‑1 (MQ=255)
            ggTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGc                                                                                    >  2:1046535/1‑100 (MQ=255)
            ggTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTACGTGCTGGTAAAGTGAGCATCATCGTCAGc                                                                                    >  1:93480/1‑100 (MQ=255)
                tATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCAttt                                                                                >  2:251258/1‑100 (MQ=255)
                  tCATTTTCCACACTTTAAGATCCTTATTATCCAGACCTGGCAGTGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTc                                                                              <  1:1034956/100‑1 (MQ=255)
                    aTGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCt                                                                            >  1:2982260/1‑100 (MQ=255)
                     tGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGTGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCtn                                                                           <  2:3021965/100‑2 (MQ=255)
                     tGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTg                                                                           >  2:215107/1‑100 (MQ=255)
                      gTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGa                                                                          <  2:857403/100‑1 (MQ=255)
                      gTACCACACGTTAAAAGCCTGATTATCAAGACCTGGCAGCGTCTGAAAGGCTTCGATCTGCGAGCTGGTAAAGTGAGCATCATCGACAGCATTTTCCTGa                                                                          >  2:1037801/1‑100 (MQ=255)
                       tACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAg                                                                         >  2:2684120/1‑100 (MQ=255)
                       tACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAg                                                                         <  1:2668563/100‑1 (MQ=255)
                         ccACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAgcg                                                                       >  1:2558149/1‑100 (MQ=255)
                         ccACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAgcg                                                                       >  2:2326790/1‑100 (MQ=255)
                         ccACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAgcg                                                                       >  1:1800496/1‑100 (MQ=255)
                           acacGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCt                                                                     <  2:2539089/100‑1 (MQ=255)
                           acacGTTAAAAAACTGATTATGCGGATCTGCCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAATGTATGCATCATTGTCAGCATTGTCCTGGGCGCt                                                                     <  2:3037717/100‑1 (MQ=255)
                            cacGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGTTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCtt                                                                    <  2:2435611/100‑1 (MQ=255)
                             acGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAGAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCttt                                                                   <  2:3018017/100‑1 (MQ=255)
                             acGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCttt                                                                   <  2:198703/100‑1 (MQ=255)
                              cGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGTCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCtttn                                                                  <  2:2591957/100‑2 (MQ=255)
                               gTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCa                                                                 <  1:325928/100‑1 (MQ=255)
                               gTTAAAAGCCTGATTAGCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCa                                                                 <  2:2982227/100‑1 (MQ=255)
                                ttAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCaa                                                                <  2:2492337/100‑1 (MQ=255)
                                  aaaaGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACa                                                              >  1:1931667/1‑100 (MQ=255)
                                  aaaaGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACa                                                              <  2:1934733/100‑1 (MQ=255)
                                  aaaaGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACa                                                              <  1:201827/100‑1 (MQ=255)
                                    aaGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGc                                                            >  2:2147582/1‑100 (MQ=255)
                                    aaGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGc                                                            >  1:953273/1‑100 (MQ=255)
                                     aGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCt                                                           <  2:953267/100‑1 (MQ=255)
                                        cTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTtct                                                        <  1:923843/100‑1 (MQ=255)
                                          gATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCg                                                      <  1:2852108/100‑1 (MQ=255)
                                          gATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTAAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCg                                                      >  2:1248696/1‑100 (MQ=255)
                                           aTTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGc                                                     <  1:1404945/100‑1 (MQ=255)
                                           aTTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGc                                                     <  1:2119408/100‑1 (MQ=255)
                                             tATCCAGACCTGGCAGCGTCTGAAAGTCTTCGTTCTGTGTGTTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTg                                                   <  1:2500261/100‑1 (MQ=255)
                                             tATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTg                                                   <  2:1823623/100‑1 (MQ=255)
                                                 cAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGc                                                >  2:1606211/1‑99 (MQ=255)
                                                 cAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCAATTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCg                                               >  1:2508193/1‑100 (MQ=255)
                                                  aGACCTGGCAGCGTCTGGAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCgg                                              <  1:2381130/100‑1 (MQ=255)
                                                  aGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCgg                                              <  2:1051790/100‑1 (MQ=255)
                                                   gACCTGGCAGCGTCTGAAAGGCTTCTTTCTGCGTTCTGGTAAAGTGAGCAGCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCgg                                              <  2:1538194/99‑1 (MQ=255)
                                                     ccTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGaa                                           <  1:2951396/100‑1 (MQ=255)
                                                      cTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGan                                           <  2:526296/99‑2 (MQ=255)
                                                       tGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATTGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGaaaa                                         <  1:793792/100‑1 (MQ=255)
                                                       tGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCCCGCTCAGCGAGaagc                                         <  1:1890871/100‑3 (MQ=255)
                                                       tGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCCCGCTCAGCGAGaagc                                         <  2:2676874/100‑3 (MQ=255)
                                                        ggCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAAt                                        >  2:369576/1‑100 (MQ=255)
                                                         gCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAAtc                                       <  2:2629672/100‑1 (MQ=255)
                                                           aGCGTCTGTAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCCCGCTCAGCGAGaagctgt                                      <  1:2028735/99‑6 (MQ=255)
                                                              gTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGa                                  <  1:1439528/100‑1 (MQ=255)
                                                              gTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGa                                  <  1:117439/100‑1 (MQ=255)
                                                              gTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCATTCAACAGCTTCTCGCTAAGCGGGAAAAtcgtc                                    >  1:2800823/1‑98 (MQ=255)
                                                                cTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATa                                <  2:36295/100‑1 (MQ=255)
                                                                cTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATa                                >  2:2797955/1‑100 (MQ=255)
                                                                cTGAAAGGCTTCGTTCTGCGAGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATa                                >  2:1855632/1‑100 (MQ=255)
                                                                 tGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATaa                               >  1:2080858/1‑100 (MQ=255)
                                                                  gAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAAc                              >  1:2192466/1‑100 (MQ=255)
                                                                   aaaGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACa                             <  1:2017673/100‑1 (MQ=255)
                                                                     aGGCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATc                           >  1:230250/1‑100 (MQ=255)
                                                                      ggCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCa                          <  2:3014690/100‑1 (MQ=255)
                                                                       gCTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCaa                         >  1:2591116/1‑100 (MQ=255)
                                                                        cTTCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAAc                        <  1:2391294/100‑1 (MQ=255)
                                                                         ttCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACg                       <  2:2417089/100‑1 (MQ=255)
                                                                          tCGTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACg                       <  2:1806975/99‑1 (MQ=255)
                                                                            gTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACt                    >  2:2409017/1‑100 (MQ=255)
                                                                            gTTCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACt                    <  2:1904539/100‑1 (MQ=255)
                                                                              tCTGCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCg                  >  2:2833508/1‑100 (MQ=255)
                                                                                tGCGTGCTGGTAAAGTGAGAATCATCGTCAGCATTTTCCTGAGCGGTTTCAAAAGCTTCTCGCTTAGCGGAAAAATCGTCTATAACATCAACGACTCgg                 >  2:901810/1‑99 (MQ=255)
                                                                                 gCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGgcg               >  1:1965983/1‑100 (MQ=255)
                                                                                 gCGTGCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGgcg               >  1:303261/1‑100 (MQ=255)
                                                                                     gCTGGTAAAGTGATCAGCATCGTCAGCATTTTACTGAGCGCGTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGAGCTgg           <  2:858354/100‑1 (MQ=255)
                                                                                     gCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTgg           >  1:2118402/1‑100 (MQ=255)
                                                                                     gCTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTgg           <  1:1167769/100‑1 (MQ=255)
                                                                                      cTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGc          >  2:1431602/1‑100 (MQ=255)
                                                                                      cTGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGc          >  2:2124438/1‑100 (MQ=255)
                                                                                       tGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGcg         <  1:1248708/100‑1 (MQ=255)
                                                                                       tGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGcg         <  2:2751727/100‑1 (MQ=255)
                                                                                       tGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGcg         <  2:255332/100‑1 (MQ=255)
                                                                                       tGGTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGcg         <  2:1346731/100‑1 (MQ=255)
                                                                                        ggTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGcgt        >  2:270928/1‑100 (MQ=255)
                                                                                        ggTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGcgt        >  2:1303171/1‑100 (MQ=255)
                                                                                         gTAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACCACTCGgcgc              >  2:2437698/1‑93 (MQ=255)
                                                                                          tAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGcgtcc      <  1:1106112/100‑1 (MQ=255)
                                                                                          tAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGcgtcc      >  2:2787057/1‑100 (MQ=255)
                                                                                          tAAAGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATAGTCGATAACATCAACGACTCGGCGCTGGcgtcc      >  1:1760123/1‑100 (MQ=255)
                                                                                           aaaGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACTTCAACGACTCGGCGCTGGcgtccg     <  2:181189/100‑1 (MQ=255)
                                                                                           aaaGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGcgtccg     <  2:2307909/100‑1 (MQ=255)
                                                                                             aGTGAGCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACGGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGcgtccgtc   <  2:433833/100‑1 (MQ=255)
                                                                                              gTGAGCATCATAGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCa  <  2:1404598/100‑1 (MQ=255)
                                                                                                  |                                                                                               
CATGGAGTTGCCGGTCTATCATGTACCACACGTTAAAAGCCTGATTATCCAGACCTGGCAGCGTCTGAAAGGCTTCGTTCTGCGTGCTGGTAAAGTGATCATCATCGTCAGCATTTTCCTGAGCGCTTTCAACAGCTTCTCGCTGAGCGGGAAAATCGTCGATAACATCAACGACTCGGCGCTGGCGTCCGTCA  >  NC_000913/3541864‑3542057

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