Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,537,741 G→A 9.9% intergenic (+27/‑6) yfeK → / → yfeS uncharacterized protein/WGR domain protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,537,7410GA 9.9% 242.4 / 3.4 93intergenic (+27/‑6)yfeK/yfeSuncharacterized protein/WGR domain protein
Reads supporting (aligned to +/- strand):  ref base G (43/40);  new base A (2/8);  total (45/48)
Fisher's exact test for biased strand distribution p-value = 9.21e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

ATCGATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCCTTTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTAGGTAGTATTGGTAAATTCCAGACAAAAGAGTTCGATAATGAAGAACA  >  NC_000913/2537594‑2537889
                                                                                                                                                   |                                                                                                                                                    
aTCGATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCCTTTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACaggagg                                                                                                                                                     >  2:458303/1‑149 (MQ=255)
 tCGATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCCTTTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGGAATTAATCCGCGTTGGGAAGGCATTTACAGGAGGt                                                                                                                                                    >  1:1735422/1‑149 (MQ=255)
     ttACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCCTTTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACa                                                                                                                                                <  2:1263278/149‑1 (MQ=255)
      tACCGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCCTTTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACAt                                                                                                                                               <  1:707593/149‑1 (MQ=255)
             aaGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCCTTTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGaaaaaa                                                                                                                                        <  1:1816376/149‑1 (MQ=255)
                ccGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCCTTTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGc                                                                                                                                     >  1:1021950/1‑149 (MQ=255)
                   tataTTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCCTTTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCttt                                                                                                                                  >  1:871336/1‑149 (MQ=255)
                          tgAAGATCCCCGGTAAGAGAGAAAAGAAAACAAAGCCTTTTTTACATTCGTTAAATGCGAAGACGGAAAAAACGGAACCTGCGGAACGAAATTAATCCGCTTTGGGAAGGCATTTACAGAAAAAAAAAAAAAAAAACGCTTTATTTATc                                                                                                                           <  2:693370/149‑1 (MQ=255)
                           gAAGCTCCCCGGTAAGAGCGATGAGAACGCACAGCCTTTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGCTTGGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCa                                                                                                                          <  1:1710305/149‑1 (MQ=255)
                           gAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCCTTTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCAGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCa                                                                                                                          <  1:347043/149‑1 (MQ=255)
                            aaGATCCCCGGTAAGAGCGATGAGAACGCACAGCCTTTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCAc                                                                                                                         >  1:631644/1‑149 (MQ=255)
                            aaGATCCCCGGTAAGAGCGATGAGAACGCACAGCCTTTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCAc                                                                                                                         >  1:632137/1‑149 (MQ=255)
                              gATCCCCGGTAAGAGCGATGAGAACGCACAGCCTTTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGa                                                                                                                       <  1:652744/149‑1 (MQ=255)
                              gATCCCCGGTAAGAGCGATGAGAACGCACAGCCTTTTTTACATGCGTTAATTGCACAGACGGATAAAACGATGCCTGCGGAAGGAAATTAATCCGCTGTGGAGAGGCATTTACAGGAGGTAACATGAAAAATAGCTTTATTTATCACGa                                                                                                                       >  2:119350/1‑149 (MQ=255)
                               aTCCCCGGTAAGAGCGATGAGAACGCACAGCCTTTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGAt                                                                                                                      >  1:756461/1‑149 (MQ=255)
                                tCCCCGGTAAGAGCGATGAGAACGCACAGCCTTTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATg                                                                                                                     >  1:1370423/1‑149 (MQ=255)
                                       aaaGAGCGATTAGAACGCACAGGCTTTTGTACATGCAGTAAGTGCGCAGACGGAAAAAACGGTGCCTCCGGAAGGAAATTAACCCGCTTTGGGAAGGCATTTACAGGAAAAAAAATGAAAAAACGCTTTATTTATCACGATGAAAAATc                                                                                                              <  2:27053/148‑1 (MQ=255)
                                        agaagCGATGAGAACGCATAGGCTGTTCTACATTGGTTAACTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCg                                                                                                             <  2:1414264/146‑1 (MQ=255)
                                         aaaacaaTGAGAACTCACAGCCGTTTTTACAAGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAAAAAAAAAAAAAAAACGCTTTATTTATCACGATGAAAAATCGa                                                                                                            <  2:409105/143‑1 (MQ=255)
                                          gagCGATGAGAACGCACAGCCTTTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGaa                                                                                                           <  2:1531238/149‑1 (MQ=255)
                                           agCGATGAGAACGCACAGCCTTTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAAt                                                                                                          >  1:1046830/1‑149 (MQ=255)
                                             cGATGAGAACGCACAGCCTTTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGGGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCCCGATGAAAAATCGAATaa                                                                                                        >  2:1558438/1‑149 (MQ=255)
                                                 gggaACGCACAGCCTTTTTTACATGCGTTAATTGCGCAGACGGATAAACCGGTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAAttt                                                                                                    <  2:1735067/147‑1 (MQ=255)
                                                 gagaACGCACAGCCTTTTTTACATTCGTTAATTGCGCAGCAGGCTAAAACGGTGCCTGCGGAAGAAAATTGATCCGCTTTGGGAAGGCATTGACAGGAGGTAACATGAAAAAACGCTTTTTTTATCACGATGAAAAATCGAATAAAttt                                                                                                    <  2:523208/149‑1 (MQ=255)
                                                 gagaACGCACAGCCTTTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAAttt                                                                                                    <  2:1735423/149‑1 (MQ=255)
                                                 gagaACGCACAGCCTTTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGCCAATTAAATGGGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAAttt                                                                                                    <  1:800762/149‑1 (MQ=255)
                                                  agaACGCACAGCCTTTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAAtttt                                                                                                   <  2:1217648/149‑1 (MQ=255)
                                                   gaACGCACAGCCTTTTTTACATGTGTTAATTGCGCAGAAGGAAAAAACGGTTCCTGCGGAAGGTAGTTAATCTGCTTTGGGCAGGCATTTACAGGAAAAAAAAAGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTtg                                                                                                  <  2:1645146/149‑1 (MQ=255)
                                                    aaCGCACAGCCTTTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTtgg                                                                                                 <  1:590129/149‑1 (MQ=255)
                                                    aaCGCACAGCCTTTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTACGCGAGTAAATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGACTAAATTTtgg                                                                                                 >  1:799640/1‑149 (MQ=255)
                                                     aCGCACAGCCTTTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTtggt                                                                                                <  2:1735422/149‑1 (MQ=255)
                                                     aCGCACAGCCTTTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTtggt                                                                                                <  2:366544/149‑1 (MQ=255)
                                                         acaGCCTTTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCCGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGGAt                                                                                            >  1:117874/1‑149 (MQ=255)
                                                         acaGCCTTTTTTACATGCGTTAATAGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAAGCGAATAAATTTTGGTGGAt                                                                                            >  1:1344337/1‑149 (MQ=255)
                                                          caGCCTTTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGGATa                                                                                           <  1:578479/149‑1 (MQ=255)
                                                           aGCCTTTTTGACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGGATAg                                                                                          <  2:1218980/149‑1 (MQ=255)
                                                               ttttttACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGGATAGATTa                                                                                      <  1:1659499/149‑1 (MQ=255)
                                                                    tccATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGGATAGATTACGAAg                                                                                 <  2:408474/147‑1 (MQ=255)
                                                                      caggCGTTAATTGCGCAGACGGATAAAACTGTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGGATAGATTACGAAggg                                                                               <  2:1148656/146‑1 (MQ=255)
                                                                       agacGTTAATTGCGCAGACGGAAAAAACGGTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAAAAAAAAAGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGGATAGATTACGAAgggg                                                                              <  2:556690/146‑1 (MQ=255)
                                                                            ttAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGt                                                                         >  1:1694048/1‑149 (MQ=255)
                                                                            ttAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGCGGTGGGAAGGAATTTACAGGAAGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTGGGTGGATAGATTACGAAGGGGATAGt                                                                         >  1:608749/1‑149 (MQ=255)
                                                                              aaTTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGttt                                                                       >  1:1646968/1‑149 (MQ=255)
                                                                               aTTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTAGGTGGATAGATTACGAAGGGGATAGTTTa                                                                      >  2:1030517/1‑149 (MQ=255)
                                                                                ttGCCCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAg                                                                     <  2:1081425/149‑1 (MQ=255)
                                                                                    gcAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAACAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGt                                                                 <  1:1275169/149‑1 (MQ=255)
                                                                                          ggATAAAACTGGGCCTGCGGAAAGAAATTAATCCCTCTTGCGAAGGCATTTACAGAAAAAAACAAAAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTa                                                           <  2:1062441/149‑1 (MQ=255)
                                                                                          ggATAAAACGGTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTa                                                           <  2:114362/149‑1 (MQ=255)
                                                                                            aTAAAACGGTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGGAAAAACGCTTTATTTATCACGATGAAAAGTCGAATAAATTTTGGTGGATAGATTGCGAAGGGGATAGTTTAGCTGTCAACTaaa                                                         >  1:2655/1‑147 (MQ=255)
                                                                                            aTAAAACGGTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATg                                                         >  1:1647135/1‑149 (MQ=255)
                                                                                            aTAAAACGGTGCCTCGGGTAGGAAATTAATCCGCTTTGGGAAGGAATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATg                                                         <  2:586213/149‑1 (MQ=255)
                                                                                            aTAAAACGGGGCCAGCGGAAGGAAATTAACCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATg                                                         <  2:484937/149‑1 (MQ=255)
                                                                                             tAAAACGGTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGCAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGATGTCAACAATgg                                                        >  2:357480/1‑149 (MQ=255)
                                                                                                aaCGGTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTCCAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGAATCGCTTCCGAAGGGGATAGTTTAGCTGTCAACTATGGCac                                                     >  2:127398/1‑148 (MQ=255)
                                                                                                aaCGGTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCaa                                                     >  2:127341/1‑149 (MQ=255)
                                                                                                   ggTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAaggt                                                  <  1:822138/149‑1 (MQ=255)
                                                                                                     ggCCAGCGGCAAGAAATTACTCCACCATGGGAAAACCTATCCAGGAAACAACGTAGAAAGACGGTTTAGTTATCACGATGAAGAATCGAATTAATTTTGGTGGATAGATTGCGAAGGGGATAGTTTAGCTGTCAACTATGGCAaggtag                                                <  2:1055549/148‑1 (MQ=255)
                                                                                                        cTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAggtaggta                                             <  2:1588295/149‑1 (MQ=255)
                                                                                                        cTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAggtaggta                                             >  1:400954/1‑149 (MQ=255)
                                                                                                         tGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGgaaggtag                                            >  2:492756/1‑149 (MQ=255)
                                                                                                            ggaaggaaATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTAGGTAGTAt                                         <  1:184455/149‑1 (MQ=255)
                                                                                                              aaggaaATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTAGGTAGTATTg                                       <  1:837734/149‑1 (MQ=255)
                                                                                                                ggaaATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTAGGTAGTATTGGt                                     <  1:820788/149‑1 (MQ=255)
                                                                                                                 gaaATTAATCCGCTTTGGGAAGGCATTTACAGAAGGTAACATGAAAAAACACTTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGGATAGAGTTCGAAGGCGATAGTTTAGCTGTCAACTCTGGAAAGGTGAGTAGTATTGGTa                                    >  1:166685/1‑149 (MQ=255)
                                                                                                                    gttAATCCGCTTTGGGAAGGCATTTACAGGAAGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTAGGTAGTATTGGTAAAt                                 >  2:1326672/2‑149 (MQ=255)
                                                                                                                    aTTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTAGGTAGTATTGGTAAAt                                 <  1:1677788/149‑1 (MQ=255)
                                                                                                                    aTTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTAGGTAGTATTGGTAAAt                                 >  1:1323457/1‑149 (MQ=255)
                                                                                                                      tAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTAGGTAGTATTGGTAAATTc                               <  1:939462/149‑1 (MQ=255)
                                                                                                                        tttcgcgTTGGGGTGGCATTTACAGGAGGTAACATGAGAAAGCGCTTTATTTATCACGCTGAAAAATCCAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTAGGTAGTATTGGTAAATTccc                             <  2:1243903/146‑2 (MQ=255)
                                                                                                                        aTCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTAGGTAGTATTGGTAAATTCCa                             <  1:571484/149‑1 (MQ=255)
                                                                                                                        aTCCGCTTTGGAAAGGCATTTACCGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAACCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTAGGTAGTATTGGTAAATTCCa                             <  2:1243628/149‑1 (MQ=255)
                                                                                                                         tCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGGATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTTGGGAGTATTGGTAAATTCCAg                            >  2:862558/1‑149 (MQ=255)
                                                                                                                         tCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTAGGTAGTATTGGTAAATTCCAg                            >  1:10302/1‑149 (MQ=255)
                                                                                                                         tCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTTGGGGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTAGGGAGTATTAGTAAGTTTCAg                            >  1:928425/1‑149 (MQ=255)
                                                                                                                         tCCGCTTTGGGAAGGCATTTACACGAGGGAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGACAACTATGGCAAGGTAGGTAGTATTGGTAAATTcccg                            >  1:1183254/1‑147 (MQ=255)
                                                                                                                          ccGCTTTGGGAAGGCATTTACAGGAGGTAACATGGAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCGAGGTAGGTAGTATTGGTAAATTACAGa                           >  1:368253/1‑149 (MQ=255)
                                                                                                                          ccGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTAGGCAGTATTGGTAAATTCCAGa                           >  2:447572/1‑149 (MQ=255)
                                                                                                                            gCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTAGGTAGTATTGGTAAATTCCAGACa                         <  2:1370423/149‑1 (MQ=255)
                                                                                                                            gCTTGGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATGTTGGTGGATGGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGGAGGTAGTATAGGTAAATTCCAGACa                         >  2:1342245/1‑149 (MQ=255)
                                                                                                                             gtgtGGGAAGGGATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTAGGTAGTATTGGTAAATTCCAGACaa                        <  1:1715488/146‑1 (MQ=255)
                                                                                                                              tttGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCGAGGGAGGTAGTATTGGTAAATTCCAGACaaa                       >  2:22195/1‑149 (MQ=255)
                                                                                                                                 gggAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTAGGTTGTATTGGTAAATTCCAGACAAAaga                    <  1:1009663/149‑1 (MQ=255)
                                                                                                                                  ggAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTAGGTAGTATTGGTAAATTCCAGACAAAagag                   <  1:1561753/149‑1 (MQ=255)
                                                                                                                                   gAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTCAACACGAAGAAAAATCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTAGGTAGTATTGGTAAATTCCAGACAAAAGAGt                  >  1:709974/1‑149 (MQ=255)
                                                                                                                                    aaGGCATTTAGAGGAGGACACGAGTATAAACGCTTTATTTATCACGATGAAAACCCGAATGAATGTTGGGGGACAGGATAAGACGGGGGTAGGTCGGGTAGAAGATAAGGCAGGGTAGGGAGTAGTGGCAAAGTCCAGCAAACagaggt                 >  2:785146/1‑147 (MQ=255)
                                                                                                                                    aaGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTAGGTAGTATTGGTAAATTCCAGACAACAAAGtt                 >  2:1079592/1‑149 (MQ=255)
                                                                                                                                    aaGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTAGGTAGTATTGGTAAATTCCAGACAAAAGAGtt                 >  1:1332090/1‑149 (MQ=255)
                                                                                                                                       gCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTAGGTAGTATTGGTAAATTCCAGACAAAAGAGTTCGa              >  1:1460200/1‑149 (MQ=255)
                                                                                                                                       gCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTAGGTAGTATTGGTAAATTCCAGACAAAAGAGTTCGa              >  1:836655/1‑149 (MQ=255)
                                                                                                                                       gCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTAGGTAGTATTGGGAAATTCCAGCCCAAAGGGTTCGa              >  1:1367792/1‑149 (MQ=255)
                                                                                                                                       gCATTTACAGGAGGGAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAGATAACAGTAGATAGATCACGTTGTGGATAGTTTAGCTGTCAACTATGGAAAGGTAGGTAGTATTGGTAAATTCCAGACAAAAGAGGTCGa              >  1:427823/1‑149 (MQ=255)
                                                                                                                                        cATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTAGGTAGTATTGGTAAATTCCAGACAAAAGAGTTCGAt             >  1:532696/1‑149 (MQ=255)
                                                                                                                                          tttACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTAGGTAGTATTGGTAAATTCCAGACAAAAGAGTTCGATaa           >  2:1161227/1‑149 (MQ=255)
                                                                                                                                            aaCAGGAAGAAAAAAGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTAGGTAGTATTGGTAAATTCCAGACAAAAGAGTTCGATAATg         <  2:1536448/148‑1 (MQ=255)
                                                                                                                                                   ggTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTAGGTAGTATTGGTAAATTCCAGACAAAAGAGTTCGATAATGAAGAACa  <  2:1262771/149‑1 (MQ=255)
                                                                                                                                                   ggTAACATGAAAAAACGCGTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTAGGTAGTATTGGTAAATTCCAGACAAAAGAGTTCGATAATGAAGAACa  <  1:686651/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ATCGATTACTGGGAAGCCGTATATTGTGAAGATCCCCGGTAAGAGCGATGAGAACGCACAGCCTTTTTTACATGCGTTAATTGCGCAGACGGATAAAACGGTGCCTGCGGAAGGAAATTAATCCGCTTTGGGAAGGCATTTACAGGAGGTAACATGAAAAAACGCTTTATTTATCACGATGAAAAATCGAATAAATTTTGGTGGATAGATTACGAAGGGGATAGTTTAGCTGTCAACTATGGCAAGGTAGGTAGTATTGGTAAATTCCAGACAAAAGAGTTCGATAATGAAGAACA  >  NC_000913/2537594‑2537889

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: