Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,581,266 G→T 100% A44D (GCC→GAC)  ttgB ← efflux pump membrane protein TtgB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,581,2660GT99.8% 230.2 / ‑7.4 73A44D (GCC→GAC) ttgBefflux pump membrane protein TtgB
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base T (38/32);  total (39/34)
Fisher's exact test for biased strand distribution p-value = 2.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GCGGTAATGGTCATGCTGCCGTCGGAGTTACTTTCCGACGACACATAACGCAGGTTGTCGATACC‑GTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGGCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCACCAAGGCGATCACCCAGGCGAAGATCGGGCGATCGATAAAGAACTTCGACATGGTTTA  >  NC_002947/1581122‑1581402
                                                                                                                                                  |                                                                                                                                        
gCGGGAATGGTCATGCTGCCGTCGGAGTTACTTTCCGACGACACATAACGCAGGTTGTCGATACC‑GTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACggcgatgtcgat                                                                                                                                      <  1:327676/149‑1 (MQ=255)
 cGGTAATGGTCATGCTGCCGTCGGAGTTACTTTCCGACGACACATAACGCAGGTTGTCGATACC‑GTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGgcgatgtcgatg                                                                                                                                     <  7:264555/149‑1 (MQ=255)
  ggTAATGGTCATGCTGCCGTCGGAGTTACTTTCCGACGACACATAACGCAGGTTGTCGATACC‑GTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGcgatgtcgatgg                                                                                                                                    >  4:362345/1‑149 (MQ=255)
      aTGGTCATGCTGCCGTCGGAGTTACTTTCCGACGACACATAACGCAGGTTGTCGATACC‑GTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGAtggctgg                                                                                                                                >  4:80350/1‑149 (MQ=255)
        ggTCATGCTGCCGTCGGAGTTACTTTCCGACGACACATAACGCAGGTTGTCGATACC‑GTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTggcg                                                                                                                              >  3:65204/1‑149 (MQ=255)
        ggTCATGCTGCCGTCGGAGTTACTTTCCGACGACACATAACGCAGGTTGTCGATACC‑GTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTggcg                                                                                                                              >  6:377181/1‑149 (MQ=255)
         gTCATGCTGCCGTCGGAGTTACTTTCCGACGACACATAACGCAGGTTGTCGATACC‑GTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTggcgg                                                                                                                             >  1:226652/1‑149 (MQ=255)
         gTCATGCTGCCGTCGGAGTTACTTTCCGACGACACATAACGCAGGTTGTCGATACC‑GTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTggcgg                                                                                                                             >  3:325666/1‑149 (MQ=255)
                   cGTCGGAGTTACTTTCCGACGACACATAACGCAGGTTGTCGATACC‑GTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTg                                                                                                                   <  2:226652/149‑1 (MQ=255)
                     tCGGAGTTACTTTCCGACGACACATAACGCAGGTTGTCGATACC‑GTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTggg                                                                                                                 <  7:354517/149‑1 (MQ=255)
                     tCGGAGTTACTTTCCGACGACACATAACGCAGGTTGTCGATACC‑GTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTggg                                                                                                                 <  7:387250/149‑1 (MQ=255)
                         aGTTACTTTCCGACGACACATAACGCAGGTTGTCGATACC‑GTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTAc                                                                                                             <  3:80350/149‑1 (MQ=255)
                            tACTTTCCGACGACACATAACGCAGGTTGTCGATACC‑GTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTgg                                                                                                          >  3:118102/1‑149 (MQ=255)
                             aCTTTCCGACGACACATAACGCAGGTTGTCGATACC‑GTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGt                                                                                                         >  2:167490/1‑149 (MQ=255)
                              cTTTCCGACGACACATAACGCAGGTTGTCGATACC‑GTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTT‑GCGCCGAGGCGCCCTGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGt                                                                                                         <  2:244069/148‑1 (MQ=255)
                               tttCCGACGACACATAACGCAGGTTGTCAATACC‑GTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGtt                                                                                                        >  8:56421/1‑148 (MQ=255)
                                ttCGGAAGACTCATAACGCTGGTTGTCGATACC‑GTTGTGCTGCTGCTGGATCTCCTGCTCCACGGTGTCCTGGACGGTTT‑GCGCCGAGGCGCCTGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTAGTGGTTg                                                                                                       <  5:335150/148‑1 (MQ=255)
                                ttCCGACGACACATAACGCAGGTTGTCGATACC‑GTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGGTGGCTGGCGGCGCGATACTGGGGTACTGGTTGa                                                                                                      >  5:243559/1‑149 (MQ=255)
                                 tcagacgacACATAACGCAGGTTGTCGATACCTGTTGAGATGCTGCTCGATCACCTCCAACACGGTGTCCTGAACGGTTT‑CCGCCGAGCCGCCCGGGTAGGTCACGGCGATGACGATGGCTGGCGCCCCGATACTGGGGTACTGGTTGa                                                                                                      >  1:53557/4‑149 (MQ=255)
                                    gacgacACATAACGCAGGTTGTCGATACC‑GTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCAATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCgg                                                                                                  >  7:176556/1‑149 (MQ=255)
                                     acgacACATAACGCAGGTTGTCGATACC‑GTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGc                                                                                                 <  4:65204/149‑1 (MQ=255)
                                     acgacACATAACGCAGGTTGTCGATACC‑GTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGc                                                                                                 <  6:143118/149‑1 (MQ=255)
                                          acaTAACGCAGGTTGTCGATACC‑GTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCtt                                                                                            >  1:239919/1‑149 (MQ=255)
                                            aTAACGCAGGTTGTCGATACC‑GTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCa                                                                                          <  2:239919/149‑1 (MQ=255)
                                                cGCAGGTTGTCGATACC‑GTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGAt                                                                                      <  1:167490/149‑1 (MQ=255)
                                                cGCAGGTTGTCGATACC‑GTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGAt                                                                                      <  3:24501/149‑1 (MQ=255)
                                                 gCAGGTTGTCGATACC‑GTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATg                                                                                     <  2:9135/149‑1 (MQ=255)
                                                 gCAGGTTGTCGATACC‑GTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATg                                                                                     <  3:362345/149‑1 (MQ=255)
                                                   aGGTTGTCGATACC‑GTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGa                                                                                   <  2:316307/149‑1 (MQ=255)
                                                     gTTGTCGATACC‑GTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGAc                                                                                  >  2:157225/1‑148 (MQ=255)
                                                             tACC‑GTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGcccc                                                                          >  4:149839/1‑148 (MQ=255)
                                                             tACC‑GTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGCCAGCTTCAGGATGGACAAGGCCCCg                                                                         >  4:424649/1‑149 (MQ=255)
                                                                       gctgctgctCGATCACCTGCACCTCGGTGTCCTGCACGGTTT‑GCTCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGa                                                                                       <  7:403742/126‑1 (MQ=255)
                                                                       gctgctgctCGATCACCTGCACCACGGAGTCCTGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGCTCACGGCGATGTCGATGCCAGGCGGAGCGATACTGGGGTACTGGTTGCTCGGCAGCTTCAGGa                                                                                       >  8:403742/1‑126 (MQ=255)
                                                                       gctgatgctCGATCACCTGCACCACGTTGTCCTGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGATGGCGGCGCGATACGGGGGTACTGGTGGATCGGCCGCTTCAGGATTGACAAGGCACCGCCCCGCAt                                                                 >  3:247877/1‑148 (MQ=255)
                                                                          gctgctCGATCACCTGCACCACGGTGTCCTGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCgg                                                                                                  >  3:262494/1‑112 (MQ=255)
                                                                          gctgctCGATCACCTGCACCACGGTGTCCTGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCgg                                                                                                  <  4:262494/112‑1 (MQ=255)
                                                                          gctgctCGATCACCTGCACCACGGTGTCCTGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGAt                                                              >  8:37472/1‑148 (MQ=255)
                                                                           ctggtCGTTCACCTGCACCACGGTGTCCTGCTCGGTTT‑GCGACGAGGCTCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTAATGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATc                                                             <  2:155376/148‑1 (MQ=255)
                                                                           ctgctCGATCACCTGCACCACGGTGTCCTGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATc                                                             <  5:113867/148‑1 (MQ=255)
                                                                             gctCGATCACCTGCACCACGGTGTCCTGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAACATGATCAc                                                           >  5:141135/1‑148 (MQ=255)
                                                                              ctCGATCACCTGCACCACGGTGTCCTGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCACCa                                                         <  2:266765/149‑1 (MQ=255)
                                                                                cGATCACCTGCACCACGGTGTCCTGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCACCAAg                                                       >  4:201793/1‑149 (MQ=255)
                                                                                cGATCACCTGCACCACGGTGTCCTGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCACCAAg                                                       >  8:250732/1‑149 (MQ=255)
                                                                                 gacaaCCTGCACCACGGTGTCCTGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCACCAAgg                                                      <  6:426120/145‑1 (MQ=255)
                                                                                   tCACCTGCACCACGGTGTCCTGCACGGTTTGGCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCACCAAGGc                                                     >  1:46255/1‑149 (MQ=255)
                                                                                   tCACCTGCACCACGGTGTCCTGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCACCAAGGCg                                                    <  7:108417/149‑1 (MQ=255)
                                                                                      ccTGCACCACGGTGTCCTGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCACCAAGGCGAt                                                  >  2:121769/1‑148 (MQ=255)
                                                                                      ccTGCACCACGGTGTCCTGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCACCAAGGCGATc                                                 <  7:37472/149‑1 (MQ=255)
                                                                                      ccTGAACCACGGTGTCCTGCACGGTTT‑GCGCCGAGGCGCCCGCGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCACCAAGGCGAt                                                  >  5:300950/1‑148 (MQ=255)
                                                                                       cttctccacGTTGTCCTGCTGTGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGGCGATGGCTGGCGGCTCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCACCAAGGCGATc                                                 <  6:300950/143‑1 (MQ=255)
                                                                                          caccacGGTGTCCTGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCACCAAGGCGATCAcc                                              <  1:121769/148‑1 (MQ=255)
                                                                                           agcacGGTGTCCTGCACGGTTT‑GCGTCGAGGCGCCCGGGTAGGTCACGGGGATGTCGATGGCTTGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCGCCGACCAGCATGATCACCAAGGCGATCAccc                                             <  2:53557/146‑1 (MQ=255)
                                                                                            ccacGGTGTCCTTCAGGGTTT‑GCGCCGAGGCGCGCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTGAGGATGGACAAGGCCCCGACCAGCATGATCACCAAGGCGATCACCCAg                                           <  3:26986/149‑1 (MQ=255)
                                                                                             cacGGTGTCCTGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCACCAAGGCGATCACCCAgg                                          >  7:227507/1‑149 (MQ=255)
                                                                                               cGGTGTCCTGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCACCAAGGCGATCACCCAGGc                                         >  2:298684/1‑148 (MQ=255)
                                                                                                  tgtCCTGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCACCAAGGCGATCACCCAGGCg                                        >  2:52793/1‑146 (MQ=255)
                                                                                                    tCCTGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCACCAAGGCGATCACCCAGGCg                                        >  4:385944/1‑144 (MQ=255)
                                                                                                      cTGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCACCAAGGCGATCACCCAGGCGAAGATCg                                 >  8:320927/1‑149 (MQ=255)
                                                                                                       tGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCACCAAGGCGATCACCCAGGCGAAGATCgg                                >  4:80038/1‑149 (MQ=255)
                                                                                                       tGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCACCAAGGCGATCACCCAGGCGAAGATCgg                                >  1:209194/1‑149 (MQ=255)
                                                                                                        gCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCACCAAGGCGATCACCCAGGCGAAGATCggg                               >  2:106465/1‑149 (MQ=255)
                                                                                                           cGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCACCAAGGCGATCACCCAGGCGAAGATCGGGcga                            >  4:128537/1‑149 (MQ=255)
                                                                                                            ggTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCACCAAGGCGATCACCCAGGCGAAGATCGGGcgat                           >  8:40481/1‑149 (MQ=255)
                                                                                                            ggTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCACCAAGGCGATCACCCAGGCGAAGATCGGGcgat                           <  8:227507/149‑1 (MQ=255)
                                                                                                              ttt‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCACCAAGGCGATCACCCAGGCGAAGATCGGGcgatcg                         <  3:201793/149‑1 (MQ=255)
                                                                                                                t‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCACCAAGGCGATCACCCAGGCGAAGATCGGGcgatcga                        >  7:131427/1‑148 (MQ=255)
                                                                                                                  gcgcCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCACCAAGGCGATCACCCAGGCGAAGATCGGGcgatcgat                       <  5:377181/148‑1 (MQ=255)
                                                                                                                    gcCGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCACCAAGGCGATCACCCAGGCGAAGATCGGGCGATCGATa                      >  1:114810/1‑147 (MQ=255)
                                                                                                                     gcGAGGCGCCCGGGTAGGGCACGGCGATGGCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCACCAAGGCGATCACCCAGGCGAAGATCGGGCGATCGATAAAg                   <  2:142622/148‑1 (MQ=255)
                                                                                                                     ccGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATc                                                                                                    <  3:186204/68‑1 (MQ=255)
                                                                                                                     ccGAGGCGCCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATc                                                                                                    >  4:186204/1‑68 (MQ=255)
                                                                                                                           cgcCCGGGTAGGTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGATGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCACCAAGGCGATCACCCAGGCGAAGATCGGGCGATCGATAAAGAACTTc             <  3:424649/149‑1 (MQ=255)
                                                                                                                                     ggTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCACCAAGGCGATCACCCAGGCGAAGATCGGGCGATCGATAAAGAACTTCGACATGGtt    <  2:209194/148‑1 (MQ=255)
                                                                                                                                     ggTCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCACCAAGGCGATCACCCAGGCGAAGATCGGGCGATCGATAAAGAACTTCGACATGGtt    <  3:385944/148‑1 (MQ=255)
                                                                                                                                       tCACGGCGATGTCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCACCAACGCGATCACCCAGGCGAAGATCCGGCGATCGATAAAGAACTACGACATCGttaa  >  7:400117/1‑146 (MQ=255)
                                                                                                                                                  |                                                                                                                                        
GCGGTAATGGTCATGCTGCCGTCGGAGTTACTTTCCGACGACACATAACGCAGGTTGTCGATACC‑GTTGAGCTGCTGCTCGATCACCTGCACCACGGTGTCCTGCACGGTTT‑GCGCCGAGGCGCCCGGGTAGGTCACGGCGATGGCGATGGCTGGCGGCGCGATACTGGGGTACTGGTTGATCGGCAGCTTCAGGATGGACAAGGCCCCGACCAGCATGATCACCAAGGCGATCACCCAGGCGAAGATCGGGCGATCGATAAAGAACTTCGACATGGTTTA  >  NC_002947/1581122‑1581402

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: ‑