Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,703,251 A→T 9.7% intergenic (+17/+17) add → / ← ydgJ adenosine deaminase/putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,703,2510AT 9.7% 163.6 / 7.3 59intergenic (+17/+17)add/ydgJadenosine deaminase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (23/30);  new base T (2/4);  total (25/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.21e-01

TACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCATCACGAATAGCCG  >  NC_000913/1703105‑1703398
                                                                                                                                                  |                                                                                                                                                   
taCCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATTGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAAt                                                                                                                                                   <  8:188386/149‑1 (MQ=255)
taCCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCaaa                                                                                                                                                    <  1:321226/148‑1 (MQ=255)
        ccGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCATTTGCAccatc                                                                                                                                           <  3:363959/149‑1 (MQ=255)
        ccGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCATTTGCAccatc                                                                                                                                           <  2:347885/149‑1 (MQ=255)
             ccAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCAttt                                                                                                                                      >  6:194986/1‑149 (MQ=255)
                  tgctgGGTTATCCCGCGAGCAAATCCGCCAGGCAGAGAGTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGCGGAAAACCGCGCACTGCGCGATAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCAt                                                                                                                                 >  4:106833/1‑149 (MQ=255)
                      gggTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCaa                                                                                                                             <  5:194986/149‑1 (MQ=255)
                       ggTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGGGCAAATGCACCATCCATTTTTCATGCaa                                                                                                                             <  7:284206/148‑1 (MQ=255)
                       ggTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCGGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAg                                                                                                                            >  1:238941/1‑149 (MQ=255)
                              ccGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGcaca                                                                                                                       >  2:167767/1‑147 (MQ=255)
                                 cgAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTc                                                                                                                  >  2:31328/1‑149 (MQ=255)
                                    gCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTcgcg                                                                                                               >  6:146025/1‑149 (MQ=255)
                                    gCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTcgcg                                                                                                               >  3:313615/1‑149 (MQ=255)
                                     cAAAACCGCCAGGCACAGTTTAATGGGCTGGAAATGGCTTTCATCTGCGCTGAGGAAAAACGCGCACTTCTAGATAAAGTCGACGCGTTGTAACAAACTGGATGGTGCATATTCACCATCCATTTTTCATGCTAGGCACATATTcgcgc                                                                                                              <  8:353286/149‑1 (MQ=255)
                                            gCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGttt                                                                                                       >  1:150977/1‑149 (MQ=255)
                                               aGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGc                                                                                                    <  1:167767/149‑1 (MQ=255)
                                               aGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGc                                                                                                    <  2:150977/149‑1 (MQ=255)
                                                ggCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCg                                                                                                   >  3:16445/1‑149 (MQ=255)
                                                        ttAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGAt                                                                                           <  5:146025/149‑1 (MQ=255)
                                                        ttAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGAt                                                                                           <  6:146332/149‑1 (MQ=255)
                                                         tAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCa                                                                                                                                  <  4:123931/109‑1 (MQ=255)
                                                         tAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCa                                                                                                                                  >  3:123931/1‑109 (MQ=255)
                                                                 tGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCt                                                                                   >  2:345192/1‑148 (MQ=255)
                                                                   gAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCaa                                                                                                                                                     <  1:46832/80‑1 (MQ=255)
                                                                   gAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCaa                                                                                                                                                     >  2:46832/1‑80 (MQ=255)
                                                                   gAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCt                                                                                   <  5:85574/146‑1 (MQ=255)
                                                                    aaaTGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCtcaat                                                                               <  1:345192/149‑1 (MQ=255)
                                                                    aaaTGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCtcaat                                                                               >  2:185809/1‑149 (MQ=255)
                                                                     aaTGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCtcaatc                                                                              <  6:135423/149‑1 (MQ=255)
                                                                         gCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAAc                                                                           <  1:31328/148‑1 (MQ=255)
                                                                           tttCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTc                                                                         >  2:71323/1‑148 (MQ=255)
                                                                           tttCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATATTTCATGCAAGGCACAAGGTCGCGCGAGGTTTGCCGGATTCGATCACCCGCTCAAGCACCTc                                                                         >  6:265777/1‑148 (MQ=255)
                                                                             tCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCa                                                                       <  7:270756/148‑1 (MQ=255)
                                                                                tCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGc                                                                                                                               >  6:14322/1‑89 (MQ=255)
                                                                                tCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGc                                                                                                                               <  5:14322/89‑1 (MQ=255)
                                                                                  aGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTAcc                                                                  <  2:171566/148‑1 (MQ=255)
                                                                                  aGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCt                                                                 <  4:16445/149‑1 (MQ=255)
                                                                                      cTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGGGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGAt                                                             <  7:389048/149‑1 (MQ=255)
                                                                                          ggAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCAt                                                                                                                                 <  1:194366/77‑1 (MQ=255)
                                                                                          ggAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCAt                                                                                                                                 >  2:194366/1‑77 (MQ=255)
                                                                                           gAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCt                                                        >  6:222600/1‑149 (MQ=255)
                                                                                             aaaaCGCGCACTGCGATAAAAAGTCGCAGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAATTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTgg                                                      <  7:104471/149‑1 (MQ=255)
                                                                                                 cgcgCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCTCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTg                                                  >  7:391204/1‑149 (MQ=255)
                                                                                                       cTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGaa                                            <  8:54847/149‑1 (MQ=255)
                                                                                                              aaaaaGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGAtt                                     <  2:91555/149‑1 (MQ=255)
                                                                                                                aaaGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGcc                                                <  1:9395/136‑1 (MQ=255)
                                                                                                                aaaGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGcc                                                >  2:9395/1‑136 (MQ=255)
                                                                                                                  aGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCa                                 <  6:75268/149‑1 (MQ=255)
                                                                                                                       ccgcGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATTCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATc                            <  5:265777/149‑1 (MQ=255)
                                                                                                                       ccgcGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATc                            <  1:71323/149‑1 (MQ=255)
                                                                                                                        cgcgAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCAGTACCTGGAGTGCCTGGCTTGCCGGCACCGGATTTTCACCATCg                           >  6:373328/1‑149 (MQ=255)
                                                                                                                           gAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCa                        >  2:327603/1‑149 (MQ=255)
                                                                                                                            aaGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCAt                       <  5:222600/149‑1 (MQ=255)
                                                                                                                                aaCAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCAttt                     >  4:239056/1‑147 (MQ=255)
                                                                                                                                aaCAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCAttt                     <  3:239056/147‑1 (MQ=255)
                                                                                                                                       tggatggTGCATTTGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCATCa            >  7:110087/1‑149 (MQ=255)
                                                                                                                                         gatggTGCATTTGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAAc                  >  2:204315/1‑141 (MQ=255)
                                                                                                                                         gatggTGCATTTGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAAc                  <  1:204315/141‑1 (MQ=255)
                                                                                                                                             gTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGGCGCGCTATGTTTGGCGGATTCGAGTCCCAGCTCAATCAACTCCATTACCTGGATGGCCTGGCTTGCCGGAACCGGATTTTCCCCATCGCCATTTAACGCATCACGAAt       <  2:238941/148‑1 (MQ=255)
                                                                                                                                              tGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCATCACGAATAg     >  7:99068/1‑149 (MQ=255)
                                                                                                                                                cAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCGACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCATCACGAATAGcc   <  4:313615/149‑1 (MQ=255)
                                                                                                                                                cAAATGAACCATCCATTTTTCATGCAAGGCACAAAGGCGCGCTATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATATTCACCATCGCCATTTAACGCATCACGAATAGc    <  5:306325/148‑1 (MQ=255)
                                                                                                                                                 aaaTGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCATCACGAATAGCCg  >  8:154425/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TACCGTTGCCGCGCCAGCTGCTGGGTTATCCCGCGAGCAAATCCGCCAGGCACAGATTAATGGTCTGGAAATGGCTTTCCTCAGCGCTGAGGAAAAACGCGCACTGCGAGAAAAAGTCGCCGCGAAGTAACAAAATGGATGGTGCAAATGCACCATCCATTTTTCATGCAAGGCACAAAGTCGCGCGATGTTTGGCGGATTCGATGCCCAGCTCAATCAACTCCATTACCTGGATTGCCTGGCTTGCCGGAACCGGATTTTCACCATCGCCATTTAACGCATCACGAATAGCCG  >  NC_000913/1703105‑1703398

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