Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 387,798 G→C 7.1% G313A (GGC→GCC)  glyA‑I → serine hydroxymethyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947387,7980GC 7.1% 260.4 / 9.0 84G313A (GGC→GCC) glyA‑Iserine hydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base G (44/34);  new base C (3/3);  total (47/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 = 8.89e-01

GCTGATGCACGTGATCGCTGCCAAGGCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGACCCGCAGTCGCCGTTCGTCACTTCGG  >  NC_002947/387652‑387941
                                                                                                                                                  |                                                                                                                                                 
gCTGATGCACGTGATCGCTGCCAAGGCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCa                                                                                                                                                >  6:140098/1‑149 (MQ=255)
     tGCACGTGATCGCTGTCAAGGCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGAc                                                                                                                                           <  7:79546/149‑1 (MQ=255)
        aCGTGATCGCTGCCAAGGCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAac                                                                                                                                        >  1:177688/1‑149 (MQ=255)
        aCGTGATCGCTGCCAAGGCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAac                                                                                                                                        <  6:233347/149‑1 (MQ=255)
        aCGTGATCGCTGCCAAGGCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAac                                                                                                                                        <  7:236163/149‑1 (MQ=255)
              tCGCTGCCAAGGCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTg                                                                                                                                 <  5:140098/149‑1 (MQ=255)
               cGCAGCCAAGGCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCt                                                                                                                                  <  4:310656/147‑1 (MQ=255)
               cGCAGCCAAGGCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCt                                                                                                                                  >  3:310656/1‑147 (MQ=255)
                gCTGCCAAGGCAGTGTTCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATAGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTTGTTTCCGGTGGCACCGACAACCACCTGt                                                                                                                                <  7:317073/148‑1 (MQ=255)
                 cTGCCAAGGCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTc                                                                                                                              >  7:66537/1‑149 (MQ=255)
                     cAAGGCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTgg                                                                                                                          >  6:63607/1‑149 (MQ=255)
                       aGGCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTc                                                                                                                        >  6:129590/1‑149 (MQ=255)
                        ggCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGt                                                                                                                         >  6:320472/1‑147 (MQ=255)
                                 cTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATcc                                                                                                              >  6:162230/1‑149 (MQ=255)
                                  ttCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCATCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACAGTGGCTAAGATTTGGTTTCCGGTGGCACAGAAAACAAAATGTTACTGGTCAGCCTGATACg                                                                                                             >  5:201815/1‑149 (MQ=255)
                                  ttCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCg                                                                                                             >  4:226331/1‑149 (MQ=255)
                                  ttCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCg                                                                                                             >  4:53547/1‑149 (MQ=255)
                                  ttCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCg                                                                                                             >  1:163042/1‑149 (MQ=255)
                                  ttCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCg                                                                                                             >  1:183728/1‑149 (MQ=255)
                                  ttCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCg                                                                                                             >  2:260916/1‑149 (MQ=255)
                                   tCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGt                                                                                                            >  5:296934/1‑149 (MQ=255)
                                    cAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGtt                                                                                                                               <  7:160946/129‑1 (MQ=255)
                                    cAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGtt                                                                                                                               >  8:160946/1‑129 (MQ=255)
                                       ggAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAgg                                                                                                        <  1:271713/149‑1 (MQ=255)
                                       ggAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTGTCGGTGCCACCGGAAaccacatc                                                                                                                                  >  3:354012/1‑121 (MQ=39)
                                       ggAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTGTCGGTGCCACCGGAAaccacatc                                                                                                                                  <  4:354012/124‑4 (MQ=39)
                                         aaGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGc                                                                                                      <  6:296934/149‑1 (MQ=255)
                                             gcTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGa                                                                                                  <  5:320472/149‑1 (MQ=255)
                                                 gAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTgg                                                                                              <  5:63607/149‑1 (MQ=255)
                                                     cTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCaa                                                                                           >  6:75825/1‑148 (MQ=255)
                                                      tGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACACCCACCTGTTCCTGGTCAGCCTGATCCGTCCGGGCCTGACTGGCa                                                                                            >  3:224283/1‑146 (MQ=255)
                                                      tGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCa                                                                                            >  8:345761/1‑146 (MQ=255)
                                                      tGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCAGGTCAGCCTGATCCGTCAGGGCCTGACTGGCa                                                                                            >  5:280763/1‑146 (MQ=255)
                                                          ttCAAGGCCTACCAGCAGCAAGTGATCGTAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAgacgc                                                                                     >  1:24846/1‑149 (MQ=255)
                                                          ttCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAgacgc                                                                                     <  2:24846/149‑1 (MQ=255)
                                                          ttCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAgacgc                                                                                     <  2:177688/149‑1 (MQ=255)
                                                          ttCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAgacgc                                                                                     >  3:166592/1‑149 (MQ=255)
                                                           tCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAATCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAgacgcc                                                                                    >  1:194321/1‑149 (MQ=255)
                                                           tCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCATGGCCTGACTGGCAAAgaagca                                                                                    >  7:235842/1‑148 (MQ=255)
                                                                 ccTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAgacgccgacgcc                                                                              >  2:86799/1‑149 (MQ=255)
                                                                     ccAGCAGCAAGTGATCGAAAACGCCCAGTCAATGGCCCAGGTGTTCATCTACAGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCTGTCAGGGCCTGACGGGCAAAGACGCCGACGCCGccc                                                                          <  4:224283/149‑1 (MQ=255)
                                                                     ccAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAATCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGccc                                                                          <  2:194321/149‑1 (MQ=255)
                                                                        gcagcaAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTgg                                                                       >  2:359441/1‑149 (MQ=255)
                                                                           gcaAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTc                                                                    <  3:226331/149‑1 (MQ=255)
                                                                              aGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGTCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGt                                                                  >  1:330429/1‑149 (MQ=255)
                                                                               gTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTgc                                                                <  1:359441/149‑1 (MQ=255)
                                                                               gTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTgc                                                                <  4:166592/149‑1 (MQ=255)
                                                                                      aaaaCGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTgg                                                                       <  8:273891/135‑1 (MQ=255)
                                                                                      aaaaCGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTgg                                                                       >  7:273891/1‑135 (MQ=255)
                                                                                      aaaaCGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATc                                                         <  6:10318/149‑1 (MQ=255)
                                                                                         aCGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCAcc                                                      >  7:297011/1‑149 (MQ=255)
                                                                                           gCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGt                                                    >  3:351631/1‑149 (MQ=255)
                                                                                               aGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTcaa                                                 >  7:84344/1‑148 (MQ=255)
                                                                                               aGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTcaa                                                 >  3:381012/1‑148 (MQ=255)
                                                                                               aGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGAGCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTcaac                                                <  5:129590/149‑1 (MQ=255)
                                                                                                 gCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCATGGTCAGCCTGATCAGTAAGGGCATGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATCTCACCGTcaacaa                                              >  6:289444/1‑149 (MQ=255)
                                                                                                    aTGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGaa                                           <  8:297011/149‑1 (MQ=255)
                                                                                                      ggCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACg                                         >  4:193751/1‑149 (MQ=255)
                                                                                                          cAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGATGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGt                                     <  8:66537/149‑1 (MQ=255)
                                                                                                                ttCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGaa                               <  5:75825/149‑1 (MQ=255)
                                                                                                                ttCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGaa                               >  7:107587/1‑149 (MQ=255)
                                                                                                                   aTCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACGCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAAc                              >  1:83106/1‑148 (MQ=255)
                                                                                                                   aTCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGa                            >  5:232533/1‑149 (MQ=255)
                                                                                                                        ccGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTcaacaa                                              <  6:72336/126‑1 (MQ=255)
                                                                                                                        ccGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTcaacaa                                              >  5:72336/1‑126 (MQ=255)
                                                                                                                         cGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGACCCGCa                      <  2:372151/149‑1 (MQ=255)
                                                                                                                           tGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGACCCGCAGt                    <  2:83106/149‑1 (MQ=255)
                                                                                                                           tGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGACCCGCAGt                    <  2:183728/149‑1 (MQ=255)
                                                                                                                                 cGATGTTGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCGTGACTGGCAAAGACGCCGACGCCGCCCTGGGTGGTGCGCATATCACCGTCAACAAGAAGGCCGTGCCGAACGTCCCGGAGTCGCCGt              <  2:334568/149‑1 (MQ=255)
                                                                                                                                  gATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCAcc                                                      >  7:274592/1‑108 (MQ=255)
                                                                                                                                  gATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCAcc                                                      <  8:274592/108‑1 (MQ=255)
                                                                                                                                  gATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGACCCGCa                      <  7:143328/140‑1 (MQ=255)
                                                                                                                                  gATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGACCCGCa                      >  8:143328/1‑140 (MQ=255)
                                                                                                                                  gATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGACCCGCAGTCGCCGtt             >  3:259855/1‑149 (MQ=255)
                                                                                                                                     gtgGTTGTCGGTGCCACCGGAAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGAc                           >  6:220771/1‑132 (MQ=255)
                                                                                                                                     gtgGTTGTCGGTGCCACCGGAAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGAc                           <  5:220771/132‑1 (MQ=255)
                                                                                                                                     gtgGTTGTCGGTGCCACCGGAAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGAc                           <  5:226423/132‑1 (MQ=255)
                                                                                                                                     gtgGTTGTCGGTGCCACCGGAAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGAc                           >  6:226423/1‑132 (MQ=255)
                                                                                                                                       ggTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGAcc                          <  6:382708/131‑1 (MQ=255)
                                                                                                                                       ggTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGAcc                          >  5:382708/1‑131 (MQ=255)
                                                                                                                                       ggTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGACCCGCAGTCGCCGTTCGTCa        >  7:70959/1‑149 (MQ=255)
                                                                                                                                        gTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGACCCGCAGTCGCCGTTCGTCAc       >  6:375231/1‑149 (MQ=255)
                                                                                                                                            ccGGTGGCACCGACAAGCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGACCCGCAGTCGCCGTTCGTCACTTc    <  1:86799/148‑1 (MQ=255)
                                                                                                                                             cGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGACCCGCAGTCGCCGTTCGTCACTTCgg  <  3:193751/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GCTGATGCACGTGATCGCTGCCAAGGCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGACCCGCAGTCGCCGTTCGTCACTTCGG  >  NC_002947/387652‑387941

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