Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,258,992 C→T A147T (GCG→ACG)  fruA ← PTS system, fructose‑specific transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,258,9920CT99.6% 880.4 / NA 230A147T (GCG→ACG) fruAPTS system, fructose‑specific transport protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (128/101);  total (128/102)
Fisher's exact test for biased strand distribution p-value = 4.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 246 total aligned reads displayed.

CAGCAAATTTCGCCAGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGCGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGGTGTGTG  >  NC_000913/2258902‑2259081
                                                                                          |                                                                                         
cAGCAAATTTCGCCAGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCgtgc                                                                                         <  1:743039/93‑1 (MQ=255)
cAGCAAATTTCGCCAGATCCACTTCGATATCTGCCGCCACAAGCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCgtgc                                                                                         >  2:4518082/1‑93 (MQ=255)
 aGCAAATTTCGCCAGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGcc                                                                                        >  2:1473936/1‑93 (MQ=255)
    aaaTTTCGCCAGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAAc                                                                                     >  1:5185120/1‑93 (MQ=255)
        ttCGCCAGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGaa                                                                                 <  2:1528257/93‑1 (MQ=255)
          cGCCAGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAAcc                                                                               <  1:2750920/93‑1 (MQ=255)
            ccAGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCAc                                                                             >  1:393202/1‑93 (MQ=255)
             cAGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACg                                                                            >  2:6217299/1‑93 (MQ=255)
             cAGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACg                                                                            >  2:2062464/1‑93 (MQ=255)
             cAGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACg                                                                            >  1:6380723/1‑93 (MQ=255)
              aGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACgg                                                                           <  2:5178646/93‑1 (MQ=255)
              aGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACgg                                                                           <  2:5218257/93‑1 (MQ=255)
              aGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACgg                                                                           <  1:4890926/93‑1 (MQ=255)
                aTCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGt                                                                         <  1:2112978/93‑1 (MQ=255)
                aTCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCAAGGGt                                                                         >  1:684071/1‑93 (MQ=255)
                 tCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGtt                                                                        >  1:77775/1‑93 (MQ=255)
                 tCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGtt                                                                        <  2:1219070/93‑1 (MQ=255)
                  ccACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGttt                                                                       >  2:2387446/1‑93 (MQ=255)
                  ccACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGttt                                                                       <  1:2739107/93‑1 (MQ=255)
                  ccACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGttt                                                                       <  2:5798918/93‑1 (MQ=255)
                   cACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTc                                                                      >  1:2177895/1‑93 (MQ=255)
                    aCTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCa                                                                     >  2:3293686/1‑93 (MQ=255)
                    aCTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCa                                                                     >  2:2715949/1‑93 (MQ=255)
                    aCTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCa                                                                     >  1:1720143/1‑93 (MQ=255)
                      ttCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAAc                                                                   <  1:4773428/93‑1 (MQ=255)
                       tCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACt                                                                  >  2:3773750/1‑93 (MQ=255)
                          atatCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTc                                                               <  2:229671/93‑1 (MQ=255)
                           tatCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTca                                                              <  1:3194620/93‑1 (MQ=255)
                           tatCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGAAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTcc                                                              >  1:4295373/1‑92 (MQ=255)
                              cTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTcaccc                                                           <  1:5736001/93‑1 (MQ=255)
                              cTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGATTCAACTTTcaccc                                                           >  2:3737365/1‑93 (MQ=255)
                               cgccgccACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTcaccca                                                          <  1:3543476/92‑1 (MQ=255)
                                   gccACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAg                                                      >  2:1408750/1‑93 (MQ=255)
                                   gccACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAg                                                      >  2:2325258/1‑93 (MQ=255)
                                    ccACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGc                                                     >  2:5362272/1‑93 (MQ=255)
                                     cacaATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGcc                                                    <  1:2960279/93‑1 (MQ=255)
                                     cacaATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGcc                                                    >  1:179112/1‑93 (MQ=255)
                                      acaATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCa                                                   <  1:2602274/93‑1 (MQ=255)
                                       caATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCAc                                                  >  2:1320129/1‑93 (MQ=255)
                                       caATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCAc                                                  >  2:6008945/1‑93 (MQ=255)
                                       caATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCAc                                                  <  2:6085861/93‑1 (MQ=255)
                                         aTCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGt                                                >  2:4331552/1‑93 (MQ=255)
                                         aTCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGt                                                <  2:673947/93‑1 (MQ=255)
                                         aTCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGt                                                <  1:3794452/93‑1 (MQ=255)
                                         aTCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGt                                                >  1:4220609/1‑93 (MQ=255)
                                         aTCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGt                                                <  1:5312088/93‑1 (MQ=255)
                                          tCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGtt                                               <  2:5749914/93‑1 (MQ=255)
                                            aCCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGtttt                                             <  1:3267424/93‑1 (MQ=255)
                                             ccAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGttttt                                            >  1:5018299/1‑93 (MQ=255)
                                             ccAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGttttt                                            <  1:2181499/93‑1 (MQ=255)
                                             ccAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGttttt                                            <  1:2259623/93‑1 (MQ=255)
                                              cAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGtttttt                                           <  2:6425106/93‑1 (MQ=255)
                                               aGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTc                                          <  2:3196543/93‑1 (MQ=255)
                                               aGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTc                                          >  2:2936082/1‑93 (MQ=255)
                                               aGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTc                                          >  2:4704102/1‑93 (MQ=255)
                                                gATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCg                                         <  1:5579449/93‑1 (MQ=255)
                                                 aTCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGc                                        <  2:5275151/93‑1 (MQ=255)
                                                   ccGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCtt                                      >  1:1959402/1‑93 (MQ=255)
                                                   ccGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCtt                                      >  1:1150341/1‑93 (MQ=255)
                                                    cGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTc                                     >  2:1840344/1‑93 (MQ=255)
                                                     gctgctGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCg                                    >  1:633337/1‑93 (MQ=255)
                                                      ctgctgCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCgg                                   >  1:6209178/1‑93 (MQ=255)
                                                       tgctgcGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGt                                  <  2:5183993/93‑1 (MQ=255)
                                                        gctgcGACTTCTTCGGGAGTGATTGCATTGCCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGtt                                 <  1:4667721/93‑1 (MQ=255)
                                                        gctgcGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGtt                                 <  2:4194570/93‑1 (MQ=255)
                                                        gctgcGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGtt                                 >  2:4721497/1‑93 (MQ=255)
                                                          tgcGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTc                               >  1:1929857/1‑93 (MQ=255)
                                                           gcGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCa                              >  2:3172563/1‑93 (MQ=255)
                                                           gcGACTGCTTCGGGAGTGATTCCATCACCCGCGACAACAGAACCACGGGTTTCAACTTTCACCCCCCAGCCACGTTTTTTCGCTTCGGTTTCa                              <  1:3462625/93‑1 (MQ=255)
                                                             gACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAAt                            <  2:4044648/93‑1 (MQ=255)
                                                             caCTTCTTCGGGAGTGATTGCATTCCCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAAt                            <  1:1420702/92‑1 (MQ=255)
                                                              aCTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATg                           >  2:2592456/1‑93 (MQ=255)
                                                              aCTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATg                           >  1:4124122/1‑93 (MQ=255)
                                                               cttcttCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATgg                          >  1:4881004/1‑93 (MQ=255)
                                                                  cttcGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCGCCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCtt                       <  1:4306927/93‑1 (MQ=255)
                                                                  cttcGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCtt                       >  2:398985/1‑93 (MQ=255)
                                                                  cttcGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCtt                       <  1:3628727/93‑1 (MQ=255)
                                                                    tcGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCa                     <  1:4517366/93‑1 (MQ=255)
                                                                      gggAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAgc                   >  2:5059633/1‑93 (MQ=255)
                                                                      cggAGTGATTGCATTCCCCGTGCCAACAGAACCACGGTTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAgc                   <  1:4245047/92‑1 (MQ=255)
                                                                       ggAGTGATTGCATTACCCGTGCCACCAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAgcc                  <  1:1548622/93‑1 (MQ=255)
                                                                       ggAGTGATTGCATTACCCGTGCCAACAGACCCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAgcc                  <  2:4766689/93‑1 (MQ=255)
                                                                       ggAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTNANCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATNGCTTCAgcc                  <  2:5977840/93‑1 (MQ=255)
                                                                       ggAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAgcc                  <  1:4422310/93‑1 (MQ=255)
                                                                        gAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAgccg                 >  1:5365380/1‑93 (MQ=255)
                                                                         aGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAgccgc                <  1:6010957/93‑1 (MQ=255)
                                                                          gTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAgccgcc               >  2:4246826/1‑93 (MQ=255)
                                                                           tGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCa              >  1:1124312/1‑93 (MQ=255)
                                                                           tGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCa              >  1:4043224/1‑93 (MQ=255)
                                                                            gATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCAt             >  1:2070502/1‑93 (MQ=255)
                                                                             aTTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATa            >  1:1847024/1‑93 (MQ=255)
                                                                               tgaaTTCCCCGTGCCAACAGACCCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATaaa          <  1:3111080/90‑1 (MQ=255)
                                                                               tGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATaaa          <  1:5765465/93‑1 (MQ=255)
                                                                               tGCATTACCCGTCCCAACAGACCTACGGTTTTCATCTTTCCCCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATaaa          <  2:3659735/93‑1 (MQ=255)
                                                                                gCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAg         <  2:1150341/93‑1 (MQ=255)
                                                                                   ttACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCATGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGgtg      >  2:4820058/1‑93 (MQ=255)
                                                                                    tACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGgtgt     <  2:4106846/93‑1 (MQ=255)
                                                                                    tACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGgtgt     >  1:2444033/1‑93 (MQ=255)
                                                                                     ccccGTGCCAACAGACCCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGgtgtg    <  2:387697/92‑1 (MQ=255)
                                                                                      cccGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGgtgtgt   >  2:2646160/1‑93 (MQ=255)
                                                                                          |                                                                                         
CAGCAAATTTCGCCAGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGCGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGGTGTGTG  >  NC_000913/2258902‑2259081

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