Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,063,315 A→G 13.6% S357G (AGC→GGC)  yihN → MFS transporter family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,063,3150AG13.6% 2151.0 / 261.4 804S357G (AGC→GGC) yihNMFS transporter family protein
Reads supporting (aligned to +/- strand):  ref base A (297/397);  new base G (54/56);  total (351/453)
Fisher's exact test for biased strand distribution p-value = 2.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 820 total aligned reads displayed.

GCCGCACGAAAGTATGCCGGTCTATTTAGGGATGGCATGTACGCTGGGCTTTGGCGCGATAGTCTTTACACAGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGATTTAAATCCGGGGATTATTGGCTACAAAATCGTGTTTGGCATTATGGCCTGCTTCGCATTCAGTGGTGCGGTGGTTTCCGTAATGCTGGT  >  NC_000913/4063167‑4063463
                                                                                                                                                    |                                                                                                                                                    
gcCGCACGAAAGTATGCCGGTCTATTTAGGGATGGCATGTACGCTGGGCTTTGGCGCGATAGTCTTTACACAGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTa                                                                                                                                                      <  4:3060716/149‑1 (MQ=255)
 ccGCACGAAAGTATGCCGGTCTATTTAGGGATGGCATGTACGCTGGGCTTTGGCGCGATAGTCTTTACACAGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAg                                                                                                                                                     >  2:1843889/1‑149 (MQ=255)
 ccGCACGAAAGTATGCCGGTCTATTTAGGGATGGCATGTACGCTGGGCTTTGGCGCGATAGTCTTTACACAGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAg                                                                                                                                                     <  2:7081992/149‑1 (MQ=255)
       gAAAGTATGCCGGTCTATTTAGGGATGGCATGTACGCTGGGCTTTGGCGCGATAGTCTTTACACAGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTGGCTTTAt                                                                                                                                               >  3:3389654/1‑149 (MQ=255)
       gAAAGTATGCCGGTCTATTTAGGGATGGCATGTACGCTGGGCTTTGGCGCGATAGTCTTTACACAGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTAt                                                                                                                                               >  3:827959/1‑149 (MQ=255)
       gAAAGTATGCCGGTCTATTTAGGGATGGCATGTACGCTGGGCTTTGGCGCGATAGTCTTTACACAGCGAGCCGTAGTTTTTGCACCTATCGTCGAAGCAAAAATTGCGGAAAATAAAACAGGTGCGGCGATGGCGTGGGGTAGcagtat                                                                                                                                               >  2:1973066/1‑144 (MQ=255)
        aaaGTATGCCGGTCTATTTAGGGATGGCATGTACGCTGGGCTTTGGCGCGATAGTCTTTACACAGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTAtt                                                                                                                                              <  1:7172388/149‑1 (MQ=255)
         aaGTATGCCGGTCTATTTAGGGATGGCATGTACGCTGGGCTTTGGCGCGATAGTCTTTACACAGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTg                                                                                                                                             <  1:4106893/149‑1 (MQ=255)
          aGTATGCCGGTCTATTTAGGGATGGCATGTACGCTGGGCTTTGGCGCGATAGTCTTTACACAGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTgg                                                                                                                                            <  4:1289748/149‑1 (MQ=255)
          aGTATGCCGGTCTATTTAGGGATGGCATGTACGCTGGGCTTTGGCGCGATAGTCTTTACACAGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTgg                                                                                                                                            <  3:3667757/149‑1 (MQ=255)
          aGTATGCCGGTCTATTTAGGGATGGCATGTACGCTGGGCTTTGGCGCGATAGTCTTTACACAGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTgg                                                                                                                                            <  1:606916/149‑1 (MQ=255)
          aGTATGCCGGTCTATTTAGGGATGGCATGTACGCTGGGCTTTGGCGCGATAGTCTTTACACAGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTgg                                                                                                                                            <  3:7017099/149‑1 (MQ=255)
            tATGCCGGTCTTTTTAGGGATGGCATGTACTCTGGGCTTTGGCGCGATAGTCTTTACACAGCGAGCCGTATTTTTTGCACCGACCGGCGAAGCAAAAATTGCTGAAAATACAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGtt                                                                                                                                          <  2:112281/149‑1 (MQ=255)
            tATGCCGGTCTATTTAGGGATGGCATGTACGCTGGGCTTTGGCGCGATAGTCTTTACACAGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGtt                                                                                                                                          <  3:1815633/149‑1 (MQ=255)
             aTGCCGGTCTATTTAGGGATGGCATGTACGCTGGGCTTTGGCGCGATAGTCTTTACACAGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTa                                                                                                                                         <  1:1775687/149‑1 (MQ=255)
                  ggTCTATTTAGGGATGGCATGTACGCTGGGCTTTGGCGCGATAGTCTTTACACAGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAGAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTc                                                                                                                                    >  4:6111624/1‑149 (MQ=255)
                    tCTATTTAGGGATGGCATGTACGCTGGGCTTTGGCGCGATAGTCTTTACACAGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTGGCTTTATTGGTTACGCTCCg                                                                                                                                  >  2:5998742/1‑149 (MQ=255)
                      tATTTAGGGATGGCATGTACGCTGGGCTTTGGCGCGATAGTCTTTACACAGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGt                                                                                                                                           <  1:3657264/138‑1 (MQ=255)
                      tATTTAGGGATGGCATGTACGCTGGGCTTTGGCGCGATAGTCTTTACACAGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGt                                                                                                                                           >  2:3657264/1‑138 (MQ=255)
                       aTTTAGGGATGGCATGTACGCTGGGCTTTGGCGCGATAGTCTTTACACAGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCg                                                                                                                               <  3:2163242/149‑1 (MQ=255)
                        tttAGGGATGGCATGTACGCTGGGCTTTGGCGCGATAGTCTTTACACAGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGa                                                                                                                              >  4:3552173/1‑149 (MQ=255)
                             ggATGGCATGTACGCTGGGCTTTGGCGCGATAGTCTTTACACAGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTc                                                                                                                         <  1:3206168/149‑1 (MQ=255)
                             ggATGGCATGTACGCTGGGCTTTGGCGCGATAGTCTTTACACAGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTc                                                                                                                         >  3:2702655/1‑149 (MQ=255)
                              gATGGCATGTACGCTGGGCTTTGGCGCGATAGTCTTTACACAGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCt                                                                                                                        <  3:5401172/149‑1 (MQ=255)
                                tGGCATGTACGCTGGGCTTTGGCGCGATAGTCTTTACACAGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGc                                                                                                                      <  3:1540067/149‑1 (MQ=255)
                                    aTGTACGCTGGGCTTTGGCGCGATAGTCTTTACACAGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCa                                                                                                                  >  2:143244/1‑149 (MQ=255)
                                    aTGTACGCTGGGCTTTGGCGCGATAGTCTTTACACAGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCa                                                                                                                  <  4:4177394/149‑1 (MQ=255)
                                          gCTGGGCTTTGGCGCGATAGTCTTTACACAGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGt                                                                                                            >  2:4131228/1‑149 (MQ=255)
                                           cTGGGCTTTGGCGCGATAGTCTTTACACAGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTa                                                                                                           >  1:5798617/1‑149 (MQ=255)
                                           cTGGGCTTTGGCGCGATAGTCTTTACACAGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTa                                                                                                           <  3:4705343/149‑1 (MQ=255)
                                           cTGGGCTTGGGCGCGATAGTCTTTACACAGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTGGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTa                                                                                                           <  4:4876954/149‑1 (MQ=255)
                                            tGGGCTTTGGCGCGATAGTCTTTACACAGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTAt                                                                                                          <  4:5628094/149‑1 (MQ=255)
                                              ggCTTTGGCGCGATAGTCTTTACACAGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATgg                                                                                                        >  2:5833744/1‑149 (MQ=255)
                                                      cgcgATAGTCTTTACACAGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTc                                                                                                <  3:6050918/149‑1 (MQ=255)
                                                         gATAGTCTTTACACAGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTa                                                                                                                                                      <  1:345766/92‑1 (MQ=255)
                                                          aTAGTCTTTACACAGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGa                                                                                            <  3:2933525/149‑1 (MQ=255)
                                                          aTAGTCTTTACACAGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGa                                                                                            <  1:6519057/149‑1 (MQ=255)
                                                          aTAGTCTTTACACAGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGa                                                                                            <  1:1246313/149‑1 (MQ=255)
                                                          aTAGTCTTTACAAAGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGa                                                                                            <  2:6635431/149‑1 (MQ=255)
                                                           tAGTCTTTACACAGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAGTTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGAt                                                                                           <  1:2258412/149‑1 (MQ=255)
                                                           tAGTCTTTACACAGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGAt                                                                                           <  1:2702515/149‑1 (MQ=255)
                                                            aGTCTTTACACAGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGAtt                                                                                          >  2:5067850/1‑149 (MQ=255)
                                                             gTCTTTACACAGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGAttt                                                                                         <  4:3366590/149‑1 (MQ=255)
                                                              tCTTTACACAGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTGGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGATTTa                                                                                        >  1:5421013/1‑149 (MQ=255)
                                                              tCTTTACACAGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGATTTa                                                                                        >  1:1001358/1‑149 (MQ=255)
                                                                tttACACAGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTAt                                                                                                                                               >  3:2372919/1‑92 (MQ=255)
                                                                   acacAGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGATTTAAATcc                                                                                   >  4:4469429/1‑149 (MQ=255)
                                                                     acaGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTg                                                                                                                                             >  3:5832071/1‑89 (MQ=255)
                                                                     acaGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGATTTAAATCCgg                                                                                 <  3:4956754/149‑1 (MQ=255)
                                                                     acaGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGATTTAAATCCgg                                                                                 >  2:3895467/1‑149 (MQ=255)
                                                                     acaGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGATTTAAATCCgg                                                                                 >  2:2728977/1‑149 (MQ=255)
                                                                     acaGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGAGATGGCGTTGGGTAGCTTTATTGGTTACGCTACGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGATTTAAATCCgg                                                                                 >  2:965704/1‑149 (MQ=255)
                                                                      caGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGATTTAAATCCggg                                                                                <  3:1245081/149‑1 (MQ=255)
                                                                      caGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGATTTAAATCCggg                                                                                <  1:6760361/149‑1 (MQ=255)
                                                                         cGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCtt                                                                                                                                                  <  2:4863055/80‑1 (MQ=255)
                                                                            gCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTAt                                                                                                                                               <  2:5724625/80‑1 (MQ=255)
                                                                             ccGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGATTTAAATCCGGGGattatt                                                                         <  2:820935/149‑1 (MQ=255)
                                                                              cGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTGGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGATTTAAATCCGGGGATTATTg                                                                        >  3:4235877/1‑149 (MQ=255)
                                                                              cGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGATTTAAATCCGGGGATTATTg                                                                        <  1:289874/149‑1 (MQ=255)
                                                                               gTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTGGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGATTTAAATCCGGGGATTATTgg                                                                       >  2:6348041/1‑149 (MQ=255)
                                                                               gTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGtt                                                                                                                                          >  2:6237006/1‑82 (MQ=255)
                                                                                  ttttttGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGATTTAAATCCGGGGATTATTGGCTa                                                                    <  4:4672289/149‑1 (MQ=255)
                                                                                       tGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGATTTAAATCCGGGGATTATTGGCTACaaaa                                                               <  2:2997808/149‑1 (MQ=255)
                                                                                       tGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGATTTAAATCCGGGGATTATTGGCTACaaaa                                                               >  1:835027/1‑149 (MQ=255)
                                                                                         cACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTAGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGATTTAAATCCGGGGATTATCGGCTACAAAATc                                                             <  4:6586537/149‑1 (MQ=255)
                                                                                          aCCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGATTTAAATCCGGGGATTATTGGCTACAAAATc                                                             <  4:2708032/148‑1 (MQ=255)
                                                                                          aCCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGATTTAAATCCGGGGATTATTGGCTACAAAATCg                                                            >  3:3154766/1‑149 (MQ=255)
                                                                                            cTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTa                                                                                                           >  1:5003464/1‑100 (MQ=255)
                                                                                            cTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGATTTAAATCCGGGGATTATTGGCTACAAAATCgtg                                                          >  1:5206870/1‑149 (MQ=255)
                                                                                            cTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGATTTAAATCCGGGGATTATTGGCTACAAAATCgtg                                                          >  1:4335984/1‑149 (MQ=255)
                                                                                             tATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGATTTAAATCCGGGGATTATTGGCTACAAAATCgtgt                                                         >  2:2512214/1‑149 (MQ=255)
                                                                                               tCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGATTTAAATCCGGGGATTATTGGCTACAAAATCGTGttt                                                       <  1:4493092/149‑1 (MQ=255)
                                                                                                cGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGATTTAAATCCGGGGATTATTGGCTACAAAATCGTGTTTg                                                      <  1:1935050/149‑1 (MQ=255)
                                                                                                  gCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGAt                                                                                           <  3:6721736/110‑1 (MQ=255)
                                                                                                  gCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGATTTAAATCCGGGGATTATTGGCTACAAAATCGTGTTTg                                                      <  4:3163093/147‑1 (MQ=255)
                                                                                                  gCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGATTTAAATCCGGGGATTATTGGCTACAAAATCGTGTTTGGc                                                    >  4:6977645/1‑149 (MQ=255)
                                                                                                     aaGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGATTTACATCCGGGGATTATTGGCTACAAAATCGTGTTTGGCAtt                                                 >  4:2368050/1‑149 (MQ=255)
                                                                                                      aGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGATTTAAATCCGGGGATTATTGGCTACAAAATCGTGTTTGGCATTa                                                <  2:4665357/149‑1 (MQ=255)
                                                                                                       gCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGAt                                                                                                                             <  3:4103372/71‑1 (MQ=255)
                                                                                                        cAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGATTTAAATCCGGGGATTATTGGCTACAAAATCGTGTTTGGCAtt                                                 >  2:5456779/1‑146 (MQ=255)
                                                                                                           aaaTTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGATTTAAATCCGGGGATTATTGGCTACAAAATCGTGTTTGGCATTATGGcc                                           <  4:4083228/149‑1 (MQ=255)
                                                                                                                  tGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGATTTAAATCCGGGGATTATTGGCTACAAAATCg                                                            >  2:4385423/1‑125 (MQ=255)
                                                                                                                  tGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGATTTAAATCCGGGGATTATTGGCTACAAAATCGTGTTTGGCATTATGGCCTGCTTCg                                    >  2:3390180/1‑149 (MQ=255)
                                                                                                                   gAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGTATTTAAATCCGGGGATTATTGGCTACAAAATCGTGTTTGGCATTATGGCCTGCTTCGc                                   <  1:6230065/149‑1 (MQ=255)
                                                                                                                    aaaaTAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGATTTAAATCCGGGGATTATTGGCTACAAAATCGTGTTTGGCATTATGGCCTGCTTCGCa                                  <  1:4925868/149‑1 (MQ=255)
                                                                                                                      aaTAAAACAGGTGCGGCGATGGCGTTGGGTGGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGATTTAAATCCGGGGATTATTGGCTACAAAATCGTGTTTGGCATTATGGCCTGCTTCGCAtt                                >  1:1493159/1‑149 (MQ=255)
                                                                                                                      aaTAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGATTTAAATCCGGGGATTATTGGCTACAAAATCGTGTTTGGCATTATGGCCTGCTTCGCAtt                                >  4:635281/1‑149 (MQ=255)
                                                                                                                       aTAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGATTTAAATCCGGGGATTATTGGCTACAAAATCGTGTTTGGCATTATGGCCTGCTTCGCATTc                               <  3:308019/149‑1 (MQ=255)
                                                                                                                        tAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGATTTAAATCCGGGGATTATTGGCTACAAAATCGTGTTTGGCATTATGGCCTGCTTCGCATTCa                              >  3:6535554/1‑149 (MQ=255)
                                                                                                                         aaaaCAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGTCTACATTCTGGATTTAAATCCGGGGTTCATTGGCTACAAAATCGTGATTGGCATTATGGCCTGCTTCGCATTCAg                             >  3:4475033/1‑149 (MQ=255)
                                                                                                                          aaaCAGGTGCGGCGATGGCGTTGGGTAGCTGTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGATTTAAATCCGGGGATTATTGGCTACAAAATCGTGTTTGGCATTATGGCCTGCTTCGCATTCAgt                            <  2:2154329/149‑1 (MQ=255)
                                                                                                                           aaCAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGATTTAAATCCGGGGATTATTGGCTACAAAATCGTGTTTGGCATTATGGCCTGCTTCGCATTCAgtg                           <  2:6898310/149‑1 (MQ=255)
                                                                                                                           aaCAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGATTTAAATCCGGGGATTATTGGCTACAAAATCGTGTTTGGCATTATGGCCTGCTTCGCATTCAgtg                           <  1:2128673/149‑1 (MQ=255)
                                                                                                                            aCAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGATTTAAATCCGGGGATTAGTGGCTACAAAATCGTGTTTGGCATTATGGCCTGCTTCGCATTCAgtgg                          >  2:5242169/1‑149 (MQ=255)
                                                                                                                              aGGTGCGGCGATGGCGTTGGGTGGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGAATTTTAATCCGGGGATTATTGGCTACAAAATCGTGTTTGGCATTATGGCCTGc                                        >  3:5170476/1‑133 (MQ=255)
                                                                                                                                  gcggcgATGGCGTTGGGTGGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGATTTAAATCCGGGGATTATTGGCTACAAAATCGTGTTTGGCATTATGGCCTGCTTCGCATTCAGTGGTGCggt                    <  4:5033283/149‑1 (MQ=255)
                                                                                                                                     gcgATGGCGTTGGGTAGCTTTATTGGTTACGCTCAGGTGATGTACTGATTCAGTCTTTATGGCTACATTCTGGATTTAAATCCGGGGATTATTGGCTACAAAATCGTGTTTGGCATTGTGGCCTGCTTCGCATTCAGTGCCGCggtggt                 <  4:2041619/149‑1 (MQ=255)
                                                                                                                                       gATGGCGTTGGGTGGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGATTTAAATCCGGGGATTATTGGCTACAAAATCGTGTTTGGCATTATGGCCTGCTTCGCATTCAGTGGTGCGGTGGttt               <  1:7079856/149‑1 (MQ=255)
                                                                                                                                             gTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGATTTAAATCCGGGGATTATTGGCTACAAAATCGTGTTTGGCATTATGGCCTGCTTCGCATTCAGTGGTGCGGTGGTTTCCGTaa         >  4:1745725/1‑149 (MQ=255)
                                                                                                                                                    aGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGATTTAAATCCGGGGATTATTGGCTACAAAATCGTGTTTGGCATTATGGCCTGCTTCGCATTCAGTGGTGCGGTGGTTTCCGTAATGCTGGt  >  4:5250622/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCCGCACGAAAGTATGCCGGTCTATTTAGGGATGGCATGTACGCTGGGCTTTGGCGCGATAGTCTTTACACAGCGAGCCGTATTTTTTGCACCTATCGGCGAAGCAAAAATTGCTGAAAATAAAACAGGTGCGGCGATGGCGTTGGGTAGCTTTATTGGTTACGCTCCGGCGATGTTCTGCTTCAGTCTGTATGGCTACATTCTGGATTTAAATCCGGGGATTATTGGCTACAAAATCGTGTTTGGCATTATGGCCTGCTTCGCATTCAGTGGTGCGGTGGTTTCCGTAATGCTGGT  >  NC_000913/4063167‑4063463

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: