Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,933,398 A→T 17.2% H97Q (CAT→CAA 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,3980AT17.2% 269.6 / 50.9 117H97Q (CAT→CAAfucAL‑fuculose‑1‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base A (49/47);  new base T (8/12);  total (57/60)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.74e-02

Only 100 of 120 total aligned reads displayed.

GTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCA  >  NC_000913/2933300‑2933497
                                                                                                  |                                                                                                   
    cATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAAc                                                                                                >  2:1303761/1‑100 (MQ=255)
     aTAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACg                                                                                               <  2:1458737/100‑1 (MQ=255)
     aTAAGGCGAGCAAGGAATAGAATTCCCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACg                                                                                               <  2:1253452/100‑1 (MQ=255)
       aaGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTTCATTGCACGGc                                                                                             <  2:170461/100‑1 (MQ=255)
       aaGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAAc                                                                                                >  2:62194/1‑97 (MQ=255)
        aGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCa                                                                                            >  1:361439/1‑100 (MQ=255)
        aGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCa                                                                                            <  2:1301921/100‑1 (MQ=255)
         ggCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTTCATTGCACGGCag                                                                                           <  1:1049359/100‑2 (MQ=255)
          gcgcgcAGGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTTCATTGCACGGCa                                                                                            <  2:942238/98‑1 (MQ=255)
           cgcgcAAGGAATAGAATTACCGCCAGCCGCCGCAATCATTTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTTCATTGCACGGCagtt                                                                                         <  1:1817216/100‑4 (MQ=255)
           cgcgcAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTTCATTGCACGGCagt                                                                                          <  2:2228243/99‑3 (MQ=255)
            gcgcAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGAt                                                                                        <  2:309551/100‑1 (MQ=255)
             cgcAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTTCATTGCACGGCagtttc                                                                                       <  1:1621575/100‑6 (MQ=255)
             cgcAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTTCATTGCACGGCagtttc                                                                                       <  2:1860508/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)
                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)
                 aGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGa                                                                                   >  2:434987/1‑100 (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:2453740/1‑100 (MQ=255)
                        gAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACcg                                                                            >  1:508413/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)
                              ccGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTg                                                                       <  2:261202/99‑1 (MQ=255)
                                gCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCa                                                                    >  2:293837/1‑100 (MQ=255)
                                 ccAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCAt                                                                   <  2:1416065/100‑1 (MQ=255)
                                  cAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATc                                                                  <  1:625485/100‑1 (MQ=255)
                                   aGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATcc                                                                 <  1:208147/100‑1 (MQ=255)
                                       gccgcAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTc                                                             >  1:53189/1‑100 (MQ=255)
                                       gccgcAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTc                                                             <  2:2412901/100‑1 (MQ=255)
                                        ccgcAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCt                                                            <  1:1051641/100‑1 (MQ=255)
                                           cAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCt                                                         <  2:979913/100‑1 (MQ=255)
                                              tgaTGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCGATGAACGGCGTGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTg                                                      <  2:2062126/98‑1 (MQ=255)
                                               cATGTAGTGAATAGCGGGGATCGATCGGTTAATAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGa                                                     <  2:1208082/100‑1 (MQ=255)
                                               cATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGa                                                     <  1:857826/100‑1 (MQ=255)
                                                  gTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAg                                                  >  1:214891/1‑100 (MQ=255)
                                                  gTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAg                                                  >  2:1768994/1‑100 (MQ=255)
                                                   tAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAgg                                                 >  1:1459474/1‑100 (MQ=255)
                                                      tGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTg                                              <  2:1383902/100‑1 (MQ=255)
                                                      tGAATAGCGGGCATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAAAGGCATGATTGTGAACAACCCCGTTGGCATCTGGACTGCTTTGATAGGc                                                >  2:2169635/1‑98 (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:1672447/100‑1 (MQ=255)
                                                             cGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATg                                       <  1:1152843/100‑1 (MQ=255)
                                                              ggggATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATgg                                      <  1:1744221/100‑1 (MQ=255)
                                                              ggggATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATgg                                      >  1:755299/1‑100 (MQ=255)
                                                              ggggATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATgg                                      >  2:2085248/1‑100 (MQ=255)
                                                              ggggATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAAAAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATgg                                      >  2:258570/1‑100 (MQ=255)
                                                              ggggATCGATCGGCTAAGAATGGAAACTGCCGTGCCATGAACGGCATGATTGTGAAAACACACGTTGTCATCCGGGCTGCGTTGATAGGCTGCCa                                           >  2:2340273/1‑95 (MQ=255)
                                                                ggATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGaa                                    >  1:1663715/1‑100 (MQ=255)
                                                                  atcgatcgGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAAc                                  <  1:251569/100‑1 (MQ=255)
                                                                   tcgatcgGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACg                                 >  1:315842/1‑100 (MQ=255)
                                                                     gatcgGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGAGAGGCTGCCATATGGAAACGc                                >  2:441083/1‑99 (MQ=255)
                                                                         ggTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTc                           >  1:612573/1‑100 (MQ=255)
                                                                         ggTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTc                           >  1:914893/1‑100 (MQ=255)
                                                                         ggTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTc                           >  1:2342631/1‑100 (MQ=255)
                                                                         ggTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTc                           >  2:1026299/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                 >  1:652296/7‑100 (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)
                                                                                     atcatGCCGTTCATTGCACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAg               >  1:1936110/5‑100 (MQ=255)
                                                                                     aaaCTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAg               >  1:769458/1‑100 (MQ=255)
                                                                                      aaCTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAg               <  2:1230967/99‑1 (MQ=255)
                                                                                        atGCCGTTCATTGCACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCtt            >  1:310356/2‑100 (MQ=255)
                                                                                           ccGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTttcc         >  2:1148729/1‑100 (MQ=255)
                                                                                            cGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTttcct        <  2:1892247/100‑1 (MQ=255)
                                                                                            cGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTttcct        >  1:531063/1‑100 (MQ=255)
                                                                                             gTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTttcctt       >  2:238413/1‑100 (MQ=255)
                                                                                             gTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTttcctt       >  1:2084299/1‑100 (MQ=255)
                                                                                              tGCAATGAACGGCATGATTGTGTACAACCGCTTTGGCCTCCGGTCTGATTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTttccttc      <  2:1732800/100‑1 (MQ=255)
                                                                                              tGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTttccttc      <  2:373502/100‑1 (MQ=255)
                                                                                                cAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTtccttcct    >  1:663353/1‑100 (MQ=255)
                                                                                                cAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTtccttcct    >  2:2520708/1‑100 (MQ=255)
                                                                                                 aaTGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTc   <  1:2061053/100‑1 (MQ=255)
                                                                                                  aTGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCa  <  2:842655/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
GTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCCTCA  >  NC_000913/2933300‑2933497

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: