Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,933,401 A→C 17.3% V96V (GTT→GTG fucA ← L‑fuculose‑1‑phosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,933,4010AC17.3% 273.6 / 52.6 120V96V (GTT→GTGfucAL‑fuculose‑1‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base A (49/50);  new base C (8/13);  total (57/63)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

Only 100 of 124 total aligned reads displayed.

GCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGT  >  NC_000913/2933303‑2933500
                                                                                                  |                                                                                                   
gCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGaa                                                                                                    >  1:782421/1‑100 (MQ=255)
 cATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAAc                                                                                                   >  2:1303761/1‑100 (MQ=255)
  aTAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACg                                                                                                  >  1:2462304/1‑100 (MQ=255)
  aTAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACg                                                                                                  <  2:1458737/100‑1 (MQ=255)
    aaGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTTCATTGCACGGc                                                                                                <  2:170461/100‑1 (MQ=255)
     aGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCa                                                                                               >  1:361439/1‑100 (MQ=255)
     aGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCa                                                                                               <  2:1301921/100‑1 (MQ=255)
      ggCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTTCATTGCACGGCag                                                                                              <  1:1049359/100‑2 (MQ=255)
        cgcgcAAGGAATAGAATTACCGCCAGCCGCCGCAATCATTTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTTCATTGCACGGCagtt                                                                                            <  1:1817216/100‑4 (MQ=255)
        cgcgcAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTTCATTGCACGGCagt                                                                                             <  2:2228243/99‑3 (MQ=255)
          cgcAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTTCATTGCACGGCagtttc                                                                                          <  2:1860508/100‑6 (MQ=255)
          cgcAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTTCATTGCACGGCagtttc                                                                                          <  2:1288384/100‑6 (MQ=255)
          cgcAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTTCATTGCACGGCagtttc                                                                                          <  1:1148594/100‑6 (MQ=255)
          cgcAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTTCATTGCACGGCagtttc                                                                                          <  1:612344/100‑6 (MQ=255)
            cAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATtgt                                                                                        >  2:2067622/1‑100 (MQ=255)
            cAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATtgt                                                                                        <  1:2162357/100‑1 (MQ=255)
             gaGGAATAGACTTACCGACAGCCGCAGAAATCATGTAGTGAATCGAGGAGCTCGATCGGTTAAGTATGGAAACTGCCGTGAAATGACCGGCATGATtgtg                                                                                       <  2:2435057/99‑1 (MQ=255)
             aaGGAATAGAATTACCGCCAGCCGCCGCAATCATGTCGGGAATCGCGGGGATCGATCTGTTAAGAATGGAAACTGCCGTTCATTGCACGGCagtttccat                                                                                       <  2:558570/100‑9 (MQ=25)
             aaGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATtgtg                                                                                       <  1:934268/100‑1 (MQ=255)
             aaGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATtgtg                                                                                       <  2:1102616/100‑1 (MQ=255)
             aaGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATtgtg                                                                                       <  1:2227921/100‑1 (MQ=255)
              aGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGa                                                                                      <  1:250507/100‑1 (MQ=255)
                gAATAGAACTACCGCCAGCCGCCGCAATTATGTAGTGACTAGCGGGGATCGATCGGTTAATAATGAAAACTGCCGTGCAATGAACGGCATGATTCTGa                                                                                      <  2:944764/98‑1 (MQ=255)
                 aaTAGAATTATCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGaaca                                                                                   >  1:2193877/1‑100 (MQ=255)
                 aaTAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGaaca                                                                                   <  2:2428870/100‑1 (MQ=255)
                   tAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGaacaac                                                                                 >  2:614999/1‑100 (MQ=255)
                     gAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACcg                                                                               >  1:508413/1‑100 (MQ=255)
                      aaTTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACcgc                                                                              >  2:1697926/1‑100 (MQ=255)
                      aaTTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACcgc                                                                              >  2:1698166/1‑100 (MQ=255)
                        ttACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGt                                                                            >  2:996864/1‑100 (MQ=255)
                          aCCGCCAGCCGCCGCAATCATTTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGt                                                                            <  1:2332114/98‑1 (MQ=255)
                          aCCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTg                                                                          <  1:2430171/100‑1 (MQ=255)
                          aCCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTg                                                                          <  1:1006547/100‑1 (MQ=255)
                           ccGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCTATGAACGGCATGATTGTGAACAACCGCGTTgg                                                                         <  1:2425755/100‑1 (MQ=255)
                           ccGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTgg                                                                         >  1:2203950/1‑100 (MQ=255)
                           ccGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTgg                                                                         >  2:321585/1‑100 (MQ=255)
                           ccGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTg                                                                          <  2:261202/99‑1 (MQ=255)
                             gCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCa                                                                       >  2:293837/1‑100 (MQ=255)
                               cAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATc                                                                     <  1:625485/100‑1 (MQ=255)
                                aGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATcc                                                                    <  1:208147/100‑1 (MQ=255)
                                    gccgcAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTc                                                                <  2:2412901/100‑1 (MQ=255)
                                     ccgcAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCt                                                               <  1:1051641/100‑1 (MQ=255)
                                        cAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCt                                                            <  2:979913/100‑1 (MQ=255)
                                          aTCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGTCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCttt                                                          <  2:547019/100‑1 (MQ=255)
                                           tgaTGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCGATGAACGGCGTGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTg                                                         <  2:2062126/98‑1 (MQ=255)
                                            cATGTAGTGAATAGCGGGGATCGATCGGTTAATAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGa                                                        <  2:1208082/100‑1 (MQ=255)
                                               gTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAg                                                     >  2:1768994/1‑100 (MQ=255)
                                                tAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAgg                                                    >  1:1459474/1‑100 (MQ=255)
                                                   tGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTg                                                 <  2:1383902/100‑1 (MQ=255)
                                                   tGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTg                                                 <  1:370021/100‑1 (MQ=255)
                                                    gAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGc                                                >  1:1608071/1‑100 (MQ=255)
                                                     aaTAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGcc                                               <  1:1490518/100‑1 (MQ=255)
                                                      aTAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCa                                              <  2:2193862/100‑1 (MQ=255)
                                                        aGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCata                                            >  2:1742553/1‑100 (MQ=255)
                                                         gCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCatat                                           >  1:1823400/1‑100 (MQ=255)
                                                         gCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCatat                                           <  1:1672447/100‑1 (MQ=255)
                                                          cGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATg                                          <  1:1152843/100‑1 (MQ=255)
                                                           ggggATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATgg                                         >  2:2085248/1‑100 (MQ=255)
                                                           ggggATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATgg                                         >  1:755299/1‑100 (MQ=255)
                                                           ggggATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAAAAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATgg                                         >  2:258570/1‑100 (MQ=255)
                                                             ggATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGaa                                       >  1:1663715/1‑100 (MQ=255)
                                                               atcgatcgGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAAc                                     <  1:251569/100‑1 (MQ=255)
                                                                  gatcgGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGcc                                  >  2:2147955/1‑100 (MQ=255)
                                                                  gatcgGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGcc                                  <  2:1773816/100‑1 (MQ=255)
                                                                  gatcgGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGcc                                  >  1:1622292/1‑100 (MQ=255)
                                                                  gatcgGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGAGAGGCTGCCATATGGAAACGc                                   >  2:441083/1‑99 (MQ=255)
                                                                      ggTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTc                              >  1:2342631/1‑100 (MQ=255)
                                                                      ggTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTc                              >  1:612573/1‑100 (MQ=255)
                                                                      ggTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTc                              >  2:1026299/1‑100 (MQ=255)
                                                                      ggTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTc                              >  1:914893/1‑100 (MQ=255)
                                                                      ggTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCGGCCATATGGAAACGCCATTc                              >  2:738282/1‑100 (MQ=255)
                                                                        ttAAGAATGGAGACTGCCGTGTAATGAACTGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACg                                    <  2:2245346/92‑1 (MQ=255)
                                                                           aGAATGGAAACTGACGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAGCGCCATTCGCTTg                         <  2:2393817/100‑1 (MQ=255)
                                                                             aaTGGAACCTGCCGTGCAATGATCGGCAGGATTGTGAACAACCGCTTTGGCAGCCGGTCTGCTTTGATAGGCTGACATATGGAAACGCCATTCGCTTGan                       <  2:394040/100‑2 (MQ=255)
                                                                             aaTGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAg                       >  2:383451/1‑100 (MQ=255)
                                                                              atgccaaCTGCCGTGCACTGCCCGGCCTGATTGTGAACAACCGCGTTGGCATCCGATCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAgg                      <  1:1790790/95‑1 (MQ=255)
                                                                              aTGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAgg                      <  1:1236820/100‑1 (MQ=255)
                                                                                caatcatGCCGTTCATTGCACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAgggg                    >  2:2312656/7‑100 (MQ=255)
                                                                                  atcatGCCGTTCATTGCACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAgggg                    >  1:1936064/5‑98 (MQ=255)
                                                                                  atcatGCCGTTCATTGCACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGa                   >  2:91029/5‑99 (MQ=255)
                                                                                  atcatGCCGTTCATTGCACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAg                  >  1:966448/5‑100 (MQ=255)
                                                                                  aaaCTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAg                  >  1:769458/1‑100 (MQ=255)
                                                                                   aaCTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAg                  <  2:1230967/99‑1 (MQ=255)
                                                                                     cTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCtt               <  1:2147969/100‑1 (MQ=255)
                                                                                     atGCCGTTCATTGCACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCtt               >  1:310356/2‑100 (MQ=255)
                                                                                        ccGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTttcc            >  2:1148729/1‑100 (MQ=255)
                                                                                         cGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTttcct           >  1:531063/1‑100 (MQ=255)
                                                                                          gTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTttcctt          >  2:238413/1‑100 (MQ=255)
                                                                                          gTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTttcctt          >  1:2084299/1‑100 (MQ=255)
                                                                                           tGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTttccttc         <  2:373502/100‑1 (MQ=255)
                                                                                            gCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTttccttcc        >  2:855028/1‑100 (MQ=255)
                                                                                             cAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTtccttcct       >  2:2520708/1‑100 (MQ=255)
                                                                                             cAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTtccttcct       >  1:663353/1‑100 (MQ=255)
                                                                                              aaTGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTc      <  1:2061053/100‑1 (MQ=255)
                                                                                               aTGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCa     >  1:616184/1‑100 (MQ=255)
                                                                                               aTGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCa     <  2:842655/100‑1 (MQ=255)
                                                                                                tGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGTCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCAt    <  1:235556/100‑1 (MQ=255)
                                                                                                 gAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATg   <  2:1135836/100‑1 (MQ=255)
                                                                                                  aaCGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGt  >  2:970916/1‑100 (MQ=255)
                                                                                                  aaCGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTGTCCTTCCTCATGt  <  1:1742565/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
GCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGT  >  NC_000913/2933303‑2933500

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