Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,039,008 G→A 10.3% L158L (CTG→TTG)  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,0080GA10.3% 184.7 / 9.3 68L158L (CTG→TTG) rffEUDP‑N‑acetylglucosamine‑2‑epimerase
Reads supporting (aligned to +/- strand):  ref base G (28/33);  new base A (4/3);  total (32/36)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

CGGCGCTCGGCTGCCAGGTAGCTGAAGTGGCTGGCGAAACGGGCCTTCAGCGCGGGGCTTTGCAGCTTTTGCTCCACCTCCAGCAAGGCGTCGATCACGGTGTTACCGGTGGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGCCCCGGTCAGGCGCCGGTTGCCTTCTTCGGGCCAGGGTGAATACAGGTTGTTGGTGCGCAGGCCGGCTTCTACATGGGCAATGGGGATTTGCTGG  >  NC_002947/2038860‑2039150
                                                                                                                                                    |                                                                                                                                              
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:140170/1‑149 (MQ=255)
                                                                          cACCTCCAGCAAGGCGTCGATCACGGTGTTACCGGTGGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGCCCCGGTCAGGCGCCGGTTGCcttctt                                                                      <  2:294781/149‑1 (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                                                                    <  1:283301/132‑1 (MQ=255)
                                                                                             aTCACGGTGTTACCGGTGGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGCCCCGGTCAGGCGCCGGTTGCCTTCTTCg                                                                    >  2:283301/1‑132 (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)
                                                                                                         ccGGT‑GTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGCCCCGGTCAGGCGCCGGTTGCCTTCTTCGGGCCAGGGTGAATACAGGTTGTTGGTgcg                                       >  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                    <  5:89100/149‑1 (MQ=255)
                                                                                                                            ccGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGCCCCGGTCAGGCGCCGGTTGCCTTCTTCGGGCCAGGGTGAATACAGGTTGTTGGTGCGCAGGCCGGCTTCTACATgg                    >  6:83421/1‑149 (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)
                                                                                                                                                    |                                                                                                                                              
CGGCGCTCGGCTGCCAGGTAGCTGAAGTGGCTGGCGAAACGGGCCTTCAGCGCGGGGCTTTGCAGCTTTTGCTCCACCTCCAGCAAGGCGTCGATCACGGTGTTACCGGTGGTGAAGATGTGCGCCGGATCGGTGCCTTCGCGCAACAGGTTGTCGCGCGAGGTGGCAGTGGGCGCGAAGTGCAGCGCGGCAAGGGCCCCGGTCAGGCGCCGGTTGCCTTCTTCGGGCCAGGGTGAATACAGGTTGTTGGTGCGCAGGCCGGCTTCTACATGGGCAATGGGGATTTGCTGG  >  NC_002947/2038860‑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: ‑