Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,933,401 A→C 15.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,4010AC15.3% 281.8 / 47.5 118V96V (GTT→GTGfucAL‑fuculose‑1‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base A (42/58);  new base C (2/16);  total (44/74)
Fisher's exact test for biased strand distribution p-value = 1.57e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.24e-01

Only 100 of 122 total aligned reads displayed.

CGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGT  >  NC_000913/2933302‑2933500
                                                                                                   |                                                                                                   
cGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGa                                                                                                     >  2:1535823/1‑100 (MQ=255)
 gCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGaa                                                                                                    >  2:1215328/1‑100 (MQ=255)
   aTAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACg                                                                                                  <  2:1628990/100‑1 (MQ=255)
     aaGGCGCGCAAGGAATAGAATTACCTCCAGCCGCCGTAATCATGTAGTTAATAGCGGGGATCGATCGGTTAAGAATGGAATCTGCCGTTCATTGCACGGc                                                                                                <  1:587561/100‑1 (MQ=255)
     aaGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTTCATTGCACGGc                                                                                                <  2:2105689/100‑1 (MQ=255)
      aGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTTCATTGCACGGc                                                                                                <  1:587259/99‑1 (MQ=255)
      aGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCa                                                                                               >  2:2838298/1‑100 (MQ=255)
       ggCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTTCATTGCACGGCag                                                                                              <  2:1790030/100‑2 (MQ=255)
       ggCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTTCATTGCACGGCag                                                                                              <  1:146398/100‑2 (MQ=255)
       ggCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCAt                                                                                              >  1:995099/1‑100 (MQ=255)
       ggCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGACCGATCGGTTAAGAATGGAAACTGCCGTTCATTGCACGGCag                                                                                              <  1:1594383/100‑2 (MQ=255)
        gctcgcAAGGAATAGAATTACCGTCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTTCATTGCACGGCagt                                                                                             <  1:933045/97‑3 (MQ=255)
        gcgcgcAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTTCATTGCACGGCag                                                                                              <  1:1337929/99‑2 (MQ=255)
         ctcgcAAGGTACAGAGTTCGCTCCAGCCGCCGCAATCTTGTAGTGAATAGCGGGGATCGATCCGTTAAGAATGGAAACTTCCGTTCATTGCACGGCa                                                                                               <  1:1473280/95‑1 (MQ=255)
         cgcgcAAGGAATAGAATTACTGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATAGATCGGTTAAGAATGGACACTGCGGTGCCATGAACGACATg                                                                                             >  1:1210142/1‑99 (MQ=255)
         cgcgcAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTTCATTGCACGGCagtt                                                                                            <  2:2516173/100‑4 (MQ=255)
         cgcgcAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGa                                                                                            >  2:58202/1‑100 (MQ=255)
         cgcgcAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGa                                                                                            >  1:2083049/1‑100 (MQ=255)
         cgcgcAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGa                                                                                            >  2:2195572/1‑100 (MQ=255)
            gcAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTTCATTGCACGGCagtttcc                                                                                         <  2:2474074/100‑7 (MQ=255)
            gcAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATtg                                                                                         >  2:1782970/1‑100 (MQ=255)
            gcAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATtg                                                                                         <  1:1193478/100‑1 (MQ=255)
            gcAAGGAATAGAACTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTTCATTGCACGGCagtttcc                                                                                         <  2:2740532/100‑7 (MQ=255)
             cAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTGAGAATGGAAACTGCTGTGCAATGAACGGCATGATtgt                                                                                        >  1:1970533/1‑100 (MQ=255)
              aaGGAATAGAACTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATtgtg                                                                                       <  1:2813545/100‑1 (MQ=255)
                 gAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGaac                                                                                    <  2:1749267/100‑1 (MQ=255)
                 gAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGaac                                                                                    <  1:2838279/100‑1 (MQ=255)
                 gAATAGAAGTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGaac                                                                                    <  2:2752614/100‑1 (MQ=255)
                  aaTAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGaaca                                                                                   >  2:1941984/1‑100 (MQ=255)
                   aTAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGaacaa                                                                                  >  2:468027/1‑100 (MQ=255)
                   aTAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGaacaa                                                                                  >  2:1511707/1‑100 (MQ=255)
                   aTAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGCGaacaa                                                                                  >  2:1667527/1‑100 (MQ=255)
                     aGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAAcc                                                                                <  2:169815/100‑1 (MQ=255)
                     aGAATTACCGCCAGCCGCCGAAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGaacaac                                                                                 <  1:144460/99‑1 (MQ=255)
                      gAATTACCGCGAGCCGCCGCAATCATGTAGTTAATAGCGGGGATCTATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACcg                                                                               <  2:2475314/100‑1 (MQ=255)
                        aTTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACcgcg                                                                             <  2:1448185/100‑1 (MQ=255)
                        aTTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACcgcg                                                                             >  1:1798783/1‑100 (MQ=255)
                          tACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGtt                                                                           >  2:1831676/1‑100 (MQ=255)
                           aCCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTg                                                                          >  1:2372692/1‑100 (MQ=255)
                           aCCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTg                                                                          >  2:286744/1‑100 (MQ=255)
                           aCCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAAGGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTg                                                                          <  2:818453/100‑1 (MQ=255)
                            ccGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTgg                                                                         >  1:1603662/1‑100 (MQ=255)
                                 aGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATcc                                                                    <  2:2045231/100‑1 (MQ=255)
                                   ccgccgCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCgg                                                                  >  2:2651219/1‑100 (MQ=255)
                                      ccgcAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCt                                                               >  2:195379/1‑100 (MQ=255)
                                        gcAATCATGTAGTGAATAGCGGTGATCGATCCGTTAATAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGc                                                             <  1:1965556/100‑1 (MQ=255)
                                        gcAATCATGTAGTGAATAGCGGTGATCGATCCGTTAACAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGc                                                             <  1:1965030/100‑1 (MQ=255)
                                            tCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTg                                                         <  1:137152/100‑1 (MQ=255)
                                             cATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGa                                                        <  1:205311/100‑1 (MQ=255)
                                             cATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGa                                                        <  2:250626/100‑1 (MQ=255)
                                              aTGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGAt                                                       >  2:765765/1‑100 (MQ=255)
                                               tGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATa                                                      >  1:1496259/1‑100 (MQ=255)
                                                 tAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAgg                                                    >  1:424914/1‑100 (MQ=255)
                                                  aGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATTATTGTGAACACCCGCGTTGGCATCCGGTCTGCTTTGATAgac                                                   <  2:1970544/100‑3 (MQ=255)
                                                    tGAGTAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTg                                                 <  1:43482/100‑1 (MQ=255)
                                                    tGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTg                                                 <  2:1338541/100‑1 (MQ=255)
                                                     gAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGc                                                <  2:225347/100‑1 (MQ=255)
                                                     gAACAGCGGGGCTCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGc                                                <  2:2734006/100‑1 (MQ=255)
                                                       aTAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGTCATCCGGTCTGCTTTGATAGGCTGCCa                                              <  1:2454141/100‑1 (MQ=255)
                                                           cGGGGATCGATCGGTTAAGAAAGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATg                                          >  1:590682/1‑100 (MQ=255)
                                                            ggggATCGATCGGTTAAGAATGGAGACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATgg                                         <  2:833465/100‑1 (MQ=255)
                                                             gggATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGa                                        <  1:2387201/100‑1 (MQ=255)
                                                                     tcgGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCAt                                <  1:373772/100‑1 (MQ=255)
                                                                      cgGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCAtt                               >  2:1891080/1‑100 (MQ=255)
                                                                      cgGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGAAAACGCCAtt                               >  1:1444245/1‑100 (MQ=255)
                                                                      cgGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAAAAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGcc                                  >  1:1612306/1‑97 (MQ=255)
                                                                       ggTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTc                              >  2:225410/1‑100 (MQ=255)
                                                                         ttAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGGGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGc                            <  1:58201/100‑1 (MQ=255)
                                                                         ttAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGc                            <  1:2195559/100‑1 (MQ=255)
                                                                          tAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCt                           <  2:1644843/100‑1 (MQ=255)
                                                                          tAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCt                           >  1:454399/1‑100 (MQ=255)
                                                                          gatctaTGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCt                           <  1:1007416/95‑1 (MQ=255)
                                                                           aaGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCtt                          <  2:2121858/100‑1 (MQ=255)
                                                                            agaatccacaCTGCCGTGCCATGAACGTCATGTTTGTGCACAACCGCGTTGGCACTCGGTCTGCTTTGATAGGCTGCCATATAGAAACGCCATTCGCTTg                         <  1:2696878/91‑1 (MQ=255)
                                                                            aGAATGGAAACTGCCGTGCAGTGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTg                         >  1:1776409/1‑100 (MQ=255)
                                                                              aaTGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAg                       <  2:28114/100‑1 (MQ=255)
                                                                               aTGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAgg                      >  2:393907/1‑100 (MQ=255)
                                                                                tGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAggg                     <  1:1763633/100‑1 (MQ=255)
                                                                                 ggAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAgggg                    >  2:1233112/1‑100 (MQ=255)
                                                                                  gAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGa                   <  1:225409/100‑1 (MQ=255)
                                                                                   atcatGCCGTTCATTGCACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAg                  >  2:533769/5‑100 (MQ=255)
                                                                                    aaCTGCCGTGCAATGAACGGCATGATTGTGAACATCCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGc                 <  1:1535817/100‑1 (MQ=255)
                                                                                    aaCTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGc                 <  2:424916/100‑1 (MQ=255)
                                                                                    aaCTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGc                 <  2:1798791/100‑1 (MQ=255)
                                                                                    aaCTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGc                 <  2:2894100/100‑1 (MQ=255)
                                                                                      cTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCtt               <  2:710061/100‑1 (MQ=255)
                                                                                        gCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTttc             <  1:592969/100‑1 (MQ=255)
                                                                                        gCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTttc             <  2:587335/100‑1 (MQ=255)
                                                                                         ccGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTttcc            >  2:1633159/1‑100 (MQ=255)
                                                                                             gCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTttccttcc        >  1:909058/1‑100 (MQ=255)
                                                                                              cAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTtccttcct       >  1:321346/1‑100 (MQ=255)
                                                                                               aaTGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTc      <  2:1017437/100‑1 (MQ=255)
                                                                                               aaTGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTc      <  1:2205633/100‑1 (MQ=255)
                                                                                               aaTGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTc      <  1:2584156/100‑1 (MQ=255)
                                                                                                aTGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCa     <  1:827620/100‑1 (MQ=255)
                                                                                                 tGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCAt    >  2:1665558/1‑100 (MQ=255)
                                                                                                 tGAACGGCATGATTGTGAACAACCGCGTTGCCAGCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTGCCTTCCTCAt    >  1:1922488/1‑100 (MQ=255)
                                                                                                  gAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATg   <  1:2338373/100‑1 (MQ=255)
                                                                                                   aaCGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGt  <  1:1664971/100‑1 (MQ=255)
                                                                                                   aaCGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGt  <  1:1665551/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
CGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCATGT  >  NC_000913/2933302‑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: