Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,509,337 A→G 6.7% A30A (GCT→GCC)  glk ← glucokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,509,3370AG 6.7% 181.0 / 3.5 60A30A (GCT→GCC) glkglucokinase
Reads supporting (aligned to +/- strand):  ref base A (33/23);  new base G (2/2);  total (35/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.10e-01

CCCAGGTATGGTTGGTCATCGCCACCCAGTCACCGGTAATTGGGCAAGCGATGGCAATACAGCCGTCTTTCACCTCGACCTTATGTTCTTCAAGATAAACGCGAATGACCGCTTCGAGGCTGGGGTAATCAAGCCCTGAATAGGTCTTAGCCTGCGAGATTTCACCACTGGCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTGTAAATACCTGGGG  >  NC_000913/2509189‑2509484
                                                                                                                                                    |                                                                                                                                                   
cccAGGTATGGTTGGTCATCGCCACCCAGTCACCGGTAATTGGGCAAGCGATGGCAATACAGCCGTCTTTCACCTCGACCTTATGTTCTTCAAGATAAACGCGAATGACCGCTTCGAGGCTGGGGTAATCAAGCCCTGAATAGGTCTTa                                                                                                                                                     <  5:287458/149‑1 (MQ=255)
cccAGGTATGGTTGGTCATCGCCACCCAGTCACCGGTAATTGGGCAAGCGATGGCAATACAGCCGTCTTTCACCTCGACCTTATGTTCTTCAAGATAAACGCGAATGACCGCTTCGAGGCTGGGGTAATCAAGCCCTGAATAGGTCTTa                                                                                                                                                     <  5:221708/149‑1 (MQ=255)
    ggTATGGTTGGTCATCGCCACCCAGTCACCGGTAATTGGGCAAGCGATGGCAATACAGCCGTCTTTCACCTCGACCTTATGTTCTTCAAGATAAACGCGAATGACCGCTTCGAGGCTGGGGTAATCAAGCCCTGAATAGGTCTTAGCCt                                                                                                                                                 <  1:336316/149‑1 (MQ=255)
    ggTATGGTTGGTCATCGCCACCCAGTCACCGGTAATTGGGCAAGCGATGGCAATACAGCCGTCTTTCACCTCGACCTTATGTTCTTCAAGATAAACGCGAATGACCGCTTCGAGGCTGGGGTAATCAAGCCCTGAATAGGTCTTAGCCt                                                                                                                                                 >  7:231071/1‑149 (MQ=255)
     gTATGGTTGGTCATCGCCACCCAGTCACCGGTAATTGGGCAAGCGATGGCAATACAGCCGTCTTTCACCTCGACCTTATGTTCTTCAAGATAAACGCGAATGACCGCTTCGAGGCTGGGGTAATCAAGCCCTGAATAGGTCTTAGCCTg                                                                                                                                                <  2:287133/149‑1 (MQ=255)
          gttggtCATCGCCACCCAGTCACCGGTAATTGGGCAAGCGATGGCAATACAGCCGTCTTTCACCTCGACCTTATGTTCTTCAAGATAAACGCGAATGACCGCTTCGAGGCTGGGGTAATCAAGCCCTGAATAGGTCTTAGCCTGc                                                                                                                                               >  7:44971/1‑145 (MQ=255)
              gtCATCGCCACCCAGTCACCGGTAATTGGGCAAGCGATGGCAATACAGCCGTCTTTCACCTCGACCTTATGTTCTTCAAGATAAACGCGAATGACCGCTTCGAGGCTGGGGTAATCAAGCCCTGAATAGGTCTTAGCCTGCGAGAttt                                                                                                                                        <  2:60751/148‑1 (MQ=255)
                 aTCGCCACCCAGTCACCGGTAATTGGGCAAGCGATGGCAATACAGCCGTCTTTCACCTCGACCTTATGTTCTTCAAGATAAACGCGAATGACCGCTTCGAGGCTGGGGTAATCAAGCCCTGAATAGGTCTTAGCCTGCGAGATTTcac                                                                                                                                     <  7:72708/148‑1 (MQ=255)
                         ccAGTCACCGGTAATTGGGCAAGCGATGGCAATACAGCCGTCTTTCACCTCGACCTTATGTTCTTCAAGATAAACGCGAATGACCGCTTCGAGGCTGGGGTAATCAAGCCCTGAATAGGTCTTAGCCTGCGAGATTTCACCACTGGCaa                                                                                                                            <  6:126147/149‑1 (MQ=255)
                           aGTCACCGGTAATTGGGCAAGCGATGGCAATACAGCCGTCTTTCACCTCGACCTTATGTTCTTCAAGATAAACGCGAATGACCGCTTCGAGGCTGGGGTAATCAAGCCCTGAATAGGTCTTAGCCTGCGAGATTTCACCACTGGCAata                                                                                                                          <  7:73956/149‑1 (MQ=255)
                               aCCGGTAATTGGGCAAGCGATGGCAATACAGCCGTCTTTCACCTCGACCTTATGTTCTTCAAGATAAACGCGAATGACCGCTTCGAGGCTGGGGTAATCAAGCCCTGAATAGGTCTTAGc                                                                                                                                                   <  4:291201/120‑1 (MQ=255)
                               aCCGGTAATTGGGCAAGCGATGGCAATACAGCCGTCTTTCACCTCGACCTTATGTTCTTCAAGATAAACGCGAATGACCGCTTCGAGGCTGGGGTAATCAAGCCCTGAATAGGTCTTAGc                                                                                                                                                   >  3:291201/1‑120 (MQ=255)
                                ccGGTAATTGGGCAAGCGATGGCAATACAGCCGTCTTTCACCTCGACCTTATGTTCTTCAAGATAAACGCGAATGACCGCTTCGAGGCTGGGGTAATCAAGCCCTGAATAGGTCTTAGCCTGCGAGATTTCACCACTGGCAATATcaca                                                                                                                     >  7:276061/1‑149 (MQ=255)
                                  ggTAATTGGGCAAGCGATGGCAATACAGCCGTCTTTCACCTCGACCTTATGTTCTTCAAGATAAACGCGAATGACCGCTTCGAGGCTGGGGTTATCAAGCCCTGAATAGGTCTTAGCCTGCGAGTTTTCACCACTGGCAATATcaca                                                                                                                     >  2:292299/1‑147 (MQ=255)
                                  ggTAATTGGGCAAGCGATGGCAATACAGCCGTCTTTCAACTCGACCTTATGTGATTCAAGATAAACGCGAATGACCGCTTCGAGGCTGGGGTAATCAAGCCCTGAATAGGTCTTAGCCTGCGAGATTTCACCACTGGCAATATcacaca                                                                                                                   <  6:41298/149‑1 (MQ=255)
                                       ttGGGCAAGCGATGGCAATACAGCCGTCTTTCACCTCGACCTTATGTTCTTCAAGATAAACGCGAATGACCGCTTCGAGGCTGGGGTAATCAAGCCCTGAATAGGTCTTAGCCTGCGAGATTTCACCACTGGCAATATCACACAGAGc                                                                                                               >  6:119743/1‑148 (MQ=255)
                                          ggCAAGCGATGTCAATACAGCCGGCTTTCACCTCGACCTTATGTTCTTCAAGATAAACGCGAATGACCGCTTCGAGGCTGGGGTAATCAAGCCCTGAATAGGTCTTAGCCTGCGAGATTTCACCACTGGCAATATCACACAGAGCAAGa                                                                                                           <  5:155654/149‑1 (MQ=255)
                                              aGCGATGGCAATACAGCCGTCTTTCACCTCGACCTTATGTTCTTCAAGATAAACGCGAATGACCGCTTCGAGGCTGGGGTAATCAAGCCCTGAATAGGTCTTGGCCTGCGAGATTTCACCACTGGCAATATCACACAGAGCAAGACGTg                                                                                                       >  1:8445/1‑149 (MQ=255)
                                                  aTGGCAATACAGCCGTCTTTCACCTCGACCTTATGTTCTTCAAGATAAACGCGAATGACCGCTTCGAGGCTGGGGTAATCAAGCCCTGAATAGGTCTTGGCCTGCGAGATTTCACCACTGGCAATATCACACAGAGCAAGACGTGCGtt                                                                                                   >  6:15274/1‑149 (MQ=255)
                                                     gCAATACAGCCGTCTTTCACCTCGACCTTATGTTCTTCAAGATAAACGCGAATGACCGCTTCGAGGCTGGGGTAATCAAGCCCTGAATAGGTCTTAGCCTGCGAGATTTCACCACTGGCAATATCACACAGAGCAAGACGTGCGTTGGt                                                                                                <  6:100592/149‑1 (MQ=255)
                                                      cAATACAGCCGTCTTTCACCTCGACCTTATGTTCTTCAAGATAAACGCGAATGACCGCTTCGAGGCTGGGGTAATCAAGCCCTGAATAGGTCTTAGCCTGCGAGCTTTCACCACTGGCAATATCACACAGAGCAAGACGTGCGTTGGTg                                                                                               >  2:168432/1‑149 (MQ=255)
                                                              ccGTCTTTCACCTCGACCTTATGTTCTTCAAGATAAACGCGAATGACCGCTTCGAGGCTGGGGTAATCAAGCCCTGAATAGGTCTTAGCCTGCGAGATTTCACCACTGGCAATATcaca                                                                                                                     >  6:338277/1‑119 (MQ=255)
                                                              ccGTCTTTCACCTCGACCTTATGTTCTTCAAGATAAACGCGAATGACCGCTTCGAGGCTGGGGTAATCAAGCCCTGAATAGGTCTTAGCCTGCGAGATTTCACCACTGGCAATATcaca                                                                                                                     <  5:338277/119‑1 (MQ=255)
                                                              ccGTCTTTCACCTCGACCTTATGTTCTTCAAGATAAACGCGAATGACCGCTTCGAGGCTGGGGTAATCAAGCCCGGAATAGGTCTTAGCCTGCGAGATTTCACCACTGGCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCca                                                                                        >  8:306587/1‑148 (MQ=255)
                                                               cGTCTTTCACCTCGACCTTATGTTCTTCAAGATAAACGCGAATGACCGCTTCGAGGCTGGGGTAATCAAGCCCTGAATAGGTCTTAGCCTGCGAGATTTCACCACTGGCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCcaca                                                                                      >  7:95250/1‑149 (MQ=255)
                                                                gTCTTTCACCTCGACCTTATGTTCTTCAAGATAAACGCGAATGACCGCTTCGAGGCTGGGGTAATCAAGCCCTGAATAGGTCTTAGCCTGCGAGATTTCACCACTGGCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACAt                                                                                     >  7:277021/1‑149 (MQ=255)
                                                                   tttCACCTCGACCTTATGTTCTTCAAGATAAACGCGAATGACCGCTTCGAGGCTGGGGTAATCAAGCCCTGAATAGGTCTTAGCCTGCGAGATTTCACCACTGGCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACATCAc                                                                                  >  8:173071/1‑149 (MQ=255)
                                                                            gACCTTATGTTCTTCAAGATAAACGCGAATGACCGCTTCGAGGCTGGGGTAATCAAGCCCTGAATAGGTCTTAGCCTGCGAGATTTCACCACTGGCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAAt                                                                          <  8:127382/148‑1 (MQ=255)
                                                                              ggttATGTTCTTCAAGATAAACGCGAATGACCGGTTCGAGGCTGGGGTAATCAATCAATGAATAGGTCTTAGCCTGCGAGCTTACCCCACTGGCAATATCACACAGCGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCa                                                                       >  8:255920/3‑149 (MQ=255)
                                                                                          cAAGATAAACGCGAATGACCGCTTCGAGGCTGGGGTAATCAAGCCCTGAATAGGTCTTAGCCTGCGAGATTTCACCACTGGCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCAtt                                                           >  1:69906/1‑149 (MQ=255)
                                                                                           aaGATAAACGCGAATGACCGCTTCGAGGCTGGGGTAATCAAGCCCTGAATAGGTCTTAGCCTGCGAGATTTCACCACTGGCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCAttc                                                          <  5:119743/149‑1 (MQ=255)
                                                                                             gATAAACGCGAATGACCGCTTCGAGGCTGGGGTAATCAAGCCCTGAATAGGTCTTGGCCTGCGAGATTTCACCACTGGCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCAttctt                                                        <  2:8445/149‑1 (MQ=255)
                                                                                                  aCGCGAATGACCGCTTCGCGGCTGGGGTAATCAAGCCCTGAATAGGTCTTAGCCTGCGAGATTTCACCACTGGCAATATCACACAGCGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACt                                                   >  5:156347/1‑149 (MQ=255)
                                                                                                  aCGCGAATGACCGCTTCGAGGCTGGGGTAATCAAGCCCTGAATAGGTCTTAGCCTGCGAGATTTCACCACTGGCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACAt                                                                                     <  2:181292/115‑1 (MQ=255)
                                                                                                  aCGCGAATGACCGCTTCGAGGCTGGGGTAATCAAGCCCTGAATAGGTCTTAGCCTGCGAGATTTCACCACTGGCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACAt                                                                                     >  1:181292/1‑115 (MQ=255)
                                                                                                    gcgAATGACCGCTTCGAGGCTGGGGTAATCAAGCCCTGAATAGGTCTTAGCCTGCGAGATTTCACCACTGGCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTg                                                                >  2:95822/1‑134 (MQ=255)
                                                                                                    gcgAATGACCGCTTCGAGGCTGGGGTAATCAAGCCCTGAATAGGTCTTAGCCTGCGAGATTTCACCACTGGCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTg                                                                <  1:95822/134‑1 (MQ=255)
                                                                                                    gcgAATGACCGCTTCGAGGCTGGGGTAATCAAGCCCTGAATAGGTCTTAGACTGCGAGATTTCACCACTGGCAATATCACACAGCGCAAGACGTTCTTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCAGTCTTCAACTGc                                                 >  1:226908/1‑149 (MQ=255)
                                                                                                      gAATGACCGCTTCGAGGCTGGGGTAATCAAGCCCTGAATAGGTCTTAGCCTGCGAGATTTCACCACTGGCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTc                                               >  4:111151/1‑149 (MQ=255)
                                                                                                       aaTGACCGCTTCGAGGCTGGGGTAATCAAGCCCTGAATAGGTCTTGGCCTGCGAGATTTCACCACTGGCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTcc                                              <  5:15274/149‑1 (MQ=255)
                                                                                                            ccGCTTCGAGGCTGGGGTAATCAAGCCCTGAATAGGTCTTAGACTGCGAGATTTCACCACTGGCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACATCAACGCCTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTa                                          >  4:13932/1‑148 (MQ=255)
                                                                                                                ttCGAGGCTGGGGTAATCAAGCCCTGAATAGGTCTTAGCCTGCGAGATTTCACCACTGGCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGt                                      >  8:340672/1‑148 (MQ=255)
                                                                                                                   gAGGCTGGGGTAATCAAGCCCTGAATAGGTCTTAGCCTGCGAGATTTCACCACTGGCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCaa                                   >  5:241305/1‑148 (MQ=255)
                                                                                                                          gggTAATCAAGCCCTGAATAGGTCTTAGCCTGCGAGATTTCACCACTGGCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAAttc                           >  6:187633/1‑149 (MQ=255)
                                                                                                                          gggTAATCAAGCCCTGAATAGGTCTTAGCCTGCGAGATTTCACCACTGGCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTATGTCATTCTTCAACTGCACCGCAAAAGTCAAAATAAttc                           >  8:175046/1‑149 (MQ=255)
                                                                                                                           ggTAATCAAGCCCTGAATAGGTCTTAGCCTGCGAGATTTCACCACTGGCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAAttct                          >  8:298205/1‑149 (MQ=255)
                                                                                                                            gTAATCAAGCCCTGAATAGGTCTTAGCCTGCGAGATTTCACCACTGGCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAAttctt                         >  8:280906/1‑149 (MQ=255)
                                                                                                                             tAATCAAGCCCTGAATAGGTCTTAGCCTGCGAGATTTCACCACTGGCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAAttcttt                        <  7:173071/149‑1 (MQ=255)
                                                                                                                                tCAAGCCCTGAATAGGTCTTAGCCTGCGAGATTTCACCACTGGCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAAttctttct                      >  4:94067/1‑148 (MQ=255)
                                                                                                                                  aaGCCCTGAATAGGTCTTAGCCTGCGAGATTTCACCACTGGCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTca                    >  2:262649/1‑148 (MQ=255)
                                                                                                                                   aGCCCTGAATAGGTCTTAGCCTGCGAGATTTCACCACTGGCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTcaca                  <  8:277021/149‑1 (MQ=255)
                                                                                                                                   aGCCCTGAATAGGTCTTAGCCTGCGAGATTTCACCACTGGCAATATCACACAGAGCAAGACGTGCGTTGGTGACGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTcaca                  <  2:33906/149‑1 (MQ=255)
                                                                                                                                   aGCCCTGAATAGGTCTTAGCATGCGAGATTTCACCACTGGCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAAGTCTTTCTcaca                  >  4:181846/1‑149 (MQ=255)
                                                                                                                                     cccTGAATAGGTCTTAGCCTGCGAGATTTCACAACTGGCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACt                >  1:34194/1‑149 (MQ=255)
                                                                                                                                        tGAATAGGTCTTAGCCTGCGAGATTTCACCACTGGCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTGt              >  2:17204/1‑148 (MQ=255)
                                                                                                                                          aaTAGGTCTTAGCCTGCGAGATTTCACCACTGGCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTGTaaa           <  7:280906/149‑1 (MQ=255)
                                                                                                                                              ggTCTTAGCCTGCGAGATTTCACCACTGGCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTGTAAATAc        <  8:231071/148‑1 (MQ=255)
                                                                                                                                               gTCTTAGCCTGCGAGATTTCACCACTGGCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTGTAAATACCt      <  5:187633/149‑1 (MQ=255)
                                                                                                                                                tCTTAGCCTGCGAGATTTCACCACTGGCAATATCACACAGAGCAAGCCGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTGTAAATACCTg     <  7:306587/149‑1 (MQ=255)
                                                                                                                                                tCTTAGCCTGCGAGATTTCACCACTGGCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTGTAAATACCTg     <  1:238621/149‑1 (MQ=255)
                                                                                                                                                 cTTAGCCTGCGAGATTTCACCACTGGCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTGTAAATACCTgg    >  4:311606/1‑149 (MQ=255)
                                                                                                                                                  ttAGCCTGCGAGATTTCACCACTGGCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTGTAAATACCTggg   >  2:229714/1‑149 (MQ=255)
                                                                                                                                                    aGCCTGCGAGATTTCACCACTGGCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTGTAAATACCTgggg  <  1:229714/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCCAGGTATGGTTGGTCATCGCCACCCAGTCACCGGTAATTGGGCAAGCGATGGCAATACAGCCGTCTTTCACCTCGACCTTATGTTCTTCAAGATAAACGCGAATGACCGCTTCGAGGCTGGGGTAATCAAGCCCTGAATAGGTCTTAGCCTGCGAGATTTCACCACTGGCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTGTAAATACCTGGGG  >  NC_000913/2509189‑2509484

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