Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 644,038 G→A 5.2% intergenic (‑71/+159) rnk ← / ← rna regulator of nucleoside diphosphate kinase/ribonuclease I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913644,0380GA 5.2% 715.4 / 17.3 211intergenic (‑71/+159)rnk/rnaregulator of nucleoside diphosphate kinase/ribonuclease I
Reads supporting (aligned to +/- strand):  ref base G (94/106);  new base A (6/5);  total (100/111)
Fisher's exact test for biased strand distribution p-value = 7.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.22e-01

Only 100 of 217 total aligned reads displayed.

AGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGCGTGTTTTGTA  >  NC_000913/643890‑644186
                                                                                                                                                    |                                                                                                                                                    
aGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAgg                                                                                                                                                      <  4:692683/149‑1 (MQ=255)
aGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAgg                                                                                                                                                      <  2:82662/149‑1 (MQ=255)
      ggCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGa                                                                                                                                                <  4:847928/149‑1 (MQ=255)
      ggCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGa                                                                                                                                                <  2:1228342/149‑1 (MQ=255)
        cGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAgg                                                                                                                                              <  3:735108/149‑1 (MQ=255)
        cGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAgg                                                                                                                                              <  3:487320/149‑1 (MQ=255)
          ggTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGt                                                                                                                                            <  2:646058/149‑1 (MQ=255)
            ttGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGttt                                                                                                                                          >  2:1259945/1‑149 (MQ=255)
              gCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTaa                                                                                                                                        <  2:748844/149‑1 (MQ=255)
                tCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACt                                                                                                                                      >  2:1180703/1‑149 (MQ=255)
                tCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACt                                                                                                                                      >  2:173593/1‑149 (MQ=255)
                tCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACt                                                                                                                                      >  2:514918/1‑149 (MQ=255)
                  cagcagAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTca                                                                                                                                    <  2:107411/149‑1 (MQ=255)
                    gcagAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCAGCGAGGGTTTAACTcaca                                                                                                                                  <  3:466989/149‑1 (MQ=255)
                        aaTATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCGCCGCCTGACGGTGAGGGt                                                                                                                                            >  2:806091/1‑135 (MQ=255)
                        aaTATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCGCCGCCTGACGGTGAGGGt                                                                                                                                            <  1:806091/135‑1 (MQ=255)
                        aaTATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACatgat                                                                                                                              <  3:142194/149‑1 (MQ=255)
                        aaTATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACatgat                                                                                                                              <  3:922990/149‑1 (MQ=255)
                               aTGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTATCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATactg                                                                                                                       >  1:1000417/1‑149 (MQ=255)
                                tGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATactga                                                                                                                      <  1:1491719/149‑1 (MQ=255)
                                   gTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATactgactg                                                                                                                   >  1:597898/1‑149 (MQ=255)
                                     tCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTAGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGtt                                                                                                                 <  3:1257143/149‑1 (MQ=255)
                                       ggCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGc                                                                                                               <  1:677158/149‑1 (MQ=255)
                                              aGGTCGTTAATGATGCTAGTTGGTGTGGACATTTTGACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGGTAAACCTACCCAAAAGTAAATCCTCACCGTCAGGCGGCGAGGGCTTAACTGACGAGATGATACTGACTGTTGCTCActct                                                                                                        >  4:918087/1‑149 (MQ=255)
                                              aGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACatgatg                                                                                                                             >  1:1219129/1‑128 (MQ=255)
                                              aGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCActct                                                                                                        >  3:390663/1‑149 (MQ=255)
                                                gTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCttt                                                                                                      <  4:1062853/149‑1 (MQ=255)
                                                 tCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTg                                                                                                     <  4:177376/149‑1 (MQ=255)
                                                  cGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGa                                                                                                    <  1:231471/149‑1 (MQ=255)
                                                   gTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGaa                                                                                                   <  3:176515/149‑1 (MQ=255)
                                                   gTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGaa                                                                                                   <  1:1457517/149‑1 (MQ=255)
                                                   gTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGaa                                                                                                   >  3:1155467/1‑149 (MQ=255)
                                                    ttAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAg                                                                                                  >  4:1335537/1‑149 (MQ=255)
                                                    ttAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAg                                                                                                  >  2:395142/1‑149 (MQ=255)
                                                    ttAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAg                                                                                                  <  4:390663/149‑1 (MQ=255)
                                                    ttAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAg                                                                                                  >  3:423693/1‑149 (MQ=255)
                                                    ttAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAg                                                                                                  <  3:878781/149‑1 (MQ=255)
                                                    ttAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAg                                                                                                  >  3:684319/1‑149 (MQ=255)
                                                      aaTGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTg                                                                                                >  2:919810/1‑149 (MQ=255)
                                                         gatgatAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGAtt                                                                                             >  2:1365505/1‑149 (MQ=255)
                                                            gatAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGttt                                                                                                                                          <  3:1129208/101‑1 (MQ=255)
                                                            gatAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGttt                                                                                                                                          >  4:1129208/1‑101 (MQ=255)
                                                             atAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCg                                                                                         <  2:249657/149‑1 (MQ=255)
                                                               aGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTc                                                                                       >  2:770648/1‑149 (MQ=255)
                                                                 ttGGTCTGGACCTTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTcac                                                                                     >  1:664566/1‑149 (MQ=255)
                                                                 ttGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCATGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTcac                                                                                     <  4:1288197/149‑1 (MQ=255)
                                                                 ttGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTcac                                                                                     <  3:1335537/149‑1 (MQ=255)
                                                                     tCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCGCCGCCTGACGGTGAggg                                                                                                                                             >  1:786921/1‑89 (MQ=37)
                                                                     tCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCGCCGCCTGACGGTGAggg                                                                                                                                             <  2:786921/89‑1 (MQ=38)
                                                                     tCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTc                                                                                 >  4:1365003/1‑149 (MQ=255)
                                                                       tGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAg                                                                               <  4:237927/149‑1 (MQ=255)
                                                                        ggACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAgg                                                                              <  1:491940/149‑1 (MQ=255)
                                                                         gACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCttt                                                                                                      >  4:627931/1‑124 (MQ=255)
                                                                         gACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAg                                                                               >  2:1129963/1‑147 (MQ=255)
                                                                            aTTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGTCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAAt                                                                          >  2:36739/1‑149 (MQ=255)
                                                                              ttttACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTc                                                                        <  4:576504/149‑1 (MQ=255)
                                                                               tttACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTcc                                                                       <  2:48768/149‑1 (MQ=255)
                                                                                 tACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTc                                                                     <  3:1481843/149‑1 (MQ=255)
                                                                                   cTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGa                                                                            <  4:1000562/140‑1 (MQ=255)
                                                                                    tCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATa                                                                                                                          <  2:1155562/93‑1 (MQ=255)
                                                                                    tCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATa                                                                                                                          >  1:1155562/1‑93 (MQ=255)
                                                                                    tCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAAt                                                                  >  1:1377894/1‑149 (MQ=255)
                                                                                     ccATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATg                                                                 >  3:246433/1‑149 (MQ=255)
                                                                                     ccATGTCGTCGGTGCTGCGAGTGTCACAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATg                                                                 >  4:324391/1‑149 (MQ=255)
                                                                                      cATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGc                                                                <  1:395142/149‑1 (MQ=255)
                                                                                      cATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGc                                                                <  3:259997/149‑1 (MQ=255)
                                                                                       aTGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCa                                                               <  2:1393774/149‑1 (MQ=255)
                                                                                       aTGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCa                                                               <  1:1181642/149‑1 (MQ=255)
                                                                                         gtcgtcGGTGCTGCGAGAGTCGCAGGTAAACATACCCAAAAGAAAACCCTGACCGCCAGGCGTCGAGGGTTTAACTCACATCATGATCCTGACTGTTGCTCACCCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCGTCAATGCAgc                                                             <  2:1341883/149‑1 (MQ=255)
                                                                                                gTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTa                                                       >  1:337083/1‑148 (MQ=255)
                                                                                                  gctgcGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGt                                                    <  1:1424318/149‑1 (MQ=255)
                                                                                                    tgcGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGtat                                                  >  1:451748/1‑149 (MQ=255)
                                                                                                            tCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTa                                          >  1:1060200/1‑149 (MQ=255)
                                                                                                             cGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTaa                                         >  2:644024/1‑149 (MQ=255)
                                                                                                                      aCATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTgg                                >  1:1129831/1‑149 (MQ=255)
                                                                                                                      aCATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTgg                                >  4:159545/1‑149 (MQ=255)
                                                                                                                       cATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGt                               <  2:432043/149‑1 (MQ=255)
                                                                                                                           cccAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTAt                           <  4:508623/149‑1 (MQ=255)
                                                                                                                            ccAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATg                          >  2:132023/1‑149 (MQ=255)
                                                                                                                            ccAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATg                          >  3:333963/1‑149 (MQ=255)
                                                                                                                            ccAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATg                          >  2:593182/1‑149 (MQ=255)
                                                                                                                             cAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCActc                                                                                                         <  3:1104101/69‑1 (MQ=255)
                                                                                                                             cAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCActc                                                                                                         >  4:1104101/1‑69 (MQ=255)
                                                                                                                             cAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATgg                         >  1:1328132/1‑149 (MQ=255)
                                                                                                                             cAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATgg                         >  4:631809/1‑149 (MQ=255)
                                                                                                                             cAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATgg                         >  1:829998/1‑149 (MQ=255)
                                                                                                                              aaaaGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGa                        >  1:633728/1‑149 (MQ=255)
                                                                                                                              aaaaGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGa                        <  4:142332/149‑1 (MQ=255)
                                                                                                                                 aGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGa                     <  3:751290/149‑1 (MQ=255)
                                                                                                                                  ttaaaCCCTCGCCGCCTGACGGTGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCttt                                                                                                      <  1:983195/65‑1 (MQ=18)
                                                                                                                                  ttaaaCCCTCGCCGCCTGACGGTGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCttt                                                                                                      >  2:983195/3‑67 (MQ=18)
                                                                                                                                  gAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATg                          >  2:30453/1‑143 (MQ=255)
                                                                                                                                  gAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGAc                    <  3:613777/149‑1 (MQ=255)
                                                                                                                                     aaCCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGaa                 >  1:732115/1‑149 (MQ=255)
                                                                                                                                      aCCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAAt                <  4:423693/149‑1 (MQ=255)
                                                                                                                                      aCCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAAt                >  3:1089141/1‑149 (MQ=255)
                                                                                                                                       cccTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCttt                                                                                                      >  4:449523/1‑62 (MQ=255)
                                                                                                                                         cTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCttt                                                                                                      >  4:1095971/1‑60 (MQ=255)
                                                                                                                                            aCCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGCgt          <  2:1000417/149‑1 (MQ=255)
                                                                                                                                                    gcggcgAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGCGTGTTTTGTa  >  3:96212/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGCGTGTTTTGTA  >  NC_000913/643890‑644186

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