Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 253,361 A→G 5.8% intergenic (‑146/+25) PP_0202 ← / ← PP_0203 CBS domain‑containing protein/dipeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947253,3610AG 5.8% 201.8 / 2.2 69intergenic (‑146/+25)PP_0202/PP_0203CBS domain‑containing protein/dipeptidase
Reads supporting (aligned to +/- strand):  ref base A (31/34);  new base G (2/2);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.61e-01

TACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCA‑TGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGCCCGGTGTACTTCACCTCCGGCAT  >  NC_002947/253215‑253508
                                                                                                                                                   |                                                                                                                                                     
tACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCA‑TGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAAcc                                                                                                                                                     <  3:66485/149‑1 (MQ=255)
      ccTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCA‑TGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAag                                                                                                                                               <  2:285636/149‑1 (MQ=255)
         ctGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAAGGTCA‑TGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAgg                                                                                                                                            >  4:142436/1‑149 (MQ=255)
          tGCGGGTGTAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCA‑TGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAagag                                                                                                                                             <  1:377873/147‑1 (MQ=255)
           gCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCA‑TGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGTTTTAGCCGCGAaaag                                                                                                                                             >  1:75835/1‑146 (MQ=255)
           gCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCA‑TGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGTTTTAGCCGCGAaaag                                                                                                                                             <  2:75835/146‑1 (MQ=255)
            cGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCACTGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGc                                                                                                                                           <  1:279792/148‑1 (MQ=255)
            cGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCA‑TGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCg                                                                                                                                         <  8:74807/149‑1 (MQ=255)
                tGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCA‑TGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGc                                                                                                                                           >  4:215885/1‑143 (MQ=255)
                tGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCA‑TGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGc                                                                                                                                           <  3:215885/143‑1 (MQ=255)
                         cACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCA‑TGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATc                                                                                                                           >  1:128450/1‑148 (MQ=255)
                         cACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCA‑TGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATc                                                                                                                           >  7:164249/1‑148 (MQ=255)
                         cACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCA‑TGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATc                                                                                                                           >  4:140930/1‑148 (MQ=255)
                              tctcGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCA‑TGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCt                                                                                                                     >  2:302541/1‑149 (MQ=255)
                                tcGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCA‑TGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAACGCTTc                                                                                                                   >  4:125014/1‑149 (MQ=255)
                                  ggAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCA‑TGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTc                                                                                                                   <  2:325272/147‑1 (MQ=255)
                                  ggAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCA‑TGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTc                                                                                                                   <  2:84400/147‑1 (MQ=255)
                                  ggAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCA‑TGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCgg                                                                                                                 <  8:164249/149‑1 (MQ=255)
                                  ggAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCA‑TGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCgg                                                                                                                 <  2:128450/149‑1 (MQ=255)
                                  ggAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCA‑TGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCg                                                                                                                  <  2:343504/148‑1 (MQ=255)
                                  ggAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCA‑TGGCTGTTCGCTTGAAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCgg                                                                                                                 >  1:351526/1‑149 (MQ=255)
                                     gTAACTAAAGCATCCTGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCA‑TGGCTGTTCGCTTGTAGCGACCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCAt                                                                                                              <  1:313386/149‑1 (MQ=255)
                                         cTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCA‑TGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGc                                                                                                          <  5:137388/149‑1 (MQ=255)
                                          tAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCA‑TGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGcc                                                                                                         <  3:101765/149‑1 (MQ=255)
                                          tAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCA‑TGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACAGCTAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGcc                                                                                                         >  3:178755/1‑149 (MQ=255)
                                            aaGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCA‑TGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGccc                                                                                                        >  5:23059/1‑148 (MQ=255)
                                              gCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCA‑TGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGAt                                                                                                     <  1:302541/149‑1 (MQ=255)
                                               cATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCA‑TGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATc                                                                                                    <  8:322121/149‑1 (MQ=255)
                                                   cGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCA‑TGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCcgcg                                                                                                >  2:112890/1‑149 (MQ=255)
                                                         ggCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCA‑TGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAagag                                                                                                                                             >  6:160238/1‑100 (MQ=255)
                                                          gcgcTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCA‑TGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCAt                                                                                         >  8:65076/1‑149 (MQ=255)
                                                            gcTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCA‑TGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGCTCCGCGCCACCAttt                                                                                       >  5:66768/1‑149 (MQ=255)
                                                                   aCGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCA‑TGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCAt                                                                                                              >  6:229050/1‑119 (MQ=255)
                                                                   aCGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCA‑TGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCAt                                                                                                              <  5:229050/119‑1 (MQ=255)
                                                                   aCGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCA‑TGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTcac                                                                                >  7:337615/1‑149 (MQ=255)
                                                                            cAGGAAACGGCATGGAGGCGATTGAAAAACGTCA‑TGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCaa                                                                       <  8:271405/149‑1 (MQ=255)
                                                                               gAAACGGCATGGAGGCGATTGAAAAACGTCA‑TGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGATTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTc                                                                    >  5:47211/1‑149 (MQ=255)
                                                                                aaaCGGCATGGAGGCGATTGAAAAACGTCA‑TGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGCCTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAg                                                                      <  1:112890/148‑1 (MQ=255)
                                                                                   cGGCATGGAGGCGATTGAAAAACGTCA‑TGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGc                                                                <  6:109759/149‑1 (MQ=255)
                                                                                     gCATGGAGGCGATTGAAAAACGTCA‑TGGCTGTTCGCTTGTACCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCAt                                                              <  3:142436/149‑1 (MQ=255)
                                                                                        tGGAGGCGATTGAAAAACGTCA‑TGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATGTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGaa                                                           <  2:351526/149‑1 (MQ=255)
                                                                                           aGGCGATTGAAAAACGTCA‑TGGCTCTTCGCATGTAGCGGCCTATTCGCGGCTAAAACCGCGAAGACGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCCATCGGCGCCACCATTTCGCTCACCACCTCGAAGACCAGCATGAACTg                                                        >  7:253141/1‑149 (MQ=255)
                                                                                            ggCGATTGAAAAACGTCA‑TGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGc                                                       <  6:47211/149‑1 (MQ=255)
                                                                                                aTTGAAAAACGTCA‑TGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCa                                                   >  2:247516/1‑149 (MQ=255)
                                                                                                    aaaaaCGTCA‑TGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTcac                                                                                <  3:210715/116‑1 (MQ=255)
                                                                                                    aaaaaCGTCA‑TGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTcac                                                                                >  4:210715/1‑116 (MQ=255)
                                                                                                    aaaaaCGTCA‑TGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACgg                                                >  7:103448/1‑148 (MQ=255)
                                                                                                    aaaaaCGTCA‑TGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGt                                               >  3:328754/1‑149 (MQ=255)
                                                                                                      aaaCGTCA‑TGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGAATCTGCCCGATCCGCGCTACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCt                                             >  6:122741/1‑149 (MQ=255)
                                                                                                       aaCGTCA‑TGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCtt                                            <  7:65076/149‑1 (MQ=255)
                                                                                                          gTCA‑TGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGaa                                         <  8:103448/149‑1 (MQ=255)
                                                                                                               tGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCTAAGAGGCCGG‑‑TACTTCGCATCAAAGCTTCGGCATCTGCCCGATCAGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACtcg                                     <  4:178755/149‑1 (MQ=255)
                                                                                                                ggCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACtcgt                                    >  8:28321/1‑149 (MQ=255)
                                                                                                                 gCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACtcgtt                                   <  6:301563/149‑1 (MQ=255)
                                                                                                                     ttCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACtcgttgtc                                <  7:28321/148‑1 (MQ=255)
                                                                                                                        gCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACtcgttgtcgttg                            >  2:58480/1‑149 (MQ=255)
                                                                                                                         cTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGCTCCGCGACAACATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTcgttgtccttgt                           >  7:249543/1‑149 (MQ=255)
                                                                                                                         cTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGGATGATCTGCTCCACGGTCGTGTACTggttgtcgttgt                           >  5:67138/1‑149 (MQ=255)
                                                                                                                         cTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTcgttgtcgttgt                           <  4:328754/149‑1 (MQ=255)
                                                                                                                         cTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTcgttgtcgttgt                           >  7:155249/1‑149 (MQ=255)
                                                                                                                                cGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGCCCGgt                    <  1:247516/149‑1 (MQ=255)
                                                                                                                                  gCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTg                                                        <  3:74623/111‑1 (MQ=255)
                                                                                                                                  gCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTg                                                        >  4:74623/1‑111 (MQ=255)
                                                                                                                                  gCCTCTTCGCGGTTTTAGCCGCGAAAAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTcgttgtcgttgt                           <  7:102687/140‑1 (MQ=255)
                                                                                                                                  gCCTCTTCGCGGTTTTAGCCGCGAAAAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTcgttgtcgttgt                           >  8:102687/1‑140 (MQ=255)
                                                                                                                                        tCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGCCCGGTGTACTTc             <  8:337615/148‑1 (MQ=255)
                                                                                                                                            ggCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGc                         <  7:113809/132‑1 (MQ=255)
                                                                                                                                            ggCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGc                         >  8:113809/1‑132 (MQ=255)
                                                                                                                                             gCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGCCCGGTGTACTTCACCTcc       >  3:222436/1‑149 (MQ=255)
                                                                                                                                                 aaaCCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGCCCGGTGTACTTCACCTCCgg     <  6:23059/147‑1 (MQ=255)
                                                                                                                                                 aaaCCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGCCCGGTGTACTTCACCTCCg      <  4:159176/146‑1 (MQ=255)
                                                                                                                                                 aaaCCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGCCCGGTGTACTTCACCTCCg      >  3:159176/1‑146 (MQ=255)
                                                                                                                                                  aaCCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGCCCGGTGTACTTCACCTCCGGCAt  >  6:186027/1‑149 (MQ=255)
                                                                                                                                                   aCCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGCTCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCCCGGTCTTGAACTCGGTGTCGTTGTGCCCGGGGTAc                >  6:381411/1‑134 (MQ=255)
                                                                                                                                                   aCCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGCCCGGTGTAc                <  5:381411/134‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCA‑TGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGCCCGGTGTACTTCACCTCCGGCAT  >  NC_002947/253215‑253508

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