Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 961,757 G→A 6.0% R186H (CGT→CAT)  cmk → cytidylate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913961,7570GA 6.0% 213.3 / 3.0 68R186H (CGT→CAT) cmkcytidylate kinase
Reads supporting (aligned to +/- strand):  ref base G (28/35);  new base A (2/2);  total (31/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

CCCTGATGCACCAGTGAAAATTTTCCTTGACGCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCAG  >  NC_000913/961614‑961900
                                                                                                                                               |                                                                                                                                               
cccTGATGCACCAGTGAAAATTTTCCTTGACGCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGAGCGTAAc                                                                                                                                             <  6:332971/148‑1 (MQ=255)
cccTGATGCACCAGTGAAAATTTTCATTGACGCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTTGCCGAGAACAAAGAACGCGACTACCGAGAACCTAAcc                                                                                                                                            >  6:335043/1‑149 (MQ=255)
 ccTGATGCACCAGTGAAAATTTTCCTTGACGCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCATAAcc                                                                                                                                            >  3:1843/1‑148 (MQ=255)
    gATGCACCAGTGAAAATTTTCCTTGACGCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGc                                                                                                                                        >  6:280962/1‑149 (MQ=255)
     aTGCACCAGTGAAAATTTTCCTTGACGCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAg                                                                                                                                         <  1:36447/147‑1 (MQ=255)
                       tCCTTGACGCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTcc                                                                                                                      <  4:102730/148‑1 (MQ=255)
                       tCCTTGACGCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCg                                                                                                                     >  2:256140/1‑149 (MQ=255)
                        ccTTGACGCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCg                                                                                                                     <  2:8491/148‑1 (MQ=255)
                         cTTGACGCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCgg                                                                                                                    <  5:280962/148‑1 (MQ=255)
                            gACGCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCa                                                                                                                  <  2:107429/147‑1 (MQ=255)
                            gACGCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGc                                                                                                                >  5:117922/1‑149 (MQ=255)
                               gCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGa                                                                                                             >  8:211023/1‑149 (MQ=255)
                               gCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGa                                                                                                             <  1:256140/149‑1 (MQ=255)
                               gCCTCCTAGGAAGAACGTGCGCAACGCCGCATGCTACAGTTGCAGTAGAAGGGATTTAGTGTTAAAATTGAGCGCCTATTGGCCGAGATCAAAGAACGCGACGCCCGCGATCGTAACCTAGAGGTAGCGCCACTGGTTCCGGCAGCCGa                                                                                                             >  8:145506/1‑149 (MQ=255)
                                  tcctcGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGc                                                                                                          >  4:325543/1‑149 (MQ=255)
                                    ctcGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGACCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCtt                                                                                                        <  3:186012/149‑1 (MQ=255)
                                     tcGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCtt                                                                                                        >  8:322651/1‑148 (MQ=255)
                                       ggAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTa                                                                                                      <  6:230266/148‑1 (MQ=255)
                                       ggAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTa                                                                                                      <  3:313156/148‑1 (MQ=255)
                                       ggAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTa                                                                                                      <  8:3923/148‑1 (MQ=255)
                                        gaagaaCGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAgt                                                                                                    >  7:30303/1‑149 (MQ=255)
                                           gaaCGTGCTAATCGCCGTATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATAAAATAACGCGACGACCGGGTTGGTAACAGAGCGGTAGCGCCACTGGTTCCGGCAGCCGTTGCTTTAGTGtt                                                                                                 <  7:145506/149‑1 (MQ=255)
                                                 gcgcATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAgcgc                                                                                                                                >  6:80982/1‑112 (MQ=255)
                                                 gcgcATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAgcgc                                                                                                                                <  5:80982/112‑1 (MQ=255)
                                                  cgcATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTcc                                                                                          <  8:301610/149‑1 (MQ=255)
                                                  cgcATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCATAACCGAGCGGTAGCGCCACTGGTTCCGGCATCCGATGCTTTAGTGTTGGATTc                                                                                           <  4:1843/148‑1 (MQ=255)
                                                   gcATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGACCGCGACGACCGCGATCGTACCCGAGCGGTAGCGCCACTTGTTCCGGCAGCCGATGCTTTAGTGTTGGATTcca                                                                                         >  2:290656/1‑149 (MQ=255)
                                                   gcATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTcca                                                                                         <  2:156722/149‑1 (MQ=255)
                                                     aTCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTccac                                                                                        >  8:115115/1‑148 (MQ=255)
                                                      tCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGCGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTccacca                                                                                      <  2:104890/149‑1 (MQ=255)
                                                          cgcATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGc                                                                                                                >  7:62931/1‑119 (MQ=255)
                                                          cgcATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGc                                                                                                                <  8:62931/119‑1 (MQ=255)
                                                            cATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTaa                                                                                <  5:53740/149‑1 (MQ=255)
                                                                 tACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGa                                                                                              >  6:147656/1‑130 (MQ=255)
                                                                 tACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGa                                                                                              <  5:147656/130‑1 (MQ=255)
                                                                 tACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCAtt                                                                           >  2:306493/1‑149 (MQ=255)
                                                                 tACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCATTTTGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGAGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTAAAGCAtt                                                                           >  1:317697/1‑149 (MQ=255)
                                                                 tACAGTTGCAGGAGAAGGGCTTTAGTGTTAAATTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGCCCGCGATCGTAACCGAGCGGTAGCGCCACGGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTGAAGCAtt                                                                           >  4:218604/1‑149 (MQ=255)
                                                                    aGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAg                                                                        <  8:37748/149‑1 (MQ=255)
                                                                    aGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAg                                                                        >  7:200000/1‑149 (MQ=255)
                                                                        gCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCaa                                                                     <  7:322651/148‑1 (MQ=255)
                                                                                        aGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAgcgc                                                     >  8:172484/1‑148 (MQ=255)
                                                                                        aGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAgcgc                                                     >  8:167664/1‑148 (MQ=255)
                                                                                         gtgtTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTa                                                   <  7:115115/149‑1 (MQ=255)
                                                                                          tgtTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTa                                                   >  6:221213/1‑148 (MQ=255)
                                                                                            ttAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCTGCAGCCGATGCTTTAGTGTTGGATTCAACCACCTAAAGCATTTAGCAAGTGATTGAAAAAGCTCTACa                                                 >  8:160078/1‑148 (MQ=255)
                                                                                            ttAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTACaa                                                >  3:187715/1‑149 (MQ=255)
                                                                                            ttAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTACaa                                                <  7:29222/149‑1 (MQ=255)
                                                                                              aaCTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTACAATa                                              >  7:156843/1‑149 (MQ=255)
                                                                                                   tGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTACAATAcgc                                           <  8:200000/147‑1 (MQ=255)
                                                                                                    gAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTACAATAcgcgc                                         <  6:117922/148‑1 (MQ=255)
                                                                                                              tGGCCGAGATCAAAGACCGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCGTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTACAATACGCGCGCCAGAAAtt                               <  7:332150/148‑1 (MQ=255)
                                                                                                               ggCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTg                              <  1:306493/148‑1 (MQ=255)
                                                                                                                gCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTg                                                                          >  3:19966/1‑103 (MQ=255)
                                                                                                                gCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTg                                                                          <  4:19966/103‑1 (MQ=255)
                                                                                                                    agaTCAAAGAACGCGACGACCGCGATCGTAACCGATCGGTAGCGCCACTGGTTCCGGCAGCCGATGCGTTAGTGTTGGATTCCACCACCTTAAGCAATGAGCAAGTGATTGAAAAAGCGCTACAATAAGCGCGCCAGAAATTGGCtctc                        <  1:220048/149‑1 (MQ=255)
                                                                                                                      aTCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCg                       >  6:100030/1‑148 (MQ=255)
                                                                                                                          aaGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCAt                    >  4:202261/1‑147 (MQ=255)
                                                                                                                          aaGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATaa                  <  3:325543/149‑1 (MQ=255)
                                                                                                                                gcgACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTACAATACgcgcgc                                       >  4:152454/1‑122 (MQ=255)
                                                                                                                                gcgACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTACAATACgcgcgc                                       <  3:152454/122‑1 (MQ=255)
                                                                                                                                gcgACGACCGCGATCATAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGAcc            >  8:81638/1‑149 (MQ=255)
                                                                                                                                    cgacCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAAt        <  7:160078/149‑1 (MQ=255)
                                                                                                                                     gacCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAAtt       >  5:208452/1‑149 (MQ=255)
                                                                                                                                       ccGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAATTTg     <  7:172484/149‑1 (MQ=255)
                                                                                                                                       ccGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAATTTg     <  7:167664/149‑1 (MQ=255)
                                                                                                                                         gcgATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCa   <  7:211023/149‑1 (MQ=255)
                                                                                                                                          cgATCATAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCAg  <  7:81638/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                               
CCCTGATGCACCAGTGAAAATTTTCCTTGACGCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCAG  >  NC_000913/961614‑961900

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