Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,789,286 G→A G101D (GGT→GAT)  glaH → glutarate dioxygenase GlaH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,789,2860GA98.0% 140.1 / NA 49G101D (GGT→GAT) glaHglutarate dioxygenase GlaH
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (24/24);  minor base C (1/0);  total (25/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGATTCACCCTCACCCCGTCGGCGCAATCCCCGCGTCTGCTGGAACTCACCTTCACCGAACAGACGACCAAACAGTTTCTTGAGCAGGTTGCCGAGTGGCCCGTGCAGGCGCTGGAGTACAAATCGTTTCTGCGTTTTCGGGTAGCCAAAATTCTTGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGGTGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGAGCGGTCAGTATTACGCGCGTTTTGTAGTGAAAAATGTCGATAACTCAGACAGCTATCTGCGCCAGCCGCACCGCGTGATGGAGCTACACAACGACGGCACGTATGTCGAAGAGATAACCGATTACGTGCTGATGATGAAAATCGACGAGCAAAACA  >  NC_000913/2789041‑2789531
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       
cGGATTCACCCTCACCCCGTCGGCGCAATCCCCGCGTCTGCTGGAACTCACCTTCACCGAACAGACGACCAAACAGTTTCTTGAGCAGGTTGCCGAGTGGCCCGTGCAGGCGCTGGAGTACAAATCGTTTCTGCGTTTTCGGGTAGCCAAAATTCTTGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCg                                                                                                                                                                                                                                                    >  1:29908/1‑251 (MQ=255)
   aTTCACCCTCTCCCCGACTGCGCAATCCCCGCGTCTGCTGGAACTCTCCTTCACTGAACCGCCGCCCAAACATTTTCTTGAGCAGGTTTCCGAGTGGCCTGTGCATGCGCTTGAGTACAAATCGTTTCTGCGTTTTCGGGTAGCCAAAATTCTTGACGATGTATGTGCGAACCAGCTGCTGCCATTTCTGGTGAAGACGTTGCTAGACCGCGCTGAAGGTGCGTTGTTGATCAATGCGTTGGATGTCgatg                                                                                                                                                                                                                                                 <  2:345228/251‑1 (MQ=255)
        ccctcacccCTGCGGCGCAATCCCCGCGGCTGCTGGAACTCACCTTCACCGAACAGACGACCAAATAGTTTCTTGGGCAGGTTGCCGAGTGGCCCGTGCAGGCGCTGGAGTATAAATCGTTTCTGCGTTTTCGGGTAGCCAAAATTCTTGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGAtgtg                                                                                                                                                                                                                                            <  2:179102/251‑1 (MQ=255)
            cacccCGTCGGCGCAATCCCCGCGTCTGCTGGAACTCACCTTCACCGAACAGACGACCAAACAGTTTCTTGAGCAGGTTGCCGAGTGGCCCGTGCAGGCGCTGGAGTACAAATCGTTTCTGCGTTTTCGGGTAGCCAAAATTCTTGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGaagc                                                                                                                                                                                                                                        <  2:175305/251‑1 (MQ=255)
                  gTCGGCGCACTCCCCGCGACTGCTGGAACTCACCTTCACCGAACAGCCGAACAAACAGTTTCTTGAGCAGGTTGCCGAGTGGCCCGCGCAGGCGCTGGAGTACACATCGTTTCTGCGTTTTCGGGTAGCCAAAATTCTTGACGATCTATGTGCGAACCAACTGCAACCATTGCTGTTGAAGGCGTTGCTAAACCGCGCTGATGGTGCGTTGTTGCTCCATGCGGTGGCTGTCGATGCTGTGGAGCAAGCgg                                                                                                                                                                                                                                  >  2:199509/1‑251 (MQ=255)
                  gTCGGCGCAATCCCCGCGTCTGCTGGAACTCACCTTCACCGAACAGACGACCAAACAGTTTCTTGAGCAGGTTGCCGAGTGGCCCGTGCAGGCGCTGGAGTACAAATCGTTTCTGCGTTTTCGGGTAGCCAAAATTCTTGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCgg                                                                                                                                                                                                                                  >  1:179330/1‑251 (MQ=255)
                        gcAATCCCCGCGTCTGCTGGAACTCACCTTCACCGAACAGACGACCAAACAGTTTCTTGAGCAGGTTGCCGAGTGGCCCGTGCAGGCGCTGGAGTACAAATCGTTTCTGCGTTTTCGGGTAGCCAAAATTCTTGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATg                                                                                                                                                                                                                               >  1:522796/1‑248 (MQ=255)
                             ccccGTGTCTGCTTGAACTCACCTTCACCGAACAGACGACCAAACAGTATCTTGAGCAGGTTTCCGAGTGGCCCGTGCAGGCGCTGGAGTACAAATCGTTTCTTCGTTTTCGGGTAGCCAAAATTCTTGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGGGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGt                                                                                                                                                                                                                        <  2:341462/250‑1 (MQ=255)
                                cgcgTCTGCTGGAACTCACCTTCACCGAACAGTCGACCAAACAGTTTCTTGAGCAGGTTGCCGAGTGGCCCGTGCAGGCGCTGGAGTACAAATCGGTTCTGCGTTTTCGGGTAGCCAAAATTCTTGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAg                                                                                                                                                                                                                    <  2:169679/251‑1 (MQ=255)
                                   gTCTTCTGGAACACACCTTCCCCGAACGGACGACCAAACAGTTTCTTGAGCAGGTTGCCGAGGGGCCCGTGCAGGCGCTGGAGTACAAAGCGGTTCTGCGTTGTCGGGTAGCCAAAATTCTGGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCTGTGGATGTCGATGATGTGAAGCATGCGGATGAGATGGTGAAGCTg                                                                                                                                                                                                                 <  2:395672/251‑1 (MQ=255)
                                       gctgGAACTCACCTTCACCGAACAGACGACCAAACAGTTTCTTGAGCAGGTTGCCGAGTGGCCCGTGCAGGCGCTGGAGTACAAATCGTTTCTGCGTTTTCGGGTAGCCAAAATTCTTGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCgg                                                                                                                                                                                                             >  1:496768/1‑250 (MQ=255)
                                        ctgGAACTCACCTTCACCGAACAGACGACCAAACAGTTTCTTGAGCAGGTTGCCGAGTGGCCCGTGCAGGCGCTGGAGTACAAATCGTTTCTGCGTTTTCGGGTAGCCAAAATTCTTGACGATCTATGTGCGCACCAGCTGCAACCATTGCTGTTGCAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCAGTGGATGTCGATGGTGTGAAGCAAGCGGGTGAGATGGTGAAGCTGGCg                                                                                                                                                                                                              >  2:15416/1‑249 (MQ=255)
                                                 aCCTTCCCCGAACAGACGACCAAACAGTTTCTTGAGCAGGCTGCCGAGTGGCCCGTGCAGGCGCTGGAGTACAAATCGTTTCTGCGTTTTCGGGTAGCCAAAATTCTTGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGgc                                                                                                                                                                                                   <  2:216949/251‑1 (MQ=255)
                                                              aGACGACCAAACAGTTTCTTGAGCAGGTTGCCGAGTGGCCCGTGCATGCGCTGGAGTACAACTCGTTTATGCGTTTTCGGGTAGCCAAAATCCTTGACGATCTATGTGCGAACCAGCTGAAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATTTCGTTGATGGGAAGCAAGCGGCTGAGATGGGGAAGCTGGCTACGGCGGTGGCGCATCTGCTTggg                                                                                                                                                                                      >  2:330535/1‑250 (MQ=255)
                                                                           gtttcgtgagatggTTGCCGAGTGTACCGTGCAGGCGCTGGAGTACAAAGCGTGTCTGCGTTTGCGTGTAGCTAAAATTCTTGACGAGCTATTGGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTTAAGCTGGCGACGGCGGTGGCGCAGCTGATTGGCCGTTCCAATTTc                                                                                                                                                                          <  2:24353/238‑1 (MQ=255)
                                                                                 tGAGCAGGTTGCCGAGTGGCCCGTGCAGGCTCTTGAGTACAAATCGTGTCGGCGTTTGCGGGTAGCCAAAATTCTTGACGATCTATGTGCGAACCAGCTGCATCCATGTCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGACGCTGGCGACGGCGGTGGCGCATCTGATGGGCCGTTGCAATTTcgatgc                                                                                                                                                                     <  2:10908/249‑1 (MQ=255)
                                                                                     cAGGTTGCCGAGTGGCCCGTGCAGGCGCTGGAGTACAAATCGTTTCTGCGTTTTCGGGTAGCCAAAATTCTTGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACCGCGGTGGGGCATCTGATTGGCCGGTCCAATTTcgatgcgatggg                                                                                                                                                               >  1:319971/1‑249 (MQ=255)
                                                                                      aGGTTGCCGAGTGGCCCGTGCAGGCGCTGGAGTACAAATCGTTTCTGCGTTTTCGGGTAGCCAAAATTCTTGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGAGc                                                                                                                                                              <  1:399896/251‑1 (MQ=255)
                                                                                           gCCGAGTGGCCCGTGCAGGCGCTGGAGTACAAATCGTGTCTGCGTTTTCGGGTAGCCAAAATTCTTGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGAGCGGTCa                                                                                                                                                         <  2:496768/251‑1 (MQ=255)
                                                                                                  ggCCCGTGCAGGCGCTGGAGTACAAATCGTTTCTGCGTTTTCGGGTAGCCAAAATTCTTGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGGGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGAGCGGTCAGTATTAc                                                                                                                                                  <  1:3948/251‑1 (MQ=255)
                                                                                                    cccGTGCAGGCGCTGGAGTACAAATCGTTTCTGCGTTTTCGGGTAGCCAAAATTCTTGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGAGGTGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGAGCGGTCAGTATTAcgc                                                                                                                                                >  1:66943/1‑251 (MQ=255)
                                                                                                     ccGTGCAGGCGCTGGAGTACAAATCGTTTCTGCGTTTTCGGGTAGCCAAAATTCTTGACGATCTATGTGCGAACCAGCTGAAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGGGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGCGATGGTGAAGCTGGCGACGGCGGGGGGGCCTCTGATGGGGCGTTACAATTTcgatgcagtgtgaggtcaagatatcgcc                                                                                                                                               >  2:413557/1‑229 (MQ=255)
                                                                                                                  ggAGTACAAATCGTTTCTGCGTTTTCGGGTAGCCAAAATTCTTGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGAGCGGTCAGTATTACGCGCGTTTTGTAGTGa                                                                                                                                 <  2:179330/251‑1 (MQ=255)
                                                                                                                  ggAGTACAAAGCGTTTCTGCGTTTTCGGGTAGCCAAAATTCTTGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGAGCGGTCAGTATTACGCGCGTTTGGTAGt                                                                                                                                    <  1:306623/249‑1 (MQ=255)
                                                                                                                         aaaTCGTTTCTGCGTTTTCGGGTAGCCAAAATTCTTGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGAGCGGTCAGTATTACGCGCGTTTTGTAGTGAAAAATGt                                                                                                                          >  1:31190/1‑251 (MQ=255)
                                                                                                                          aaTCGTTTCTGCGTTTTCGGGTCGCCAAAATTCTTGACGCTCTATGTGCGAACCCGCTGCAACCATTGCTGTTGAAGACGTTGCTACACCGCGCTGCAGGTGCGTTGTTGATCAATGCGGCGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGACGCTGGCGACGGCGGCGGCGCATCTGATTGGCCGTTCCCATTTCCATGCGTACAGCGGTCAGCATTACGCGCGTTTTGTCGTGAAAAAtctc                                                                                                                         >  2:52473/1‑248 (MQ=255)
                                                                                                                                        tttCGGGTAGCCAAAATTCTTGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGATGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGGGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGAGGCGACTCTGATTGGCCGTTCCGATTTTGATGCGATGCGCGGTCCGGATTACGCGCGTTTTGTAGTGATAAATGGCAATAACTCAggc                                                                                                             >  2:367290/1‑247 (MQ=255)
                                                                                                                                                   cAAAATTCTTGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGAGCGGTCAGTATTACGCGCGTTTTGTa                                                                                                                                      <  1:43645/214‑1 (MQ=255)
                                                                                                                                                   cAAAATTCTTGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGGGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGAGCGGGCAGTATTACGCGCGTTTTGTa                                                                                                                                      >  2:43645/1‑214 (MQ=255)
                                                                                                                                                             gACGATCTATGTGCGAACCAGCTGCAACCCTTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGATGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGAAAGCGGATGAGATGGTGAAGCTGGCGACGGCGCGGGAGCATCTGATTGGCCGATCCAATTTCGAAGCGGAGAGCGGTCCGTATTACGCGCGGTTTGTAGGGCAAAAATGTGATAACTCAGGAAGCTAACTGCGCCAGCCGCAcc                                                                                      >  2:547722/1‑250 (MQ=255)
                                                                                                                                                               cGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCCATTTCGATGCGATGAGCGGTCTGTGTTACGCGCGTTTTGTAGTGCAAAATGTCGATAAACCAGAGCGCTCTCTGCGCCAGCCGCACcgcg                                                                                   >  2:14038/1‑251 (MQ=255)
                                                                                                                                                               cGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGa                                                                                                                                                                >  1:81694/1‑176 (MQ=255)
                                                                                                                                                               cGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTCGACGCGATGa                                                                                                                                                                <  2:81694/176‑1 (MQ=255)
                                                                                                                                                                           gAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCCATTTCGATGCGATGAGCGGTCAGTATTACGCGCGTTTTGTTGTGAAAAATGTCGATAACTCAGACAGCTATCTGCGCCAGCCGccacgcggggtgggacatc                                                                       >  2:265985/1‑233 (MQ=255)
                                                                                                                                                                              ccAGCTGCAACCATTGCTGTTGAAGAAGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGAGCGGGCAGTATTACGCGCGTTTTGTAGTGAAAAATGTCGATAACTCAGACAGCTATCTGCGCCAGCCGCACCGCGGGATGGAGCTacac                                                                     >  1:141674/1‑250 (MQ=255)
                                                                                                                                                                                         cATTGCTGTTGAAGACGTTGCTAAACCGCGATGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGAGCGGTCAGTATTACGCGCGTTTTGTAGGGTAAAATGTCGATAACTCAGAACACTATCTGCGCCAGCCGCACAGCGGGATGGAGCTACACAacgacgacccg                                                         >  2:374297/1‑246 (MQ=255)
                                                                                                                                                                                           tggttgtggAAGTCGTTGCTAAACCGCGCTGAAGGTGCGGTGTGGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGAGCGGTCAGTATTACGCGCGTTTTGTAGTGAAAAATGTCGATAACTCAGACAGCTATCTGCGCCAGCCGCACCGCGTGATGGAGCTACACAACGACGGCACg                                                         <  2:319971/241‑1 (MQ=255)
                                                                                                                                                                                                gTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGCGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGAGCGGTCAGTATTACGCGCGTTTTGTAGTGAAAAATGTCGATAACTCAGACAGCTATCTGCGCCAGCCGCACCGCGTGATGGAGCTACACAACGACGGCACGTATGTc                                                   <  1:456810/250‑1 (MQ=255)
                                                                                                                                                                                                   gAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGAGCGGTCAGTATTACGCGCGTTTTGTAGTGAAAAATGTCGATAACTCAGACAGCTATCTGCGCCAGCCGCACCGCGTGATGGAGCTACACAACGACGGCACGTATGTCGAag                                               <  1:15416/251‑1 (MQ=255)
                                                                                                                                                                                                            gCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGAGCGGTCAGTATTACGCGCGTTTTGTAGTGAAAAATGTCGATAACTCAGACAGCTATCTGCGcc                                                                                               <  1:334611/194‑1 (MQ=255)
                                                                                                                                                                                                              aaaaCCGCGTGGAAGGTGCGTTGTTGATAAATGCGGTGTATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTGCGATGCGATGAGCGGTCAGTATTACGCGCGTTTTGTAGTGAAAAATGTCGATAACTCAGTCAGCTATCTGCGCCAGCCGCAACGCGTGATGGAGCTACACAACGACGGCACGTATGTCGAAGAGATATCCGAt                                    <  2:518182/250‑1 (MQ=255)
                                                                                                                                                                                                                aaCCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGAGCGGTCAGTATTACGCGCGTTTTGTAGTGAAAAATGTCGATAACTCAGACAGCTATCTGCGCCAGCCGCACCGCGTGATGGAGCTACACAACGACGGCACGTATGTCGAAGAGATAACCGATTa                                  <  2:522796/251‑1 (MQ=255)
                                                                                                                                                                                                                    gcgcTGAAGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGAGCGGTCAGTATTACGCGCGTTTTGTAGTGAAAAATGTCGATAACTCAGACAGCTATCTGCGCCAGCCGCACCGCGTGATGGAGCTACACAACGACGGCACGTATGTCGAAGAGATAACCGATTACGTg                              <  1:34754/251‑1 (MQ=255)
                                                                                                                                                                                                                           aGGTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGAGCGGTCAGTATTACGCGCGTTTTGTAGTGAAAAATGTCGATAACTCAGACAGCTATCTGCGCCAGCCGCACCGCGTGATGGAGCTACACAACGACGGCACGTATGTCGAAGAGATAACCGATTACGTGCtgatga                       <  1:199509/251‑1 (MQ=255)
                                                                                                                                                                                                                             gTGCGTTGTTGATCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGAGCGGTCAGTATTACGCGCGTTTTGTAGTGAAAAATGTCGATAACTCAGACAGCTATCTGCGCCAGCCGCACCGCGTGATGGAGCTACACAACGACGGCACGTATGTCGAAGAGATAACCGATTACGTGCtgatggtg                     >  1:294764/1‑248 (MQ=255)
                                                                                                                                                                                                                                        aTCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGAGCGGTCAGTATTACGCGCGTTTTGTAGGGAAAAATGTCGATAACTCAGACAGCTATCTGCGCCAGCCGCACCGCGTGATGGAGCTACACAACGACGGCACGTATGTCGAAGAGATAACCGATTACGTGCTGATGTTGAAAATcgacg           >  1:565668/1‑250 (MQ=255)
                                                                                                                                                                                                                                        aTCAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATATGATTGGCCGTTTCAATTTCGATGCGATGAGCGGTCAGTCTTACCCGCGTTTTGTAGTGAAAAATGTCGATAACTCAGACAGCTATCTGCGCCCGCCGGACCGCGTGGTGGAGCTTCCCAACGACGGCCCGTATGGCATAACGTTAACCGATTACGTGccgtatatgaaaatagacc           >  2:437482/1‑232 (MQ=255)
                                                                                                                                                                                                                                          cAATGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGAGCGGTCAGTATTACGCGCGTTTTGTAGTGAAAAATGTCGATAACTCAGACAGCTATCTGCGCCAGCCGCACCGCGTGATGGAGCTACACAACGACGGCACGTATGTCGAAGAGATAACCGATTACGTGCTGATGATGAAAATCGACGAGc        >  1:44426/1‑251 (MQ=255)
                                                                                                                                                                                                                                           aaTGCGGTGGATGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGAGCGGTCAGTATTACGCGCGTTTTGTAGTGAAAAATGTCGATAACTCAGACAGCTATCTGCGCCAGCCGCACCGCGTGATGGAGCTACACAACGACGGCACGTATGTCGAAGAGATAACCGATTACGTGCTGATGATGAAAATCGACGAGc        >  1:105199/1‑250 (MQ=255)
                                                                                                                                                                                                                                                ggtggatgTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGAGCGGTCAGTATTACGCGCGTTTTGTAGTGAAAAATGTCGATAACTCAGACAGCTATCTGCGCCAGCCGCACCGCGTGATGGAGCTACACAACGACGGCACGTATGTCGAAGAGATAACCGATTACGTGCTGATGATGAAAATCGACGAGCAAAACa  <  1:209714/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       
CGGATTCACCCTCACCCCGTCGGCGCAATCCCCGCGTCTGCTGGAACTCACCTTCACCGAACAGACGACCAAACAGTTTCTTGAGCAGGTTGCCGAGTGGCCCGTGCAGGCGCTGGAGTACAAATCGTTTCTGCGTTTTCGGGTAGCCAAAATTCTTGACGATCTATGTGCGAACCAGCTGCAACCATTGCTGTTGAAGACGTTGCTAAACCGCGCTGAAGGTGCGTTGTTGATCAATGCGGTGGGTGTCGATGATGTGAAGCAAGCGGATGAGATGGTGAAGCTGGCGACGGCGGTGGCGCATCTGATTGGCCGTTCCAATTTCGATGCGATGAGCGGTCAGTATTACGCGCGTTTTGTAGTGAAAAATGTCGATAACTCAGACAGCTATCTGCGCCAGCCGCACCGCGTGATGGAGCTACACAACGACGGCACGTATGTCGAAGAGATAACCGATTACGTGCTGATGATGAAAATCGACGAGCAAAACA  >  NC_000913/2789041‑2789531

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