Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 387,792 C→T 7.7% S311F (TCC→TTC) ‡ glyA‑I → serine hydroxymethyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947387,7920CT 7.7% 226.1 / 6.1 78S311F (TCC→TTC) ‡glyA‑Iserine hydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base C (33/39);  new base T (3/3);  total (36/42)
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.66e-01

TGGCCCGCTGATGCACGTGATCGCTGCCAAGGCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGACCCGCAGTCGCCGTTCGT  >  NC_002947/387646‑387933
                                                                                                                                                  |                                                                                                                                             
tGGCCCGCTGATGCACGTGATCGCTGCCAAGGCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCg                                                                                                                                             >  1:2753/1‑149 (MQ=255)
   cccGCTGATGCACGTGATCGCTGCCAAGGCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTg                                                                                                                                          >  6:285343/1‑149 (MQ=255)
             cACGTGATCGCTGCCAAGGCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACaa                                                                                                                                >  1:336842/1‑149 (MQ=255)
              aCGTGATCGCTGCCAAGGCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGCCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAac                                                                                                                               >  3:94758/1‑149 (MQ=255)
                gTGATCGCTGCCAAGGCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAacca                                                                                                                             <  2:139855/149‑1 (MQ=255)
                    tCGCTGCCAAGGCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTGTCGGTGCCACCGGAAaccacatc                                                                                                                         <  4:128314/149‑4 (MQ=255)
                      gCTGCCAAGGCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGt                                                                                                                        >  3:321517/1‑148 (MQ=255)
                      gCTGCCAAGGCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGt                                                                                                                        <  4:321517/148‑1 (MQ=255)
                       cTGCCAAGGCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTc                                                                                                                      <  2:336842/149‑1 (MQ=255)
                          ccAAGGCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTg                                                                                                                   <  5:242519/149‑1 (MQ=255)
                           cAAGGCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTgg                                                                                                                  >  3:73241/1‑149 (MQ=255)
                            aaGGCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTGTCGGTGCCACCGGAAaccacatcgta                                                                                                                      <  1:8696/144‑7 (MQ=255)
                            aaGGCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTGTCGGTGCCACCGGAAaccacatcgta                                                                                                                      >  2:8696/1‑138 (MQ=255)
                             aGGCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTc                                                                                                                <  7:11670/149‑1 (MQ=255)
                             aGGCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTc                                                                                                                <  1:125353/149‑1 (MQ=255)
                                     tgCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGa                                                                                                         <  4:304615/148‑1 (MQ=255)
                                           aaGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTACTGGTCAGCCTGATCCGt                                                                                                    >  7:234602/1‑147 (MQ=255)
                                           aaGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGCTTTCCGGTGGCACAGACAACCACCTGTTCATGGACAGCCTGATCCGt                                                                                                    >  8:217625/1‑147 (MQ=255)
                                            aGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTc                                                                                                   >  4:143694/1‑147 (MQ=255)
                                            aGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTc                                                                                                   <  3:143694/147‑1 (MQ=255)
                                             ggAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAgg                                                                                                >  8:17995/1‑149 (MQ=255)
                                             ggAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAgg                                                                                                <  7:17995/149‑1 (MQ=255)
                                              gAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAggg                                                                                               >  4:189410/1‑149 (MQ=255)
                                                    cttgAGCCTGTATTGAAGGCCTACCAGCAGCAAGTGATCGAAATCGCCCAGGCCATGGCCCAGGTGTTCATCGACGGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACATGTTCGTGGTCAGCCTGATCCGTCAGGGCCTGAc                                                                                         <  6:282982/146‑1 (MQ=255)
                                                      cgAGCCTGGAATCAAGCCCTCCCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACAGTGGCTACGATGTGGTTTCCGGCGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTg                                                                                       <  4:16702/148‑1 (MQ=255)
                                                       gAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTgg                                                                                      <  7:217625/149‑1 (MQ=255)
                                                         gCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGc                                                                                     >  5:313051/1‑148 (MQ=255)
                                                             ggATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAacca                                                                                                                             <  3:225719/104‑1 (MQ=255)
                                                             ggATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAacca                                                                                                                             >  4:225719/1‑104 (MQ=255)
                                                              gATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAgac                                                                               <  7:116460/149‑1 (MQ=255)
                                                                ttCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAgacgc                                                                             >  5:94396/1‑149 (MQ=255)
                                                                  cAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAgacgcc                                                                            <  7:168579/148‑1 (MQ=255)
                                                                   aaGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGAc                                                                                                                                  <  2:39532/93‑1 (MQ=255)
                                                                   aaGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGAc                                                                                                                                  >  1:39532/1‑93 (MQ=255)
                                                                     ggCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAgacgccgacg                                                                        >  5:314434/1‑149 (MQ=255)
                                                                       ccTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTc                                                                                                                >  3:147623/1‑107 (MQ=255)
                                                                       ccTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTc                                                                                                                <  4:147623/107‑1 (MQ=255)
                                                                             agcagcAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTg                                                                <  6:341337/149‑1 (MQ=255)
                                                                                 gcaAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTc                                                            >  8:282070/1‑149 (MQ=255)
                                                                                  caAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTg                                                                <  6:175970/144‑1 (MQ=255)
                                                                                  caAGTGATAGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTg                                                                >  5:175970/1‑144 (MQ=255)
                                                                                    aGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTg                                                         <  2:2753/149‑1 (MQ=255)
                                                                                       gATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTgcg                                                       >  1:151100/1‑148 (MQ=255)
                                                                                        aTCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCTGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTTGGTCGTGCGCa                                                     <  4:94758/149‑1 (MQ=255)
                                                                                           gAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCata                                                   >  1:110393/1‑148 (MQ=255)
                                                                                            aaaaCGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCatat                                                  >  5:251911/1‑148 (MQ=255)
                                                                                             aaaCGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCa                                                >  3:64549/1‑149 (MQ=255)
                                                                                                 gCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCg                                             >  7:299032/1‑148 (MQ=255)
                                                                                                 gCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGt                                            >  6:248448/1‑149 (MQ=255)
                                                                                                  cccAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTc                                           <  3:189410/149‑1 (MQ=255)
                                                                                                   ccAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCaaa                                                                                  <  4:304215/109‑1 (MQ=255)
                                                                                                   ccAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTca                                          <  5:2624/149‑1 (MQ=255)
                                                                                                    cAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTcaa                                         >  5:253703/1‑149 (MQ=255)
                                                                                                          aTGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGaa                                   <  6:313051/149‑1 (MQ=255)
                                                                                                               ccAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAggg                                                                                               <  6:242788/84‑1 (MQ=255)
                                                                                                               ccAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAggg                                                                                               >  5:242788/1‑84 (MQ=255)
                                                                                                                cAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGt                             <  8:299032/149‑1 (MQ=255)
                                                                                                                cAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGt                             >  1:72554/1‑149 (MQ=255)
                                                                                                                 aGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTg                                                                                                                   <  1:90653/62‑1 (MQ=255)
                                                                                                                 aGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTg                                                                                                                   >  2:90653/1‑62 (MQ=255)
                                                                                                                  ggTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCTTGc                           <  6:314434/149‑1 (MQ=255)
                                                                                                                  ggTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGc                           <  5:248448/149‑1 (MQ=255)
                                                                                                                  ggTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGc                           <  4:73241/149‑1 (MQ=255)
                                                                                                                         aTCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTggg                                                              <  6:169719/107‑1 (MQ=255)
                                                                                                                         aTCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATGACCGTCAACAAGAAc                                  >  5:307248/1‑135 (MQ=255)
                                                                                                                         aTCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAAc                                  <  6:307248/135‑1 (MQ=255)
                                                                                                                             aCCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGACCCg                <  4:64549/149‑1 (MQ=255)
                                                                                                                             aCCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGACCCg                <  2:110393/149‑1 (MQ=255)
                                                                                                                               cGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGACCCGCa              <  6:253703/149‑1 (MQ=255)
                                                                                                                                 tGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGaa                       <  7:12507/138‑1 (MQ=255)
                                                                                                                                 tGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGaa                       >  8:12507/1‑138 (MQ=255)
                                                                                                                                    cTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGACCCGCAGTCGc         <  5:285343/149‑1 (MQ=255)
                                                                                                                                      aacaggtgGTTGTCGGTGCCACCGGAAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGACCCGCAGTCGCCg       >  6:82916/6‑149 (MQ=255)
                                                                                                                                        caggtgGTTGTCGGTGCCACCGGAAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGACCCGCAGTCGCCGtt     >  8:284030/4‑149 (MQ=255)
                                                                                                                                          tgtgGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCg                                             >  1:116062/1‑107 (MQ=255)
                                                                                                                                          tgtgGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCg                                             <  2:116062/107‑1 (MQ=255)
                                                                                                                                           gtgGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGACCCGCAGTCGCCGTTCGt  >  4:90447/1‑149 (MQ=255)
                                                                                                                                           gtgGTTGTCGGTGCCACCGGAAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCAcc                                              <  6:267281/105‑1 (MQ=39)
                                                                                                                                                  ccGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTcaac                                        <  7:305274/104‑1 (MQ=255)
                                                                                                                                                  ccGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTcaac                                        >  8:305274/1‑104 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
TGGCCCGCTGATGCACGTGATCGCTGCCAAGGCAGTGTGCTTCAAGGAAGCGCTGGAGCCTGGATTCAAGGCCTACCAGCAGCAAGTGATCGAAAACGCCCAGGCCATGGCCCAGGTGTTCATCGACCGTGGCTACGATGTGGTTTCCGGTGGCACCGACAACCACCTGTTCCTGGTCAGCCTGATCCGTCAGGGCCTGACTGGCAAAGACGCCGACGCCGCCCTGGGTCGTGCGCATATCACCGTCAACAAGAACGCCGTGCCGAACGACCCGCAGTCGCCGTTCGT  >  NC_002947/387646‑387933

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: ‑