Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,032,336 G→C 6.8% intergenic (‑19/+48) dcm ← / ← yedJ DNA cytosine methyltransferase/putative HD superfamily phosphohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,032,3360GC 6.8% 659.4 / 28.4 190intergenic (‑19/+48)dcm/yedJDNA cytosine methyltransferase/putative HD superfamily phosphohydrolase
Reads supporting (aligned to +/- strand):  ref base G (76/101);  new base C (8/5);  total (84/106)
Fisher's exact test for biased strand distribution p-value = 2.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

Only 100 of 194 total aligned reads displayed.

GGATGATGTGCCGGTGGTTTGGGTAATAACGTTTGCAGATGGGCGAACTCTTTCTCACTTAAACGGTGCCATGCCGAGTCATTCGCCAGCGCACGTTTTAAAATTGCCGCGCTCCAGTGATTCTCACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAA  >  NC_000913/2032070‑2032608
                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                
       gtgCCGGTGGTTTGGGTAATAACGTTTGCAGATGGGCGAACTCTTTCTCACTTAAACGGTGCCATGCCGAGTCATTCGCCAGCGCACGTTTTAAAATTGCCGCGCTCCAGTGATTCTCACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCCGTTTCTCCAGCATTTCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCAGAGATTTACCCGGCCATCTGAGATTGCCGACAgg                                                                                                                                                                                                                                                                    >  2:243257/1‑274 (MQ=255)
                                 tGCAGATGGGCGAACTCTTTCTCACTTAAACGGTGCCATGCCGAGTCATTCGCCAGCGCACGTTTTAAAATTGCCGCGCTCCAGTGATTCTCACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGCTTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGTCATTTAGGGc                                                                                                                                                                                                                                          >  2:274/1‑274 (MQ=255)
                                                                  tGCCATGCCGAGTCATTCGCCAGCGCACGTTTTAAAATTGCCGCGCTCCAGTGATTCTCACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATgg                                                                                                                                                                                                                                                                            >  2:365925/1‑207 (MQ=255)
                                                                      atgccggtggaTTCGCCAGTGCACGGTTGGAAATGGCGACGTTCTAGTGATTCTCACCTAAACCATTTAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCATCTGTGTCTCCAGCATTGCCTGTGCGGCATTACGGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGGGTTGGCCGACAGGTTAACATAATTACAGCAATTTAGGGCGACGCGCTCAGCCAGCAGATGAAAACGCATCTATCAcc                                                                                                                                                                                                    <  4:74059/265‑1 (MQ=255)
                                                                            aGTCATTCGCCAGCGCACGTTTTAAAATTGCCGCGCTCCAGTGATTCTCACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTAAAGCTATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTtgtg                                                                                                                                                                                               >  2:196478/1‑274 (MQ=255)
                                                                                                              gcTCCAGAGATTCTCACCTAAACCATTAAGCAGGGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGGGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCa                                                                                                                                                             <  2:139739/274‑1 (MQ=255)
                                                                                                                      gATTCTCACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTg                                                                                                                                                                                                                                                                                  <  5:814593/149‑1 (MQ=255)
                                                                                                                         tctcACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTgaga                                                                                                                                                                                                                                                                               <  8:542424/149‑1 (MQ=255)
                                                                                                                           tcACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATg                                                                                                                                                                                                                                                                             >  8:1367356/1‑149 (MQ=255)
                                                                                                                           tcACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATg                                                                                                                                                                                                                                                                             <  5:285221/149‑1 (MQ=255)
                                                                                                                               cTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCg                                                                                                                                                                                                                                                                         >  5:241656/1‑149 (MQ=255)
                                                                                                                                tACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCgg                                                                                                                                                                                                                                                                        <  6:725680/149‑2 (MQ=255)
                                                                                                                                  cacCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATgg                                                                                                                                                                                                                                                                            >  4:347253/1‑143 (MQ=255)
                                                                                                                                     cATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGt                                                                                                                                                                                                                                                                   <  8:108467/149‑1 (MQ=255)
                                                                                                                                     cATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGt                                                                                                                                                                                                                                                                   <  6:554882/149‑1 (MQ=255)
                                                                                                                                     cATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGt                                                                                                                                                                                                                                                                   <  5:1184258/149‑1 (MQ=255)
                                                                                                                                     cATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGt                                                                                                                                                                                                                                                                   >  7:57927/1‑149 (MQ=255)
                                                                                                                                      aTTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGtt                                                                                                                                                                                                                                                                  >  5:188225/1‑149 (MQ=255)
                                                                                                                                       ttAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTa                                                                                                                                                                                                                                                                 >  7:1066648/1‑149 (MQ=255)
                                                                                                                                         aaGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAAc                                                                                                                                                                                                                                                               <  5:1057597/149‑1 (MQ=255)
                                                                                                                                           gCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACAt                                                                                                                                                                                                                                                             <  7:650156/149‑1 (MQ=255)
                                                                                                                                            cTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTACCTGTGCGGCATTCCCGATGCTGTATGAATCGGTTACTGTTATATTTTCCTCCATAGATTTCACCGTCCATCGCAGATGGCCGACAGGTTAACATa                                                                                                                                                                                                                                                            >  8:156103/1‑149 (MQ=255)
                                                                                                                                             tGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATaa                                                                                                                                                                                                                                                           <  8:1035904/149‑1 (MQ=255)
                                                                                                                                              gcgcCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAAt                                                                                                                                                                                                                                                          <  8:878428/149‑1 (MQ=255)
                                                                                                                                               cgcCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATgg                                                                                                                                                                                                                                                                            >  4:292467/1‑130 (MQ=255)
                                                                                                                                               cgcCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACGCCTTCATTCTCCCCCGCCAGTTCGGCCCTGAGCTTCGCCATAAACTCGACTAAAAAATgc                                                                                                                            >  2:125920/1‑274 (MQ=255)
                                                                                                                                                  caccaACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACa                                                                                                                                                                                                                                                      <  6:228464/149‑1 (MQ=255)
                                                                                                                                                  caccaACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATGGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACa                                                                                                                                                                                                                                                      <  7:720228/149‑1 (MQ=255)
                                                                                                                                                   accaACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAg                                                                                                                                                                                                                                                     <  4:208330/149‑1 (MQ=255)
                                                                                                                                                      aaCGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGt                                                                                                                                                                                                                                                                   <  5:461789/132‑1 (MQ=255)
                                                                                                                                                      aaCGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCa                                                                                                                                                                                                                                                  >  7:507540/1‑149 (MQ=255)
                                                                                                                                                           tttAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAg                                                                                                                                                                                                                                             <  5:704680/149‑1 (MQ=255)
                                                                                                                                                            ttAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAgg                                                                                                                                                                                                                                            <  8:952392/149‑1 (MQ=255)
                                                                                                                                                            ttAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAgg                                                                                                                                                                                                                                            <  8:629086/149‑1 (MQ=255)
                                                                                                                                                            ttAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAgg                                                                                                                                                                                                                                            <  6:1010696/149‑1 (MQ=255)
                                                                                                                                                            ttAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAgg                                                                                                                                                                                                                                            <  8:326841/149‑1 (MQ=255)
                                                                                                                                                              aaCATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGc                                                                                                                                                                                                                                          >  7:1285043/1‑149 (MQ=255)
                                                                                                                                                                catcatAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCa                                                                                                                                                                                                                                        <  6:776253/149‑1 (MQ=255)
                                                                                                                                                                 atcatAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCAc                                                                                                                                                                                                                                       >  7:352809/1‑149 (MQ=255)
                                                                                                                                                                   catAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCCTTTTGCTGCGCCAgctgc                                                                                                       >  2:219405/1‑275 (MQ=255)
                                                                                                                                                                    ctAAATTTGCAGCGGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTGCCTGCATAGATTTCCCCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCAcgcg                                                                                                                                                                                                                                    >  8:1264325/2‑149 (MQ=255)
                                                                                                                                                                    atAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCAcgcg                                                                                                                                                                                                                                    >  6:1302044/1‑149 (MQ=255)
                                                                                                                                                                         tttGCAGCAGTGTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAAc                                                                                                                                                                                                                                                               >  5:772277/1‑117 (MQ=255)
                                                                                                                                                                             cagcagTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAg                                                                                                                                                                                                                                             >  5:84788/1‑131 (MQ=255)
                                                                                                                                                                             cagcagTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAg                                                                                                                                                                                                                                             <  6:84788/131‑1 (MQ=255)
                                                                                                                                                                              gtgagtctatcCAGCAGTGCCTGTGGGGGAGTCAGCGGCGTGTATGCTACGGTAAGTGATACAGTTTCCGAGATAGATGTCACGGGCGTTCTGGGATGGCCGACAGGTTAACATAATACCAGACATGTAGGGCCACGCGCTCAGCACGCGGAAGAAAACGCAGCTATCACCTGGTGATCGAAACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGc                                                                                            <  4:243257/266‑1 (MQ=255)
                                                                                                                                                                                 agTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGa                                                                                                                                                                                                                       <  7:1341087/149‑1 (MQ=255)
                                                                                                                                                                                     tctcCAGCATTGCATGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCTGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGc                                                                                     <  1:6768/275‑1 (MQ=255)
                                                                                                                                                                                       tcCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGaaaa                                                                                                                                                                                                                 >  6:306062/1‑149 (MQ=255)
                                                                                                                                                                                       tcCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCggt                                                                                                                                                                                                                                                                       >  8:228874/1‑93 (MQ=255)
                                                                                                                                                                                       tcCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTCAGATGGCCggt                                                                                                                                                                                                                                                                       <  7:228874/95‑3 (MQ=255)
                                                                                                                                                                                       tcCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCGGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGaaaa                                                                                                                                                                                                                 <  8:276374/149‑1 (MQ=255)
                                                                                                                                                                                        ccAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAAc                                                                                                                                                                                                                <  6:543009/149‑1 (MQ=255)
                                                                                                                                                                                         cAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACg                                                                                                                                                                                                               >  6:893385/1‑149 (MQ=255)
                                                                                                                                                                                                 ccTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATc                                                                                                                                                                                                       <  7:294141/149‑1 (MQ=255)
                                                                                                                                                                                                 ccTGTGCGGCATTCCCGGTGCTGTATGAATCGATTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTCTCACCTTGTGCTCGGCACCTTCATTCTCCACCGCCAGTTCGGCCCTGAGCTTCGCCATAAACTCGACTAAAAAAGGCGCTTTGTCCTGCGCCCCACGCTTCCCCCTTGCGGCTTCACTGGTTTGc                                                                            >  4:98495/1‑272 (MQ=255)
                                                                                                                                                                                                      gCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCtt                                                                                                                                                                                                  <  5:622690/149‑1 (MQ=255)
                                                                                                                                                                                                       cGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTtg                                                                                                                                                                                                 <  7:1431275/149‑1 (MQ=255)
                                                                                                                                                                                                          cATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGa                                                                                                                                                                                              <  7:838044/149‑1 (MQ=255)
                                                                                                                                                                                                          cATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCAGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGa                                                                                                                                                                                              >  7:449578/1‑149 (MQ=255)
                                                                                                                                                                                                            ttCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATc                                                                                                                                                                                            <  5:219523/149‑1 (MQ=255)
                                                                                                                                                                                                                ccGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGac                                                                                                                                                                                         >  8:1239136/1‑149 (MQ=255)
                                                                                                                                                                                                                  ggTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGacac                                                                                                                                                                                       <  7:1385537/149‑1 (MQ=255)
                                                                                                                                                                                                                  ggTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGacac                                                                                                                                                                                       <  8:507540/149‑1 (MQ=255)
                                                                                                                                                                                                                    tGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGAATTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCt                                                                                                                                                                                     <  8:193655/149‑1 (MQ=255)
                                                                                                                                                                                                                       tGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCa                                                                                                                                                                                  <  6:247084/149‑1 (MQ=255)
                                                                                                                                                                                                                       tGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCa                                                                                                                                                                                  <  8:1427153/149‑1 (MQ=255)
                                                                                                                                                                                                                        gTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAATGCATCTATCACCTTGTGATCGACACCTTCAt                                                                                                                                                                                 <  6:1077098/149‑1 (MQ=255)
                                                                                                                                                                                                                        gTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCa                                                                                                                                                                                  <  6:767788/148‑1 (MQ=255)
                                                                                                                                                                                                                            gAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATtctc                                                                                                                                                                             <  8:783265/149‑1 (MQ=255)
                                                                                                                                                                                                                              aTCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTGAGTCTGGAAATGa                                              >  2:304495/1‑274 (MQ=255)
                                                                                                                                                                                                                               tCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTcagccagc                                                                                                                                                                                                                          <  5:80389/101‑1 (MQ=255)
                                                                                                                                                                                                                               tCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTcccc                                                                                                                                                                          >  7:1288891/1‑149 (MQ=255)
                                                                                                                                                                                                                                cGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTccccc                                                                                                                                                                         <  7:121859/149‑1 (MQ=255)
                                                                                                                                                                                                                                   ttACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGa                                                                                                                                                                                              >  4:136957/1‑125 (MQ=255)
                                                                                                                                                                                                                                   ttACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAACAAATGCGCGTTGTTCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGACAATAATTCAgg                                        >  2:69388/1‑274 (MQ=255)
                                                                                                                                                                                                                                      cTGATATATTTTCCTGCATAGATTTCACCCGCCATCTGAGATGGCCGACAGGTTAACATAATTACa                                                                                                                                                                                                                                                      <  8:1232810/66‑1 (MQ=255)
                                                                                                                                                                                                                                       tGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGtt                                                                                                                                                                  >  6:197370/1‑149 (MQ=255)
                                                                                                                                                                                                                                         atatatTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTaa                                                                                                                                                                                                                                                                <  4:44896/53‑1 (MQ=255)
                                                                                                                                                                                                                                           atatTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGc                                                                                                                                                              >  5:109730/1‑149 (MQ=255)
                                                                                                                                                                                                                                             atTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGTCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCAc                                                                                                                                                            >  6:387540/1‑149 (MQ=255)
                                                                                                                                                                                                                                             atTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCgg                                                                                                                                                               <  5:1334491/146‑1 (MQ=255)
                                                                                                                                                                                                                                             atTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCgg                                                                                                                                                               >  6:1334491/1‑146 (MQ=255)
                                                                                                                                                                                                                                              ttttCCGGCATAGAGTTCCCCTGGCACTTGAGAGGGCAGACAGGTTAACAAAATTACCGCCTTTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGAATCTATCACCTTGTGAGCGAAACCTCCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCAAAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGttt                             <  4:274/275‑1 (MQ=255)
                                                                                                                                                                                                                                                 tCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTcag                                                                                                                                                                                                                               >  1:63704/1‑78 (MQ=255)
                                                                                                                                                                                                                                                   cTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCt                                                                                                                                                      <  5:1414518/149‑1 (MQ=255)
                                                                                                                                                                                                                                                   cTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGtt                                                                                        >  2:335121/1‑211 (MQ=255)
                                                                                                                                                                                                                                                     gCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTcagccag                                                                                                                                                                                                                           <  3:63704/78‑1 (MQ=255)
                                                                                                                                                                                                                                                     gCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTc                                                                                                                                                    <  7:958392/149‑1 (MQ=255)
                                                                                                                                                                                                                                                       aTAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGc                                                                                                                                                                       <  8:117465/128‑1 (MQ=255)
                                                                                                                                                                                                                                                       aTAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGc                                                                                                                                                                       >  7:117465/1‑128 (MQ=255)
                                                                                                                                                                                                                                                       aTAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGc                                                                                                                                                  >  5:483843/1‑149 (MQ=255)
                                                                                                                                                                                                                                                       aTAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCa                                                                                    <  4:335121/211‑1 (MQ=255)
                                                                                                                                                                                                                                                        tAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGcc                                                                                                                                                 <  7:1238201/149‑1 (MQ=255)
                                                                                                                                                                                                                                                           aTTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATa                                                                                                                                              >  7:109243/1‑149 (MQ=255)
                                                                                                                                                                                                                                                            tttCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATaa                                                                                                                                             <  8:1288891/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                   ggCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGAc                                                                                                                                      <  5:1302044/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                   ggCCATCTCAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGTCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGAc                                                                                                                                      >  7:997826/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                   ggCCATCTCAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCa                                                                                                                                                                                                             <  7:1431517/78‑1 (MQ=255)
                                                                                                                                                                                                                                                                        tCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTaaaa                                                                                                                                 >  8:713498/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                
GGATGATGTGCCGGTGGTTTGGGTAATAACGTTTGCAGATGGGCGAACTCTTTCTCACTTAAACGGTGCCATGCCGAGTCATTCGCCAGCGCACGTTTTAAAATTGCCGCGCTCCAGTGATTCTCACCTACACCATTAAGCTGCGCCACCAACGTTTTAACATCATAAATTTGCAGCAGTTTCTCCAGCATTGCCTGTGCGGCATTCCCGGTGCTGTATGAATCGGTTACTGATATATTTTCCTGCATAGATTTCACCGGCCATCTGAGATGGCCGACAGGTTAACATAATTACAGCCATTTAGGGCCACGCGCTCAGCCAGCGGATGAAAACGCATCTATCACCTTGTGATCGACACCTTCATTCTCCCCCGCCAGTTCGGCACTGAGCTTCGCCATAAACTCGACTAAAAAATGCGCGTTGTGCTGCGCCAGCTGCTTGCCTCTTGCGGTTTGCATGGTTTGCGGTAGCTTCAGCAACTTAGTCTGGAAATGATCCAGCGCGTAGCGTTTATCATCAAGCGGGCGATGCTGTGCGAA  >  NC_000913/2032070‑2032608

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