Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,580,567 A→T 20.6% V277D (GTC→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,580,5670AT20.6% 176.2 / 42.6 97V277D (GTC→GAC) ttgBefflux pump membrane protein TtgB
Reads supporting (aligned to +/- strand):  ref base A (35/40);  new base T (9/11);  total (45/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.03e-01

GGAATACCGCCTTGACTCCCGGCGGGAAGAACGGTTCCAGGCCTTTGATGGTCTCGCGCAGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCG‑‑TCTTTGTTGACCTTGAGCAGGATACTCTCGAACTGCTCGGCGGTTTGCAGGCGGGTCTTGCCGATGATGGTGGCGTTCAGCTGG  >  NC_002947/1580421‑1580713
                                                                                                                                                   |                                                                                                                                                    
ggAATACCGCCTTGACTCCCGGCGGGAAGAACGGTTCCAGGCCTTTGATGGTCTCGCGCAGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGTCg                                                                                                                                                    <  5:98762/149‑1 (MQ=255)
 gAATACCGCCTTGACTCCCGGCGGGAAGAACGGTTCCAGGCCTTTGATGGTCTCGCGCAGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGCCgg                                                                                                                                                   >  4:213534/1‑149 (MQ=255)
 gAATACCGCCTTGACTCCCGGCGGGAAGAACGGTTCCAGGCCTTTGATGGTCTCGCGCAGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACgt                                                                                                                                                   >  6:203953/1‑148 (MQ=255)
 gAATACCGCCTTGACTCCCGGCGGGAAGAACGGTTCCAGGCCTTTGATGGTCTCGCGCAGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACgg                                                                                                                                                   >  5:91865/1‑149 (MQ=255)
      ccGCCTTGACTCCCGGCGGGAAGAACGGTTCCAGGCCTTTGATGGTCTCGCGCAGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGTCGTAg                                                                                                                                              >  5:347962/1‑149 (MQ=255)
      ccGCCTTGACTCCCGGCGGGAAGAACGGTTCCAGGCCTTTGATGGTCTCGCGCAGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGTCGTAg                                                                                                                                              <  3:246827/149‑1 (MQ=255)
      ccGCCTTGACTCCCGGCGGGAAGAACGGTTCCAGGCCTTTGATGGTCTCGCGCAGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGGCGTAg                                                                                                                                              >  1:69413/1‑149 (MQ=255)
      ccGCCTTGACTCCCGGCGGGAAGAACGGTTCCAGGCCTTTGATGGTCTCGCGCAGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCGTGCCGTTGAACTGGGCGCTGACGTCGTAg                                                                                                                                              >  5:342038/1‑149 (MQ=255)
         ccTTGACTCCCGGCGGGAAGAACGGTTCCAGGCCTTTGATGGTCTCGCGCAGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGGCGTAGttt                                                                                                                                           <  6:91865/149‑1 (MQ=255)
            tGACTCCCGGCGGGAAGAACGGTTCCAGGCCTTTGATGGTCTCGCGCAGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGGCGTAGTTTTCa                                                                                                                                        >  4:239482/1‑149 (MQ=255)
             gACTCCCGGCGGGAAGAACGGTTCCAGGCCTTTGAAGGTCTCGCGCAGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGTCGGCGTAGTTTTCa                                                                                                                                        <  1:255125/148‑1 (MQ=255)
              aCTCCCGGCGGGAAGAACGGTTCCAGGCCTTTGATGGTCTCGCGCAGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGTCGTAGTTTTCAc                                                                                                                                       <  8:126675/148‑1 (MQ=255)
                tCCCGGCGGGAAGAACGGTTCCAGGCCTTTGATGGTCTCGCGCAGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGGCGTAGTTTTCAcccc                                                                                                                                    >  6:348653/1‑149 (MQ=255)
                 cccGGGGGGACGCCCGGGGCCCGGCCGTTGATGGGCTCGCGCAGGGCCTGGGCGGTGCCCAGCGCGGGGGCGCCGTTTGCCAGTTTTCCCGCCCGGCCGGAAGCCGGCTTGCCGTTGACCTGGGCGCTGCCGGCGTAGTTTTCAccccc                                                                                                                                   <  8:107316/149‑1 (MQ=255)
                   cggcggGAAGAACGGTTCCAGGCCTTTGATGGTCTCGCGCAGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGTCGTAGtt                                                                                                                                            >  4:329536/1‑138 (MQ=255)
                   cggcggGAAGAACGGTTCCAGGCCTTTGATGGTCTCGCGCAGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGTCGTAGtt                                                                                                                                            <  3:329536/138‑1 (MQ=255)
                    ggcggGAAGAACGGTTCCAGGCCTTTGATGGTCTCGCGCAGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGTCGGCGTAGTTTTCACCCCCCAg                                                                                                                                >  5:280684/1‑149 (MQ=255)
                    ggcggGAAGAACGGTTCCAGGCCTTTGATGGTCTCGCGCAGTGCCTTGGCGGTGGCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGGCGTAGTTTTCACCCCCCAg                                                                                                                                <  7:255196/149‑1 (MQ=255)
                            gaaCGGTTCCAGGCCTTTGATGGTCTCGCGCAGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGTCGTAGTTTTCACCCCCCAGGCCTACtt                                                                                                                        >  5:318226/1‑149 (MQ=255)
                             aaCGGTTCCAGGCCTTTGATGGTCTCGCGCAGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGTCGTAGTTTTCACCCCCCAGGCCTACTTg                                                                                                                       >  8:281514/1‑149 (MQ=255)
                             aaCGGTTCCAGGCCTTTGATGGTCTCGCGCAGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGTCGTAGTTTTCACCCCCCAGGCCTACTTg                                                                                                                       >  1:321748/1‑149 (MQ=255)
                              aCGGTTCCAGGCCTTTGATGGTCTCGCGCAGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGGCGTAGTTTTCACCCCCCAGGCCTACTTgg                                                                                                                      <  1:63480/149‑1 (MQ=255)
                                ggTTCCAGGCCTTTGATGGTCTCGCGCAGTGCGTTGGCGGTGTCGAGCGCGTTGGCGCCGGTTGCGAGTTTTTCCGGCATGCCGGAAGCCGGGTTGGCGTTGAACTGGGCGCTGTCGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGc                                                                                                                    <  5:4206/149‑1 (MQ=255)
                                ggTTCCAGGCCTTTGATGGTCTCGCGCAGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGc                                                                                                                    <  3:239482/149‑1 (MQ=255)
                                ggTTCCAGGCCTTTGATGGTCTCGCGCAGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGc                                                                                                                    <  5:348653/149‑1 (MQ=255)
                                      aGGCCTTTGATGGTCTCGCGCAGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATc                                                                                                              <  5:387720/149‑1 (MQ=255)
                                       ggCCTTTGATGGTCTCGCGCAGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCg                                                                                                             <  1:146940/149‑1 (MQ=255)
                                          cTTTGATGGTCTCGCGCAGTGCCTTGGCGGGGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCGAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGccc                                                                                                          <  6:318245/149‑1 (MQ=255)
                                           tttGATGGTCTCGCGCAGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGGCGTAGTTTTCACCCCCCAGGCCTCCTTGGGCGACATCGcc                                                                                                           >  6:244109/1‑147 (MQ=255)
                                               aTGGTCTCGCGCAGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGTCGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGgcg                                                                                                     <  5:38897/149‑1 (MQ=255)
                                               aTGGTCTCGCGCAGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGTCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGATGTCGCCCaatta                                                                                                     <  6:48940/149‑5 (MQ=255)
                                               aTGGTCTCGCGCAGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGgcg                                                                                                     >  4:120269/1‑149 (MQ=255)
                                               aTGGTCTCCCGGAGTGCCTTGGGGGTGTCCAGCGCGTTGGCGCCGGGTGCCATTTTTACCGCCAGGCCGGAAGCCGGCTTGCTGGTTGAACTGGGCGCTGACGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGATGTCGCCCaagt                                                                                                      <  1:311280/149‑4 (MQ=255)
                                                 ggactcGCGCAGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGGCGTAGTTTTCACCCCCCAGGCCAACTTGGGCGACAACGCCCAGGCGCa                                                                                                   >  2:341721/4‑149 (MQ=255)
                                                  gTCTCGCGCAGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGTCGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCAc                                                                                                  <  6:280684/149‑1 (MQ=255)
                                                   tctcGCGCAGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGTCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCAcc                                                                                                 >  3:313809/1‑149 (MQ=255)
                                                   tctcGCGCAGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGTCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCAcc                                                                                                 >  3:95854/1‑149 (MQ=255)
                                                    ctcGCGCAGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCt                                                                                                <  3:213534/149‑1 (MQ=255)
                                                     tcGCGCAGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCtg                                                                                               >  7:189713/1‑149 (MQ=255)
                                                       gcgcAGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGTCGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCtgt                                                                                              >  3:141817/1‑148 (MQ=255)
                                                       gcgcAGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGTCGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCtgt                                                                                              >  5:253582/1‑148 (MQ=255)
                                                       gcgcAGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGTCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCtgt                                                                                              >  5:40933/1‑148 (MQ=255)
                                                       gcgcAGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCtgt                                                                                              >  2:257206/1‑148 (MQ=255)
                                                       gcgcAGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCtgt                                                                                              <  1:257206/148‑1 (MQ=255)
                                                          cAGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGTCGTTGAACTGGGCGCTGACGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAAc                                                                                          >  7:133324/1‑149 (MQ=255)
                                                           aGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGTCGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAAcc                                                                                         >  6:391476/1‑149 (MQ=255)
                                                           aGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGCAcc                                                                                         >  3:259087/1‑149 (MQ=255)
                                                                cTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCG‑‑TCtt                                                                                  <  2:69413/149‑1 (MQ=255)
                                                                   ggCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCG‑‑TCttt                                                                                 <  2:172095/147‑1 (MQ=255)
                                                                   ggCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCG‑‑TCTttg                                                                                <  3:120269/148‑1 (MQ=255)
                                                                    gCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACGGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGTCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCG‑‑TCTttgtt                                                                              <  4:95854/149‑1 (MQ=255)
                                                                    gCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGTCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCG‑‑TCTttgtt                                                                              <  4:313809/149‑1 (MQ=255)
                                                                     cGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCG‑‑TCTttgttg                                                                             <  8:189713/149‑1 (MQ=255)
                                                                      ggTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGTCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCG‑‑TCTTTGTTGa                                                                            <  5:56633/149‑1 (MQ=255)
                                                                        tgtCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGTCGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCG‑‑TCTTTGTTGAcc                                                                          <  5:391476/149‑1 (MQ=255)
                                                                            cAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCGGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCG‑‑TCTTTGTTGACCtt                                                                        >  6:289061/1‑147 (MQ=255)
                                                                            cAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGTCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCG‑‑TCTTTGTTGACCtt                                                                        >  4:56494/1‑147 (MQ=255)
                                                                            cAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGTCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCG‑‑TCTTTGTTGACCTTGa                                                                      <  7:281514/149‑1 (MQ=255)
                                                                            cAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCG‑‑TCTTTGTTGACCTTg                                                                       <  2:258163/148‑1 (MQ=255)
                                                                            cAGCGCGTTGGCGCCGGTTGCCAGTTTTAACGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGCGCGCTGACGGCGTAGTTTTCACCCCCCAGACCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCG‑‑TCTTTGTTGACCtt                                                                        >  2:31014/1‑147 (MQ=255)
                                                                               cgcgTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGTCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCG‑‑TCTTTGTTGACCTTGAGCa                                                                   <  2:321748/149‑1 (MQ=255)
                                                                                     ggCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGTCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCG‑‑TCTTTGTTGACCTTGAGCAGGATa                                                              <  5:293201/148‑1 (MQ=255)
                                                                                      gcgcCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGTCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCG‑‑TCTTTGTTGACCTTGAGCAGGATAct                                                            >  5:7117/1‑149 (MQ=255)
                                                                                      gcgaCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAAATGGGCGCTGTCGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCG‑‑TCTTTGTTGACCTTGAGCAGGATAct                                                            >  2:138466/1‑149 (MQ=255)
                                                                                       cgcCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGTCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCG‑‑TCTTTGTTGACCTTGAGCAGGATActc                                                           <  6:40933/149‑1 (MQ=255)
                                                                                            gTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGTCGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGgc                                                                                                      >  2:170151/1‑103 (MQ=255)
                                                                                            gTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGTCGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGgc                                                                                                      <  1:170151/103‑1 (MQ=255)
                                                                                                 cAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCG‑‑TCTTTGTTGACCTTGAGCAGGATACTCTCGAACTGCt                                                 >  1:293647/1‑149 (MQ=255)
                                                                                                      ttACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCG‑‑TCTTTGTTGACCTTGAGc                                                                    >  4:26525/1‑125 (MQ=255)
                                                                                                      ttACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCG‑‑TCTTTGTTGACCTTGAGc                                                                    <  3:26525/125‑1 (MQ=255)
                                                                                                        aCCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGTCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCG‑‑TCTTTGTTGACCTTGAGCAGGATACTCTCGAACTGCTcggcgg                                           <  5:18315/148‑1 (MQ=255)
                                                                                                        aCCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCG‑‑TCTTTGTTGACCTTGAGCAGGATACTCTCGAACTGCTCGGCGGt                                          >  2:236593/1‑149 (MQ=255)
                                                                                                         ccGCCAGGCCGGAAGCGGGCTTGCCGTTGAACTGGGCGCTGACGGCGTAGTTTTCTCCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCG‑‑TCTTTGTTGACCTTGAGCAGGATACTCTCGAACTGCTCGGCGGtt                                         <  1:341721/149‑1 (MQ=255)
                                                                                                          cGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGTCGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCG‑‑TCTTTGTTGACCTTGAGCAGGATACTCTCGAACTGCTCGGCGGttt                                        >  1:145828/1‑149 (MQ=255)
                                                                                                          cGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCG‑‑TCTTTGTTGACCTTGAGCAGGATACTCTCGAACTGCTCGGCGGttt                                        >  2:92211/1‑149 (MQ=255)
                                                                                                           gCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGTCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCG‑‑TCTTTGTTGACCTTGAGAAGGATACTCTCGAACTGCTCGGCGGTTTg                                       >  5:33634/1‑149 (MQ=255)
                                                                                                            ccAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCG‑‑TCTTTGTTGACGTTGAGCGGGATAATCGCGAAATGGGCGGCGGTTTGc                                      >  4:142143/1‑149 (MQ=255)
                                                                                                            ccAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCG‑‑TCTTTGTTGACCTTGAGCAGGATACTCTCGAACTGCTCGGCGGTTTGc                                      <  1:92211/149‑1 (MQ=255)
                                                                                                              aGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGTCGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCG‑‑TCTTTGTTGACCTTGAGCAGGATACTCTCGAACTGCTCGGCGGTTTGCAg                                    <  3:59336/149‑1 (MQ=255)
                                                                                                                  cGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGTCGGCGTAGTTTTCACCCCCCAGGCCTACTTggg                                                                                                                     >  8:239770/1‑66 (MQ=255)
                                                                                                                  cGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGTCGGCGTAGTTTTCACCCCCCAGGCCTACTTggg                                                                                                                     <  7:239770/66‑1 (MQ=255)
                                                                                                                  cGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGTCGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCG‑‑TCTTTGTTGACCTTGAGCAGGATACTCTCGAACTGCTCGGCGGTTTGCAGGCgg                                <  6:253582/149‑1 (MQ=255)
                                                                                                                  cGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGTCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCG‑‑TCTTTGTTGACCTTGAGCAGGATACTCTCGAACTGCTCGGCGGTTTGCAGGCgg                                <  5:203953/149‑1 (MQ=255)
                                                                                                                        ccGGCTTGCCGTTGAACTGGGCGCTGACGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCG‑‑TCTTTGTTGACCTTGAGCAGGATACTCTCGAACTGCTCGGCGGTTTGCAGGCGGGTCTTg                          >  2:43642/1‑149 (MQ=255)
                                                                                                                         cGGCTTGTCGTTGAACTGGGCGCTGACGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCG‑‑TCTTTGTTGACCTTGAGCAGGATACTCTCGAACTGCTCGGCGGTTTGCAGGCGGGTCTTGc                         <  8:133324/149‑1 (MQ=255)
                                                                                                                                ccGTTGAACTGGGCGCTGACGTCGTAGTTTTCACCCCCCAGTCCTACTTGGGCGACTTCGCCCAGGCGCACCTGTGAACCG‑‑TCTTTGTTGACCTTGAGCAGGATACTCTCGAACTGCTCGGCGGTTTGCAGGCGGGTCTTGCCgatga                   <  6:347962/148‑1 (MQ=255)
                                                                                                                                ccGTTGAACTGGGCGCTGACGTCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCG‑‑TCTTTGTTGACCTTGAGCAGGATACTCTCGAACTGCTCGGCGGTTTGCAGGCGGGTCTTGCCgatga                   <  6:342038/148‑1 (MQ=255)
                                                                                                                                  cGTTGAACTGGGCGCTGACGTCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCG‑‑TCTTTGTTGACCTTGAGCAGGATACTCTCGAACTGCTCGGCGGTTTGCAGGCGGGTCTTGCCGatgatg                 <  6:7117/149‑1 (MQ=255)
                                                                                                                                   gTTGAACTGGGCGCTGACGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCG‑‑TCTTTGTTGACCTTGAGCAGGATACTCTCGAACTGCTCGGCGGttt                                        >  2:43735/1‑125 (MQ=255)
                                                                                                                                   gTTGAACTGGGCGCTGACGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCG‑‑TCTTTGTTGACCTTGAGCAGGATACTCTCGAACTGCTCGGCGGttt                                        <  1:43735/125‑1 (MQ=255)
                                                                                                                                      gAACTGGGCGCTGTCGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCG‑‑TCTTTGTTGACCTTGAGCAGGATACTCTCGAACTGCTCGGCGGTTTGCAGGCGGGTCTTGCCGATGAtggtg              <  4:141817/148‑1 (MQ=255)
                                                                                                                                      gAACTGGGCGCTGACGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCG‑‑TCTTTGTTGACCTTGAGCAGGATACTCTCGAACTGCTCGGCGGTTTGCAGGCGGGTCTTGCCGATGAtggtg              <  2:293647/148‑1 (MQ=255)
                                                                                                                                       aaCTGGGCGCTGACGTCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCG‑‑TCTTTGTTGACCTTGAGCAGGATACTCTCGAACTGCTCGGCGGTTTGCAGGCGGGTCTTGCCGATGATGGTGGc            >  1:62416/1‑149 (MQ=255)
                                                                                                                                        aCTGGGCGCTGACGGCGTAGTTTTCACCCCCCAGACCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCG‑‑TCTTTGTTGACCTTGAGCAGGATACTCTCGAACTGCTCGGCGGTNTGCAGGCGGGTCTTGCCGATGATGGTGGCg           <  1:31014/149‑1 (MQ=255)
                                                                                                                                            ggCTCTGACGGCGTAGTTTTCACCCCCCAGGCCTACATGGGCGACAACGCCCAGGCGCACCTGTGAACCG‑‑TCTTTGTTGACCTTGAGCAGGATACTCTCGAACTGCTCGGCGGTTTGCAGGCGGGTCTTGCCGATGATGGTGGCGtt         <  6:33634/147‑1 (MQ=255)
                                                                                                                                              cgcTGTCGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCG‑‑TCTTTGTTGACCTTGAGCAGGATACTCTCGAACTGCTCGGCGGTTTGCAGGCGGGTCTTGCCGATGATGGTGGCGTTCAGc     >  6:82902/1‑149 (MQ=255)
                                                                                                                                               gcTGACGGCGTATTTTTCACCACCCAGGCCTACTTGGGCGAGATCGCCCAGGCGCACCTGTGAACCGACGCTTTGTTGACCTTGAGCAGGATACTCTCGAACTGCTCGGCGGTTTGCAGGCGGGTCTTGCCGATGATGGTGGAGTTCAg      >  1:315937/1‑149 (MQ=255)
                                                                                                                                                   aCGTCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCG‑‑TCTTTGTTGACCTTGAGCAGGATACTCTCGAACTGCTCGGCGGTTTGCAGGCGGGTCTTGCCGATGATGGTGGCGTTCAGCTgg  <  3:56494/147‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GGAATACCGCCTTGACTCCCGGCGGGAAGAACGGTTCCAGGCCTTTGATGGTCTCGCGCAGTGCCTTGGCGGTGTCCAGCGCGTTGGCGCCGGTTGCCAGTTTTACCGCCAGGCCGGAAGCCGGCTTGCCGTTGAACTGGGCGCTGACGGCGTAGTTTTCACCCCCCAGGCCTACTTGGGCGACATCGCCCAGGCGCACCTGTGAACCG‑‑TCTTTGTTGACCTTGAGCAGGATACTCTCGAACTGCTCGGCGGTTTGCAGGCGGGTCTTGCCGATGATGGTGGCGTTCAGCTGG  >  NC_002947/1580421‑1580713

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: