Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,066,699 C→T A372V (GCC→GTC)  agp → glucose‑1‑phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,066,6990CT100.0% 111.9 / NA 40A372V (GCC→GTC) agpglucose‑1‑phosphatase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/19);  total (21/19)

ATATCGACAAAGCTCTGGTCACCGATCGCACCAGCGCACCGAAAATTACAGTGTTGGTTGGGCACGACTCCAACATTGCCTCTCTGTTAACGGCGCTGGATTTCAAACCGTATCAGTTGCATGACCAGAACGAACGCACGCCGATTGGCGGCAAAATCGTTTTCCAGCGTTGGCATGACAGCAAAGCCAATCGCGATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGCCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCGTGTTGAATGAAGCGGTGAAATAACA‑GAAAACTCCCCCGCGAGAAGCGGGGGAGTCGCTGGTTAAACGTTTTTACGTTCGATGGTCTGTTCGCCCCAAAAAAGCGAATCTTTATCGGTCTTAGCAAAGGCTTTGACT  >  NC_000913/1066453‑1066939
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 
atatCGACAAAGCTCTGGTCCCCGATCGCACCAGCGCACCGAAAATTACAGTGTTGCTTGGGACCGACTCCAACATTGCCTCTCTGTTAACGGCGCTGGATTTCAAACAGTATCAGTTGCATGACAAGAACGAACGCACGCCGATTGGCGGCAAAATAGTTTTCCCGCGTTGGCATGACAGGAACTCCAATCGCGATTTGATGCACATTTAATATGAGTATAAGAGTGAGGACCAGTTAATTAATGTCGAt                                                                                                                                                                                                                                               >  2:569658/1‑251 (MQ=255)
  atCGACAAAGCTCTGGTCACCGATCGCACCAGCGCACCGAAAATTACAGTGTTGGTTGGGCACGACTCCAACATTGCCTCTCTGTTAACGGCGCTGGATTTCAAACCGTATCAGTTGCATGACCAGAACGAACGCACGCCGATTGGCGGCAAAATCGTTTTCCAGCGTTGGCATGACAGCAAAGCCAATCGCGATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATg                                                                                                                                                                                                                                              <  1:578992/250‑1 (MQ=255)
       cAAATCTCTGGTCACCGATCGCACCAGCGCACCGAAAATTACAGTGTTAGTTGGGCACGACTCCAACATTGCCTCTCTGTTAACGGCGCTGGCTTTCAAACCGTATCAGTTGCATGAACAGAACGACAGCACGCCGATTGGCGGCAAAATCGATTTCCAGCGGTGGCATGACCGCAAAGCCCATCGCGATTTGATGAAAATTGAATTTGTGTATCCGAGTGCGGAACCGTTAAGTAAAGTCGCTGCGTTa                                                                                                                                                                                                                                         >  2:426411/1‑250 (MQ=255)
                      ccaccgccccAGCGCACCTAAAATTACTGTGTTGGTTGGGCACCCCCCCACCATTGCCTCTCTTGTAACGGCGCTTGATTTCAAACCGTATCAGTTGCATGACCAGAACGAACGCACGCCGATTGGCGGCAAAATCGTTTTCCAGCGTTGGCATGACAGCAAAGCCAATCTCGATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGTTGCGTTAACCCTGCATGCACCTg                                                                                                                                                                                                                         <  2:542139/243‑1 (MQ=255)
                       gctCGCACCAGCGCACCGAAAATATCGGTGTTGGTTGGGCACGCCTACAACATTGCCTCTCTGTTAACGGCGCTGGATTTCAAACCGTATCAGTTGCATGACCAGAACGAACTCACGCAGCTTGGCGGCAAAATCGTTTTCCAGCGTTGGCCTCACAGCAAAGCCAATCGCGGTTTGATGAAAATTGCATATGTGTATCAGAGTGCGGAACAGTTACGTACTTTCGATGAGTTAACCATGCAGGCACCTgg                                                                                                                                                                                                                        >  2:132677/3‑250 (MQ=255)
                           gCACCAGCGCACCGAAAATTACAGTGTTGGTTGGGCACGACTCCAACATTGCCTCTCTGTTAACGGCGCTGGATTTCAAACCGTATCAGTTGCATGACCAGAACGAACGCACGCCGATTGGCGGCAAAATCGTTTTCCAGCGTTGGCATGACAGCAAAGCCAATCGCGATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAg                                                                                                                                                                                                                    <  1:299581/251‑1 (MQ=255)
                              ccAGCGCCCCGAAAATTACAGCGTTGGTTTGGCAAGACTCCAACATTGCCTCGCTGTTAACGGCGCTGGATTTCAAACCGTACCAGTTGCATGACCAGAACGAACGCACGCCGATTGGCGGCAAAATCGTTTTCCAGCGTTGGCATGACAGCAAAGCCAATCGCGATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGCgt                                                                                                                                                                                                                 <  2:530293/251‑1 (MQ=255)
                                     aCCGAAAATTACAGTGTTGGTTGGGCACGACTCCAACATTGCCTCTCTGTTAACGGCGCTGGATTTCAAACCGTATCAGTTGCATGACCAGAACGAACGCACGCCGATTGGCGGCCAAATCGTTTTCCAGCGTTGGCATGACAGCAAAGCCAATCGCGATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGc                                                                                                                                                                                                          >  1:393294/1‑251 (MQ=255)
                                                    tttggttggGCACGACTCCAACATTGCCTCCCTGTAAACGGCGCTGGATTTCACACCGTATCCGTTGCATGACCAGAACGAACGCACGCCGATTGGCGGCAAAATCGTTTTCCAGCGATGGCATGACAGCTAAGACGATCGCGCGTTGGTGGACATTGACGATGTGTCTCAGAGAGCGGAACGGGTAAATAATGTTGGTGCATTAAATCTGCAGGCAACTGCGCAGAgtggggcacggaaatttagcggtt                                                                                                                                                                                           >  2:438502/2‑230 (MQ=255)
                                                           gggCACGACTCCAACATTGCCTCGCTTTTAACGGCGCTGGATTTCAAACCGTATCAGTTGCATGACCAGAACGAACGCACGCCGATTGGCGGCAAAATCGTTTTCCAGCGTTGGCATGACAGCAAAGCCAATCGCGATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGAt                                                                                                                                                                                    <  2:104089/251‑1 (MQ=255)
                                                                  aCTCCAACATTGCCTCTCTGTTAACGGCGCTGGATTTCAAACCGTATCAGTTGCATGACCAGAACGAACGCACGCCGATTGGCGGCAAAATCGTTTTCCAGCGTTGGCATGACAGCAAAGCCAATCGCGATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGACGcg                                                                                                                                                                             >  2:367857/1‑250 (MQ=255)
                                                                     ccAACATTGCCTCTCTGTTAACGGCGCTGGATTTCAAACCGTATCAGTTGCATGACCAGAACGAACGCACGCCGATTGGCGGCAAAATCGTTTTCCAGCGTTGGCATGACAGCAAAGCCAATCGCGATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGACGCTGAt                                                                                                                                                                          >  1:563692/1‑251 (MQ=255)
                                                                          aTTGCCTCTCTGTTAACGGCGCTGGATTTCAAACCGTATCAGTTGCATGACCAGAACGAACGCACGCCGATTGGCGGCAAAATCGTTTTCCAGCGTTGGCATGACAGCAAAGCCAATCGCGATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGACGCTGATGGttt                                                                                                                                                                     >  1:85014/1‑251 (MQ=255)
                                                                                 ctctGTTAACGGCGCTGGATTTCAAACCGTATCAGTTGCATGACCAGAACGAACGCACGCCGATTGGCGGCAAAATCGTTTTCCAGCGTTGGCATGACAGCAAAGCCAATCGCGATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGACGCTGATGGTTTCTGcc                                                                                                                                                                <  1:271888/249‑1 (MQ=255)
                                                                                       tAACGGCGCTTGGTTGCAAACCTTATCAGTTGAATGACCAGGACGAACGCACGCCGATTGGCGGCAAAATCGTTTTCCAGCGGTGGCATGACAGCAAAGCCAATCGCGATTTGATGAAAAATGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGTCACCTGCGCAGCGTGTGACGCTGGAATTAAGCTGTTGCCCGATAGACGCTGATGGTTGCTGCCCGATGGa                                                                                                                                                         <  2:619719/250‑1 (MQ=255)
                                                                                          cGGCGCTGGATTTCAAACCGTATCAGTTGCATGACCAGAACGAACGCACGCCGATTGGCGGCAAAATCGTTTTCCAGCGTTGGCATGACAGCAAAGCCAATCGCGATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAg                                                                                                                                                     <  1:367857/251‑1 (MQ=255)
                                                                                                      tCAAACCGTATCAGTTGCATGACCAGAACGAACGCACGCCGATTGGCGGCAAAATCGTTTTCCAGCGTTGGCATGACAGCAAAGCCAATCGCGATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCgt                                                                                                                                          <  1:260632/250‑1 (MQ=255)
                                                                                                                      gCATGACCAGATCGAACGCACGCCGATTGGCGGCAAAATCGGTTGCCAGCGGTGGCATGACAGCAAAGCCAATCGCGATTTGATGAAAATTGAATATGTGTATCATAGTGCGGAACAGTTACGTAATGTCGATGCGGTAACCCTGCAGTCACCTGCTCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGACGCTGATGGTTTCTGCTCGATGGATAAGTTTGAGAGCGTGTTGAATGAAGCGGTg                                                                                                                          <  2:393294/250‑1 (MQ=255)
                                                                                                                          gACCAGAACGAACGCACGCCGATTGGCGGCAAAATCGTTTTCCAGCGTTGGCATGACAGCAAAGCCAATCGCGATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCGTGTTGAATGAAGCGGTGAAATa                                                                                                                     <  1:19193/251‑1 (MQ=255)
                                                                                                                               gaacgaacGCACGCCGATTGTCGGCAAAATCGTTTACCAGCGTTGGCATGACAGCAAAGCCAATCGCGATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGTCGCTGGAATTAAGCGGTGGCCCGATAGACGCTGATGGTGTCTGCCCGATGGATAAGTTTTATAGCGTGTTGAATGAAGCGGTGAAATAACA‑Ga                                                                                                               <  1:438502/251‑1 (MQ=255)
                                                                                                                                          cgcCGATTGGCGGCAAAATCGGTTTCCAGCGTTGGCATGACAGCAAAGCCAATCGCGATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGAGGCGTTATCCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCTGTTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCGTGTTGAATGAATCGGTGAAATAACA‑GAAAACTCCCCcg                                                                                                    <  2:513105/251‑1 (MQ=255)
                                                                                                                                            ccGATTGGCGGCAAAATCGTTTTCCAGCGTTGGCATGACAGCAAAGCCAATCGCGATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCGTGTTGAATGAAGCGGTGAAATAACA‑GAAAACTCCCCcgcg                                                                                                  <  1:208458/251‑1 (MQ=255)
                                                                                                                                                   gcggcAAAATCGTTTTCCAGCGTTGGCATGACAGCAAAGCCAATCGCGATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCG‑TAAACGGTGGTGGTTTCTGCCCCAAGGATTAGTTTGATAACGTTTTGAATTAAGCGGTGAAATAACAAGAAACCTccccgccagaagggg                                                                                           >  1:86130/1‑240 (MQ=255)
                                                                                                                                                      gaaaaaTCGTTGTCCAGCGGCGGCATGACAGCAAAGCCATGCGCGATTTGATGAAAATTGAATATGTGTACCAGAGGGCTGAACAGATACGTAATGTCGATGCGGTATCCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCTGATAGACGCTGAGGGTTTCGGCCCGAGGGACAAGTTTGATAGAGTGTTGAATTAAGCGGTGAAATAAGA‑GAGAACCCCACCGCGAGAAGCgggg                                                                                        <  2:297928/249‑1 (MQ=255)
                                                                                                                                                        aaaaTCGTTTTCCAGCGTTGGCATGACAGCAAAGCCAATCGCGATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGACGCTGATGGATTCTGCCCGATGGATAAGTTTGATAGCGTGTTGAATGAAGCGGTGAAATAACA‑GAAACCTCCCCCGCGAGAAGCggggg                                                                                       >  2:491265/1‑250 (MQ=255)
                                                                                                                                                          aaTCGTTTTCCAGCGTGGGCATGACAGCAAAGCCAATCGCGATTTGATGAAAATTGCATATGTGTATCAGAGTGTGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGATGCCCGATAGCCGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCGTGTTGAATGAAGCGGTGACATAACC‑GCAAACTCCCCCGCGAGAAGCGGGGGAGt                                                                                    >  2:383291/1‑251 (MQ=255)
                                                                                                                                                              gTTTTCCAGCGTTGGCATGACAGCAAAGCCAATCGCGATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGCCGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCGTGTTGAATGAAGCGGTGAAATAACA‑GAAAACTCCCCCGCGAGAAGCGGGGGAGTCgtg                                                                                >  2:133207/1‑249 (MQ=255)
                                                                                                                                                                   ccAGCGTTGGCATGAAAGCAAAGCCAATCGAGATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGAGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGCCGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCGTGTTGAATGAAGCGGTGAAATAACA‑GAAAACTCCCCCGCGAGAAGCGGGGGAGTCGCTGGtt                                                                            >  1:258727/1‑250 (MQ=255)
                                                                                                                                                                    cAGCGTTGGCATGACAGCAAAGCCAATCGCGATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCGTGTTGAATGAAGCGGTGAAATAACA‑GAAAACTCCCCCGCGAGAAGCGGGGGAGTCGCTGGTTa                                                                           >  1:373578/1‑250 (MQ=255)
                                                                                                                                                                          tGGCATGACAGCAAAGCCAATCGCGATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCGTGTTGAATGAAGCGGTGAAATAACA‑GAAAACTCCCCCGCGAGAAGCGGGGGAGTCGCTGGTTAAACGttt                                                                    >  1:168011/1‑251 (MQ=255)
                                                                                                                                                                           ggCATGACAGCAAAGCCAATCGCGATCTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCTCCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCGTGTTGAATGAAGCGGTGAAATAACA‑GAAAACTCCCCCGCGAGAAGCGGGGGAGGCGCTGGTTAAACGttt                                                                    <  2:458892/250‑1 (MQ=255)
                                                                                                                                                                                      aaaGCCAATCGCGATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCGTGTTGAATGAAGCGGTGAAATAAAA‑GAAAACTCCCCCGCGAGAAGCGGGGGAGTCGCTGGTTAAACGTTTTTACGTTCGATg                                                        >  2:551982/1‑251 (MQ=255)
                                                                                                                                                                                           cAATCGCGATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCGTGTTGAATGAAGCGGTGAAATAACA‑GAAAACTCCCCCGCGAGAAGCGGGGGAGTCGCTGGTTAAACGTTTTTACGTTCGATGGTCTg                                                   >  1:123474/1‑251 (MQ=255)
                                                                                                                                                                                                  gctttGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCGTGTTGAATGAAGCGGTGAAATAACA‑GAAAACTCCCCCGCGAGAAGCGGGGGAGTCGCTGGTTAAACGTTTTTACGTTCGATGGTCTGTTCGccc                                            >  2:67453/3‑251 (MQ=255)
                                                                                                                                                                                                  gATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCGTGTTGAATGAAGCGGTGAAATAACA‑GAAAACTCCCCCGCGAGAAGCGGGGGAGTCGCTGGTTAAACGTTTTTACGTTCGATGGTCTGTTCGccc                                            <  1:551982/251‑1 (MQ=255)
                                                                                                                                                                                                      tgatgaAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGCATTAAGCGGTTGCCCGCTAGCCGCTGATGGATTCTGCCCGATGGATAAGTTTGATAGCGTGTTGAATGCAGCGGTGAAATAACA‑GAAAACTCCCCCGTGCGAAGCGGGGGAGTCACTGGTTAAACGTTTTTACGTTCGATGGTCTGTTCGCCCCaa                                         >  1:648622/1‑250 (MQ=255)
                                                                                                                                                                                                        atgaAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCGTGTTGAATGAAGCGGTGAAATAACA‑GAAAACTCCCCCGCGAGAAGCGGGGGAGTCGCTGGTTAAACGTTTTTACGTTCGATGGTCTGTTCGCCCCaaaaa                                      <  1:382851/251‑1 (MQ=255)
                                                                                                                                                                                                                     atGTGGATCAGAGTGCGGAACAGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCGTGTTGAATGAAGCGGTGAAATAACA‑GAAAACTCCCCCGCGAGAAGCGGGGGAGTCGCTGGTTAAACGTTTTTACGTTCGATGGTCTGTTCGCCCCAAAAAAGCGAATCTTTAt                         <  2:532635/251‑1 (MQ=255)
                                                                                                                                                                                                                                          aGTTACGTAATGTCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCGTGTTGAATGAAGCGGTGAAATAACA‑GAAAACTCCCCCGCGAGAAGCGGGGGAGTCGCTGGTTAAACGTTTTTACGTTCGATGGTCTGTTCGCCCCAAAAAAGCGAATCTTTATCGGTCTTAGCAAAGGCTTTg     >  1:106614/1‑250 (MQ=255)
                                                                                                                                                                                                                                             tACGTAATGTCGATGCGATAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCGTGTTGAATGAAGCGGTGAAATAACA‑GAAAACTCCCCCGCGAGACGCGGGTGAATCGCTGGTTAACAGTGTTTACGTTCGATGGACTGTTAGCCCCAAAAAAGCGAATCTTTATCGTTCTTAGCAATGGCTTTTACt  >  2:649961/1‑250 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 
ATATCGACAAAGCTCTGGTCACCGATCGCACCAGCGCACCGAAAATTACAGTGTTGGTTGGGCACGACTCCAACATTGCCTCTCTGTTAACGGCGCTGGATTTCAAACCGTATCAGTTGCATGACCAGAACGAACGCACGCCGATTGGCGGCAAAATCGTTTTCCAGCGTTGGCATGACAGCAAAGCCAATCGCGATTTGATGAAAATTGAATATGTGTATCAGAGTGCGGAACAGTTACGTAATGCCGATGCGTTAACCCTGCAGGCACCTGCGCAGCGTGTGACGCTGGAATTAAGCGGTTGCCCGATAGACGCTGATGGTTTCTGCCCGATGGATAAGTTTGATAGCGTGTTGAATGAAGCGGTGAAATAACA‑GAAAACTCCCCCGCGAGAAGCGGGGGAGTCGCTGGTTAAACGTTTTTACGTTCGATGGTCTGTTCGCCCCAAAAAAGCGAATCTTTATCGGTCTTAGCAAAGGCTTTGACT  >  NC_000913/1066453‑1066939

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