Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,933,395 G→T 17.2% C98* (TGC→TGA 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,3950GT17.2% 290.5 / 46.3 115C98* (TGC→TGAfucAL‑fuculose‑1‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base G (46/49);  new base T (8/12);  total (54/61)
Fisher's exact test for biased strand distribution p-value = 6.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.14e-01

Only 100 of 116 total aligned reads displayed.

AAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCC  >  NC_000913/2933297‑2933494
                                                                                                  |                                                                                                   
aaGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGc                                                                                                    <  2:1044720/100‑1 (MQ=255)
   gTCGCATAAGGCGCGCAAGGAATAGAATTACCGGCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAAt                                                                                                 <  2:2326976/100‑1 (MQ=255)
        aTAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACg                                                                                            >  1:2462304/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                                                                                    <  2:1860508/100‑6 (MQ=255)
                cgcAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTTCATTGCACGGCagtttc                                                                                    <  1:1621575/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                                                                                  <  1:2162357/100‑1 (MQ=255)
                   gaGGAATAGACTTACCGACAGCCGCAGAAATCATGTAGTGAATCGAGGAGCTCGATCGGTTAAGTATGGAAACTGCCGTGAAATGACCGGCATGATtgtg                                                                                 <  2:2435057/99‑1 (MQ=255)
                   aaGGAATAGAATTACCGCCAGCCGCCGCAATCATGTCGGGAATCGCGGGGATCGATCTGTTAAGAATGGAAACTGCCGTTCATTGCACGGCagtttccat                                                                                 <  2:558570/100‑9 (MQ=25)
                   aaGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATtgtg                                                                                 <  2:1102616/100‑1 (MQ=255)
                   aaGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATtgtg                                                                                 <  1:934268/100‑1 (MQ=255)
                   aaGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATtgtg                                                                                 <  1:2227921/100‑1 (MQ=255)
                    aGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGa                                                                                >  2:434987/1‑100 (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:2453740/1‑100 (MQ=255)
                           gAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACcg                                                                         >  1:508413/1‑100 (MQ=255)
                            aaTTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACcgc                                                                        >  2:1698166/1‑100 (MQ=255)
                            aaTTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACcgc                                                                        >  2:1697926/1‑100 (MQ=255)
                              ttACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGt                                                                      >  2:996864/1‑100 (MQ=255)
                                aCCGCCAGCCGCCGCAATCATTTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGt                                                                      <  1:2332114/98‑1 (MQ=255)
                                aCCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTg                                                                    <  1:1006547/100‑1 (MQ=255)
                                aCCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTg                                                                    <  1:2430171/100‑1 (MQ=255)
                                 ccGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCTATGAACGGCATGATTGTGAACAACCGCGTTgg                                                                   <  1:2425755/100‑1 (MQ=255)
                                 ccGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTgg                                                                   >  2:321585/1‑100 (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                                                          <  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)
                                                  cATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGa                                                  <  1:857826/100‑1 (MQ=255)
                                                     gTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAg                                               >  2:1768994/1‑100 (MQ=255)
                                                     gTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAg                                               >  1:214891/1‑100 (MQ=255)
                                                         tGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTg                                           <  2:1383902/100‑1 (MQ=255)
                                                         tGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTg                                           <  1:370021/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)
                                                               gCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCatat                                     >  1:1823400/1‑100 (MQ=255)
                                                                cGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATg                                    <  1:1152843/100‑1 (MQ=255)
                                                                 ggggATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATgg                                   >  1:755299/1‑100 (MQ=255)
                                                                 ggggATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATgg                                   >  2:2085248/1‑100 (MQ=255)
                                                                 ggggATCGATCGGCTAAGAATGGAAACTGCCGTGCCATGAACGGCATGATTGTGAAAACACACGTTGTCATCCGGGCTGCGTTGATAGGCTGCCa                                        >  2:2340273/1‑95 (MQ=255)
                                                                   ggATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGaa                                 >  1:1663715/1‑100 (MQ=255)
                                                                      tcgatcgGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACg                              >  1:315842/1‑100 (MQ=255)
                                                                        gatcgGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGcc                            >  2:2147955/1‑100 (MQ=255)
                                                                        gatcgGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGcc                            >  1:1622292/1‑100 (MQ=255)
                                                                        gatcgGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGcc                            >  1:118946/1‑100 (MQ=255)
                                                                        gatcgGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGcc                            <  2:1773816/100‑1 (MQ=255)
                                                                        gatcgGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGAGAGGCTGCCATATGGAAACGc                             >  2:441083/1‑99 (MQ=255)
                                                                            ggTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTc                        >  1:612573/1‑100 (MQ=255)
                                                                            ggTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTc                        >  1:914893/1‑100 (MQ=255)
                                                                            ggTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTc                        >  2:1026299/1‑100 (MQ=255)
                                                                            ggTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTc                        >  1:2342631/1‑100 (MQ=255)
                                                                              ttAAGAATGGAGACTGCCGTGTAATGAACTGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACg                              <  2:2245346/92‑1 (MQ=255)
                                                                                 aGAATGGAAACTGACGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAGCGCCATTCGCTTg                   <  2:2393817/100‑1 (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)
                                                                                        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)
                                                                                           cTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCtt         <  1:2147969/100‑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)
                                                                                                 tGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTttccttc   <  2:373502/100‑1 (MQ=255)
                                                                                                  gCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTttccttcc  >  2:855028/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
AAGGTCGCATAAGGCGCGCAAGGAATAGAATTACCGCCAGCCGCCGCAATCATGTAGTGAATAGCGGGGATCGATCGGTTAAGAATGGAAACTGCCGTGCAATGAACGGCATGATTGTGAACAACCGCGTTGGCATCCGGTCTGCTTTGATAGGCTGCCATATGGAAACGCCATTCGCTTGAGGGGAGCTTTCCTTCC  >  NC_000913/2933297‑2933494

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: