Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 253,361 A→G 7.5% 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 7.5% 221.1 / 6.8 80intergenic (‑146/+25)PP_0202/PP_0203CBS domain‑containing protein/dipeptidase
Reads supporting (aligned to +/- strand):  ref base A (37/36);  new base G (3/3);  total (40/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

GGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGCCCGGTGTACTTCACCTCCGGCATG  >  NC_002947/253218‑253509
                                                                                                                                               |                                                                                                                                                    
gggCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACcgc                                                                                                                                                  >  7:347582/1‑148 (MQ=255)
    ctctGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGaa                                                                                                                                               >  1:79187/1‑147 (MQ=255)
    ctctGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAaga                                                                                                                                             <  2:240074/149‑1 (MQ=255)
      ctGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGaa                                                                                                                                               >  6:168382/1‑145 (MQ=255)
      ctGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGaa                                                                                                                                               <  5:168382/145‑1 (MQ=255)
       tGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTACCGGCCTCTTCGCGGTTTTAGCCGCGAaaag                                                                                                                                            <  3:123201/147‑1 (MQ=255)
        gCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGcc                                                                                                                                         >  1:48746/1‑149 (MQ=255)
         ctggTGAAAAAAGCACTGTCTCGGAGTAACAAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGAGTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGTTTTAGCCGCGAAAAGGcc                                                                                                                                         <  6:264711/146‑1 (MQ=255)
         cGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCg                                                                                                                                        <  3:248322/149‑1 (MQ=255)
               aaaaaaGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTg                                                                                                                                  >  6:129360/1‑149 (MQ=255)
                aaaaaGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTg                                                                                                                                  <  8:347582/148‑1 (MQ=255)
                      cACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATc                                                                                                                            >  2:72249/1‑148 (MQ=255)
                        cTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACTGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCaaa                                                                                                                         <  8:26458/149‑1 (MQ=255)
                        cTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCaaa                                                                                                                         <  7:183757/149‑1 (MQ=255)
                        cTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCa                                                                                                                           >  4:22105/1‑147 (MQ=255)
                           tctcGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAg                                                                                                                        <  5:129360/147‑1 (MQ=255)
                              cGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCg                                                                                                                   >  4:315340/1‑149 (MQ=255)
                               ggAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCaa                                                                                                                          <  6:221403/141‑1 (MQ=255)
                               ggAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCaa                                                                                                                          >  5:221403/1‑141 (MQ=255)
                               ggAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTc                                                                                                                    <  1:345200/147‑1 (MQ=255)
                               ggAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTc                                                                                                                    <  1:363146/147‑1 (MQ=255)
                                gAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGc                                                                                                                 <  1:72249/149‑1 (MQ=255)
                                     aCTAACGGATCCGGTACGTCGATCTTCACGGAAATGCATGAAAGGGCATTGTGGCGTTTGAAATACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTg                                                                                                            <  7:71484/149‑1 (MQ=255)
                                     aCTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCt                                                                                                             <  4:299736/148‑1 (MQ=255)
                                      cTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGc                                                                                                           <  2:48746/149‑1 (MQ=255)
                                           ccATCCGTTACGGCGCTCTGCACGGCAATGAAGTAAACGGCAGGGAGGCTATTTAAAAACGTCATGGCTGTGCGCTTGTAGCGGCCTTTTCGCGGCTAAACCCGCTAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGAt                                                                                                      <  2:182598/148‑1 (MQ=255)
                                            cATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATc                                                                                                     >  5:293514/1‑149 (MQ=255)
                                             aTCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATcc                                                                                                    >  1:149067/1‑149 (MQ=255)
                                                cGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCcgc                                                                                                  <  6:293514/148‑1 (MQ=255)
                                                   tACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGcca                                                                                              <  2:79187/149‑1 (MQ=255)
                                                   tACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGcca                                                                                              >  8:122053/1‑149 (MQ=255)
                                                         gcTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCAttt                                                                                        <  2:149067/149‑1 (MQ=255)
                                                             tGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCg                                                                                      <  5:224622/147‑1 (MQ=255)
                                                               cACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCt                                                                                    <  5:320307/147‑1 (MQ=255)
                                                               cACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTca                                                                                  >  3:66328/1‑149 (MQ=255)
                                                                     aaTGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCt                                                                            <  3:315340/149‑1 (MQ=255)
                                                                        gCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGc                                                                          >  8:92307/1‑148 (MQ=255)
                                                                          aGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAg                                                                       <  3:316794/149‑1 (MQ=255)
                                                                            gAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTAGGCATCTGAGATATCCGCGCCACCATTTCGATCACCCGCTGCAAGTc                                                                     >  6:127350/1‑149 (MQ=255)
                                                                              aaCGGCATGGATGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCa                                                                   <  7:122053/149‑1 (MQ=255)
                                                                              aaCGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTc                                                                     >  5:199787/1‑147 (MQ=255)
                                                                                 ggCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTcc                                                                    <  2:273323/145‑1 (MQ=255)
                                                                                 ggCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTcc                                                                    >  1:273323/1‑145 (MQ=255)
                                                                                        aGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTg                                                         >  1:57623/1‑149 (MQ=255)
                                                                                            gATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTc                                                      >  3:124666/1‑148 (MQ=255)
                                                                                                 aaaaaCGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCCCCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGt                                                >  2:357201/1‑149 (MQ=255)
                                                                                                 aaaaaCGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCg                                                                                      <  6:67241/111‑1 (MQ=255)
                                                                                                 aaaaaCGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCg                                                                                      >  5:67241/1‑111 (MQ=255)
                                                                                                    accggcATGGCTGTTCGGTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCtt                                             <  8:319537/144‑1 (MQ=255)
                                                                                                       gTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGaa                                          <  3:22105/149‑1 (MQ=255)
                                                                                                        tCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAAc                                         >  3:363814/1‑149 (MQ=255)
                                                                                                          aTGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACtc                                       <  4:255570/149‑1 (MQ=255)
                                                                                                          aTGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACtc                                       >  3:255570/1‑149 (MQ=255)
                                                                                                           tGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACtcg                                      >  3:325302/1‑149 (MQ=255)
                                                                                                            ggCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACtcgt                                     >  7:156213/1‑149 (MQ=255)
                                                                                                            ggCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACtcgt                                     >  3:364277/1‑149 (MQ=255)
                                                                                                             gCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACtcgtt                                    <  8:156213/149‑1 (MQ=255)
                                                                                                                   cGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAAc                                         >  1:320391/1‑138 (MQ=255)
                                                                                                                   cGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAAc                                         <  2:320391/138‑1 (MQ=255)
                                                                                                                   cGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAAc                                         <  2:325707/138‑1 (MQ=255)
                                                                                                                   cGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCAACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAAc                                         >  1:325707/1‑138 (MQ=255)
                                                                                                                    gCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACtcgttgtcgtt                              <  4:66328/148‑1 (MQ=255)
                                                                                                                     cTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCCCCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTcgttgtcgttgt                            >  3:132419/1‑149 (MQ=255)
                                                                                                                      ttGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTtgtg                           >  7:171474/1‑149 (MQ=255)
                                                                                                                        gTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTtgtg                           <  4:325302/147‑1 (MQ=255)
                                                                                                                          aGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCACGATCCGCGCAACCATTTCGCTAACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTATTGAACTCGTTTGCGTTGTGCCCg                       >  6:339312/1‑149 (MQ=255)
                                                                                                                          aGCGGCCTTTTCGCGGCTAAAACAGCGAAGAGGCCGGTAATGCGCATCAAAGCATCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCAAGTTCTTGAACTCGTGGTGGTTGTGcgcg                       >  2:355710/1‑149 (MQ=255)
                                                                                                                            cGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGCCCGgt                     <  4:364277/149‑1 (MQ=255)
                                                                                                                             ggCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGACATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTcgttgtcgttgt                            <  8:333129/141‑1 (MQ=255)
                                                                                                                             ggCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGACATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTcgttgtcgttgt                            >  7:333129/1‑141 (MQ=255)
                                                                                                                             ggCCTCTTCGCGGTTTTAGCCGCGAAAAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAAc                                                           <  3:256218/110‑1 (MQ=38)
                                                                                                                             ggCCTCTTCGCGGTTTTAGCCGCGAAAAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAAc                                                           >  4:256218/1‑110 (MQ=38)
                                                                                                                             ggCCTCTTCGCGGTTTTAGCCGCGAAAAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGCCCGgtg                    >  5:185550/1‑149 (MQ=255)
                                                                                                                               ccTCTTCGCGGTTTTAGCCGCGAAAAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGCCCGGTGTa                  >  4:142380/1‑149 (MQ=255)
                                                                                                                                cTTTTCGTGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGCCCGGTGTAc                 >  8:47783/1‑149 (MQ=255)
                                                                                                                                cTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGCCCGGTGTAc                 >  3:303756/1‑149 (MQ=255)
                                                                                                                                    tCGCGTCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGCCCGGTGTACTTc              <  4:132419/148‑1 (MQ=255)
                                                                                                                                    tCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGCCCGGTGTACTTc              <  6:199787/148‑1 (MQ=255)
                                                                                                                                      gcgGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGCCCGGTGTACTTCAc            >  1:155136/1‑148 (MQ=255)
                                                                                                                                       cgGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGCCCGGTGTACTTCACCt          >  8:57469/1‑149 (MQ=255)
                                                                                                                                              aaCCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGCCCGGGGTACTTCACCTCCTGcct   >  3:253640/1‑147 (MQ=255)
                                                                                                                                               aCCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCATGAACTCGTTGTCGTTGTGCCCGGTGTACTTCACCTCCGGCATg  >  3:361031/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
GGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGCCCGGTGTACTTCACCTCCGGCATG  >  NC_002947/253218‑253509

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: