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,5030GA98.8% 289.7 / ‑5.9 86W699* (TGG→TGAfeoBferrous iron transporter protein B and GTP‑binding protein; membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (49/36);  minor base C (0/1);  total (49/37)
Fisher's exact test for biased strand distribution p-value = 4.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TTTTCGTCCTGCTGTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGGGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGGTGTGCATTCTGGCGGTT  >  NC_000913/3542412‑3542599
                                                                                           |                                                                                                
ttttCGTCCTGCTGTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGaa                                                                                          <  2:2271903/100‑1 (MQ=255)
ttttCGTCCTGCTGTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGaa                                                                                          <  2:1025203/100‑1 (MQ=255)
ttttCGTCCTGCTGTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGaa                                                                                          >  1:2879199/1‑100 (MQ=255)
  ttCGTCCTGCTGTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAata                                                                                        >  1:1375545/1‑100 (MQ=255)
  ttCGTCCTGCTGTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAata                                                                                        >  2:594632/1‑100 (MQ=255)
  ttCGTCCTGCTGTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAata                                                                                        >  2:594631/1‑100 (MQ=255)
   tCGTCCTGCTGTATGTACCATGTATCTCGGGGATGGGGGCTATCGCCCGTTAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTTAat                                                                                         <  1:966929/98‑1 (MQ=255)
    cGTCCTGCTGTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATc                                                                                      <  1:522012/100‑1 (MQ=255)
     gTCCTGCTGTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCg                                                                                     >  2:2015972/1‑100 (MQ=255)
     gTCCTGCTGTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCg                                                                                     <  1:584823/100‑1 (MQ=255)
      tCCTGCTGTATGTACCATGTATCTCGGTGATGGGGGCTATCGCACGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGc                                                                                    >  1:701952/1‑100 (MQ=255)
       ccTGCTGTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGc                                                                                    <  2:202562/99‑1 (MQ=255)
        ctgctgTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCtt                                                                                  >  1:1184743/1‑100 (MQ=255)
        ctgctgTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCt                                                                                   <  1:2429556/99‑1 (MQ=255)
        atgctgTATGTACCATGTATCTCGGAGATGGGGGCTATCGCGCGTGAATAAAGACGTGGATGGATGGGCTTATCCATCCTGTGAGGGCTGAata                                                                                        >  1:1000001/2‑94 (MQ=255)
          gctgTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGCTGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTAc                                                                                >  1:2517997/1‑100 (MQ=255)
          gctgTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTAc                                                                                >  1:2819825/1‑100 (MQ=255)
          gctgTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGAGCTGAATATCGCTTAc                                                                                >  1:1703201/1‑100 (MQ=255)
               atgtaCCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACt                                                                           >  2:698860/1‑100 (MQ=255)
                tgtaCCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTg                                                                          >  1:2732362/1‑100 (MQ=255)
                 gtaCCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTgg                                                                         >  2:1207264/1‑100 (MQ=255)
                 gtaCCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTgg                                                                         >  2:2446747/1‑100 (MQ=255)
                  taCCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGc                                                                        >  2:682198/1‑100 (MQ=255)
                  taCCATGTATCTCGGTGATGGGGGCTATCGCACGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGc                                                                        >  1:564158/1‑100 (MQ=255)
                     cATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATAGCTTACTCACTGGCAAc                                                                     >  1:1821323/1‑100 (MQ=255)
                       tGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACAt                                                                   >  2:1408669/1‑100 (MQ=255)
                        gTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACAtt                                                                  >  1:1934872/1‑100 (MQ=255)
                         tATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGCATATCGCTTACTCACTGGCAACATTg                                                                 <  2:1772328/100‑1 (MQ=255)
                          aTCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGt                                                                >  1:2108896/1‑100 (MQ=255)
                            ctcGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTc                                                              <  2:296244/100‑1 (MQ=255)
                               ggTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTAt                                                           >  1:1875662/1‑100 (MQ=255)
                                gTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATc                                                          <  1:557419/100‑1 (MQ=255)
                                gTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATc                                                          <  1:948793/100‑1 (MQ=255)
                                gTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATc                                                          <  1:2311282/100‑1 (MQ=255)
                                  ttcacgggCGATAGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAg                                                        >  2:1219225/6‑100 (MQ=255)
                                  gATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAg                                                        >  1:2833699/1‑100 (MQ=255)
                                    tGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGt                                                      <  1:98197/100‑1 (MQ=255)
                                    tGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGt                                                      >  2:1695116/1‑100 (MQ=255)
                                    tGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGt                                                      >  1:1163474/1‑100 (MQ=255)
                                         gCTATCGCCCGTGAATCACGCCGTGGCTGGATGGGCTCCCCCATCCTGTGCGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCCGGTCGCCa                                                 <  2:2466876/100‑1 (MQ=255)
                                             tCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTa                                             >  2:1037135/1‑100 (MQ=255)
                                              cGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTAc                                            <  2:1613835/100‑1 (MQ=255)
                                              cGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTAc                                            >  2:1484432/1‑100 (MQ=255)
                                              cGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTAc                                            <  2:121865/100‑1 (MQ=255)
                                              cGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTAc                                            <  2:22930/100‑1 (MQ=255)
                                              cGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTAc                                            >  1:843210/1‑100 (MQ=255)
                                               gCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGTAACATTGTTCTATCAGGTCGCCAGCTACa                                           >  2:2687732/1‑100 (MQ=255)
                                               gCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACa                                           <  2:1629006/100‑1 (MQ=255)
                                               gCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACa                                           >  1:1190956/1‑100 (MQ=255)
                                                cccGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAg                                          <  1:868510/100‑1 (MQ=255)
                                                cccGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAg                                          >  2:2871856/1‑100 (MQ=255)
                                                 ccGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATGTCGCGTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACa                                           <  1:2477833/98‑1 (MQ=255)
                                                    tGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAg                                      >  2:1319528/1‑100 (MQ=255)
                                                     gAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGc                                     <  1:980535/100‑1 (MQ=255)
                                                     gAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGc                                     >  1:1840855/1‑100 (MQ=255)
                                                      aaTCAAGCCGTGGCTGGATGGGCTTCTTCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCa                                    <  1:1312697/100‑1 (MQ=255)
                                                      aaTCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCa                                    >  2:858383/1‑100 (MQ=255)
                                                       aTCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCAt                                   >  2:69407/1‑100 (MQ=255)
                                                       aTCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCAt                                   >  2:26821/1‑100 (MQ=255)
                                                       aTCAAGCCGAGGCTGGACGGGCTTCTCCATCCTGTTAGGGCTGAATATCGCTTACTCACTGGCAACATAGTTCTATCAGGTCGCCAGCTACAGTCAGc                                     >  1:390943/1‑98 (MQ=255)
                                                         cAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTCTCAGGTCGCCAGCTACAGTCAGCATc                                  >  1:2830877/1‑99 (MQ=255)
                                                         cAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATcc                                 >  2:1206278/1‑100 (MQ=255)
                                                          aaGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCa                                >  1:849521/1‑100 (MQ=255)
                                                          aaGCCGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCa                                >  1:818950/1‑100 (MQ=255)
                                                              cGTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACtt                            <  2:2734086/100‑1 (MQ=255)
                                                               gTGGCTGGATGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTa                           >  1:2842048/1‑100 (MQ=255)
                                                                       aTGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTgg                   <  2:800537/100‑1 (MQ=255)
                                                                        tGGGCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTgg                   <  1:1984647/99‑1 (MQ=255)
                                                                         gggCTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGgtg                 <  1:698854/100‑1 (MQ=255)
                                                                            cTTCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGGTGTGc              >  1:753986/1‑100 (MQ=255)
                                                                            cTTCTCCATCCTGTGAAGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGGTGTGc              <  1:1695090/100‑1 (MQ=255)
                                                                             ttCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGGTGTGCa             <  2:1163492/100‑1 (MQ=255)
                                                                             ttCTCCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGGTGTGCa             <  2:2816851/100‑1 (MQ=255)
                                                                                tcCATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGGTGTGCATTc          <  2:2819863/100‑1 (MQ=255)
                                                                                 ccATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGGTGTGCATTCt         >  1:87400/1‑100 (MQ=255)
                                                                                 ccATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGGTGTGCATTCt         >  2:588265/1‑100 (MQ=255)
                                                                                 ccATCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGGTGTGCATTCt         >  1:589290/1‑100 (MQ=255)
                                                                                   aTCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGGTGTGCATTCTgg       <  2:115742/100‑1 (MQ=255)
                                                                                   aTCCTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGGTGTGCATTCTgg       <  2:1603990/100‑1 (MQ=255)
                                                                                     ccTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGGTGTGCATTCTGGCg     <  1:2181224/100‑1 (MQ=255)
                                                                                     ccTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGGTGTGCATTCTGGCg     <  2:839038/100‑1 (MQ=255)
                                                                                      cTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGGTGTGCATTCTGGCgn    <  1:26820/100‑2 (MQ=255)
                                                                                      cTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGGTGTGCATTCTGGCgg    <  2:2833737/100‑1 (MQ=255)
                                                                                      cTGTGAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGGTGTGCATTCTGGCgg    <  2:689696/100‑1 (MQ=255)
                                                                                       tgtgAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGGTGTGCATTCTGGCGGt   >  2:2216890/1‑100 (MQ=255)
                                                                                        gtgAGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGGTGTGCATTCTGGCGGtt  >  1:1527875/1‑100 (MQ=255)
                                                                                           |                                                                                                
TTTTCGTCCTGCTGTATGTACCATGTATCTCGGTGATGGGGGCTATCGCCCGTGAATCAAGCCGTGGCTGGATGGGCTTCTCCATCCTGTGGGGGCTGAATATCGCTTACTCACTGGCAACATTGTTCTATCAGGTCGCCAGCTACAGTCAGCATCCAACTTACAGCCTGGTGTGCATTCTGGCGGTT  >  NC_000913/3542412‑3542599

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