Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,331,979 T→C 15.4% T253A (ACG→GCG)  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,9790TC15.4% 471.2 / 71.0 200T253A (ACG→GCG) yhbEEamA family inner membrane putative transporter
Reads supporting (aligned to +/- strand):  ref base T (97/72);  new base C (20/11);  total (117/83)
Fisher's exact test for biased strand distribution p-value = 5.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.22e-01

Only 100 of 214 total aligned reads displayed.

GAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGC  >  NC_000913/3331881‑3332078
                                                                                                  |                                                                                                   
 agaTAACCTAAAAGGTTTAACATCGGCCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGtga                                                                                                   <  2:1559072/100‑1 (MQ=255)
  gaTAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGtgat                                                                                                  <  1:3424615/100‑1 (MQ=255)
    tAACCTAAAAGGTTTAACATCGGTCTGGCGAAGGAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGtgatga                                                                                                >  2:2511352/1‑100 (MQ=255)
    tAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGtgatga                                                                                                >  1:789969/1‑100 (MQ=255)
    tAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGtgatga                                                                                                >  1:1285109/1‑100 (MQ=255)
    tAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGtgatga                                                                                                >  2:507846/1‑100 (MQ=255)
     aaCCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTgatgat                                                                                               <  1:2602310/100‑1 (MQ=255)
      aCCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATc                                                                                              >  2:3327262/1‑100 (MQ=255)
        cTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATcgc                                                                                            >  2:3435738/1‑100 (MQ=255)
           aaaGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCt                                                                                         <  2:1802144/100‑1 (MQ=255)
              ggTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCAc                                                                                      >  1:3742575/1‑100 (MQ=255)
              ggTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGCGATCATCACGCTCAc                                                                                      <  1:271475/100‑1 (MQ=255)
              ggTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGCGATCATCACGCTCAc                                                                                      <  1:1714344/100‑1 (MQ=255)
               gTTTACCATCGATCTAGCGAATAAGTCGGGCCAGGCCGGTGATAAAAAAAAAAAAAAAAACAGCGTAAACAGTGGGGTGAGCGCGATCATCACGCTCAcc                                                                                     <  2:1704155/100‑1 (MQ=255)
               gTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCAcc                                                                                     <  1:1829221/100‑1 (MQ=255)
               gTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGCGATCATCACGCTCAcc                                                                                     <  1:714254/100‑1 (MQ=255)
                 ttAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTg                                                                                   <  2:1169338/100‑1 (MQ=255)
                   aataTCGGTCTGGCGAAGAAATCGGGCCAGGCTAGTGATAAAAGATCTTAAAGAAACAGCATACACCGTGGGGTGAGCGCGATCATCACGCTCAccc                                                                                    <  2:1438030/94‑2 (MQ=25)
                   aaCATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTgcg                                                                                 >  1:867535/1‑100 (MQ=255)
                     cATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGct                                                                               <  1:3927019/100‑1 (MQ=255)
                     cATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGct                                                                               >  2:1575818/1‑100 (MQ=255)
                     cATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGct                                                                               >  1:1476843/1‑100 (MQ=255)
                      aTCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGctg                                                                              >  2:3716980/1‑100 (MQ=255)
                       tCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGctgc                                                                             <  1:2293915/100‑1 (MQ=255)
                           tctcGCGAAGAAATCGGGCGATGCCAGTGATAATAGATCTTAAAACAACAGCGTAAACAGTGGGGTGAGCATGATGATCGCGCTCACCTGCGctgcctg                                                                          <  2:1289498/99‑1 (MQ=255)
                            cTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCtgcctgcc                                                                        <  1:3018259/100‑1 (MQ=255)
                            cTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCtgcctgcc                                                                        >  1:1118654/1‑100 (MQ=255)
                            cTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCtgcctgcc                                                                        >  2:3220179/1‑100 (MQ=255)
                            cTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGNGAGCGTGATGATCGCGCTCACCTGCGCtgcctgcc                                                                        >  1:2404607/1‑100 (MQ=255)
                              ggCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAg                                                                      >  2:1211130/1‑100 (MQ=255)
                                cGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCg                                                                    <  2:3272421/100‑1 (MQ=255)
                                cGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCg                                                                    >  1:3497382/1‑100 (MQ=255)
                                cGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATCATCGCGCCCACCTGCGctgcct                                                                           >  2:3489115/1‑93 (MQ=255)
                                  aagaaATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAg                                                                  <  2:590832/100‑1 (MQ=255)
                                  aagaaATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAg                                                                  <  1:38804/100‑1 (MQ=255)
                                   agaaATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGc                                                                 <  2:1664654/100‑1 (MQ=255)
                                   agaaATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGc                                                                 <  1:1178162/100‑1 (MQ=255)
                                    gaaATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGGCAGGCAGCgcagg                                                                <  2:4140991/100‑5 (MQ=255)
                                     aaaTCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCa                                                               >  2:2702817/1‑100 (MQ=255)
                                     aaaTCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCa                                                               <  1:3024931/100‑1 (MQ=255)
                                      aaTCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCAt                                                              <  2:1929006/100‑1 (MQ=255)
                                      aaTCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCAt                                                              >  2:932375/1‑100 (MQ=255)
                                       aTCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATc                                                             >  1:255069/1‑100 (MQ=255)
                                        tCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCg                                                            >  1:4026853/1‑100 (MQ=255)
                                         cGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGc                                                           >  2:2816379/1‑100 (MQ=255)
                                          gggCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCt                                                          >  2:3691318/1‑100 (MQ=255)
                                           ggccaggccaGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCtt                                                         >  1:1796878/1‑100 (MQ=255)
                                             ccaggccagTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCTCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTcc                                                       <  1:2948953/100‑1 (MQ=255)
                                                ggccagTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGcc                                                    >  1:3360577/1‑100 (MQ=255)
                                                 gccagTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCa                                                   >  2:3392834/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)
                                                         cttccgATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCa                                           <  2:1308230/95‑1 (MQ=255)
                                                              gATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATcc                                      <  1:484804/100‑1 (MQ=255)
                                                                  tGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTAcc                                  >  1:1112978/1‑100 (MQ=255)
                                                                  tGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTAcc                                  >  2:1254685/1‑100 (MQ=255)
                                                                  tGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTAcc                                  >  2:4031661/1‑100 (MQ=255)
                                                                   gagcaaaaCAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCCGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCa                                 >  2:2123863/5‑100 (MQ=255)
                                                                   gAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCaagatcggaag                                 >  1:1308231/1‑90 (MQ=255)
                                                                   gAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCa                                 >  1:1269400/1‑100 (MQ=255)
                                                                   gAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCa                                 >  1:2349137/1‑100 (MQ=255)
                                                                   caaaaaaaCAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTAcc                                  <  2:2548972/98‑1 (MQ=255)
                                                                    aaaaaaaCGGCGGAGAGAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCaa                                >  2:2886482/1‑100 (MQ=255)
                                                                      aaaaaCAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAgg                              <  2:3288359/100‑1 (MQ=255)
                                                                      aaaaaCAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAgg                              <  1:702296/100‑1 (MQ=255)
                                                                      aaaaaCAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAgg                              <  1:1545698/100‑1 (MQ=255)
                                                                        aaaCAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCATACCAAGGTa                            >  2:2746951/1‑100 (MQ=255)
                                                                        aaaCAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTCAGCTGCATGCCAGCGAGCCATCGCGTCCGCCAGTGCTCCATATCCTACGaa                                >  2:2183935/1‑96 (MQ=255)
                                                                             gCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAg                       >  2:147316/1‑100 (MQ=255)
                                                                             gCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAg                       <  2:431111/100‑1 (MQ=255)
                                                                              cGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAgatcggaaga                      >  1:3179394/1‑91 (MQ=255)
                                                                                tAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTcc                    <  2:730001/100‑1 (MQ=255)
                                                                                 aaaCAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCg                   <  2:2917586/100‑1 (MQ=255)
                                                                                    cAGTGGGGCGAGCGCGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGCGGCATCGCTTCCGCCAGGGCGCCCTATCCTACCACAGTATTCAGTccccca                >  2:3643580/1‑96 (MQ=255)
                                                                                      gTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCaaaa              <  1:2680182/100‑1 (MQ=255)
                                                                                      gTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCaaaa              >  1:3049660/1‑100 (MQ=255)
                                                                                       tGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCaaaaa             >  2:3972624/1‑100 (MQ=255)
                                                                                       tGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCaaaaa             <  2:3127421/100‑1 (MQ=255)
                                                                                       tGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATACAGTCCGCaaaaa             <  1:2818692/100‑1 (MQ=255)
                                                                                        caggTGAGCGCGATCATCACGCTCACCTGCGCTGCCTGCCAGTGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCaaaaaa            >  1:2476438/3‑100 (MQ=38)
                                                                                        caggTGAGCGCGATCATCACGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCaaaaaa            >  1:1860812/3‑100 (MQ=38)
                                                                                        caggTGAGCGCGATCATCACGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCaaaaaa            >  1:965249/3‑100 (MQ=39)
                                                                                        caggTGAGCGCGATCATCACGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCaaaaaa            >  1:3890288/3‑100 (MQ=38)
                                                                                          ggTGAGCGCGATCATCACGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAAtt          >  1:3446011/1‑100 (MQ=255)
                                                                                          ggTGAGCGCGATCATCACGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAAtt          >  1:3194070/1‑100 (MQ=255)
                                                                                          ggTGAGCGCGATCATCACGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAAtt          >  1:4060426/1‑100 (MQ=255)
                                                                                           gTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTa         >  2:2364996/1‑100 (MQ=255)
                                                                                           gTGAGCGCGATCATCACGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTa         >  1:825712/1‑100 (MQ=255)
                                                                                            tGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTaa        >  2:850019/1‑100 (MQ=255)
                                                                                             gAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCACCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTaaa       <  1:4005196/100‑1 (MQ=255)
                                                                                              aGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAAc      <  2:3049647/100‑1 (MQ=255)
                                                                                              aGCGCGATCATCACGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAAc      >  2:4050687/1‑100 (MQ=255)
                                                                                               gCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACa     >  1:3606645/1‑100 (MQ=255)
                                                                                               gCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACa     >  1:1905245/1‑100 (MQ=255)
                                                                                                cGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACAt    >  1:919331/1‑100 (MQ=255)
                                                                                                cGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACAt    >  2:3742793/1‑100 (MQ=255)
                                                                                                 gTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATg   <  1:1427737/100‑1 (MQ=255)
                                                                                                  tgatgaTCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGc  <  1:2259042/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
GAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGC  >  NC_000913/3331881‑3332078

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: