Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,039,004 A→G 9.7% L159S (TTG→TCG)  rffE ← UDP‑N‑acetylglucosamine‑2‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,039,0040AG 9.7% 196.9 / 10.1 72L159S (TTG→TCG) rffEUDP‑N‑acetylglucosamine‑2‑epimerase
Reads supporting (aligned to +/- strand):  ref base A (30/35);  new base G (4/3);  total (34/38)
Fisher's exact test for biased strand distribution p-value = 7.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.35e-01

GCCGGCGCTCGGCTGCCAGGTAGCTGAAGTGGCTGGCGAAACGGGCCTTCAGCGCGGGGCTTTGCAGCTTTTGCTCCACCTCCAGCAAGGCGTCGATCACGGTGTTACCGGTGGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGCCCCGGTCAGGCGCCGGTTGCCTTCTTCGGGCCAGGGTGAATACAGGTTGTTGGTGCGCAGGCCGGCTTCTACATGGGCAATGGGGATTTGCTGG  >  NC_002947/2038858‑2039150
                                                                                                                                                  |                                                                                                                                                  
gCCGGCGCTCGGCTGCCAGGTAGCTGAAGTGGCTGGCGAAACGGGCCTTCAGCGCGGGGCTTTGCAGCTTTTGCTCCACCTCCAGCAAGGCGTCGATCACGGTGTTACCGGTGGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAAc                                                                                                                                                  <  4:312733/149‑1 (MQ=255)
  cGGCGCTCGGCTGCCAGGTAGCTGAAGTGGCTGGCGAAACGGGCCTTCAGCGCGGGGCTTTGCAGCTTTTGCTCCACCTCCAGCAAGGCGTCGATCACGGTGTTACCGGTGGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAg                                                                                                                                                >  4:58175/1‑149 (MQ=255)
    ccgcTCGGCTGCCAGGTAGCTGAAGTGGCTGGCGAAACGGGCCTTCAGCGCGGGGCTTTGCAGCTTTTGCTCCACCTCCAGCAAGGCGTCGATCACGGAGTTACCGGAGGTGAAGATGTGCGCCGCATCGGTGCCTTCCCCCAACAGGt                                                                                                                                              >  8:122531/2‑149 (MQ=255)
        tCGGCTGCCAGGTAGCTGAAGTGGCTGGCGAAACGGGCCTTCAGCGCGGGGCTTTGCAGCTTTTGCTCCACCTCCAGCAAGGCGTCGATCACGGTGTTACCGGTGGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAg                                                                                                                                                <  2:338761/143‑1 (MQ=255)
        tCGGCTGCCAGGTAGCTGAAGTGGCTGGCGAAACGGGCCTTCAGCGCGGGGCTTTGCAGCTTTTGCTCCACCTCCAGCAAGGCGTCGATCACGGTGTTACCGGTGGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAg                                                                                                                                                >  1:338761/1‑143 (MQ=255)
            gtGCCAGTTATCTTAAGTTGCTGGCGAAACGGGACTTCAGCGCGGGGCTTTGCAGCTTTTGCTCCACCTCCATCAAGGCGTCGATCACGGTGTTACCGGTGGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTcgcgc                                                                                                                                      <  7:183172/148‑1 (MQ=255)
            cTGCCAGGTAGCTGAAGTGGCTGGCGAAACGGGCCTTCAGCGCGGGGCTTTGCAGCTTTTGCTCCACCTCCAGCAAGGCGTCGATCACGGTGTTACCGGTGGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTcgcgc                                                                                                                                      <  8:296128/149‑1 (MQ=255)
            cTGCCAGGTAGCTGAAGTGGCTGGCGAAACGGGCCTTCAGCGCGGGGCTTTGCAGCTTTTGCTCCACCTCCAGCAAGGCGTCGATCACGGTGTTACCGGTGGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTcgcgc                                                                                                                                      >  7:296128/1‑149 (MQ=255)
             tGCCAGGTAGCTGAAGTGGCTGGCGAAACGGGCCTTCAGCGCGGGGCTTTGCAGCTTTTGCTCCACCTCCAGCAAGGCGTCGATCACGGTGTTACCGGTGGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTcgcgc                                                                                                                                      >  6:138504/1‑148 (MQ=255)
               ccAGGTAGCTGAAGTGGCTGGCGAAACGGGCCTTCAGCGCGGGGCTTTGCAGCTTTTGCTCCACCTCCAGCAAGGCGTCGATCACGGTGTTACCGGTGGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAg                                                                                                                                   <  1:139368/149‑1 (MQ=255)
                cAGGTAGCTGAAGTGGCTGGCGAAACGGGCCTTCAGCGCGGGGCTTTGCAGCTTTTGCTCCACCTCCAGCAAGGCGTCGATCACGGTGTTACCGGTGGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCGACAACCTGTTGCGCGaag                                                                                                                                  <  6:92051/149‑3 (MQ=255)
                 aGGTAGCTGAAGTGGCTGGCGAAACGGGCCTTCAGCGCGGGGCTTTGCAGCTTTTGCTCCACCTCCAGCAAGGCGTCGATCACGGTGTTACCGGTGGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGt                                                                                                                                 <  2:35232/149‑1 (MQ=255)
                       cTGAAGTGGCTGGCGAAACGGGCCTTCAGCGCGGGGCTTTGCAGCTTTTGCTCCACCTCCAGCAAGGCGTCGATCACGGTGTTACCGGTGGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGt                                                                                                                           <  4:17895/149‑1 (MQ=255)
                        tGAAGTGGCTGGCGAAACGGGCCTTCAGCGCGGGGCTTTGCAGCTTTTGCTCCACCTCCAGCAAGGCGTCGATCACGGTGTTACCGGTGGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTg                                                                                                                          >  3:46361/1‑149 (MQ=255)
                        tGAAGTGGCTGGCGAAACGGGCCTTCAGCGCGGGGCTTTGCAGCTTTTGCTCCACCTCCAGCAAGGCGTCGATCACGGTGTTACCGGTGGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTg                                                                                                                          >  6:290738/1‑149 (MQ=255)
                         gAAGTGGCTGGCGAAACGGGCCTTCAGCGCGGGGCTTTGCAGCTTTTGCTCCACCTCCAGCAAGGCGTCGATCACGGTGTTACCGGTGGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTgg                                                                                                                         >  1:281742/1‑149 (MQ=255)
                          aaGTGGCTGGCGAAACGGGCCTTCAGCGCGGGGCTTTGCAGCTTTTGCTCCACCTCCAGCAAGGCGTCGATCACGGTGTTACCGGTGGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTg                                                                                                                          <  5:186584/147‑1 (MQ=255)
                          aaGTGGCTGGCGAAACGGGCCTTCAGCGCGGGGCTCTGCAGCTTTTGCTCCACCTCCAGCAAGGCGTCGATCACGGTGTTACCGGTTGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGt                                                                                                                           <  7:46032/146‑1 (MQ=255)
                              ggctggcGAAACGGGCCTTCAGCGCGGGGCTTTGCAGCTTTTGCTCCACCTCCAGCAAGGCGTCGATCACGGTGTTACCGGTGGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTGGgcgc                                                                                                                     >  1:147974/1‑148 (MQ=255)
                                ctggcGAAACGGGCCTTCAGCGCGGGGCTTTGCAGCTTTTGCTCCACCTCCAGCAAGGCGTCGATCACGGTGTTACCGGTGGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTGGGcgcg                                                                                                                    >  3:173738/1‑147 (MQ=255)
                                ctggcGAAACGGGCCTTCAGCGCGGGGCTTTGCAGCTTTTGCTCCACCTCCAGCAAGGCGTCGATCACGGTGTTACCGGTGGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTGGGCGCGaa                                                                                                                  <  2:239334/149‑1 (MQ=255)
                                   gcGAAACGGGCCTTCAGCGCGGGGCTTTGCAGCTTTTGCTCCACCTCCAGCAAGGCGTCGATCACGGTGTTACCGGTGGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTGGGCGCGAAGTg                                                                                                               <  3:58175/149‑1 (MQ=255)
                                               ttCAGCGCGGGGCTTTGCAGCTTTTGCTCCACCTCCAGCAAGGCGTCGATCACGGTGTTACCGGTGGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAg                                                                                                   >  7:239691/1‑149 (MQ=255)
                                                  aGCGCGGGGCTTTGCAGCTTTTGCTCCACCTCCAGCAAGGCGTCGATCACGGTGTTACCGGTGGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGc                                                                                                <  4:173738/149‑1 (MQ=255)
                                                            tttGCAGCTTTTGCTCCACCTCCAGCAAGGCGTCGATCACGGTGTTACCGGTGGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGCCCCGGTCAgg                                                                                      >  7:204515/1‑149 (MQ=255)
                                                              tGCAGCTTTTGCTCCACCTCCAGCAAGGCGTCGATCACGGTGTTACCGGTGGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGCCCCGGTCAGgcg                                                                                    >  1:310471/1‑149 (MQ=255)
                                                                      ttGCTCCACCTCCAGCAAGGCGTCGATCACGGTGTTACCGGTGGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGCCCCGGTCAGGCGCCGGtt                                                                              <  7:310946/147‑1 (MQ=255)
                                                                      ttGCTCCACCTCCAGCAAGGCGTCGATCACGGTGTTACCGGTGGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGCCCCGGTCAGGCGCCGGTTg                                                                             <  2:281742/148‑1 (MQ=255)
                                                                        gCTCCACCTCCAGCAAGGCGTCGATCACGGTGTTACCGGTGGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGCCCCGGTCAGGCGCCGGTTGCct                                                                          >  3:100734/1‑149 (MQ=255)
                                                                            cACCTCCAGCAAGGCGTCGATCACGGTGTTACCGGTGGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCTGTGGGCGCGATGTGCAGCGCGGCAAGGGCCCCGGTCAGGCGCCGGTTGCcttctt                                                                      >  2:139217/1‑149 (MQ=255)
                                                                            cACCTCCAGCAAGGCGTCGATCACGGTGTTACCGGTGGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGCCCCGGTCAGGCGCCGGTTGCcttctt                                                                      <  2:294781/149‑1 (MQ=255)
                                                                            cACCTCCAGCAAGGCGTCGATCACGGTGTTACCGGTGGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGCCCCGGTCAGGCGCCGGTTGCcttctt                                                                      >  2:140170/1‑149 (MQ=255)
                                                                             aCCTCCAGCAAGGCGTCGATCACGGTGTTACCGGTGGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGCCCCGGTCAGGCGCCGGTTGCCttcttc                                                                     >  6:228776/1‑149 (MQ=255)
                                                                               cTCCAGCAAGGCGTCGATCACGGTGTTACCGGTGGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGCCCCGGTCAGGCGCCGGTTGCCTTCTTCgg                                                                   >  5:228678/1‑149 (MQ=255)
                                                                                  cAGCAAGGCGTCGATCACGGTGTTACCGGTGGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGCCCCGGTCAGGCGCCGGTTGCCTTCTTCggg                                                                  <  4:46361/147‑1 (MQ=255)
                                                                                    gCAAGGCGTCGATCACGGTGTTACCGGTGGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGCCCCGGTCAGGCGCCGGTTGCCTTCTTCGGGCCAg                                                              <  2:310471/149‑1 (MQ=255)
                                                                                       aGGCGTCGATCACGGTGTTACCGGTGGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGCCCCGGTCAGGCGCCGGTTGCCTTCTTCGGGCCAGGGt                                                           <  5:138504/149‑1 (MQ=255)
                                                                                        ggCGTCGATCACGGTGTTACCGGTGGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGCCCCGGTCAGGCGCCGGTTGCCTTCTTCGGGCCAGGGTg                                                          <  2:147974/149‑1 (MQ=255)
                                                                                          cGTCGATCACGGTGTTACCGGTGGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGCCCCGGTCAGGCGCCGGTTGCCTTCTTCGGGCCAGGGTGaa                                                        <  5:228776/149‑1 (MQ=255)
                                                                                               aTCACGGTGTTACCGGTGGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGCCCCGGTCAGGCGCCGGTTGCCTTCTTCg                                                                    >  2:283301/1‑132 (MQ=255)
                                                                                               aTCACGGTGTTACCGGTGGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGCCCCGGTCAGGCGCCGGTTGCCTTCTTCg                                                                    <  1:283301/132‑1 (MQ=255)
                                                                                                      tgtTACCGGTGGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTGGGCGCGAAGTGCTGCGCGTCAAGGGCCCCGGTCAGGCGCCGGTTGCCTTCTTCGGGCCAGGGTGAATACAGGttgttg                                            <  5:277304/149‑1 (MQ=255)
                                                                                                      tgtTACCGGTGGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGCCCCGGTCAGGCGCCGGTTGCCTTCTTCGGGCCAGGGTGAATACAGGttgttg                                            <  5:290738/149‑1 (MQ=255)
                                                                                                      tgtTACCGGTGGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGCCCCGGTCAGGCGCCGGTTGCCTTCTTCGGGCCAGGGTGAATACAGGttgttg                                            <  2:319545/149‑1 (MQ=255)
                                                                                                      tgtTACCGGTGGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGCCCCGGTCAGGCGCCGGTTGCCTTCTTCCGGCCAGGGTGAATACAGGttgttg                                            <  7:48249/149‑1 (MQ=255)
                                                                                                        ttACCGGTGGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGCCCCGGTCAGGCGCCGGTTGCCTTCTTCGGGCCAGGGTGAATACAGGTTGTTGGt                                          >  1:300054/1‑149 (MQ=255)
                                                                                                           ccGGTGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGCCCCGGTCAGGCGCCGGTTGCCTTCTTCGGGCCAGGGTGAATACAGGTTGTTGGTgcg                                       >  2:212811/1‑148 (MQ=255)
                                                                                                            cGGTGGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGCCCCGGTCAGGCGCCGGTTGCCTTCTTCGGGCCAGGGTGAATACAGGTTGTTGGTgcg                                       >  3:246391/1‑148 (MQ=255)
                                                                                                                 gtgAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGCCCCGGTCAGGCGCCGGTTGCCTTCTTCGGGCCAGGGTGAATACAGGTTGTTGGTGCGCAGGcc                                 >  6:333993/1‑149 (MQ=255)
                                                                                                                  tgAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGAGCGGCAAGGGCCCCGGTCAGGCGCCGGTAGCCTTCTTCGGGCCAGGGTGAATACAGGTTGTTGGTGCGCAGGCCg                                >  8:131919/1‑149 (MQ=255)
                                                                                                                   gAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGCCCCGGTCAGGCGCCGGTTGCCTTCTTCGGGCCAGGGTGAATACAGGTTGTTGGTGCGCAGGCCgg                               >  7:268891/1‑149 (MQ=255)
                                                                                                                    aaGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGCCCCGGTCAGGCGCCGGTTGCCTTCTTCGGGCCAGGGTGAATACAGGTTGTTGGTGCGCAGGCCGGc                              <  2:114638/149‑1 (MQ=255)
                                                                                                                    aaGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGCCCCGGTCAGGCGCCGGTTGCCTGCTTCGGGCCAGGGTGAATACAGGTTGTTGGTGCGCAGGCCGGc                              >  5:41669/1‑149 (MQ=255)
                                                                                                                      gttgtgCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTCGGGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGCCCCGGTCAGGCGCCGGTTGCCTTCTTCGGGCCAGGGTGAATACAGGTTGTTGGTGCGCAGGCCGGCtt                            <  7:122531/147‑1 (MQ=255)
                                                                                                                           gcgcCGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGCCCCGGTCAGGCGCCGGTTGCCTTCTTCGGGCCAGGGTGAATACAGGTTGTTGGTGCGCAGGCCGGCTTCTACa                       >  6:89100/1‑149 (MQ=255)
                                                                                                                           gcgcCGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGCCCCGGTCAGGCGCCGGTTGCCATCTTCGGGCCAGGGTGAATACAGGTTGTTGGTGCGCAGGCCGGCTTCTACa                       >  6:286160/1‑149 (MQ=255)
                                                                                                                              ccGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGCCCCGGTCAGGCGCCGGTTGCCTTCTTCGGGCCAGGGTGAATACAGGTTGTTGGTGCGCAGGCCGGCTTCTACATgg                    >  6:83421/1‑149 (MQ=255)
                                                                                                                              ccGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGCCCCGGTCAGGCGCCGGTTGCCTTCTTCGGGCCAGGGTGAATACAGGTTGTTGGTGCGCAGGCCGGCTTCTACATgg                    <  5:89100/149‑1 (MQ=255)
                                                                                                                              ccGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGCCCCGGTCAGGCGCCGGTTGCCTTCTTCGGGCCAGGGTGAATACAGGTTGTTGGTGCGCAGGCCGGCTTCTACATgg                    <  8:204515/149‑1 (MQ=255)
                                                                                                                              ccGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGCCCCGGTCAGGCGCCGGTTGCCTTCTTCGGGCCAGGGTGAATACAGGTTGTTGGTGCGCAGGCCGGCTTCTACATgg                    <  1:212811/149‑1 (MQ=255)
                                                                                                                              ccGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGCCCCGGTCAGGCGCCGGTTGCCATCTTCGGGCCAGGGTGAATACAGGTTGTTGGTGCGCAGGCCGGCTTCTACATgg                    <  5:286160/149‑1 (MQ=255)
                                                                                                                                     ggTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGCCCCGGTCAGGCGCCGGTTGCCTTCTTCGGGCCAGGGTGAATACAGGTTGTTGGTGCGCAGGCCGGCTTCTACATGGGCAATgg             <  8:239691/149‑1 (MQ=255)
                                                                                                                                     ggTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGCCCCGGTCAGGCGCCGGTTGCCTTCTTCGGGCCAGGGTGAATACAGGTTGTTGGTGCGCAGGCCGGCGTCTACATGGGCAATg              <  4:100734/148‑1 (MQ=255)
                                                                                                                                       ccacctCGCGCGACAACCTGTTGCGCGAAGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGCCCCGGTCAGGCGCCGGTTGCCTTCTTCGGGCCAGGGTGAATACAGGTTGTTGGTGCGCAGGCCGGCTTCTACATGGGCAATggg            >  8:17073/6‑148 (MQ=255)
                                                                                                                                        cacctCGCGCGACAACCTGTTGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGCCCCGGTCAGGCGCCGGTTGCCTTCTTCGGGCCAGGGTGAATACAGGTTGTTGGTGCGCAGGCCGGCTTCTACATGGGCAATGGGGa          >  7:66348/5‑149 (MQ=255)
                                                                                                                                         acctCGCGCGACAACCTGTTGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCTGCAAGGGCCCCGGTCAGGCGCCGGtt                                                                              <  2:220519/77‑1 (MQ=25)
                                                                                                                                         acctCGCGCGACAACCTGTTGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGCCCCGGTCAGGCGCCGGTTGCCTTCTTCGGGCCAGGGTGAATACAg                                                   <  1:303455/104‑1 (MQ=255)
                                                                                                                                         acctCGCGCGACAACCTGTTGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGCCCCGGTCAGGCGCCGGTTGCCTTCTTCGGGCCAGGGTGAATACAg                                                   >  2:303455/4‑107 (MQ=255)
                                                                                                                                           ctCGCGCGACAACCTGTTGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGCCCCGGTCAGGCGCCGGTTGCCTTCTTCGGGCCAGGGTGAATACAGGTTGTTGGTGCGCAGGCCGGCTTCTACATGGGCAATGGGGAttt       >  1:63821/2‑149 (MQ=255)
                                                                                                                                              gcgcAACAGGTTGTCGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGCCCCGGTCAGGCGCCGGTTGCCTTCTTCGGGCCAGGGTGAATACAGGTTGTTGGTGCGCAGGCCGGCTTCTACATGGGCAATGGGGATTTGCt    <  5:333993/149‑1 (MQ=255)
                                                                                                                                                gcAACAGGTTGTCGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGCCCCGGTCAGGCGCCGGTTGCCTTCTTCGGGCCAGGGTGAATACAGGTTGTTGGTGCGCAGGCCGGCTTCTACATGGGCAATGGGGATTTGCTgg  <  5:83421/149‑1 (MQ=255)
                                                                                                                                                 cAACAGGTTGTCGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGCCCCGGTCAGGCGCCGGTTGCCTTCTTCGGGCCAGGGTGAATACAGGTTGTTGGTGCGCAGGCCGGCTTCTACATGGGCAATGGGGATTTGCTgg  >  4:125795/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GCCGGCGCTCGGCTGCCAGGTAGCTGAAGTGGCTGGCGAAACGGGCCTTCAGCGCGGGGCTTTGCAGCTTTTGCTCCACCTCCAGCAAGGCGTCGATCACGGTGTTACCGGTGGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGCCCCGGTCAGGCGCCGGTTGCCTTCTTCGGGCCAGGGTGAATACAGGTTGTTGGTGCGCAGGCCGGCTTCTACATGGGCAATGGGGATTTGCTGG  >  NC_002947/2038858‑2039150

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