Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 253,359 A→T 5.8% intergenic (‑144/+27) 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,3590AT 5.8% 207.6 / 2.4 69intergenic (‑144/+27)PP_0202/PP_0203CBS domain‑containing protein/dipeptidase
Reads supporting (aligned to +/- strand):  ref base A (31/34);  new base T (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 = 7.05e-01

TACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGCCCGGTGTACTTCACCTCCGG  >  NC_002947/253215‑253505
                                                                                                                                                 |                                                                                                                                                    
tACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAAcc                                                                                                                                                  <  3:66485/149‑1 (MQ=255)
      ccTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAag                                                                                                                                            <  2:285636/149‑1 (MQ=255)
         ctGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAAGGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAgg                                                                                                                                         >  4:142436/1‑149 (MQ=255)
          tGCGGGTGTAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAagag                                                                                                                                          <  1:377873/147‑1 (MQ=255)
           gCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGTTTTAGCCGCGAaaag                                                                                                                                          <  2:75835/146‑1 (MQ=255)
           gCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGTTTTAGCCGCGAaaag                                                                                                                                          >  1:75835/1‑146 (MQ=255)
            cGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCACTGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGc                                                                                                                                        <  1:279792/148‑1 (MQ=255)
            cGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCg                                                                                                                                      <  8:74807/149‑1 (MQ=255)
                tGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGc                                                                                                                                        >  4:215885/1‑143 (MQ=255)
                tGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGc                                                                                                                                        <  3:215885/143‑1 (MQ=255)
                         cACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATc                                                                                                                        >  1:128450/1‑148 (MQ=255)
                         cACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATc                                                                                                                        >  7:164249/1‑148 (MQ=255)
                         cACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATc                                                                                                                        >  4:140930/1‑148 (MQ=255)
                              tctcGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCt                                                                                                                  >  2:302541/1‑149 (MQ=255)
                                tcGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAACGCTTc                                                                                                                >  4:125014/1‑149 (MQ=255)
                                  ggAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTc                                                                                                                <  2:325272/147‑1 (MQ=255)
                                  ggAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTc                                                                                                                <  2:84400/147‑1 (MQ=255)
                                  ggAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCgg                                                                                                              <  8:164249/149‑1 (MQ=255)
                                  ggAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCgg                                                                                                              <  2:128450/149‑1 (MQ=255)
                                  ggAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCg                                                                                                               <  2:343504/148‑1 (MQ=255)
                                  ggAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGAAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCgg                                                                                                              >  1:351526/1‑149 (MQ=255)
                                     gTAACTAAAGCATCCTGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGACCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCAt                                                                                                           <  1:313386/149‑1 (MQ=255)
                                         cTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGc                                                                                                       <  5:137388/149‑1 (MQ=255)
                                          tAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGcc                                                                                                      <  3:101765/149‑1 (MQ=255)
                                          tAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACAGCTAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGcc                                                                                                      >  3:178755/1‑149 (MQ=255)
                                            aaGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGccc                                                                                                     >  5:23059/1‑148 (MQ=255)
                                              gCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGAt                                                                                                  <  1:302541/149‑1 (MQ=255)
                                               cATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATc                                                                                                 <  8:322121/149‑1 (MQ=255)
                                                   cGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCcgcg                                                                                             >  2:112890/1‑149 (MQ=255)
                                                         ggCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAagag                                                                                                                                          >  6:160238/1‑100 (MQ=255)
                                                          gcgcTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCAt                                                                                      >  8:65076/1‑149 (MQ=255)
                                                            gcTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGCTCCGCGCCACCAttt                                                                                    >  5:66768/1‑149 (MQ=255)
                                                                   aCGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCAt                                                                                                           >  6:229050/1‑119 (MQ=255)
                                                                   aCGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCAt                                                                                                           <  5:229050/119‑1 (MQ=255)
                                                                   aCGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTcac                                                                             >  7:337615/1‑149 (MQ=255)
                                                                            cAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCaa                                                                    <  8:271405/149‑1 (MQ=255)
                                                                               gAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGATTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTc                                                                 >  5:47211/1‑149 (MQ=255)
                                                                                aaaCGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGCCTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAg                                                                   <  1:112890/148‑1 (MQ=255)
                                                                                   cGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGc                                                             <  6:109759/149‑1 (MQ=255)
                                                                                     gCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTACCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCAt                                                           <  3:142436/149‑1 (MQ=255)
                                                                                        tGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATGTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGaa                                                        <  2:351526/149‑1 (MQ=255)
                                                                                           aGGCGATTGAAAAACGTCATGGCTCTTCGCATGTAGCGGCCTATTCGCGGCTAAAACCGCGAAGACGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCCATCGGCGCCACCATTTCGCTCACCACCTCGAAGACCAGCATGAACTg                                                     >  7:253141/1‑149 (MQ=255)
                                                                                            ggCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGc                                                    <  6:47211/149‑1 (MQ=255)
                                                                                                aTTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCa                                                >  2:247516/1‑149 (MQ=255)
                                                                                                    aaaaaCGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTcac                                                                             >  4:210715/1‑116 (MQ=255)
                                                                                                    aaaaaCGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTcac                                                                             <  3:210715/116‑1 (MQ=255)
                                                                                                    aaaaaCGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACgg                                             >  7:103448/1‑148 (MQ=255)
                                                                                                    aaaaaCGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGt                                            >  3:328754/1‑149 (MQ=255)
                                                                                                      aaaCGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGAATCTGCCCGATCCGCGCTACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCt                                          >  6:122741/1‑149 (MQ=255)
                                                                                                       aaCGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCtt                                         <  7:65076/149‑1 (MQ=255)
                                                                                                          gTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑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                                                     >  4:74623/1‑111 (MQ=255)
                                                                                                                                  gCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTg                                                     <  3:74623/111‑1 (MQ=255)
                                                                                                                                  gCCTCTTCGCGGTTTTAGCCGCGAAAAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTcgttgtcgttgt                        >  8:102687/1‑140 (MQ=255)
                                                                                                                                  gCCTCTTCGCGGTTTTAGCCGCGAAAAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTcgttgtcgttgt                        <  7:102687/140‑1 (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   >  3:159176/1‑146 (MQ=255)
                                                                                                                                                 aaaCCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGCCCGGTGTACTTCACCTCCg   <  4:159176/146‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGG‑‑TACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGCCCGGTGTACTTCACCTCCGG  >  NC_002947/253215‑253505

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: