Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,331,983 G→C 15.4% I251M (ATC→ATG yhbE ← EamA family inner membrane putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,331,9830GC15.4% 551.5 / 82.4 199I251M (ATC→ATGyhbEEamA family inner membrane putative transporter
Reads supporting (aligned to +/- strand):  ref base G (97/71);  new base C (20/11);  total (117/82)
Fisher's exact test for biased strand distribution p-value = 5.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 216 total aligned reads displayed.

TAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAG  >  NC_000913/3331885‑3332081
                                                                                                  |                                                                                                  
tAACCTAAAAGGTTTAACATCGGTCTGGCGAAGGAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGtgatga                                                                                                   >  2:2511352/1‑100 (MQ=255)
tAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGtgatga                                                                                                   <  1:1867635/100‑1 (MQ=255)
 aaCCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTgatgat                                                                                                  <  1:2602310/100‑1 (MQ=255)
   ccTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATcg                                                                                                <  1:3441415/100‑1 (MQ=255)
   ccTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATcg                                                                                                <  1:2694694/100‑1 (MQ=255)
     tAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATcgcg                                                                                              >  2:3262934/1‑100 (MQ=255)
       aaaGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCt                                                                                            <  2:1802144/100‑1 (MQ=255)
          ggTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCAc                                                                                         >  1:3742575/1‑100 (MQ=255)
          ggTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGCGATCATCACGCTCAc                                                                                         <  1:1714344/100‑1 (MQ=255)
           gTTTACCATCGATCTAGCGAATAAGTCGGGCCAGGCCGGTGATAAAAAAAAAAAAAAAAACAGCGTAAACAGTGGGGTGAGCGCGATCATCACGCTCAcc                                                                                        <  2:1704155/100‑1 (MQ=255)
           gTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCAcc                                                                                        >  1:2291817/1‑100 (MQ=255)
           gTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCAcc                                                                                        >  1:2232294/1‑100 (MQ=255)
           gTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGCGATCATCACGCTCAcc                                                                                        <  1:579522/100‑1 (MQ=255)
            tttAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCt                                                                                       <  2:2621695/100‑1 (MQ=255)
             ttAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTg                                                                                      <  2:1169338/100‑1 (MQ=255)
               aataTCGGTCTGGCGAAGAAATCGGGCCAGGCTAGTGATAAAAGATCTTAAAGAAACAGCATACACCGTGGGGTGAGCGCGATCATCACGCTCAccc                                                                                       <  2:1438030/94‑2 (MQ=25)
               aaCATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTgcg                                                                                    <  1:3568067/100‑1 (MQ=255)
                    cGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGctgcc                                                                               <  1:1494103/100‑1 (MQ=255)
                     ggTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGctgcct                                                                              >  1:143552/1‑100 (MQ=255)
                      gTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGctgcctg                                                                             <  1:712932/100‑1 (MQ=255)
                        cTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGctgcctg                                                                             >  2:1635920/1‑98 (MQ=255)
                        cTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCtgcctgcc                                                                           >  1:1118654/1‑100 (MQ=255)
                        cTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCtgcctgcc                                                                           <  1:3018259/100‑1 (MQ=255)
                        cTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGNGAGCGTGATGATCGCGCTCACCTGCGCtgcctgcc                                                                           >  1:2404607/1‑100 (MQ=255)
                         tGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCa                                                                          <  1:318403/100‑1 (MQ=255)
                          ggCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAg                                                                         >  2:1211130/1‑100 (MQ=255)
                            cGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCg                                                                       >  1:2320824/1‑100 (MQ=255)
                               agaaATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGc                                                                    <  1:3693776/100‑1 (MQ=255)
                                gaaATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGGCAGGCAGCgcagg                                                                   <  2:4140991/100‑5 (MQ=255)
                                 aaaTCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCa                                                                  >  1:730002/1‑100 (MQ=255)
                                 aaaTCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCa                                                                  <  1:3024931/100‑1 (MQ=255)
                                  aaTCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCAt                                                                 <  2:3165337/100‑1 (MQ=255)
                                  aaTCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCAt                                                                 >  2:932375/1‑100 (MQ=255)
                                  aaTCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCAt                                                                 <  2:1963419/100‑1 (MQ=255)
                                   aTCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATc                                                                >  1:255069/1‑100 (MQ=255)
                                      gggCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCt                                                             >  2:3691318/1‑100 (MQ=255)
                                      gggCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCt                                                             <  1:4020789/100‑1 (MQ=255)
                                       ggccaggccaGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCtt                                                            >  1:1796878/1‑100 (MQ=255)
                                       ggccaggccaGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCt                                                             >  1:1238163/1‑99 (MQ=255)
                                         ccaggccagTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCTCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTcc                                                          <  1:2948953/100‑1 (MQ=255)
                                           aggccagTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGc                                                        >  2:1308205/1‑100 (MQ=255)
                                            ggccagTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGcc                                                       >  1:3360577/1‑100 (MQ=255)
                                            ggccagTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGcc                                                       >  1:1301751/1‑100 (MQ=255)
                                               cagTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAgg                                                    >  1:2116110/1‑100 (MQ=255)
                                               cagTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAgg                                                    >  2:2356416/1‑100 (MQ=255)
                                               cagTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACATGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAgg                                                    >  1:1582509/1‑100 (MQ=255)
                                                    aTAAAAGATCTGAAAAAAACAGCGTAACCAGTGGGGAGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGcc                                               >  2:2610287/1‑100 (MQ=255)
                                                    aTAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGcc                                               >  2:1083789/1‑100 (MQ=255)
                                                     tAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCa                                              <  1:519441/100‑1 (MQ=255)
                                                     tAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCa                                              >  1:127916/1‑100 (MQ=255)
                                                     tAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATAGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCa                                              >  2:3435065/1‑100 (MQ=255)
                                                      aaaaGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGccn                                              <  2:4044052/99‑2 (MQ=255)
                                                         aGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATc                                          <  2:3497365/100‑1 (MQ=255)
                                                          gATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATcc                                         <  1:484804/100‑1 (MQ=255)
                                                            tCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTa                                       >  2:4025930/1‑100 (MQ=255)
                                                             cTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTAc                                      >  1:897150/1‑100 (MQ=255)
                                                             cTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTAc                                      <  1:3918299/100‑1 (MQ=255)
                                                              tGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTAcc                                     <  2:1891115/100‑1 (MQ=255)
                                                              tGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTAcc                                     <  1:1843962/100‑1 (MQ=255)
                                                              tGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTAcc                                     >  1:1112978/1‑100 (MQ=255)
                                                              tGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTAcc                                     >  2:1254685/1‑100 (MQ=255)
                                                              tGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTAcc                                     >  1:3883534/1‑100 (MQ=255)
                                                              tGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTAACCTGCGCTGGCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTa                                       >  2:1569335/1‑98 (MQ=255)
                                                               gAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCaagatcggaag                                    >  1:1308231/1‑90 (MQ=255)
                                                               caaaaaaaCAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTAcc                                     <  2:2548972/98‑1 (MQ=255)
                                                                aaaaaaaCGGCGGAGAGAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCaa                                   >  2:2886482/1‑100 (MQ=255)
                                                                 aaaaaaCAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAg                                  >  1:3101734/1‑100 (MQ=255)
                                                                  aaaaaCAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAgg                                 <  1:3961349/100‑1 (MQ=255)
                                                                  aaaaaCAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAgg                                 <  1:1545698/100‑1 (MQ=255)
                                                                  aaaaaCAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAgg                                 <  1:4031531/100‑1 (MQ=255)
                                                                  aaaaaCAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAgg                                 >  2:3374609/1‑100 (MQ=255)
                                                                  aaaaaCAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAgg                                 <  2:3288359/100‑1 (MQ=255)
                                                                    aaaCAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTa                               >  1:3454002/1‑100 (MQ=255)
                                                                    aaaCAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTa                               >  2:2733903/1‑100 (MQ=255)
                                                                         gCGTCAACAGTGGGGAGAGCGTGATGATCGCGCGCAACTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGATATTc                            >  2:1390845/1‑98 (MQ=255)
                                                                         gCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAg                          >  2:147316/1‑100 (MQ=255)
                                                                             aaaCAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCg                      <  2:2917586/100‑1 (MQ=255)
                                                                                cAGTGGGGCGAGCGCGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGCGGCATCGCTTCCGCCAGGGCGCCCTATCCTACCACAGTATTCAGTccccca                   >  2:3643580/1‑96 (MQ=255)
                                                                                  gTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCaaaa                 >  1:3049660/1‑100 (MQ=255)
                                                                                   tGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCaaaaa                <  2:3127421/100‑1 (MQ=255)
                                                                                   tGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATACAGTCCGCaaaaa                <  1:2818692/100‑1 (MQ=255)
                                                                                   gcaggTGAGCGCGATCATCACGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCaaaaa                >  2:4115984/4‑100 (MQ=38)
                                                                                    caggTGAGCGCGATCATCACGCTCACCTGCGCTGCCTGCCAGTGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCaaaaaa               >  1:2476438/3‑100 (MQ=38)
                                                                                    caggTGAGCGCGATCATCACGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCaaaaaa               >  1:965249/3‑100 (MQ=39)
                                                                                     aggTGAGCGCGATCATCACGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCaaaaaa               >  2:546850/2‑99 (MQ=39)
                                                                                      ggTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAAtt             <  2:2811730/100‑1 (MQ=255)
                                                                                      ggTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAAtt             >  2:476160/1‑100 (MQ=255)
                                                                                       gTGAGCGCGATCATCACGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTa            >  1:825712/1‑100 (MQ=255)
                                                                                        tGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTaa           >  2:850019/1‑100 (MQ=255)
                                                                                         gAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTaaa          <  2:3438947/100‑1 (MQ=255)
                                                                                          aGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAAc         <  2:3049647/100‑1 (MQ=255)
                                                                                            cGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACAt       >  2:3742793/1‑100 (MQ=255)
                                                                                            cGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACAt       >  2:2142841/1‑100 (MQ=255)
                                                                                            cGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACAt       >  1:919331/1‑100 (MQ=255)
                                                                                             gTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCATTCCGCAAAAAATTAAACAt       >  2:623597/1‑99 (MQ=255)
                                                                                             gTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATg      <  1:2621175/100‑1 (MQ=255)
                                                                                               gatgatCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCAATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGc     >  1:1130152/1‑99 (MQ=255)
                                                                                                 tgatCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAg  >  2:3933323/1‑100 (MQ=255)
                                                                                                 tgatCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAg  <  1:4098822/100‑1 (MQ=255)
                                                                                                  gatCGCGCTCACCTGCACTGCCTGCCAGCGAGCCATAGCTTCCGCCAGGGCGCCATATCCTACCAAAGTATTCAGTCCGCAAAAAATTAAACATGCGa   >  2:76930/1‑98 (MQ=255)
                                                                                                  |                                                                                                  
TAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAG  >  NC_000913/3331885‑3332081

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: