Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,622,354 Δ1 bp 10.0% coding (121/639 nt) purN → phosphoribosylglycinamide formyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,622,3540T.10.0% 161.1 / 14.8 60coding (121/639 nt)purNphosphoribosylglycinamide formyltransferase 1
Reads supporting (aligned to +/- strand):  ref base T (19/34);  new base . (2/4);  total (21/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.50e-02

CGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTGCTGGCTGGTTT  >  NC_000913/2622217‑2622496
                                                                                                                                          |                                                                                                                                              
cgctTGTTTATCGATTAAGGAATATTGTGGTGCTTATTTCCGCCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGc                                                                                                                                           <  2:26578/143‑1 (MQ=255)
   gTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAAATATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCtt                                                                                                                                        <  2:52818/143‑1 (MQ=255)
   gTGGTTATCGAATAATGAATATTGTGGTGCTCATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGAACGCGTCGGCCtt                                                                                                                                        <  2:12970/143‑1 (MQ=255)
     ggTTATCGAATAATGAATACTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGAACGCGTCGGCCttat                                                                                                                                      <  1:53244/143‑3 (MQ=255)
      gTTATGGAAAAATGAATATTGTGGGGTTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGaa                                                                                                                                     <  2:673235/143‑1 (MQ=255)
       ttATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGAACGCGTCGGCCttattg                                                                                                                                    <  1:311121/143‑5 (MQ=255)
        tATCGAATGACTAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAAcg                                                                                                                                   <  2:496079/143‑1 (MQ=255)
         aTCGAATAATGAATATTGGGGTGCTCATTGCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAAcgc                                                                                                                                  <  1:69159/143‑1 (MQ=255)
             aataatGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACGAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCGTGAACGagccc                                                                                                                              >  2:1012479/1‑143 (MQ=255)
             aataatAAAAATGGGGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGCCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGcgccc                                                                                                                              <  1:597879/143‑1 (MQ=255)
                   gAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGACAGTTTTCAGCCATGGGTCCGACGCGTTCGGCCTTGAACGCGCCCGCCAgg                                                                                                                        >  2:1018646/1‑143 (MQ=255)
                    aaTATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGc                                                                                                                       >  2:592908/1‑143 (MQ=255)
                      tatTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGAACGCGTCGGCCttatt                                                                                                                                     >  1:657929/1‑124 (MQ=255)
                      tatTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGAACGCGTCGGCCttatt                                                                                                                                     <  2:657929/127‑4 (MQ=255)
                      tatTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCgg                                                                                                                     <  1:791308/143‑1 (MQ=255)
                       atTGTGGTGTTTATTTCGGGCAATGGAAGTAATTTACAGGCAATTATTGATGCTTGTGCAACCAACAAAATTAAAGGCAACGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCggg                                                                                                                    <  1:768534/143‑1 (MQ=255)
                           gggtgCGTATGTCCGTCGACGGAAGGAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGTGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTAtt                                                                                                                <  1:921062/142‑1 (MQ=255)
                                cTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAAc                                                                                                           >  1:540583/1‑143 (MQ=255)
                                 ttATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAAAGCGCCCGCCAGGCGGGTGTTGCAACg                                                                                                          >  1:84791/1‑143 (MQ=255)
                                       ccGTCAACAAAAGTAGTCTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACg                                                                                                    <  2:532432/143‑1 (MQ=255)
                                       ccGGCACTGGAAGTAATTTCCAGGCAATTATTGACGCCTGCAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACg                                                                                                    <  2:91869/143‑1 (MQ=255)
                                         ggCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCGCGCCAGGCGGGTATTCCAACGCCTACGcc                                                                                                  >  2:761062/1‑142 (MQ=255)
                                         ggCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCt                                                                                                  >  1:861083/1‑143 (MQ=255)
                                            aaCGGAAGTAATGTACAGGCAATGACCGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCAt                                                                                               <  1:632935/143‑1 (MQ=255)
                                                    tAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTGGAACGCGCCCGCCAGGCGGGCATTGCAACGCATACGCACATCGCAAgcg                                                                                       >  2:153129/1‑143 (MQ=255)
                                                     aaTTTACAGGTAATTATTGACGCCTGTAAAACTAACAAAATTAAAGGATCCGTACTGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTTATCGCCAgcgc                                                                                      <  2:405346/143‑1 (MQ=255)
                                                        ttACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGtt                                                                                   <  1:594856/143‑1 (MQ=255)
                                                           cAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGa                                                                                <  1:862346/143‑1 (MQ=255)
                                                                  ttattGGCGCCGGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGt                                                                         <  2:864511/143‑1 (MQ=255)
                                                                    agtGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGa                                                                       <  2:140630/141‑1 (MQ=255)
                                                                     ttGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCCCATCGCCAGCGCGTTTGACAGTCGTGaa                                                                      >  1:599669/1‑143 (MQ=255)
                                                                     gtGACGCCTGCAAAACCAACAAAATTAAAGGCACCGTCGGGGCAGTTCTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGaa                                                                      <  1:508929/142‑1 (MQ=255)
                                                                         cGCCTGTAAAACCAACAAAATTAAAGGCACCGAGCGTTCATTGTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCt                                                                  <  1:713689/143‑1 (MQ=255)
                                                                          gCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTAGGCCTTGAACGAGCCCGCCAGGCGGGACTTGCAACGCATAGGACCATCGCCACCGCGTTTGCCAGCCGTGAAGccca                                                                 >  2:117548/1‑141 (MQ=255)
                                                                          gCCTGGAAAACGAACAAAATTAAAGGGACCGTACGGGCAGTGTGCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTa                                                                 <  2:120504/143‑1 (MQ=255)
                                                                           ccTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTGTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTAt                                                                <  1:323207/143‑1 (MQ=255)
                                                                                aaaaCCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGCTTGACAGTCGTGAAGCCTATGACCg                                                           <  1:1029573/143‑1 (MQ=255)
                                                                                  caCCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCggg                                                         <  1:735278/142‑1 (MQ=255)
                                                                                  aaCCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCggg                                                         >  1:670094/1‑143 (MQ=255)
                                                                                   aCCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGa                                                        >  1:942468/1‑143 (MQ=255)
                                                                                      aacaaAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGAAGtt                                                     >  2:779486/1‑143 (MQ=255)
                                                                                              tAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATg                                             <  1:576043/143‑1 (MQ=255)
                                                                                                 aGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGaaa                                          <  1:484505/143‑1 (MQ=255)
                                                                                                 aGGCACCGTACGGGCAGGTTTCAGCAATAAGGCCGACGCGTGCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGaaa                                          <  1:107403/143‑1 (MQ=255)
                                                                                                     aCCGTGAGGGCAGCTTTCAGCAATAGGGCCGACGCGTCCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGa                                      <  2:983003/143‑1 (MQ=255)
                                                                                                             ggCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGACGCCTATGACCGGGAGGTGATTCATGAAATCGACATGTACg                              >  2:582424/1‑143 (MQ=255)
                                                                                                                aGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCAc                           <  1:807119/143‑1 (MQ=255)
                                                                                                                     tCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGAt                      <  1:1018646/143‑1 (MQ=255)
                                                                                                                         gttcaaGGCCGAACGCGTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGACGTCTATGACCGGGAGTTGATTCATGAAATCGACGTGTACGCACCCGATGTgg                  >  2:351131/5‑143 (MQ=255)
                                                                                                                           ctagtgCCTACCCGCTCGTCCTTGAACGCACCGGCCCGGCGCTTATTGCAACGCATACGCTCATCGGCATCGCGTTTTGCATGCGTGAAGCCTATGACCGGGAGTCGATTCATGAAATCGACATGTACGCTCCCGATGTGGTc                <  2:488340/138‑1 (MQ=255)
                                                                                                                            tacggCCGACGCGTCCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCg               <  1:934981/140‑1 (MQ=255)
                                                                                                                            tAAGCCCGGCGCGTTTGGCCGTGATCGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCg               <  1:412921/143‑1 (MQ=255)
                                                                                                                                gCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGGGGGTATTGCAACGCATACGCTCATCGCCAGCGCGGTTGACAGTCGTGAAGCCTATGACCGGGAGTCGATTCATGACATCGGCCGGTACGCGCCAGCTGTGGTcatcgt           >  2:560035/1‑138 (MQ=255)
                                                                                                                                 ccGACGCGTTCGGCCTTGAACGCGCCCGCGAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgctg          <  2:455539/143‑1 (MQ=255)
                                                                                                                                  cGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGCAGCCTATGACCGGGAGTTGACTCATGAAATCGACATGTACGCACCCGATGTGGTCCTgttgg         >  1:930402/1‑143 (MQ=255)
                                                                                                                                  cGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCt                                                                  >  2:360011/1‑86 (MQ=255)
                                                                                                                                  cGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCt                                                                  <  1:360011/86‑1 (MQ=255)
                                                                                                                                  cGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgctgg         <  2:683890/143‑1 (MQ=255)
                                                                                                                                   gACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCACCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATCTACGCACCCGATGTGGTCGTgctggc        >  2:570689/1‑143 (MQ=255)
                                                                                                                                       gcgTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTGGGACAGTCGTGAAGCCTATGACCGGGAGTTGCTTCCTGAAATCGGCATGTACGCACACGATACGGGCGTGcgcgctgg     >  2:782672/1‑143 (MQ=255)
                                                                                                                                          ttCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGCCATGTACGCACCCGATGTGGTCGCGCTGGCTGGctt  >  1:1026816/1‑143 (MQ=255)
                                                                                                                                          |                                                                                                                                              
CGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTGCTGGCTGGTTT  >  NC_000913/2622217‑2622496

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