Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 644,034 T→C 10.3% intergenic (‑67/+163) 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,0340TC10.3% 380.3 / 27.2 126intergenic (‑67/+163)rnk/rnaregulator of nucleoside diphosphate kinase/ribonuclease I
Reads supporting (aligned to +/- strand):  ref base T (56/57);  new base C (7/6);  total (63/63)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.13e-01

Only 100 of 129 total aligned reads displayed.

AACCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGCGTGTT  >  NC_000913/643886‑644181
                                                                                                                                                    |                                                                                                                                                   
aaCCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGAGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATCCCCAAAAGAAAACCCTCACCGt                                                                                                                                                     >  4:1210075/1‑149 (MQ=255)
  ccAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCa                                                                                                                                                   <  3:1210075/149‑1 (MQ=255)
            cGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAgg                                                                                                                                         <  3:803492/149‑1 (MQ=255)
            cGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAgg                                                                                                                                         <  2:893340/149‑1 (MQ=255)
              ggTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATCATAGTTGGTCTTGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCTGTTAAACATACCCAAAAGAAAACCCTCGCCGCCTGACGGTGAGGGt                                                                                                                                       <  4:901974/149‑1 (MQ=255)
               gTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGa                                                                                                                                           >  1:1073779/1‑144 (MQ=255)
               gTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGa                                                                                                                                           <  2:1073779/144‑1 (MQ=255)
                ttGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAg                                                                                                                                          <  2:819855/144‑1 (MQ=255)
                    tCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCGCCGCCTGACGGTGAgg                                                                                                                                         <  1:483691/141‑1 (MQ=255)
                    tCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCGCCGCCTGACGGTGAgg                                                                                                                                         >  2:483691/1‑141 (MQ=255)
                    tCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACt                                                                                                                                 >  1:182269/1‑149 (MQ=255)
                           gAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTaa                                                                                                                                   <  4:1379218/140‑1 (MQ=255)
                           gAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTaa                                                                                                                                   >  3:1379218/1‑140 (MQ=255)
                              tatCGATGCGTTCGGCATCCAGGTCGTTAACGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACAtaatga                                                                                                                       >  4:1451730/1‑149 (MQ=255)
                                       gTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTaa                                                                                                                                   >  2:55054/1‑128 (MQ=255)
                                       gTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATacagactg                                                                                                              >  2:1131371/1‑149 (MQ=255)
                                            gCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACat                                                                                                                            >  2:843373/1‑130 (MQ=255)
                                            gCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACat                                                                                                                            <  1:843373/130‑1 (MQ=255)
                                                     tCGTTAATGATGATAGTTTGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCATGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTg                                                                                                <  2:473196/149‑1 (MQ=255)
                                                       gTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGaa                                                                                              >  4:720705/1‑149 (MQ=255)
                                                       gTTAATCATGATAGTGGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGaa                                                                                              <  1:1131371/149‑1 (MQ=255)
                                                        ttAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAg                                                                                             >  4:1441608/1‑149 (MQ=255)
                                                                gatAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGaa                                                                                              <  1:694770/140‑1 (MQ=255)
                                                                gatAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGaa                                                                                              >  2:694770/1‑140 (MQ=255)
                                                                gatAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGc                                                                                     >  1:847127/1‑149 (MQ=255)
                                                                gatAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGc                                                                                     >  3:504020/1‑149 (MQ=255)
                                                                    gTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTca                                                                                 >  4:1460863/1‑149 (MQ=255)
                                                                     ttGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGGTAATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTcac                                                                                >  1:808960/1‑149 (MQ=255)
                                                                     ttGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTcac                                                                                <  3:1460863/149‑1 (MQ=255)
                                                                      tGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACACAGAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGATCACTCTTTGAAGTGATTTGCGTcaca                                                                               <  3:1451730/149‑1 (MQ=255)
                                                                          cTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCa                                                                           <  4:504020/149‑1 (MQ=255)
                                                                           tGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAg                                                                          <  1:824454/149‑1 (MQ=255)
                                                                            ggACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAgg                                                                         <  2:808960/149‑1 (MQ=255)
                                                                             gACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAAc                                                                                                                                  <  1:642681/91‑1 (MQ=255)
                                                                             gACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTc                                                                                                                                <  4:1368655/93‑1 (MQ=255)
                                                                             gACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTc                                                                                                                                >  3:1368655/1‑93 (MQ=255)
                                                                             gACATTTTTACTCCATGTCGGCGGTGCTGCGAGTGTCGCAGATAAACATACGCAAAAGAACACCCTCACCGTCAGGCGGCGAGGGTTTAAc                                                                                                                                  >  2:642681/1‑91 (MQ=255)
                                                                                aTTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAAt                                                                     >  1:506136/1‑149 (MQ=255)
                                                                                  ttttACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTc                                                                   <  2:506136/149‑1 (MQ=255)
                                                                                       cTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGAttt                                                                                       <  2:772489/124‑1 (MQ=255)
                                                                                       cTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGAttt                                                                                       >  1:772489/1‑124 (MQ=255)
                                                                                        tCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTc                                                                >  2:805644/1‑146 (MQ=255)
                                                                                        tCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTc                                                                <  1:805644/146‑1 (MQ=255)
                                                                                          cATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCGCCGCCTGACGGTGAg                                                                                                                                          >  3:629346/1‑70 (MQ=37)
                                                                                          cATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTc                                                                   >  3:586414/1‑141 (MQ=255)
                                                                                          cATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTc                                                                   <  4:586414/141‑1 (MQ=255)
                                                                                           aTGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAg                                                                                                                                          >  4:680441/1‑69 (MQ=255)
                                                                                            tGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAg                                                         <  3:1441608/149‑1 (MQ=255)
                                                                                                 tcgGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAgg                                                                         <  3:777544/128‑1 (MQ=255)
                                                                                                   ggTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTa                                                  <  2:847127/149‑1 (MQ=255)
                                                                                                        tgcGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAgc                                                        >  4:577708/1‑138 (MQ=255)
                                                                                                        tgcGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGtat                                             >  3:1009635/1‑149 (MQ=255)
                                                                                                             gtgtCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATc                                        >  3:374686/1‑149 (MQ=255)
                                                                                                             gtgtCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATc                                        >  3:114864/1‑149 (MQ=255)
                                                                                                              tgtCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCt                                       <  4:1009635/149‑1 (MQ=255)
                                                                                                               gtCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCtt                                      <  4:114864/149‑1 (MQ=255)
                                                                                                                tCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTa                                     >  3:1250480/1‑149 (MQ=255)
                                                                                                                 cGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTaa                                    <  4:374686/149‑1 (MQ=255)
                                                                                                                   cAGATAAACATACCCAAAAGAAAACCCTCACCGTCAggcg                                                                                                                                               <  4:315387/40‑1 (MQ=255)
                                                                                                                   cAGATAAACATACCCAAAAGAAAACCCTCACCGTCAggcg                                                                                                                                               >  3:315387/1‑40 (MQ=255)
                                                                                                                      aTAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACatgatg                                                                                                                        <  2:700364/60‑1 (MQ=255)
                                                                                                                          gcATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTg                            <  4:896636/147‑1 (MQ=255)
                                                                                                                          aCATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCAtt                                                   >  1:1078283/1‑125 (MQ=255)
                                                                                                                          aCATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCAtt                                                   <  2:1078283/125‑1 (MQ=255)
                                                                                                                          aCATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTg                            >  3:896636/1‑148 (MQ=255)
                                                                                                                           cATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGt                          <  3:720705/149‑1 (MQ=255)
                                                                                                                            aTACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGtt                         <  4:1250480/149‑1 (MQ=255)
                                                                                                                                ccAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTAt                                                 <  1:771766/121‑1 (MQ=255)
                                                                                                                                ccAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGta                                              >  1:521998/1‑124 (MQ=255)
                                                                                                                                ccAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGta                                              <  2:521998/124‑1 (MQ=255)
                                                                                                                                 cAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGc                                                                                                          <  1:680132/63‑1 (MQ=255)
                                                                                                                                 cAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGc                                                                                                          >  2:680132/1‑63 (MQ=255)
                                                                                                                                 cAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTc                                                                            >  1:521152/1‑93 (MQ=255)
                                                                                                                                 cAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCa                                                                           <  4:954098/94‑1 (MQ=255)
                                                                                                                                 cAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCa                                                                           >  3:954098/1‑94 (MQ=255)
                                                                                                                                 cAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCt                                                                 >  2:281949/1‑104 (MQ=255)
                                                                                                                                  aaaaGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAg                                                                          >  3:525603/1‑94 (MQ=255)
                                                                                                                                  aaaaGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAg                                                                          <  1:1436351/94‑1 (MQ=255)
                                                                                                                                  aaaaGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAg                                                                          >  2:1436351/1‑94 (MQ=255)
                                                                                                                                  aaaaGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCaa                                                              >  1:1275248/1‑106 (MQ=255)
                                                                                                                                  aaaaGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTAt                                                 <  3:69957/119‑1 (MQ=255)
                                                                                                                                  aaaaGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTAt                                                 >  4:69957/1‑119 (MQ=255)
                                                                                                                                      ttaaaCCCTCGCCGCCTGACGGTGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTAt                      <  3:125613/140‑1 (MQ=255)
                                                                                                                                      ttaaaCCCTCGCCGCCTGACGGTGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTAt                      >  4:125613/3‑142 (MQ=255)
                                                                                                                                       aaaaCCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGaa                                                                                              >  4:1438286/1‑69 (MQ=255)
                                                                                                                                       aaaaCCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACg              >  2:1376350/1‑149 (MQ=255)
                                                                                                                                        aaaCCCTCGCCGCCTGACGGTGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGaa                                                                                              >  3:615298/1‑68 (MQ=21)
                                                                                                                                        aaaCCCTCGCCGCCTGACGGTGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGaa                                                                                              <  4:615298/68‑1 (MQ=21)
                                                                                                                                         aaCCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGaa            <  1:1376350/149‑1 (MQ=255)
                                                                                                                                          aCCCTCGCCGCCTGACGGTGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATg                     >  3:697028/1‑139 (MQ=255)
                                                                                                                                           cccTCGCCGCCTGACGGTGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAAtt          >  4:71639/1‑149 (MQ=255)
                                                                                                                                             cTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAAGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAg        >  4:1239604/1‑149 (MQ=255)
                                                                                                                                              tCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGa                                                                                          >  4:1071534/1‑66 (MQ=255)
                                                                                                                                              tCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGa                                                                                          <  3:1071534/66‑1 (MQ=255)
                                                                                                                                                aCCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGCgt     >  2:1450578/1‑149 (MQ=255)
                                                                                                                                                 ccGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGcc                               <  2:966955/122‑1 (MQ=255)
                                                                                                                                                   gTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTc                                                                   <  3:1052088/84‑1 (MQ=255)
                                                                                                                                                   gTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTc                                                                   >  4:1052088/1‑84 (MQ=255)
                                                                                                                                                   gTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGCGTGtt  >  2:1140885/1‑149 (MQ=255)
                                                                                                                                                    tCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAAt                                                             <  3:1413431/89‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AACCAGCATAGGCGGGTTGCTCCAGCAGAATATCGATGCGTTCGGCATCCAGGTCGTTAATGATGATAGTTGGTCTGGACATTTTTACTCCATGTCGTCGGTGCTGCGAGTGTCGCAGATAAACATACCCAAAAGAAAACCCTCACCGTCAGGCGGCGAGGGTTTAACTCACATGATGATACTGACTGTTGCTCACTCTTTGAAGTGATTTGCGTCACATTCAGGGAATTCCTCAATGCAGCGCATTATGTATAAATCTTAATCGCCTTGGTTTATGGAAGACGAATTAGCGTGTT  >  NC_000913/643886‑644181

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