Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 766,922 C→T 6.3% C138C (TGC→TGT)  metB → cystathionine gamma‑synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947766,9220CT 6.3% 188.1 / 2.2 66C138C (TGC→TGT) metBcystathionine gamma‑synthase
Reads supporting (aligned to +/- strand):  ref base C (32/28);  new base T (2/2);  total (36/30)
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.12e-01

TCGCGTTCAGCAGCGGCATGGCGGCCATCAGTGCCGTGCTGCACACCTTCCTGTCCAACGGCAAACGCGTGGTGTCGACCCGCGACAGCTATGGCGGCACCAACAAGATCTTCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGCCTGTGCGATA‑CCTTGGACACTGAAGCGCTGGAGCGCGAGATTGCCGCCGGTTGCGACCTGCTGTACCTGGAAACCCCCACCAACCCGACGCTGAAGGTGCTGGACATTCGTCGGCTGAGCGCCGCTGCGCACAAGGTGGGGGCGCTGG  >  NC_002947/766774‑767070
                                                                                                                                                     |                                                                                                                                                     
tCGCGTTCAGCAGCGGCATGGCGGCCATCAGTGCCGTGCTGCACACCTTCCTGTCCAACGGCAAACGCGTGGTGTCGACCCGCGACAGCTATGGCGGCACCAACAAGATCTTCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGc                                                                                                                                                       <  1:160850/149‑1 (MQ=255)
  gcgTTCAGCAGCGGCATGGCGGCCATCAGTGCCGTGCTGCACACCTTCCTGTCCAACGGCAAACGCGTGGTGTCGACCCGCGACAGCTATGGCGGCACCAACAAGATCTTCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGCCt                                                                                                                                                     >  4:6550/1‑149 (MQ=255)
       cagcagCGGCATGGCGGCCATCAGTGCCGTGCTGCACACCTTCCTGTCCAACGGCAAACGCGTGGTGTCGACCCGCGACAGCTATGGCGGCACCAACAAGATCTTCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGCCTGTGc                                                                                                                                                 <  4:124423/148‑1 (MQ=255)
          cagcGGCATGGCGGCCATCAGTGCCGTGCTGCACACCTTCCTGTCCAACGGCAAACGCGTGGTGTCGACCCGCGACAGCTATGGCGGCACCAACAAGATCTTCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGCCTGTGCGATa                                                                                                                                             >  3:215115/1‑149 (MQ=255)
            gcGGCATGGCGGCCATCAGTGCCGTGCTGCACACCTTCCTGTCCAACGGCAAACGCGTGGTGTCGACCCGCGACAGCTATGGCGGCACCAACAAGAACTTCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGCCTGTGCGATA‑cc                                                                                                                                          <  1:90934/149‑1 (MQ=255)
                    gcggcCATCAGTGCCGTGCTGCACACCTTCCTGTCCAACGGCAAACGCGTGGTGTCGACCCGCGACAGCTATGGCGGCACCAACAAGATCTTCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGTCTGTGCGATA‑CCTTGGacac                                                                                                                                  >  2:68188/1‑149 (MQ=255)
                    gcggcCATCAGTGCCGTGCTGCACACCTTCCTGTCCAACGGCAAACGCGTGGTGTCGACCCGCGACAGCTATGGCGGCACACACAAGATGTTCGAAT‑AGTTCCTGCAGGGCAGGGGCGTGCAAGGGTGGCTGTGGGAAA‑CAGTGGacaa                                                                                                                                  >  8:49071/1‑148 (MQ=255)
                      ggcCATCAGTGCCGTGCTGCACACCTTCCTGTCCAACGGCAAACGCGTGGTGTCGACCCGCGACAGCTATGGCGGCACCAACAAGATCTTCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGCCTGTGCGATA‑CCTTGGACACt                                                                                                                                 >  3:302768/1‑148 (MQ=255)
                           tCAGTGCCGTGCTGCACACCTTCCTGTCCAACGGCAAACGCGTGGTGTCGACCCGCGACAGCTATGGCGGCACCAACAGGGGGGGGGTAGAAGGGAGAGACGCGCAAAGGCGTGCAAATGTGACAGTGCGATA‑ACATGGACACTGAAga                                                                                                                            >  1:283928/1‑148 (MQ=255)
                              gTGCCGGGCTGCACACCTTCCTGGCCAACGGCAAACGCGTGGTGTCGACCCGCGACAGCTATGGCGGCACCAACAAGATCTTCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGCCTGTGCGATA‑CCTTGGACACTGAAGCGCTg                                                                                                                        <  7:295366/149‑1 (MQ=255)
                               tGCCGTGCTGCACACCTTCCTGTCCAACGGCAAACGCGTGGTGTCGACCCGCGACAGCTATGGCGGCACCAACAAGATCTTCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGCCTGTGCGa                                                                                                                                               >  5:278227/1‑126 (MQ=255)
                               tGCCGTGCTGCACACCTTCCTGTCCAACGGCAAACGCGTGGTGTCGACCCGCGACAGCTATGGCGGCACCAACAAGATCTTCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGCCTGTGCGATA‑CCTTGGACACTGAAGCGCTg                                                                                                                        >  6:203542/1‑148 (MQ=255)
                                      ctgcACACCTTCCTGTCCAACGGCAAACGCGTGGTGTCGACCCGCGACAGCTATGGCGGCACCAACAAGATCTTCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGCCTGTGCGATA‑CCTTGGACACTg                                                                                                                                >  6:245068/1‑133 (MQ=255)
                                      ctgcACACCTTCCTGTCCAACGGCAAACGCGTGGTGTCGACCCGCGACAGCTATGGCGGCACCAACAAGATCTTCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGCCTGTGCGATA‑CCTTGGACACTg                                                                                                                                <  5:245068/133‑1 (MQ=255)
                                       tgcACACCTTCCTGTCCAACGGCAAACGCGTGGTGTCGACCCGCGACAGCTATGGCGGCACCAACAAGATCTTCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGCCTGTGCGATA‑CCTTGGACACTGAAGCGCTGGAgcgc                                                                                                                  >  2:107398/1‑146 (MQ=255)
                                         cacaCCTTCCTGTCCAACGGCAAACGCGTGGTGTCGACCCGCGACAGCTATGGCGGCACCAACAAGATCTTCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGCCTGTGCGATA‑CCTTGGACACTGAAGCGCTGGAGCGCGAGAt                                                                                                             >  8:59292/1‑149 (MQ=255)
                                             ccTTCCTGTCCAACGGCAAACGCGTGGTGTCGACCCGCGACAGCTATGGCGGCACCAACAAGATCTTCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGCCTGTGCGATA‑cc                                                                                                                                          <  6:54018/116‑1 (MQ=255)
                                             ccTTCCTGTCCAACGGCAAACGCGTGGTGTCGACCCGCGACAGCTATGGCGGCACCAACAAGATCTTCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGCCTGTGCGATA‑cc                                                                                                                                          >  5:54018/1‑116 (MQ=255)
                                                   tGTCCAACGGCAAACGCGTGGTGTCGACCCGCGACAGCTATGGCGGCACCAACAAGATCTTCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGCCTGTGCGATA‑CCTTGGACACTGAAGCGCTGGAGCGCGAGATTGCCGCCGGt                                                                                                   >  5:151116/1‑149 (MQ=255)
                                                    gTCCAACGGCAAACGCGTGGTGTCGACCCGCGACAGCTATGGCGGCACCAACAAGATCTTCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGCCTGTGCGATA‑CCTTGGACACTGAAGCGCTGGAGCGCGAGATTGCCGCCGGtt                                                                                                  >  1:63145/1‑149 (MQ=255)
                                                       cAACGGCAAACGCGTGGTGTCGACCCGCGACAGCTATGGCGGCACCAACAAGATCTTCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGCCTGTGCGATA‑CCTTGGACACTGAAGCGCTGGAGCGCGAGATTGCCGCCGGTTGCg                                                                                               >  7:60022/1‑149 (MQ=255)
                                                        aaCGGCAAACGCGTGGTGTCGACCCGCGACAGCTATGGCGGCACCAACAAGATCTTCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGCCTGTGCGATA‑CCTTGGACACTGAAGCGCTGGAGCGCGAGATTGCCGCCGGTTGc                                                                                                <  5:252581/147‑1 (MQ=255)
                                                        aaCGGCAAACGCGTGGTGTCGACCCGCGACAGCTATGGCGGCACCAACAAGATCTTCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGCCTGTGCGATA‑CCTTGGACACTGAAGCGCTGGAGCGCGAGATTGCCGCCGGTTGc                                                                                                >  6:252581/1‑147 (MQ=255)
                                                            gCAAACGCGTGGTGTCGACCCGCGACAGCTATGGCGGCACCAACAAGATCTTCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGCCTGTGCGATA‑CCTTGGACACTGAAGCGCTGGAGCGCGAGATTGCCGCCGGTTGCGACctg                                                                                          <  3:6550/149‑1 (MQ=255)
                                                                        tgtCGACCCGCGACAGCTATGGCGGCACCAACAAGATCTTCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGTCTGTGCGATA‑CCTTGGACACTGAAGCGCTGGAGCGCGAGATTGCCGCCGGTTGCGACCTGCTGTACCTGGaa                                                                              <  1:68188/149‑1 (MQ=255)
                                                                         gtCGACCCGCGACAGCTATGGCGGCACCAACAAGATCTTCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGCCTGTGCGATA‑CCTTGGACACTGAAGCGCTGGAGCGCGAGATTGCCGCCGGTTGCGACct                                                                                           >  8:317048/1‑135 (MQ=255)
                                                                         gtCGACCCGCGACAGCTATGGCGGCACCAACAAGATCTTCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGCCTGTGCGATA‑CCTTGGACACTGAAGCGCTGGAGCGCGAGATTGCCGCCGGTTGCGACct                                                                                           <  7:317048/135‑1 (MQ=255)
                                                                         gtCGACCCGCGACAGCTATGGCGGCACCAACAAGATCTTCGAAG‑AGTTCCTGCCGCGCATAGGCGTGCAAGTGTGCCTGTGCGATA‑CCTTGGACACTGAAGCGCTGGAGCGCGAGATTGCCGCCGGTTGCGACCTGCTGTACCTGGaaa                                                                             >  1:272435/1‑149 (MQ=255)
                                                                          tCGACCCGCGACAGCTATGGCGGCACCAACAAGATCTTCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGCCTGTGCGATA‑CCTTGGACACTGAAGCGCTGGAGCGCGAGATTGCCGCCGGTTGCGACCTGCTGTACCTGGAAAc                                                                            <  4:215115/149‑1 (MQ=255)
                                                                            gACCCGCGACAGCTATGGCGGCACCAACAAGATCTTCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGTCTGTGCGATA‑CCTTGGACACTGAAGCGCTGGAGCGCGAGATTGCCGCCGGTTGCGACCTGCTGTACCTGGAAAcac                                                                          >  6:323684/1‑147 (MQ=255)
                                                                            gACCCGCGACAGCTATGGCGGCACCAACAAGATCTTCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGCCTGTGCGATA‑CCTTGGACACTGAAGCGCTGGAGCGCGAGATTGCCGCCGGTTGCGACCTGCTGTACCTGGAAAccc                                                                          >  3:31775/1‑149 (MQ=255)
                                                                             aCCCGCGACAGCTATGGCGGCACCAACAAGATCTTCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGCCTGTGCGATA‑CCTTGGACACTGAAGCGCTGGAGCGCGAGATTGCCGCCGGTTGCGACCTGCTGTACCTGGaa                                                                              <  1:309568/144‑1 (MQ=255)
                                                                             aCCCGCGACAGCTATGGCGGCACCAACAAGATCTTCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGCCTGTGCGATA‑CCTTGGACACTGAAGCGCTGGAGCGCGAGATTGCCGCCGGTTGCGACCTGCTGTACCTGGaa                                                                              >  2:309568/1‑144 (MQ=255)
                                                                             aCCCGCGACAGCTATGGCGGCACCAACAAGATCTTCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGCCTGTGCGATA‑CCTTGGACACTGAAGCGCTGGAGCGCGAGATTGCCGCCGGTTGCGACCTGCTGTACCTGGAAAcccc                                                                         <  6:292345/149‑1 (MQ=255)
                                                                               ccGCGACAGCTATGGCGGCACCAACAAGATCTTCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGCCTGTGCGATA‑CCTTGGACACTGAAGCGCTGGAGCGCGAGATTGCCGCCGGTTGCGACCTGCTGTACCTGGAAACCCcca                                                                       >  6:115818/1‑149 (MQ=255)
                                                                                cgcgACTGCTATGGCTGCACCAACAAGATATTCGAAT‑AGTTCCTGCCGCGCATGGGCGTGCATGTGTGCCTGTGCGATAGTCTTGGACACTGAAGCGCTGGAGCGCGAGATTGCCGCCGGTTGCGACCTGCTGTACCTGGAAACCCcca                                                                       <  1:107398/149‑1 (MQ=255)
                                                                                      aGCTATGGCGGCACCAACAAGATCTTCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGCCTGTGCGATA‑CCTTGGACACTGAAGCGCTg                                                                                                                        >  5:167020/1‑93 (MQ=255)
                                                                                      aGCTATGGCGGCACCAACAAGATCTTCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGCCTGTGCGATA‑CCTTGGACACTGAAGCGCTg                                                                                                                        <  6:167020/93‑1 (MQ=255)
                                                                                        cTATGGCGGCACCAACAAGATCTTCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGCCTGTGCGATA‑CCTTGGACACTGAAGCGCTGGAGCGCGAGATTGCCGCCGGTTGCGACCTGCTGTACCTGGAAACCCCCACCAACCCGa                                                              >  3:327949/1‑149 (MQ=255)
                                                                                         tATGGCGGCACCAACAAGATCTTCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGCCTGTGCGATA‑CCTTGGACACTGAAGCGCTGGa                                                                                                                      <  8:46697/92‑1 (MQ=255)
                                                                                         tATGGCGGCACCAACAAGATCTTCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGCCTGTGCGATA‑CCTTGGACACTGAAGCGCTGGa                                                                                                                      >  7:46697/1‑92 (MQ=255)
                                                                                         tATGGCGGCACCAACAAGATCTTCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGCCTGTGCGATA‑CCTTGGACACTGAAGCGCTGGAGCGCGAGATTGCCGCCGGTTGCGACCTGCTGTACCTGGAAACCCCCACCAACCCGAc                                                             >  5:24923/1‑149 (MQ=255)
                                                                                           tGGCGGCACCAACAAGATCTTCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGTCTGTGCGATA‑CCTTGGACACTGAAGCGCTGGAGCGCGAGATTGCCGCCGGTTGCGACCTGCTGTACCTGGAAACCCCCACCAACCCGACGc                                                           <  5:323684/149‑1 (MQ=255)
                                                                                             gcggcACCAACAAGATCTTCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGCCTGTGCGATA‑CCTTGGACACTGAAGCGCTGGAGCGCGAGATTGCCGCCGGTTGCGACCTGCTGTACCTGGAAACCCCCACCAACCCGACGCt                                                          >  3:186202/1‑148 (MQ=255)
                                                                                             gcggcACCAACAAGATCTTCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGCCTGTGCGAGA‑CCTTGGACACTGAAGCGCTGGAGCGCGAGATTGCCGCCGGTTGCGACCTGCTGTACCTGGAAACCCCCACCAACCCGACGCTg                                                         <  2:272435/149‑1 (MQ=255)
                                                                                               ggcACCAACAAGATCTTCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGCCTGTGCGATA‑CCTTGGACACTGAAGCGCTGGAGCGCGAGATTGCCGCCGGTTGCGACCTGCTGTACCTGGAAACCCCCACCAACCCGACGCTGaa                                                       <  6:151116/149‑1 (MQ=255)
                                                                                                     aacaaGATCTTCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGCCTGTGCGATA‑CCTTGGACACTGAAGCGCTGGAGCGCGAGATTGCCGCCGGTTGCGACCTGCTGTACCTGGAAACCCCCACCAACCCGACGCTGAAGGTGCt                                                 <  7:59292/149‑1 (MQ=255)
                                                                                                      acaaGATCTTCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGCCTGTGCGATA‑CCTTGGACACTGAAGCGCTGGAGCGCGAGATTGCCGCCGGTTGCGACCTGCTGTACCTGGAAACCCCCACCAACCCGACGCTGAAGGTGCTg                                                <  3:152289/149‑1 (MQ=255)
                                                                                                      acaaGATCTTCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGCCTGTGCGATA‑CCTTGGACACTGAAGCGCTGGAGCGCGAGATTGCCGCCGGTTGCGACCTGCTGTACCTGGAAACCCCCACCAACCCGACGCTGAAGGTGCTg                                                >  1:125753/1‑149 (MQ=255)
                                                                                                       caaGATCTTCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGCCTGTGCGATA‑CCTTGGACACTGAAGCGCTGGAGCGCGAGATTGCCGCCGGTTGCGACCTGCTGTACCTGGAAACCCCCACCAACCCGACGCTGAAGGTGCTg                                                >  6:223470/1‑148 (MQ=255)
                                                                                                       caaGATCTTCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGCCTGTGCGATA‑CCTTGGACACTGAAGCGCTGGAGCGCGAGATTGCCGCCGGTTGCGACCTGCTGTACCTGGAAACCCCCACCAACCCGACGCTGAAGGTGCTg                                                >  3:224459/1‑148 (MQ=255)
                                                                                                             cTTCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGCCTGTGCGATA‑CCTTg                                                                                                                                       >  1:79877/1‑55 (MQ=255)
                                                                                                             cTTCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGCCTGTGCGATA‑CCTTg                                                                                                                                       <  2:79877/55‑1 (MQ=255)
                                                                                                              ttCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGCCTGTGCGATA‑CCTTGGACACTGAAGCGCTGGAGCGCGAGATTGCCGCCGGTTGCGACCTGCTGTACCTGGAAACCCCCACCAACCCGACGCTGAAGGTGCTGGACATtcg                                        >  5:91304/1‑149 (MQ=255)
                                                                                                               tCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGCCTGTGCGATA‑CCTTGGACACTGAAGCGCTGGAGCGCGAGATTGCCGCCGGTTGCGACCTGCTGTACCTGGAAACCCCCACCAACCCGACGCTGAAGGTGCTGGACATtcgt                                       >  7:113840/1‑149 (MQ=255)
                                                                                                                cGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGCCTGTGCGATA‑CCTTGGACACTGAAGCGCTGGAGCGCGAGATTGCCGCCGGTTGCGACCTGCTGTACCTGGAAACCCCCACCAACCCGACGCTGAAGGTGCTGGACATtc                                         <  4:224459/146‑1 (MQ=255)
                                                                                                                 gAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGCCTGTGCGATA‑CCTTGGACACTGAAGCGCTGGAGCGCGAGATTGCCGCCGGTTGCGACCTGCTGTACCTGGAAACCCCCACCAACCCGACGCTGAAGGTGCTGGACATtcgtc                                      <  8:113840/148‑1 (MQ=255)
                                                                                                                  aaG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGCCTGTGCGATA‑CCTTGGACACTGAAGCGCTGGAGCGCGAGATTGCCGCCGGTTGCGACCTGCTGTACCTGGAAACCCCCACCAACCCGACGCTGAAGGTGCTGGACATTCGTCgg                                    >  4:281114/1‑149 (MQ=255)
                                                                                                                  aaG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGCCTGTGCGATA‑CCTTGGACACTGAAGCGCTGGAGCGCGAGATTGCCGCCGGTTGCGACCTGCTGTACCTGGAAACCCCCACCAACCCGACGCTGAAGGTGCTGGACATTCGTCgg                                    <  2:283928/149‑1 (MQ=255)
                                                                                                                           cTGCCGCGCATGGGCGTGCAAGTGTGCCTGTGCGATA‑CCTTGGACACTGAAGCGCTGGAGCGCGAGATTGCCGCCGGTTGCGACCTGCTGTACCTGGAAACCCCCACCAACCCGACGCTGAAGGTGCTGGACATTCGTCGGCTGAgcgc                            <  4:327949/149‑1 (MQ=255)
                                                                                                                               cgcgCATGGGCGTGCAAGTGTGCCTGTGCGATA‑CCTTGGACACTGAAGCGCTGGAGCGCGAGATTGCCGCCGGTTGCGACCTGCTGTACCTGGAAACCCCCACCAACCCGACGCTGAAGGTGCTGGACATTCGTCGGCTGAGCGCCGCt                        <  4:302768/149‑1 (MQ=255)
                                                                                                                               cgcgCATGGGCGTGCAAGTGTGCCTGTGCGATA‑CCTTGGACACTGAAGCGCTGGAGCGCGAGATTGCCGCCGGTTGCGACCTGCTGTACCTGGAAACCCCCACCAACCCGACGCTGAAGGTGCTGGACATTCGTCGGCTGAGCGCCGCt                        <  5:203542/149‑1 (MQ=255)
                                                                                                                                    aTGGGCGTGCAAGTGTGCCTGTGCGATA‑CCTTGGACACTGAAGCGCTGGAGCGCGAGATTGCCGCCGGTTGCGACCTGCTGTACATGGAAACCCCCACCAACCCGACGCTGAAGGTGCTGGACATTCGTCGGCTGAGCGCCGCTgcgaa                   >  6:225582/1‑147 (MQ=255)
                                                                                                                                             cAAGTGTGCCTGTGCGATA‑CCTTGGACACTGAAGCGCTGGAGCGCGAGATTGCCGCCGGTTGCGACCTGCTGTACCTGGAAACCCCCACCAACCCGACGCTGAAGGTGCTGGACATTCGTCGGCTGAGCGCCGCTGCGCACAAGGTgg           <  3:281114/148‑1 (MQ=255)
                                                                                                                                              aaGTGTGCCTGTGCGATA‑CCTTGGACACTGAAGCGCTGGAGCGCGAGATTGCCGCCGGTTGCGAc                                                                                             >  3:52280/1‑65 (MQ=255)
                                                                                                                                              aaGTGTGCCTGTGCGATA‑CCTTGGACACTGAAGCGCTGGAGCGCGAGATTGCCGCCGGTTGCGAc                                                                                             <  4:52280/65‑1 (MQ=255)
                                                                                                                                                   tgCCTGTGCGATA‑CCTTGGACACTGAAGCGCTGGAGCGCGAGATTGCCGCCGGTTGCGACCTGCTGTACCTGGAAACCCCCACCAACCCGACGCTGAAGGTGCTGGACATTCGTCGGCTGAGCGCCGCTGCGCACAAGGTGGGGGCGCt    <  7:49071/149‑1 (MQ=255)
                                                                                                                                                     ccTGTGCGATA‑CCTTGGACACTGAAGCGCTGGAGCGCGAGATTGCCGCCGGTTGCGACCTGCTGTACCTGGAAACCCCCACCACCCCGACGCTGAAGGTGCTGGACATTCGTCGGCTGAGCGCCGCTGCGCACAAGGGGGGGGCGCtgg  >  4:193954/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCGCGTTCAGCAGCGGCATGGCGGCCATCAGTGCCGTGCTGCACACCTTCCTGTCCAACGGCAAACGCGTGGTGTCGACCCGCGACAGCTATGGCGGCACCAACAAGATCTTCGAAG‑AGTTCCTGCCGCGCATGGGCGTGCAAGTGTGCCTGTGCGATA‑CCTTGGACACTGAAGCGCTGGAGCGCGAGATTGCCGCCGGTTGCGACCTGCTGTACCTGGAAACCCCCACCAACCCGACGCTGAAGGTGCTGGACATTCGTCGGCTGAGCGCCGCTGCGCACAAGGTGGGGGCGCTGG  >  NC_002947/766774‑767070

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