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,9920CT85.2% 698.3 / 75.9 264A147T (GCG→ACG) fruAPTS system, fructose‑specific transport protein
Reads supporting (aligned to +/- strand):  ref base C (22/17);  new base T (99/125);  total (122/142)
Fisher's exact test for biased strand distribution p-value = 1.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01

Only 100 of 268 total aligned reads displayed.

CAGCAAATTTCGCCAGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGCGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGGTGTGTGC  >  NC_000913/2258902‑2259082
                                                                                          |                                                                                          
cAGCAAATTTCGCCAGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCgtgc                                                                                          >  2:1979065/1‑93 (MQ=255)
cAGCAAATTTCGCCAGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCgtgc                                                                                          >  1:1757135/1‑93 (MQ=255)
 aGCAAATTTCGCCAGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGcc                                                                                         >  2:1140443/1‑93 (MQ=255)
 aGCAAATTTCGCCAGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGcc                                                                                         >  2:3273959/1‑93 (MQ=255)
     aaTTTCGCCAGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACa                                                                                     <  2:778731/93‑1 (MQ=255)
     aaTTTCGCCAGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACa                                                                                     <  1:901909/93‑1 (MQ=255)
       tttCGCCAGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGa                                                                                   <  1:5287894/93‑1 (MQ=255)
       tttCGCCAGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGa                                                                                   <  1:4672390/93‑1 (MQ=255)
        ttCGCCAGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGaa                                                                                  <  1:4164976/93‑1 (MQ=255)
        tccGCCAGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGaa                                                                                  <  1:5603786/91‑1 (MQ=255)
           gCCAGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCa                                                                               >  2:1268881/1‑93 (MQ=255)
            ccAGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCAc                                                                              <  2:4292667/93‑1 (MQ=255)
            ccAGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCAc                                                                              <  1:4068532/93‑1 (MQ=255)
            ccAGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCAc                                                                              <  2:4028696/93‑1 (MQ=255)
             cAGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACg                                                                             >  1:3837602/1‑93 (MQ=255)
             cAGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACg                                                                             <  2:6558460/93‑1 (MQ=255)
             cAGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGCGCCAACAGAACCACg                                                                             <  1:2456416/93‑1 (MQ=255)
              aGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGCGCCAACAGAACCACgg                                                                            >  1:3217847/1‑93 (MQ=255)
              aGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGCGCCAACAGAACCACgg                                                                            >  1:2084522/1‑93 (MQ=255)
               gATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACggg                                                                           >  2:6857593/1‑93 (MQ=255)
                aTCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAATAGAACCACGGGt                                                                          >  1:7094341/1‑93 (MQ=255)
                aTCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGt                                                                          >  2:5371165/1‑93 (MQ=255)
                aTCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGt                                                                          >  1:2715093/1‑93 (MQ=255)
                 tCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGtt                                                                         >  1:6049896/1‑93 (MQ=255)
                 tCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGtt                                                                         <  2:2045653/93‑1 (MQ=255)
                  ccACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGttt                                                                        <  1:3170264/93‑1 (MQ=255)
                  ccACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGCGCCAACAGAACCACGGGttt                                                                        >  2:6464756/1‑93 (MQ=255)
                   cACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTc                                                                       >  2:2761459/1‑93 (MQ=255)
                    aCTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCa                                                                      >  1:5745211/1‑93 (MQ=255)
                     cTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCaa                                                                     >  2:7379488/1‑93 (MQ=255)
                      ttCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAAc                                                                    <  1:7350278/93‑1 (MQ=255)
                      ttCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAAc                                                                    >  1:2927632/1‑93 (MQ=255)
                      ttCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAAc                                                                    >  2:4904886/1‑93 (MQ=255)
                      tccGATGTCTGCCGCCACAATCACCAGCTCCGCTGCTGCGTCTTCTTCGGGAGTGATTGCATTTCCCGTGCCAACAGAACCACGGGTTTCAAc                                                                    <  1:5398998/91‑1 (MQ=255)
                       tCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACt                                                                   >  2:3455753/1‑93 (MQ=255)
                         gATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACttt                                                                 >  2:4857756/1‑93 (MQ=255)
                          atatCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTc                                                                <  2:2931734/93‑1 (MQ=255)
                          atatCGCCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTc                                                                <  1:162595/93‑1 (MQ=255)
                           tatCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTca                                                               <  2:1427968/93‑1 (MQ=255)
                           tatCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTca                                                               <  1:7050409/93‑1 (MQ=255)
                           tatCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGCGCCAACAGAACCACGGGTTTCAACTTTca                                                               >  2:6155406/1‑93 (MQ=255)
                            atCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTcac                                                              >  2:4812539/1‑93 (MQ=255)
                            atCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTcac                                                              <  1:2480484/93‑1 (MQ=255)
                             tCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTcacc                                                             >  1:4856682/1‑93 (MQ=255)
                             tCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTcacc                                                             <  2:7248391/93‑1 (MQ=255)
                              cTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTcaccc                                                            <  2:5564443/93‑1 (MQ=255)
                                  cgccACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACccacca                                                        >  1:553410/1‑93 (MQ=255)
                                   gccACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAg                                                       >  1:4734195/1‑93 (MQ=255)
                                      acaATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACTCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCa                                                    >  2:377057/1‑93 (MQ=255)
                                        aaTCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACg                                                  <  1:3421356/93‑1 (MQ=255)
                                         aTCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGt                                                 <  1:5878905/93‑1 (MQ=255)
                                           cACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTCCCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCCCCAGCCACGttt                                               <  1:2196613/93‑1 (MQ=255)
                                            aCCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGtttt                                              >  1:6009270/1‑93 (MQ=255)
                                              cAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGtttttt                                            <  2:136166/93‑1 (MQ=255)
                                              cAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGCGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGtttttt                                            <  2:278303/93‑1 (MQ=255)
                                               aGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTc                                           <  2:6622323/93‑1 (MQ=255)
                                                  tCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCt                                        <  2:864625/93‑1 (MQ=255)
                                                  tCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCt                                        >  1:2905762/1‑93 (MQ=255)
                                                   ccGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCtt                                       >  1:7511412/1‑93 (MQ=255)
                                                   ccGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGCGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCtt                                       <  2:1712469/93‑1 (MQ=255)
                                                    cGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTc                                      <  2:6111494/93‑1 (MQ=255)
                                                    cGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGCGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTc                                      <  2:4047120/93‑1 (MQ=255)
                                                     gctgctGCGACTTCTTCGGGAGTGATTGCATTCCCCGCGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCg                                     <  1:437266/93‑1 (MQ=255)
                                                     gctgctGCGACTTCTTCGGGAGTGATTGCATTACCCGCGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCg                                     >  1:6681539/1‑93 (MQ=255)
                                                       tgctgcGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCCCCCACCAGCCACGTTTTTTCGCTTCGGt                                   <  2:7087700/93‑1 (MQ=255)
                                                       tgctgcGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGt                                   <  2:909467/93‑1 (MQ=255)
                                                       tgctgcGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGt                                   >  1:4344586/1‑93 (MQ=255)
                                                        gctgcGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGtt                                  <  2:2349620/93‑1 (MQ=255)
                                                        gctgcGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGtt                                  >  1:6520691/1‑93 (MQ=255)
                                                         ctgcGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGttt                                 <  2:2659982/93‑1 (MQ=255)
                                                            cGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCaa                              >  2:2534913/1‑93 (MQ=255)
                                                            cGACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCaa                              <  2:3453895/93‑1 (MQ=255)
                                                             gACTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAAt                             >  1:3706876/1‑93 (MQ=255)
                                                              aCTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATg                            >  1:1063002/1‑93 (MQ=255)
                                                              aCTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATg                            >  2:7475545/1‑93 (MQ=255)
                                                              aCTTCTTCGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATg                            <  2:5173883/93‑1 (MQ=255)
                                                               cttcttCGGGAGTGATTGCATTACCCGCGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATgg                           >  1:5477177/1‑93 (MQ=255)
                                                                ttcttcGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCGCGTTTTTTCGCTTCGGTTTCAATGGc                          <  1:5905505/93‑1 (MQ=255)
                                                                ttcttcGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGc                          >  1:2292361/1‑93 (MQ=255)
                                                                   ttcGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCCCCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTc                       <  1:5786009/93‑1 (MQ=255)
                                                                   ttcGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTc                       <  2:243474/93‑1 (MQ=255)
                                                                    tcGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCa                      <  1:1004948/93‑1 (MQ=255)
                                                                    tcGGGAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCa                      <  2:2066403/93‑1 (MQ=255)
                                                                      gggAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAgc                    <  2:4965972/93‑1 (MQ=255)
                                                                       ggAGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAgcc                   >  2:1554713/1‑93 (MQ=255)
                                                                         aGTGATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAgccgc                 <  2:3359089/93‑1 (MQ=255)
                                                                          gTGATTGCATTACCCGGGCCAACAGAACCACGGGTTTCAACTTTCACCCCCCAGCCCCGTTTTTTCCCTTCGGTTTAAATGGCTTCAgccccc                >  1:6205302/1‑93 (MQ=255)
                                                                            gATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCAt              >  1:3906910/1‑93 (MQ=255)
                                                                            gATTGCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTCTCAATGGCTTCAGCCGCCAt              >  2:1871011/1‑93 (MQ=255)
                                                                                gCATTACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAg          >  1:6759643/1‑93 (MQ=255)
                                                                                  tttccccGTGCCACCAAACCCACGGTTTTAAACTTCCCCCCCCCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGgt        <  1:301278/89‑1 (MQ=255)
                                                                                    tccccGTGCCAACAGACCCACGGTTTTCACCTTCCCCCCCCCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGgtgt      <  1:129618/91‑1 (MQ=255)
                                                                                    tACCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGgtgt      >  1:558473/1‑93 (MQ=255)
                                                                                     aCCCGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGgtgtg     >  1:5047409/1‑93 (MQ=255)
                                                                                     aCCCGCGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGgtgtg     >  2:4671654/1‑93 (MQ=255)
                                                                                      cccGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGgtgtgt    <  2:3116537/93‑1 (MQ=255)
                                                                                       ccGTGCCAACAGACCCAGGGGTTTCAACTTTCCCCCCCCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGGtgtgtg   <  1:1538567/93‑1 (MQ=255)
                                                                                       ccGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGGtgtgtg   <  2:463898/93‑1 (MQ=255)
                                                                                       ccGTGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGGtgtgtg   <  2:3602640/93‑1 (MQ=255)
                                                                                        cgcgCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGGTGTgcgc  >  1:6852915/1‑93 (MQ=255)
                                                                                          |                                                                                          
CAGCAAATTTCGCCAGATCCACTTCGATATCTGCCGCCACAATCACCAGATCCGCTGCTGCGACTTCTTCGGGAGTGATTGCATTACCCGCGCCAACAGAACCACGGGTTTCAACTTTCACCCACCAGCCACGTTTTTTCGCTTCGGTTTCAATGGCTTCAGCCGCCATAAAGGTGTGTGC  >  NC_000913/2258902‑2259082

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