Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,622,348:1 +A 10.0% coding (115/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,3481.A10.0% 203.7 / 20.4 60coding (115/639 nt)purNphosphoribosylglycinamide formyltransferase 1
Reads supporting (aligned to +/- strand):  ref base . (18/36);  new base A (2/4);  total (20/40)
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.18e-01

ATTCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTGCTGGC  >  NC_000913/2622206‑2622490
                                                                                                                                               |                                                                                                                                              
aTTCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCg                                                                                                                                                 >  2:713689/1‑143 (MQ=255)
  tCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGAc                                                                                                                                              <  2:970280/143‑1 (MQ=255)
     gAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGAcgca                                                                                                                                           >  2:252412/1‑142 (MQ=255)
           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                                                            <  1:360011/86‑1 (MQ=255)
                                                                                                                                             cGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCt                                                            >  2:360011/1‑86 (MQ=255)
                                                                                                                                             cGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgctgg   <  2:683890/143‑1 (MQ=255)
                                                                                                                                              gACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCACCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATCTACGCACCCGATGTGGTCGTgctggc  >  2:570689/1‑143 (MQ=255)
                                                                                                                                               |                                                                                                                                              
ATTCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTGCTGGC  >  NC_000913/2622206‑2622490

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: