Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 253,356 C→T 7.3% intergenic (‑141/+30) 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,3560CT 7.3% 241.6 / 5.2 81intergenic (‑141/+30)PP_0202/PP_0203CBS domain‑containing protein/dipeptidase
Reads supporting (aligned to +/- strand):  ref base C (37/38);  new base T (3/3);  total (40/41)
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.52e-01

GTCTTCATACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGCCCGGTGTACTTCACCT  >  NC_002947/253208‑253501
                                                                                                                                                    |                                                                                                                                                 
gTCTTCATACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGc                                                                                                                                                   <  4:217054/149‑1 (MQ=255)
     cATACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTaa                                                                                                                                                <  7:170666/147‑1 (MQ=255)
          gggCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACcgc                                                                                                                                          >  7:347582/1‑148 (MQ=255)
              ctctGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGaa                                                                                                                                       >  1:79187/1‑147 (MQ=255)
              ctctGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAaga                                                                                                                                     <  2:240074/149‑1 (MQ=255)
                ctGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGaa                                                                                                                                       <  5:168382/145‑1 (MQ=255)
                ctGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGaa                                                                                                                                       >  6:168382/1‑145 (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:363146/147‑1 (MQ=255)
                                         ggAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTc                                                                                                            <  1:345200/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                                                            >  1:273323/1‑145 (MQ=255)
                                                                                           ggCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTcc                                                            <  2:273323/145‑1 (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                             >  3:364277/1‑149 (MQ=255)
                                                                                                                      ggCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACtcgt                             >  7:156213/1‑149 (MQ=255)
                                                                                                                       gCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACtcgtt                            <  8:156213/149‑1 (MQ=255)
                                                                                                                             cGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAAc                                 >  1:320391/1‑138 (MQ=255)
                                                                                                                             cGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAAc                                 <  2:325707/138‑1 (MQ=255)
                                                                                                                             cGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAAc                                 <  2:320391/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                    >  7:333129/1‑141 (MQ=255)
                                                                                                                                       ggCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGACATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTcgttgtcgttgt                    <  8:333129/141‑1 (MQ=255)
                                                                                                                                       ggCCTCTTCGCGGTTTTAGCCGCGAAAAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAAc                                                   >  4:256218/1‑110 (MQ=38)
                                                                                                                                       ggCCTCTTCGCGGTTTTAGCCGCGAAAAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAAc                                                   <  3:256218/110‑1 (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)
                                                                                                                                                    |                                                                                                                                                 
GTCTTCATACGGGCCTCTGCGGGTGAAAAAAGCACTGTCTCGGAGTAACTAAAGCATCCGGTACGGCGCTCTGCACGGCAATGCAGGAAACGGCATGGAGGCGATTGAAAAACGTCATGGCTGTTCGCTTGTAGCGGCCTTTTCGCGGCTAAAACCGCGAAGAGGCCGGTACTGCGCATCAAAGCTTCGGCATCTGCCCGATCCGCGCCACCATTTCGCTCACCACCTGCAAGTCCAGCATGAACTGCTCCACGGTCTTGAACTCGTTGTCGTTGTGCCCGGTGTACTTCACCT  >  NC_002947/253208‑253501

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: