Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,723,065 G→A A113A (GCG→GCA pssA → phosphatidylserine synthase, CDP‑diacylglycerol‑serine O‑phosphatidyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,723,0650GA100.0% 223.7 / NA 69A113A (GCG→GCApssAphosphatidylserine synthase, CDP‑diacylglycerol‑serine O‑phosphatidyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (33/36);  total (33/36)

GTATGAGGCTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCGGTGCCAGCCTGAACGATGT  >  NC_000913/2722939‑2723199
                                                                                                                              |                                                                                                                                      
gTATGAGGCTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCg                                                                                                                            <  2:678114/139‑1 (MQ=255)
         tAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGAt                                                                                                                   >  1:1123862/1‑139 (MQ=255)
         tAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGAt                                                                                                                   <  1:1262516/139‑1 (MQ=255)
             aGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGAATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTc                                                                                                               <  1:30939/139‑1 (MQ=255)
                 aGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGt                                                                                                           >  2:527786/1‑139 (MQ=255)
                     tCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGATGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTAt                                                                                                       >  2:299266/1‑139 (MQ=255)
                     tCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTAt                                                                                                       >  1:687953/1‑139 (MQ=255)
                          aaCTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGtt                                                                                                 <  2:1040462/139‑1 (MQ=255)
                               gATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCcaa                                                                                             >  2:1340273/1‑139 (MQ=255)
                                 tgtgCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCcaatc                                                                                           <  1:527786/139‑1 (MQ=255)
                                   tgCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCcaatcaa                                                                                         >  2:316244/1‑139 (MQ=255)
                                    gCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCcaatcaat                                                                                        >  2:1424307/1‑139 (MQ=255)
                                    gCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCcaatcaat                                                                                        >  1:702311/1‑139 (MQ=255)
                                         tGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATgg                                                                                                     <  1:1111831/121‑1 (MQ=255)
                                         tGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATgg                                                                                                     >  2:1111831/1‑121 (MQ=255)
                                           cTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCg                                                                                                   >  2:615311/1‑121 (MQ=255)
                                           cTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCg                                                                                                   <  1:615311/121‑1 (MQ=255)
                                            tGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGa                                                                                <  1:730614/139‑1 (MQ=255)
                                             ggTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGa                                                                                <  2:1321882/138‑1 (MQ=255)
                                             ggTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGa                                                                                >  1:1321882/1‑138 (MQ=255)
                                                  aCTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCt                                                                          <  2:1215117/139‑1 (MQ=255)
                                                  aCTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCt                                                                          <  1:531920/139‑1 (MQ=255)
                                                            tGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAg                                                                                                                                     <  2:1133800/70‑1 (MQ=255)
                                                            tGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAg                                                                                                                                     >  1:1133800/1‑70 (MQ=255)
                                                              cacaACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCa                                                              <  1:305123/139‑1 (MQ=255)
                                                                   cGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTa                                                         <  2:702311/139‑1 (MQ=255)
                                                                       gACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTg                                                                >  1:1295801/1‑128 (MQ=255)
                                                                       gACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTg                                                                <  2:1295801/128‑1 (MQ=255)
                                                                        aCGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGtt                                                                                                                <  2:1238487/79‑1 (MQ=255)
                                                                        aCGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGtt                                                                                                                >  1:1238487/1‑79 (MQ=255)
                                                                         tgCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCt                                                   <  1:299266/138‑1 (MQ=255)
                                                                             ttGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTat                                               >  1:771241/1‑139 (MQ=255)
                                                                             ttGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTat                                               >  2:941724/1‑139 (MQ=255)
                                                                             ttGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTat                                               >  2:62427/1‑139 (MQ=255)
                                                                                 cgcTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTg                                                                        >  2:201596/1‑110 (MQ=255)
                                                                                 cgcTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTg                                                                        <  1:201596/110‑1 (MQ=255)
                                                                                   cTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATcga                                         <  1:8006/139‑1 (MQ=255)
                                                                                          aTCTAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGc                                  <  1:1424307/139‑1 (MQ=255)
                                                                                             tAACACTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTa                               <  2:771241/139‑1 (MQ=255)
                                                                                                cacTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTa                                                         <  1:111273/110‑1 (MQ=255)
                                                                                                cacTAACGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTa                                                         >  2:111273/1‑110 (MQ=255)
                                                                                                     aCGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAg                       >  2:903829/1‑139 (MQ=255)
                                                                                                     aCGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAg                       <  1:903829/139‑1 (MQ=255)
                                                                                                     aCGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAg                       >  1:1117884/1‑139 (MQ=255)
                                                                                                     aCGCTGACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAg                       <  1:316244/139‑1 (MQ=255)
                                                                                                          gACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATa                                                                                       >  2:15819/1‑70 (MQ=255)
                                                                                                          gACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATa                                                                                       <  1:15819/70‑1 (MQ=255)
                                                                                                          gACTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCGGTg                  <  1:941724/139‑1 (MQ=255)
                                                                                                            cTGGTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCGGTGcc                <  1:1340273/139‑1 (MQ=255)
                                                                                                               gTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTa                               <  2:148666/121‑1 (MQ=255)
                                                                                                               gTACTGCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTa                               >  1:148666/1‑121 (MQ=255)
                                                                                                                    gCCGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCGGTGCCAGCCTGaa        <  1:1419180/139‑1 (MQ=255)
                                                                                                                     ccGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCcaatcaa                                                                                         <  1:804374/57‑1 (MQ=255)
                                                                                                                     ccGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCcaatcaa                                                                                         >  2:804374/1‑57 (MQ=255)
                                                                                                                     ccGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCGGTGCCAGCCTGAAc       >  2:934605/1‑139 (MQ=255)
                                                                                                                     ccGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCGGTGCCAGCCTGAAc       >  1:126391/1‑139 (MQ=255)
                                                                                                                     ccGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCGGTGCCAGCCTGAAc       >  1:439168/1‑139 (MQ=255)
                                                                                                                     ccGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCGGTGCCAGCCTGAAc       >  2:396346/1‑139 (MQ=255)
                                                                                                                      cGCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCGGTGCCAGCCTGAACg      <  2:126391/139‑1 (MQ=255)
                                                                                                                       gCATGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCGGTGCCAGCCTGAACGa     >  1:234049/1‑139 (MQ=255)
                                                                                                                          tGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCcaat                                                                                            >  1:1106954/1‑49 (MQ=255)
                                                                                                                          tGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCcaat                                                                                            <  2:1106954/49‑1 (MQ=255)
                                                                                                                          tGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCt                                                                          <  1:772414/67‑1 (MQ=255)
                                                                                                                          tGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCt                                                                          >  2:772414/1‑67 (MQ=255)
                                                                                                                          tGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATa                                    >  2:773628/1‑105 (MQ=255)
                                                                                                                          tGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATa                                    <  1:773628/105‑1 (MQ=255)
                                                                                                                          tGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAg                       <  1:961147/118‑1 (MQ=255)
                                                                                                                          tGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAg                       >  2:961147/1‑118 (MQ=255)
                                                                                                                          tGGCACAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCGGTGCCAGCCTGAACGATGt  <  2:234049/139‑1 (MQ=255)
                                                                                                                              |                                                                                                                                      
GTATGAGGCTAAAAGGCAGCGTCCGGAACTGGATGTGCGGGTGCTGGTCGACTGGCATCGTGCACAACGTGGACGCATTGGCGCTGCGGCATCTAACACTAACGCTGACTGGTACTGCCGCATGGCGCAGGAAAATCCGGGCGTAGATGTTCCGGTTTATGGCGTTCCAATCAATACTCGTGAAGCCCTTGGTGTTCTGCACTTTAAAGGCTTTATCATCGACGATAGCGTACTTTATAGCGGTGCCAGCCTGAACGATGT  >  NC_000913/2722939‑2723199

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