Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,212,540 C→T 6.0% intergenic (+56/+38) PP_4593 → / ← PP_4594 hypothetical protein/cystathionine gamma‑synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,212,5400CT 6.0% 306.0 / 6.0 106intergenic (+56/+38)PP_4593/PP_4594hypothetical protein/cystathionine gamma‑synthase
Reads supporting (aligned to +/- strand):  ref base C (44/51);  new base T (3/3);  total (49/57)
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.54e-01

Only 100 of 109 total aligned reads displayed.

GGCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGCTGATCCCCAT  >  NC_002947/5212392‑5212682
                                                                                                                                                     |                                                                                                                                               
ggCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGCCAAGCCCGCTCCTAC‑ACcgacg                                                                                                                                                >  2:372170/1‑150 (MQ=255)
   aCCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAg                                                                                                                                             <  1:1435889/150‑1 (MQ=255)
   aCCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAg                                                                                                                                             <  2:325585/150‑1 (MQ=255)
         gacgaAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑CCCGTCGTAGGAGCgg                                                                                                                                       <  2:954509/150‑1 (MQ=255)
         gacgaAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTa                                                                                                                                              >  1:1193538/1‑143 (MQ=255)
         gacgaAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCGTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTa                                                                                                                                              <  2:1193538/143‑1 (MQ=255)
          acgaAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACGACGGT‑GTAGGAGCggg                                                                                                                                      <  2:276688/150‑1 (MQ=255)
          acgaAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑GACGGTGTAGGAGCggg                                                                                                                                      <  1:131859/150‑1 (MQ=255)
          acgaAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCggg                                                                                                                                      <  2:1028833/150‑1 (MQ=255)
           cgaAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGt                                                                                                                                     <  2:1745617/150‑1 (MQ=255)
               gCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑c                                                                                                                                <  2:1937591/150‑1 (MQ=255)
               gCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACAGCCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑c                                                                                                                                >  1:610365/1‑150 (MQ=255)
                 cGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑ccc                                                                                                                              >  2:81615/1‑150 (MQ=255)
                 cGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑ccc                                                                                                                              >  2:1128089/1‑150 (MQ=255)
                     aaCCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑c                                                                                                                                >  2:1337641/1‑144 (MQ=255)
                     aaCCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑c                                                                                                                                <  1:1337641/144‑1 (MQ=255)
                         gCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAgg                                                                                                                      <  1:467151/150‑1 (MQ=255)
                          ccTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAggg                                                                                                                     <  2:123543/150‑1 (MQ=255)
                           cTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGc                                                                                                                    <  1:1300005/150‑1 (MQ=255)
                               gcAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGc                                                                                                                >  1:958719/1‑150 (MQ=255)
                                cAAAGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCg                                                                                                               <  2:370016/150‑1 (MQ=255)
                                 aaGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCgg                                                                                                              >  2:303441/1‑150 (MQ=255)
                                  aGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAg                                                                                                                                             >  2:1731103/1‑119 (MQ=255)
                                  aGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAg                                                                                                                                             <  1:1731103/119‑1 (MQ=255)
                                    gAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGccc                                                                                                           >  2:683517/1‑150 (MQ=255)
                                    gAAAAGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGccc                                                                                                           >  1:279817/1‑150 (MQ=255)
                                           tGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTcacaca                                                                                                    <  1:1128089/150‑1 (MQ=255)
                                           tGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTcacaca                                                                                                    >  2:1587234/1‑150 (MQ=255)
                                            gCTGATCAAGGACAAGATCGAACGGATCATGCATCCGCCCAAATCCTGACAGGCATATATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGGTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAg                                                                                                   >  2:201179/1‑150 (MQ=255)
                                              tGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAGCg                                                                                                 <  2:87045/150‑1 (MQ=255)
                                                 tCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑ccc                                                                                                                              >  2:24788/1‑118 (MQ=255)
                                                 tCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑ccc                                                                                                                              <  1:24788/118‑1 (MQ=255)
                                                     ggACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGTAAACCCGCTCCTACGACGGT‑GTAGGAGCGGGCTTG‑CCcgc                                                                                                                            >  2:170140/1‑116 (MQ=21)
                                                     ggACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGTAAACCCGCTCCTACGACGGT‑GTAGGAGCGGGCTTG‑CCcgc                                                                                                                            <  1:170140/116‑1 (MQ=21)
                                                     ggACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGTAAACCCGCTCCTACGACGGT‑GTAGGAGCGGGCTTG‑CCcgc                                                                                                                            <  1:169801/116‑1 (MQ=21)
                                                     ggACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGTAAACCCGCTCCTAC‑GACGGTGTAGGAGCGGGCTTG‑CCcgc                                                                                                                            >  2:169801/1‑116 (MQ=21)
                                                      gACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAgcgc                                                                                         <  2:769233/150‑1 (MQ=255)
                                                      gACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCCTTTCGCGGGTAAACCCGCTCCTAC‑GACGGTGTAGGAGCGGGCTTG‑CCCGCGAAGAGGccggtg                                                                                                               >  1:1663152/1‑124 (MQ=14)
                                                         aaGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCttt                                                                                      >  2:378980/1‑150 (MQ=255)
                                                         aaGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCttt                                                                                      <  1:1381801/150‑1 (MQ=255)
                                                             tCGAACGGATCATTCATCCGCCCAAATACTGACAGGCATCTATTGAATCCACCGGACTCTTCGCGGGAAAGCCCGCTCCTAC‑ACCGTCGTAGAAGCGGGCTTA‑CCCGCGAAAGGGCCTGCAGCCCTCACACAGCGACCACCGCTTTCAAc                                                                                  >  2:841396/1‑150 (MQ=255)
                                                               gAACGGATCATTCATCCGCCCAAATCCTGACCGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGc                                                                                <  2:610365/150‑1 (MQ=255)
                                                                  cGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTca                                                                                                        >  2:807239/1‑123 (MQ=255)
                                                                  cGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTca                                                                                                        <  1:807239/123‑1 (MQ=255)
                                                                   ggATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAGCGAc                                                                                               <  2:1161614/131‑1 (MQ=255)
                                                                   ggATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAGCGAc                                                                                               >  1:1161614/1‑131 (MQ=255)
                                                                     aTCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATc                                                                          >  1:705535/1‑150 (MQ=255)
                                                                     aTCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATc                                                                          <  1:1390019/150‑1 (MQ=255)
                                                                        attcatCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGTAAACCCGCTCCTACGACGGT‑GTAGGAGCgggc                                                                                                                                     >  2:1811389/1‑88 (MQ=18)
                                                                        attcatCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGTAAACCCGCTCCTAC‑GACGGTGTAGGAGCgggc                                                                                                                                     <  1:1811389/89‑2 (MQ=18)
                                                                        attcatCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAg                                                                       <  1:1659359/150‑1 (MQ=255)
                                                                        attcatCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAg                                                                       <  2:1573131/150‑1 (MQ=255)
                                                                             tCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCt                                                                  >  1:967338/1‑150 (MQ=255)
                                                                             tCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCt                                                                  >  2:1023409/1‑150 (MQ=255)
                                                                             tCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCt                                                                  >  2:1227948/1‑150 (MQ=255)
                                                                               cGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGc                                                                                                                    <  1:606444/98‑1 (MQ=255)
                                                                               cGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGc                                                                                                                    >  2:606444/1‑98 (MQ=255)
                                                                                 cccAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTTGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAGc                                                                                                  <  2:1629821/114‑1 (MQ=255)
                                                                                 cccAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTTGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAGc                                                                                                  >  1:1629821/1‑114 (MQ=255)
                                                                                    aaaTCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAg                                                           <  1:1514318/150‑1 (MQ=255)
                                                                                      aTCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAg                                                           <  2:1604535/148‑1 (MQ=255)
                                                                                      aTCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAg                                                           >  1:1604535/1‑148 (MQ=255)
                                                                                          tGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACg                                                     <  1:1167247/150‑1 (MQ=255)
                                                                                              aGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCt                                                                  >  2:3131/1‑133 (MQ=255)
                                                                                              aGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCt                                                                  <  1:3131/133‑1 (MQ=255)
                                                                                                gCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATc                                               >  2:1132638/1‑150 (MQ=255)
                                                                                                gCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATc                                               >  1:739756/1‑150 (MQ=255)
                                                                                                gCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATc                                               >  1:1488002/1‑150 (MQ=255)
                                                                                                gCATCTATTGAATCCACCGGCCTCTTCGCGGACAAGGCCGCTCCTAC‑ACCGTCGTAGGATCGGGGTTA‑CGCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGTCCAGGTCAGGGTAATc                                               >  2:1647258/1‑150 (MQ=255)
                                                                                                 cATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATcc                                              >  1:1655895/1‑150 (MQ=255)
                                                                                                  aTCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑ccc                                                                                                                              >  1:218568/1‑69 (MQ=255)
                                                                                                  aTCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑ccc                                                                                                                              <  2:218568/69‑1 (MQ=255)
                                                                                                  aTCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCt                                             <  2:1024865/150‑1 (MQ=255)
                                                                                                     tATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGa                                          <  1:81615/150‑1 (MQ=255)
                                                                                                      aTTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGAt                                         >  2:287172/1‑150 (MQ=255)
                                                                                                       ttGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACTCTCTTAGCGGGTCGTACAGCGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATc                                        <  2:958719/150‑1 (MQ=25)
                                                                                                          aaTCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCt                                     >  1:384398/1‑150 (MQ=255)
                                                                                                          aaTCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCAGACTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCt                                     >  2:957628/1‑150 (MQ=255)
                                                                                                            tCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTAc                                   <  1:1023409/150‑1 (MQ=255)
                                                                                                              cACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCgg                                                                                                                                       <  1:978955/49‑1 (MQ=255)
                                                                                                                      tctTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCTGGTTTA‑CCCGCGAAAGGGCCTGCTGCCCTCACACAGCGACCATCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCgc                         >  2:882179/1‑150 (MQ=255)
                                                                                                                        ttCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTc                                                                 >  1:458933/1‑108 (MQ=255)
                                                                                                                        ttCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTc                                                                 <  2:458933/108‑1 (MQ=255)
                                                                                                                        ttCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCgcag                       >  1:451515/1‑150 (MQ=255)
                                                                                                                             gggCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAGCGAc                                                                                               >  2:1706004/1‑73 (MQ=255)
                                                                                                                             gggCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAGCGAc                                                                                               <  1:1706004/73‑1 (MQ=255)
                                                                                                                             gggCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATcc                                                                         >  2:1263535/1‑95 (MQ=255)
                                                                                                                             gggCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATcc                                                                         <  1:1263535/95‑1 (MQ=255)
                                                                                                                                   gCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAg                                                                       <  2:102496/91‑1 (MQ=255)
                                                                                                                                   gCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAg                                                                       >  1:102496/1‑91 (MQ=255)
                                                                                                                                   gCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAAcc                           <  1:452370/135‑1 (MQ=255)
                                                                                                                                   gCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAAcc                           >  2:452370/1‑135 (MQ=255)
                                                                                                                                   gCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGc            <  2:705535/150‑1 (MQ=255)
                                                                                                                                   gCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGc            >  1:1216102/1‑150 (MQ=255)
                                                                                                                                       gCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGCTGAt        >  1:213246/1‑150 (MQ=255)
                                                                                                                                         tCCTAC‑ACCGTCGTAGGGGCGGGTTTACCCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGCTGATc       <  1:287263/150‑1 (MQ=255)
                                                                                                                                          ccTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGCTGATccc     <  1:900952/150‑1 (MQ=255)
                                                                                                                                             ac‑acCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGCTGATCCCCAt  <  1:1132638/150‑1 (MQ=255)
                                                                                                                                                   gtcgtAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGa                                          >  1:1324976/1‑105 (MQ=255)
                                                                                                                                                   gtcgtAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGa                                          <  2:1324976/105‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
GGCACCAACGACGAAGCCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAC‑ACCGTCGTAGGAGCGGGTTTA‑CCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCGACAACCGCAGCAGGTTCTCGCTGATCCCCAT  >  NC_002947/5212392‑5212682

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