Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,297,193 G→T 5.1% intergenic (‑68/+26) yjcO ← / ← fdhF Sel1 family TPR‑like repeat protein/formate dehydrogenase‑H, selenopolypeptide subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,297,1930GT 5.1% 419.5 / 7.3 142intergenic (‑68/+26)yjcO/fdhFSel1 family TPR‑like repeat protein/formate dehydrogenase‑H, selenopolypeptide subunit
Reads supporting (aligned to +/- strand):  ref base G (73/61);  new base T (5/3);  total (78/64)
Fisher's exact test for biased strand distribution p-value = 7.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 142 total aligned reads displayed.

TAAATACTGGCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGG  >  NC_000913/4297045‑4297341
                                                                                                                                                    |                                                                                                                                                    
 aaaTACTGGCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGaa                                                                                                                                                    <  1:1258142/150‑1 (MQ=255)
        ggCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTATCCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGc                                                                                                                                             >  2:895751/1‑150 (MQ=255)
         gCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATTAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCt                                                                                                                                            >  1:449662/1‑150 (MQ=255)
          ttGCCAGGCTCGGAGTCGTTGGCGTTGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGTGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTg                                                                                                                                           <  1:1280343/149‑1 (MQ=255)
             ccAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTAAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGTGTAGGGAGTAACCAGTATAATGAGTTCGAATGGATAAAAAAACAGACTCCGAAAGGAGGCTGTat                                                                                                                                        >  1:1510354/1‑149 (MQ=255)
              cAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCTTTCGGAGGCTGtttt                                                                                                                                       <  2:747061/150‑4 (MQ=255)
              cAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCTTTCGGAGGCTGtttt                                                                                                                                       <  1:1441122/150‑4 (MQ=255)
                 gCTCGGAGTCGTTGGCGGGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATATTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAg                                                                                                                                    <  2:1719210/150‑1 (MQ=255)
                  cTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCTTTCGGAGGCTGttttttta                                                                                                                                   <  2:40904/150‑8 (MQ=255)
                   tCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGa                                                                                                                                  <  1:1891195/150‑1 (MQ=255)
                     ggAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTTATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACg                                                                                                                                <  1:851070/150‑1 (MQ=255)
                      gAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACgg                                                                                                                               <  1:1058012/150‑1 (MQ=255)
                      gAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACgg                                                                                                                               <  2:1085571/150‑1 (MQ=255)
                      gAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACgg                                                                                                                               <  2:1545690/150‑1 (MQ=255)
                        gTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTa                                                                                                                             <  2:1050225/150‑1 (MQ=255)
                         tCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATTATTTTTTTCATTTTGATAATACGATAGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACACCCTCCGAAAGGAGGCTGTAGAAAGTACGGTAt                                                                                                                            >  1:256278/1‑150 (MQ=255)
                         tCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTAt                                                                                                                            <  1:988775/150‑1 (MQ=255)
                          cGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTAtt                                                                                                                           >  2:1904498/1‑150 (MQ=255)
                           gTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTa                                                                                                                          <  1:694075/150‑1 (MQ=255)
                           gTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTa                                                                                                                          >  2:129095/1‑150 (MQ=255)
                            ttGGCGTGGGCAAAGAATTTCAAAAACAACATCAATTCGATAATTTTTTTCATTTTGATAATACGATTGGGGAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTAc                                                                                                                         <  2:613585/150‑1 (MQ=255)
                              ggCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGc                                                                                                                       >  1:1081838/1‑150 (MQ=255)
                                 gTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAg                                                                                                                    >  2:209995/1‑150 (MQ=255)
                                   gggCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTg                                                                                                                  <  2:1519627/150‑1 (MQ=255)
                                    ggCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGc                                                                                                                 >  1:186386/1‑150 (MQ=255)
                                    ggCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGATGTAGAAAGGACGGGATTACGCCAGTGc                                                                                                                 >  2:345278/1‑150 (MQ=255)
                                     gCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGcc                                                                                                                >  2:1807051/1‑150 (MQ=255)
                                      cAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTTATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCg                                                                                                               >  2:482735/1‑150 (MQ=255)
                                      cAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCg                                                                                                               >  1:1721518/1‑150 (MQ=255)
                                      cAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAAAACGATTGGGTAGGGAGTAACCACTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTAATACGCCAGTGccc                                                                                                               >  2:1752407/1‑149 (MQ=255)
                                       aaaGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGc                                                                                                              >  1:154625/1‑150 (MQ=255)
                                       aaaGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGACAGTGCCGc                                                                                                              >  1:736259/1‑150 (MQ=255)
                                          gAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACTGTATTACGCCAGTGCCGCtta                                                                                                           >  1:505412/1‑149 (MQ=255)
                                           aaTGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTACCCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTcg                                                                                                          >  1:939479/1‑150 (MQ=255)
                                           aaTGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTcg                                                                                                          <  1:184551/150‑1 (MQ=255)
                                            aTGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTTGAAAGGACGGTATTACGCCAGTGCCGCTTcgc                                                                                                         >  2:133108/1‑150 (MQ=255)
                                             tGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTcgcg                                                                                                        <  1:721129/150‑1 (MQ=255)
                                              gTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCgcgc                                                                                                       >  2:512590/1‑150 (MQ=255)
                                               tCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCa                                                                                                      >  2:1537836/1‑150 (MQ=255)
                                               tCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCa                                                                                                      <  1:1540710/150‑1 (MQ=255)
                                               tCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCAGAATGGAGGCTGTAGAAAGTACGGTATTACGCCAGTGCAGCTTCGCGCa                                                                                                      >  2:65171/1‑150 (MQ=255)
                                                cAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAg                                                                                                     >  2:387975/1‑150 (MQ=255)
                                                cAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAAAACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAAAAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAg                                                                                                     >  2:833157/1‑150 (MQ=255)
                                                 aaaaaCAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAgg                                                                                                    <  1:1599888/150‑1 (MQ=255)
                                                 aaaaaCAACATCAATGCGATAATTTTTTTCATTTAGATAATACGATTGGGTAGGGAGTAACCAGTATATAGAGTTCGAATGGTTAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAgg                                                                                                    <  1:319930/150‑1 (MQ=255)
                                                   aaaCAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCg                                                                                                  <  2:807760/150‑1 (MQ=255)
                                                      caacaTCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCAGCTTAGCGCAGGCGAGt                                                                                               >  2:92857/1‑150 (MQ=255)
                                                        acaTCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGttt                                                                                             <  2:556914/150‑1 (MQ=255)
                                                         caTCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGtttt                                                                                            >  2:874576/1‑150 (MQ=255)
                                                           tccaTGCGATAATTTTTATCATTATGATAATACGATTGGGTCGGGAGTAACCAGTCTACAGAGTTCGAATCGATAAAAAAACAGCCTCCGAACGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCa                                                                                          <  2:287917/147‑1 (MQ=255)
                                                           tCAATTCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCa                                                                                          <  1:1115433/150‑1 (MQ=255)
                                                            cAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCaa                                                                                         >  2:1567824/1‑150 (MQ=255)
                                                             aaTGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAAc                                                                                        <  1:2451/150‑1 (MQ=255)
                                                              aTGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACt                                                                                       <  1:1223342/150‑1 (MQ=255)
                                                                gCGATAATTTTTTTCATTTTGATAATACGATTGGGGAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACttg                                                                                     <  2:1692416/150‑1 (MQ=255)
                                                                  taTAATTTTTTTCATTTTGATAATACGATTTGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACttgtt                                                                                   <  1:375522/149‑1 (MQ=255)
                                                                   aTAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACttgttg                                                                                  <  1:345547/150‑1 (MQ=255)
                                                                    tAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATATAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTtgttgt                                                                                 >  1:1376012/1‑150 (MQ=255)
                                                                    tAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAATGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTtgttgt                                                                                 >  1:320799/1‑150 (MQ=255)
                                                                      aTTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTAc                                                                               >  2:1815284/1‑150 (MQ=255)
                                                                       tttttttCATTTTGATTATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTTCCGCTTCGCGCACGCGAGTTTTCAACTTGTTGTAca                                                                              >  1:1185200/1‑149 (MQ=255)
                                                                           tttCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGATTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCACTGCCGCTAAGCGCAGGCGAGTTTTCCACATGTTGTACTCGTc                                                                          >  2:1180669/1‑150 (MQ=255)
                                                                           tttCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTc                                                                          >  1:1625155/1‑150 (MQ=255)
                                                                              cATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAAt                                                                       >  2:818458/1‑150 (MQ=255)
                                                                                ttttGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCa                                                                     >  1:767850/1‑150 (MQ=255)
                                                                                   tGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCTCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGt                                                                  >  1:158967/1‑150 (MQ=255)
                                                                                         taCGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCt                                                            <  1:619168/150‑1 (MQ=255)
                                                                                          aCGATTGGGTAGGGAGTAACCAGTATAAAGATTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCTCGCAGGCGAGTTTTAAAATTGTTGTACTCGTCAATAACTTACTGCTc                                                           >  1:919289/1‑150 (MQ=255)
                                                                                           cGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTcg                                                          >  1:1471668/1‑150 (MQ=255)
                                                                                           cGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTcg                                                          >  1:849033/1‑150 (MQ=255)
                                                                                           cGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATAACGTACTGCTcg                                                          >  1:1163563/1‑150 (MQ=255)
                                                                                                gggTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCggcggc                                                     >  2:1443041/1‑150 (MQ=255)
                                                                                                gggTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCggcggc                                                     >  1:1670085/1‑150 (MQ=255)
                                                                                                  gTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGGCAATCACGTACTGCTCGGCGgcgc                                                   <  2:302438/150‑1 (MQ=255)
                                                                                                     gggAGTAACCTGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCt                                                <  1:1807051/150‑1 (MQ=255)
                                                                                                     gggAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCt                                                <  2:1466887/150‑1 (MQ=255)
                                                                                                           aaCCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTgatcgg                                          <  1:482735/150‑1 (MQ=255)
                                                                                                           aaCCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTgatcgg                                          <  1:536132/150‑1 (MQ=255)
                                                                                                           aaCCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTgatcgg                                          >  2:1117009/1‑150 (MQ=255)
                                                                                                           aaCCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTgatcgg                                          >  2:757740/1‑150 (MQ=255)
                                                                                                                gTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTgatcggcgatc                                     <  2:1838132/150‑1 (MQ=255)
                                                                                                                 tataAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTgatcggcgatcg                                    <  1:784501/150‑1 (MQ=255)
                                                                                                                   taAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTgatcggcgatcggc                                  <  1:209995/150‑1 (MQ=255)
                                                                                                                   taAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTAGTAGTACTCGGCAAGCACGTACTGCTCGGCGACGCGCAcatcggcgctcggc                                  >  1:929573/1‑150 (MQ=255)
                                                                                                                    aaaGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGGCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCt                                 <  2:1127241/150‑1 (MQ=255)
                                                                                                                     aaGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTc                                <  1:1154865/150‑1 (MQ=255)
                                                                                                                           tCGAATGGATATAAAAACAGCCTCAGAAATTAGTCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCTGCGGCGCGCTGATCGGCGATCGGCTCGAcgcg                          <  2:1606579/150‑1 (MQ=255)
                                                                                                                           tCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGAcgcg                          <  1:1234921/150‑1 (MQ=255)
                                                                                                                                tGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACgg                     <  1:1360270/150‑1 (MQ=255)
                                                                                                                                  gATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGATAGGCCGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGgcg                   >  2:604502/1‑150 (MQ=255)
                                                                                                                                  gATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGgcg                   >  1:169768/1‑150 (MQ=255)
                                                                                                                                   aTAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGgcgc                  >  2:1464296/1‑150 (MQ=255)
                                                                                                                                     aaaaaaaCAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAg                >  2:1017567/1‑150 (MQ=255)
                                                                                                                                      tttctaCAGCCTCCTTTCGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGt               >  1:1675576/6‑150 (MQ=255)
                                                                                                                                      aaaaaaCAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGt               >  2:1839894/1‑150 (MQ=255)
                                                                                                                                      aaaaaaCAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCACGCGAGTTTTCAACTTGATGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGt               >  1:691695/1‑150 (MQ=255)
                                                                                                                                          taCAGCCTCCTTTCGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTAttt           >  2:959614/2‑150 (MQ=255)
                                                                                                                                           aCAGCCTCCTTTCGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCTCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGTCTGCGCGCTGATCGGCGATCGGCTCGACGCGAAAGGCGCAGTAtttt          >  2:679160/1‑149 (MQ=255)
                                                                                                                                            cAGCCTCCGGTCGGAGGCTGTAGAAAGGCCGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCCACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGCTCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTgg         >  1:70439/1‑149 (MQ=255)
                                                                                                                                                   cGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCTGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCgg  <  1:752911/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TAAATACTGGCTGCCAGGCTCGGAGTCGTTGGCGTGGGCAAAGAATGTCAAAAACAACATCAATGCGATAATTTTTTTCATTTTGATAATACGATTGGGTAGGGAGTAACCAGTATAAAGAGTTCGAATGGATAAAAAAACAGCCTCCGAAAGGAGGCTGTAGAAAGGACGGTATTACGCCAGTGCCGCTTCGCGCAGGCGAGTTTTCAACTTGTTGTACTCGTCAATCACGTACTGCTCGGCGGCGCGCTGATCGGCGATCGGCTCGACGCGAACGGCGCAGTATTTGTACTCCGG  >  NC_000913/4297045‑4297341

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