Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,985,983 T→G K49T (AAA→ACA)  araF ← L‑arabinose ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,985,9830TG94.8% 485.6 / 7.8 153K49T (AAA→ACA) araFL‑arabinose ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (4/4);  new base G (73/72);  total (77/76)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 153 total aligned reads displayed.

GCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTTTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATAGCGGATTGTGACA  >  NC_000913/1985887‑1986079
                                                                                                |                                                                                                
gcgcCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTGTAt                                                                                               <  1:2073806/100‑1 (MQ=255)
 cgcCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTTTATc                                                                                              <  2:3441901/100‑1 (MQ=255)
 cgcCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTGTATc                                                                                              <  2:1260934/100‑1 (MQ=255)
  gcCACTGGCAGCCAGGCTGTCGATCGCGTTTAATGTTTTTTCGCCATCCGGCACGGCAATCTTACTAACCTCAAACCCTCAATCCTTCCCGGCTGTATCg                                                                                             >  2:2462927/1‑100 (MQ=255)
  gcCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAACCTTAATAACCCCAAACCCTAAATCCTTCCCGGCTGTATCg                                                                                             >  1:1529394/1‑100 (MQ=255)
    cACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTGTATCGGc                                                                                           >  1:619616/1‑100 (MQ=255)
     aCTGGCATCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCa                                                                                          <  2:1984812/100‑1 (MQ=255)
       tGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTTAATTACCTCCAACCCTAAATCCTTCCCGGCTGTATCGGCaaa                                                                                        <  1:3628839/100‑1 (MQ=255)
        ggCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAAc                                                                                       <  1:1455417/100‑1 (MQ=255)
        ggCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAAc                                                                                       <  1:1455446/100‑1 (MQ=255)
        ggCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAAc                                                                                       >  2:3464139/1‑100 (MQ=255)
           aGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTc                                                                                    <  2:867218/100‑1 (MQ=255)
            gCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCTGTCCGGCACGGCAATCTTAACAACCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACCTcc                                                                                   >  2:3083572/1‑100 (MQ=255)
             ccAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCa                                                                                  >  1:515259/1‑100 (MQ=255)
              cAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCAt                                                                                 >  2:1965232/1‑100 (MQ=255)
              cAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCAt                                                                                 >  1:1576800/1‑100 (MQ=255)
                ggCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTTAATAACGTCAAACCCTAAATCCTTCCCGGCTGTGTCGGCAAACATCCATTc                                                                               >  1:3752689/1‑100 (MQ=255)
                  cTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCgg                                                                             >  1:3570129/1‑100 (MQ=255)
                   tGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTTTATCGGCAAACTTCCATTCGGt                                                                            >  2:1709847/1‑100 (MQ=255)
                   tGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGt                                                                            >  1:3846069/1‑100 (MQ=255)
                    ntCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTc                                                                           >  1:3316132/2‑100 (MQ=255)
                       gATCGCGTTCAATGTTTTTTCGCCATCCGACACGGCAATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTgg                                                                        >  1:1725449/1‑100 (MQ=255)
                        aTCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGa                                                                       >  1:2755666/1‑100 (MQ=255)
                        aTCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGa                                                                       <  1:3662070/100‑1 (MQ=255)
                           gcgTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAAcc                                                                    >  2:2103742/1‑100 (MQ=255)
                             gTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCAc                                                                  >  2:1570939/1‑100 (MQ=255)
                             gTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTTAATAACCTCAAACCCGAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCAc                                                                  <  2:1260536/100‑1 (MQ=255)
                                cAATGTTTTTTCGCCATCCGGCACGGCAATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGc                                                               <  2:2378727/100‑1 (MQ=255)
                                cAATGTTTTTTCGCCATCCGGCACGGCAATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGc                                                               >  2:3550281/1‑100 (MQ=255)
                                 aaTGTTTTTTCGCCATCCGGCACGGCAATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGct                                                              >  1:287424/1‑100 (MQ=255)
                                  aTGTTTTTTCGCCATCCGGCACGGCAATCTTAATAACCTCAAACCCTATATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGctc                                                             <  2:3617066/100‑1 (MQ=255)
                                   tGTTTTTTCGCCATCCGGCACGGCAATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGctct                                                            >  1:1039511/1‑100 (MQ=255)
                                   tGTTTTTTCGCCATCCGGCACGGCAATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGctct                                                            >  2:2230231/1‑100 (MQ=255)
                                   tGTTTTTTCGCCATCCGGCACGGCAATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGctct                                                            >  2:1411741/1‑100 (MQ=255)
                                   tGTTTTTTCGCCATCCGGCACGGCAATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGctct                                                            <  1:1534657/100‑1 (MQ=255)
                                      tttttCGCCATCCGGCACGGCAATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTcc                                                         <  1:1119582/100‑1 (MQ=255)
                                      tttttCGCCATCCGGCACGGCAATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTcc                                                         <  2:586016/100‑1 (MQ=255)
                                      tttttCGCCATCCGGCACGGCAATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTcc                                                         <  2:3704789/100‑1 (MQ=255)
                                       ttttCGCCATCCGGCACGGCAATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCg                                                        <  2:606812/100‑1 (MQ=255)
                                         ttCGCCATCCGGCACGGCAATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGt                                                      <  1:2614305/100‑1 (MQ=255)
                                          tCGCCATCCGGCACGGCAATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGtt                                                     >  1:1395990/1‑100 (MQ=255)
                                           cGCCATCCGGCACGGCAATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTg                                                    >  1:3870912/1‑100 (MQ=255)
                                           cGCCATCCGGCACGGCAATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTg                                                    >  2:195991/1‑100 (MQ=255)
                                            gCCATCCGGCACGGCAATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGc                                                   <  2:3630534/100‑1 (MQ=255)
                                            gCCATCCGGCACGGCAATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGc                                                   >  2:2491493/1‑100 (MQ=255)
                                            gCCATCCGGCACGGCAATCTTAATAACCTCAAACCCTAAAACCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGc                                                   >  1:2009964/1‑100 (MQ=255)
                                             ccATCCGGCACGGCAATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCt                                                  <  1:1718392/100‑1 (MQ=255)
                                             ccATCCGGCACGGCAATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCt                                                  <  2:2866142/100‑1 (MQ=255)
                                               aTCCGGCACGGCAATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTc                                                >  2:2561680/1‑100 (MQ=255)
                                               aTCCGGCACGGCAATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTc                                                >  1:1649478/1‑100 (MQ=255)
                                                tCCGGCACGGCAATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCa                                               <  2:3199638/100‑1 (MQ=255)
                                                tCCGGCACGGCAATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCa                                               <  2:1587738/100‑1 (MQ=255)
                                                 ccGGCACGGCAATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTTTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCAc                                              >  1:3762805/1‑100 (MQ=255)
                                                 ccGGCACGGCAATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCAc                                              <  1:2264928/100‑1 (MQ=255)
                                                 ccGGCACGGCAATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCAc                                              <  2:2132164/100‑1 (MQ=255)
                                                 ccGGCACCGCAATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCAc                                              <  2:2132779/100‑1 (MQ=255)
                                                   ggcacggcaATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCa                                            <  1:217098/100‑1 (MQ=255)
                                                     cacggcaATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGa                                          <  2:1971141/100‑1 (MQ=255)
                                                       cggcaATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGaaa                                        >  1:987142/1‑100 (MQ=255)
                                                          caATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAAcc                                     <  1:2415100/100‑1 (MQ=255)
                                                            aTCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGa                                   <  1:2931688/100‑1 (MQ=255)
                                                             tCTTAATAACCTCAGACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAg                                  >  2:3759398/1‑100 (MQ=255)
                                                             tCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTTTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAg                                  <  1:2547338/100‑1 (MQ=255)
                                                                   taaCCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAg                            <  2:2679529/100‑1 (MQ=255)
                                                                   taaCCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAg                            >  2:2464810/1‑100 (MQ=255)
                                                                    acccTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAgg                           <  1:1306117/98‑1 (MQ=255)
                                                                    aaCCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAgg                           <  1:1108744/100‑1 (MQ=255)
                                                                    aaCCACAAACCCTAAATCCTTCCCGGCTGTATCAGCAAACTTCCATTCGCTCTGGAACCACGGCGCTTCCGGTTGCTTCACCAGATAACCGAGCTTCAgg                           <  2:80061/100‑1 (MQ=255)
                                                                     aCCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGt                          <  1:2206715/100‑1 (MQ=255)
                                                                     aCCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGt                          <  2:282710/100‑1 (MQ=255)
                                                                     aCCTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGt                          >  1:1841709/1‑100 (MQ=255)
                                                                       cTCAAGCGCTGATTCCTTCCCGGCTGTAGCGGCAAGATTCTATCCGGTCGGGATACGCGGCCCTTCCGGTGGCTTCACCAGAAAACCGAGCTTCAGGtcc                        <  2:3133248/100‑3 (MQ=255)
                                                                       cTCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTtc                        >  1:1647096/1‑100 (MQ=255)
                                                                        tCAAACCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTtct                       <  1:1049708/100‑1 (MQ=255)
                                                                          aaaCCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTcc                     <  1:3615889/100‑1 (MQ=255)
                                                                           aaCCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGCCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCg                    >  1:446939/1‑100 (MQ=255)
                                                                            aCCCTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGc                   >  1:824801/1‑100 (MQ=255)
                                                                             cccTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGcc                  >  1:2712520/1‑100 (MQ=255)
                                                                             cccTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGcc                  <  2:1147061/100‑1 (MQ=255)
                                                                              ccTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCa                 >  1:2598779/1‑100 (MQ=255)
                                                                              ccTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCa                 >  2:462661/1‑100 (MQ=255)
                                                                              ccTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCa                 <  2:1022441/100‑1 (MQ=255)
                                                                              ccTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCa                 <  1:3612958/100‑1 (MQ=255)
                                                                               cTAAATCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCAt                <  2:3939832/100‑1 (MQ=255)
                                                                                 aaaTCCTTCCCGTCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTACTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATAg              >  2:3844122/1‑100 (MQ=255)
                                                                                 aaaTCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTCCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATAg              >  2:3843869/1‑100 (MQ=255)
                                                                                  aaTCCTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATAGc             >  1:1548758/1‑100 (MQ=255)
                                                                                    ggcTTCGCGGCTGTACCGGCGAACTTCCATTCGGTATGGGACCACGGCGCTGCCGGTTGCTTCGCCAGAGGACCGAGCTTCAGGTTCTCCGCCATAGCgg           <  2:423837/98‑1 (MQ=255)
                                                                                     ccTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATAGCGGa          <  2:2916034/100‑1 (MQ=255)
                                                                                     ccTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATAGCGGa          >  2:1687437/1‑100 (MQ=255)
                                                                                      cTTCCCGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATAGCGGAt         >  1:74407/1‑100 (MQ=255)
                                                                                         cccGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATAGCGGATtgt      >  2:1502641/1‑100 (MQ=255)
                                                                                         cccGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATAGCGGATtgt      >  1:2330688/1‑100 (MQ=255)
                                                                                          ccGGCTTTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATAGCGGATtgtg     <  2:3858032/100‑1 (MQ=255)
                                                                                          ccGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATAGCGGATtgtg     <  2:2315992/100‑1 (MQ=255)
                                                                                           cGGCTTTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATAGCGGATTGTGa    <  1:3122343/100‑1 (MQ=255)
                                                                                           cGGCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATAGCGGATTGTGa    <  2:2090967/100‑1 (MQ=255)
                                                                                            ggCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATAGCGGATTGTGAc   <  2:1089532/100‑1 (MQ=255)
                                                                                             gCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATAGCGGATTGTGACa  <  2:2737074/100‑1 (MQ=255)
                                                                                             gCTGTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATAGCGGATTGTGACa  <  2:677058/100‑1 (MQ=255)
                                                                                                |                                                                                                
GCGCCACTGGCAGCCAGGCTGTCGATCGCGTTCAATGTTTTTTCGCCATCCGGCACGGCAATCTTAATAACCTCAAACCCTAAATCCTTCCCGGCTTTATCGGCAAACTTCCATTCGGTCTGGAACCACGGCTCTTCCGGTTGCTTCACCAGAAAACCGAGCTTCAGGTTCTCCGCCATAGCGGATTGTGACA  >  NC_000913/1985887‑1986079

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: