Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,331,983 G→C 17.2% 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,9830GC17.2% 296.7 / 52.1 116I251M (ATC→ATGyhbEEamA family inner membrane putative transporter
Reads supporting (aligned to +/- strand):  ref base G (42/53);  new base C (18/2);  total (61/55)
Fisher's exact test for biased strand distribution p-value = 1.51e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.43e-01

Only 100 of 126 total aligned reads displayed.

ATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGC  >  NC_000913/3331884‑3332082
                                                                                                   |                                                                                                   
aTAACCTAACAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGAAAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGtgatg                                                                                                     <  1:1430499/100‑1 (MQ=255)
aTAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGtgatg                                                                                                     <  1:797671/100‑1 (MQ=255)
 tAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAACTCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGtgatga                                                                                                    >  1:367734/1‑100 (MQ=255)
   aCCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATc                                                                                                  <  1:159059/100‑1 (MQ=255)
    ccTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATcg                                                                                                 >  2:1483918/1‑100 (MQ=255)
     cTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATcgc                                                                                                >  2:2169703/1‑100 (MQ=255)
       aaaaGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCgcgc                                                                                              <  1:497875/100‑1 (MQ=255)
       aaaaGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCgcgc                                                                                              >  2:2821864/1‑100 (MQ=255)
       aaaaGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCgcgc                                                                                              >  1:2143967/1‑100 (MQ=255)
       aaaaGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCgcgc                                                                                              <  2:2428182/100‑1 (MQ=255)
           ggTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGAGGTGAGCGCGATCATCACGCTCAc                                                                                          <  1:1356597/100‑1 (MQ=255)
              ttAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTg                                                                                       <  2:1413571/100‑1 (MQ=255)
                aaCATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTgcg                                                                                     >  1:2420336/1‑100 (MQ=255)
                aaCATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCACTCACCTgcg                                                                                     <  1:1176884/100‑1 (MQ=255)
                     cGGTCTGGCGAAGAAATCGGGCCAGGCCAGTTATAAACGATCTGAAAGAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGctgcc                                                                                <  1:1969932/100‑1 (MQ=255)
                     cGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGctgcc                                                                                <  2:942748/100‑1 (MQ=255)
                       gTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGctgcctg                                                                              >  2:1740636/1‑100 (MQ=255)
                           ggCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAg                                                                          >  2:2377136/1‑100 (MQ=255)
                           ggCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAg                                                                          >  1:962877/1‑100 (MQ=255)
                           ggCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAg                                                                          <  2:120952/100‑1 (MQ=255)
                           ggCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAg                                                                          >  1:2292392/1‑100 (MQ=255)
                               aagaaATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAg                                                                      <  1:1728716/100‑1 (MQ=255)
                               aagaaATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAg                                                                      <  2:1296126/100‑1 (MQ=255)
                                agaaATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGc                                                                     <  2:478146/100‑1 (MQ=255)
                                 gaaATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGcc                                                                    >  2:1992237/1‑100 (MQ=255)
                                  aaaTCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCa                                                                   <  2:2844298/100‑1 (MQ=255)
                                  aaaTAGGGCCAGGCCAGTGATAAAAGTTCTGTAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCTCCTGCGCTGGCAGGCAGCgcagg                                                                    <  1:854103/99‑5 (MQ=255)
                                   aaTCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCAt                                                                  >  2:832251/1‑100 (MQ=255)
                                       gggCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCt                                                              >  2:766468/1‑100 (MQ=255)
                                                cagTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAgg                                                     >  2:2696294/1‑100 (MQ=255)
                                                cagTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAgg                                                     >  2:1267931/1‑100 (MQ=255)
                                                  gTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGTGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGgc                                                   <  2:2050308/100‑1 (MQ=255)
                                                   tGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGgcg                                                  <  1:1639447/100‑1 (MQ=255)
                                                   tGATAAAAGATCTCAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGgcg                                                  >  2:1982167/1‑100 (MQ=255)
                                                    gATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCAACTGCGCTGCCAGACAGCGAGCCATCGCTTCCGCCAg                                                      >  1:2906263/1‑95 (MQ=255)
                                                     aTAAAAGGTCTGAAAACAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGcc                                                <  1:2802173/100‑1 (MQ=255)
                                                     aTAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGcc                                                >  2:1907835/1‑100 (MQ=255)
                                                     aTAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGcc                                                <  1:9964/100‑1 (MQ=255)
                                                     aTAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGcc                                                >  1:1503735/1‑100 (MQ=255)
                                                     aTAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGcc                                                <  2:1834832/100‑1 (MQ=255)
                                                      tAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCa                                               <  1:1811505/100‑1 (MQ=255)
                                                       aaaaGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCat                                              >  2:553202/1‑100 (MQ=255)
                                                         aaGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCatat                                            <  2:2423091/100‑1 (MQ=255)
                                                          aGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATc                                           <  1:2811273/100‑1 (MQ=255)
                                                              cTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTAc                                       <  1:325197/100‑1 (MQ=255)
                                                               tGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTAcc                                      <  2:1567942/100‑1 (MQ=255)
                                                                gAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCAGCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTAcc                                      <  1:2370510/99‑1 (MQ=255)
                                                                gAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCa                                     >  2:2182670/1‑100 (MQ=255)
                                                                gAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCa                                     >  2:1899623/1‑100 (MQ=255)
                                                                gAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCAGCGCTTCCCCCAGGGCGCCATATCCt                                         >  1:283309/1‑96 (MQ=255)
                                                                 aaaaaaaCAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTAc                                       <  1:1901134/97‑1 (MQ=255)
                                                                 aaaaaaaCAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCaa                                    >  1:628839/1‑100 (MQ=255)
                                                                  aaaaaaCAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAg                                   <  1:1483912/100‑1 (MQ=255)
                                                                   aaaaaCAGCGTGAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAgg                                  <  1:1405438/100‑1 (MQ=255)
                                                                    aaaaCAGCGTAAACAGTGGGGTGAGCGTGAGGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGt                                 <  2:1761234/100‑1 (MQ=255)
                                                                      aaCAGCGTACCCAGTGGGTTGAGCGTGATGATCGCGCTAACCCGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTAt                               <  1:2863279/100‑1 (MQ=255)
                                                                      aaCAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTAt                               <  2:1329603/100‑1 (MQ=255)
                                                                      aaCAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTAt                               >  2:941130/1‑100 (MQ=255)
                                                                       aCAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTAtt                              >  1:1851335/1‑100 (MQ=255)
                                                                       aCAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTAtt                              <  2:2638016/100‑1 (MQ=255)
                                                                       aCAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTAtt                              >  1:219281/1‑100 (MQ=255)
                                                                        cAGCGTAAACAGTGGGGTCAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTACCGCCAGGGCGCCATATCCTACCAAGGTATTc                             <  2:209207/100‑1 (MQ=255)
                                                                          gCGTAAACAGTTGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAg                           <  2:127957/100‑1 (MQ=255)
                                                                          gCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAg                           <  2:2584145/100‑1 (MQ=255)
                                                                          gCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAg                           <  1:267161/100‑1 (MQ=255)
                                                                             tAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTcc                        <  1:1893396/100‑1 (MQ=255)
                                                                             tAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTcc                        >  1:262421/1‑100 (MQ=255)
                                                                             tAAACAGTGGGGTGAGCGTGATGATCGCGCGCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTcc                        <  1:2468073/100‑1 (MQ=255)
                                                                                 cAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCaa                    >  2:144578/1‑100 (MQ=255)
                                                                                 agcgcaggTGAGCGCGATCATCACGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCaa                    >  1:2082572/7‑100 (MQ=38)
                                                                                  gcgcaggTGAGCGCGATCATCACGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCaaa                   >  2:844806/6‑100 (MQ=38)
                                                                                  gcgcaggTGAGCGCGATCATCACGCTCACCCGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCaaa                   >  2:2243623/6‑100 (MQ=37)
                                                                                  aGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCaaa                   <  2:2245987/100‑1 (MQ=255)
                                                                                   ggggggTAGGCGTGTCGATCTCGCTCACCTGCGCTGCCGGCTAGCGAGGCATCGACTCCGCGAGGGAGCCACATCCTAGCAAGTTATTCAGTcc                        <  1:2191674/92‑1 (MQ=255)
                                                                                   cgcaggTGAGCGCGATCATCACGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCaaaa                  >  1:803190/5‑100 (MQ=38)
                                                                                   cgcaggTGAGCGCGATCATCACGCTCACCCGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCaaaa                  >  1:166262/5‑100 (MQ=37)
                                                                                    gcaggagAGCGCGATCATCACGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCGTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCacga                  >  1:2248335/7‑96 (MQ=38)
                                                                                    gcaggTGAGCGCGATCATCACGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTGTTCAGTCCGCaaaa                  >  2:2867605/4‑99 (MQ=38)
                                                                                     ggggTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGTCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCaaaaaa                >  1:2648350/1‑100 (MQ=255)
                                                                                     caggTGAGCGCGATCATCACGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCaaaaaa                >  2:1181959/3‑100 (MQ=38)
                                                                                     caggTGAGCGCGATCATCACGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCaaaaaa                >  1:1495260/3‑100 (MQ=39)
                                                                                     caggTGAGCGCGATCATCACGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCaaaaaa                >  1:2089989/3‑100 (MQ=38)
                                                                                      gggTGAGCGTGATGATCGCCCTATCCTTCGCTTCCTGCCAGCGAGCCACCCCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCaaaaaa                <  1:1794818/99‑1 (MQ=255)
                                                                                      aggTGAGCGCGATCATCACGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAAt               >  2:613578/2‑100 (MQ=39)
                                                                                      aggTGAGCGCGATCATCACGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAAt               >  2:1537353/2‑100 (MQ=39)
                                                                                       ggTGAGCGTGATCATCGCGCTCAGCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAAtt              >  1:1214573/1‑100 (MQ=255)
                                                                                       ggTGAGCGCGATCATCACGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAAtt              >  2:2492212/1‑100 (MQ=255)
                                                                                       ggTGAGCGCGATCATCACGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAAtt              >  1:1164598/1‑100 (MQ=255)
                                                                                       ggTGAGCGCGATCATCACGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAAtt              >  2:1748974/1‑100 (MQ=255)
                                                                                       ggTGAGCGCGATCATCACGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAAt               >  2:1698092/1‑99 (MQ=39)
                                                                                          gAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTaaa           <  1:1015042/100‑1 (MQ=255)
                                                                                           aGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAAc          <  1:1267926/100‑1 (MQ=255)
                                                                                           aGCGTGATGATCGCGCTCACCTGCCCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAAc          <  2:1892477/100‑1 (MQ=255)
                                                                                             cGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACAt        >  2:391172/1‑100 (MQ=255)
                                                                                                gatgatCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCg     <  2:2601896/100‑1 (MQ=255)
                                                                                                  tgatCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAg   <  2:417849/100‑1 (MQ=255)
                                                                                                  tgatCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAg   >  1:269169/1‑100 (MQ=255)
                                                                                                  tgatCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAg   >  2:2667759/1‑100 (MQ=255)
                                                                                                   gatCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAg   >  1:338286/1‑99 (MQ=255)
                                                                                                   gatCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGc  >  2:1737963/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
ATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGC  >  NC_000913/3331884‑3332082

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: