Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,279,422 A→G 8.8% A113A (GCT→GCC bcr ← bicyclomycin/cysteine/sulfonamide efflux transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,279,4220AG 8.8% 181.2 / 6.9 69A113A (GCT→GCCbcrbicyclomycin/cysteine/sulfonamide efflux transporter
Reads supporting (aligned to +/- strand):  ref base A (34/27);  new base G (4/2);  total (40/29)
Fisher's exact test for biased strand distribution p-value = 6.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.95e-01

CTCAGCCACACCAGCACCCAGCCGCCAACTATCGGTGCCATCAGCGGTGCAATGGTTGTCACCAGCATGACAAACGACATCATCCGCGAGAACTCTTCTTTCGGGTAAATATCGCGCATCAGGGCGTTAATGACCACGCTGGCCGCAGCCGCAGCCAGCCCGTGGAAGAAACGCATCACAATCAGCTGATCGATGGTGTTTGCCAACGCACACGCCACCGCGGCGGCGGCAAACACCAGCGTACCGCCGAGCACCACCGGCTTACGCCCGAAGCTGTCTGCCATCGGCCCGT  >  NC_000913/2279276‑2279567
                                                                                                                                                  |                                                                                                                                                 
cTCAGCCACACCAGCACCCAGCCGCCAACTATCGGTGCCATCAGCGGTGCAATGGTTGTCACCAGCATGACAAACGACATCATCCGCGAGAACTCTTCTTTCGGGTAAATATCGCGCATCAGGGCGTTAATGACCACGCTGgccgcagcc                                                                                                                                                >  1:51368/1‑150 (MQ=255)
        cacCAGCACCCAGCCGCCAACTATCGGTGCCATCAGCGGTGCAATGGTTGTCACCAGCATGACAAACGACATCATCCGCGAGAACTCTTCTTTCGGGTAAATATCGCGCATCAGGGCGTTAATGACCACGCTGGCTGCGGCTGcggccag                                                                                                                                        <  2:1400481/150‑1 (MQ=255)
          ccAGCACCCAGCCGCCAACTATCGGTGCCATCAGCGGTGCAATGGTTGTCACCAGCATGACAAACGACATCATCCGCGAGAACTCTTCTTTCGGGTAAATATCGCGCATCAGGGCGTTAATGACCACGCTGGCTGCGGCTGcggccagcg                                                                                                                                      <  1:479302/150‑2 (MQ=255)
          ccAGCACCCAGCCGCCAACTATCGGTGCCATCAGCGGTGCAATGGTTGTCACCAGCATGACAAACGACATCATCCGCGAGAACTCTTCTTTCGGGTAAATATCGCGCATCAGGGCGTTAATGACCACGCTGGCCGCAGCCGCagccagcc                                                                                                                                      >  1:1467846/1‑150 (MQ=255)
             gCACCCAGCCGCCAACTATCGGTGCCATCAGCGGTGCAATGGTTGTCACCAGCATGACAAACGACATCATCCGCGAGAACTCTTCTTTCGGGTAAATATCGCGCATCAGGGCGTTAATGACCACGCTGGCCGCAGCCGCAGCCAGCCCGt                                                                                                                                   >  1:408434/1‑150 (MQ=255)
               aCCCAGCCGCCAACTATCGGTGCCATCAGCGGTGCAATGGTTGTCACCAGCATGACAAACGACATCATCCGCGAGAACTCTTCTTTCGGGTAAATATCGCGCATCAGGGCGTTAATGACCACGCTGGCCGCAGCCGCAGCCAGCCCGTgg                                                                                                                                 <  1:1472289/150‑1 (MQ=255)
                cccAGCCGCCAACTATCGGTGCCATCAGCGGTGCAATGGTTGTCACCAGCATGACAAACGACATCATCCGCGAGTACTCTTCTTTCGGGTAAATATCGCGCATCAGGGCGTTAATGACCACGCTGGCCGCAGCCGCAGCCAGCCCGTGga                                                                                                                                <  1:854527/150‑1 (MQ=255)
                         cAACTATCGGTGCCATCAGCGGTGCAATGGTTGTCACCAGCATGACAAACGACATCATCCGCGAGAACTCTTCTTTCGGGTAAATATCGCGCATCAGGGCGTTAATGACCACGCTGGCCGCAGCCGCAGCCAGCCCGTGGAAGAAACGCa                                                                                                                       >  2:623375/1‑150 (MQ=255)
                           aCTATCGGTGCCATCAGCGGTGCAATGGTTGTCACCAGCATGACAAACGACATCATCCGCGAGAACTCTTCTTTCGGGTAAATATCGCGCATCAGGGCGTTAATGACCACGCTGGCCGCAGCCGCAGCCAGCCCGTGGAAGAAACGCATc                                                                                                                     <  1:1566990/150‑1 (MQ=255)
                            cTATCGGTGCCATCAGCGGTGCAATGGTTGTCACCAGCATGACAAACGACATCATCCGCGAGAACTCTTCTTTCGGGTAAATATCGCGCATCAGGGCGTTAATGACCACGCTGGCCGCAGCCGCAGCCAGCCCGTGGAAGAAACGCATca                                                                                                                    <  1:1186178/150‑1 (MQ=255)
                                 ggTGCCATCAGCGGTGCAATGGTTGTCACCAGCAGGACAAACGACATCATCCGCGAGAACTCTTCTTTCGGGTAAATATCGAGCATCAGGGCGGTAATGACCACGCTGGCCGCGGCCGCAGCCAGCACGTGGAAGAATAGCATCACAAtg                                                                                                               >  1:230548/1‑149 (MQ=255)
                                    gCCATCAGCGGTGCAATGGTTGTCACCAGCATGACAAACGACATCCTCCGCGAGAACTCTTCTTTCGGGTAAATATCGCGCATCAGGGCGTTAATGACCACGCTGGCCGCAGCCGGAGCCAGCCCCTGGAAGAAACGCATCAAAAACAGc                                                                                                            >  1:303274/1‑150 (MQ=255)
                                          atcGGTGCAATGGTTGTCACCAGCATGACAAACGACATCATCCGAGATAACTCTTCTTTCGGGTAAATATCGCGCATCAGGGCGTTAATGACCACGCTGGCCGCAGCCGCAGCCAGCCGGTGGAAGAAACGCATCACAATCAGCTGATCg                                                                                                      <  2:35051/148‑1 (MQ=255)
                                            cGGTGCAATGGTTGTCACCAGCATGACAAACGACATCATCCGCGAGAACTCTTCTTTCGGGTAAATATCGCGCATCAGGGCGTTAATGACCACGCTGGCCGCAGCCGCAGCCAGCCCGTGGAAGAAACGCATCACAATCAGCTGATCGAt                                                                                                    <  1:279371/150‑1 (MQ=255)
                                            cGGTGCAATGGTTGTCACCAGCATGACAAACGACATCATCCGCGAGAACTCTTCTTTCGGGTAAATATCGCGCATCAGGGCGTTAATGACCACGCTGGCCGCAGCCGCAGCCAGCCCGTGGAAGAAACGCATCACAATCAGCTGATCGAt                                                                                                    <  1:166187/150‑1 (MQ=255)
                                             ggTGCAATGGTTGTCACCAGCATGACAAACGACATCATCCGCGAGAACTCTTCTTTCGGGTAAATATCGCGCATCAGGGCGTTAATGACCACGCTGGCCGCAGCCGCAGCCAGCCCGTGGAAGAAACGCATCACAATCAGCTGATCGATg                                                                                                   <  2:1654787/150‑1 (MQ=255)
                                             ggTGCAATGGTTGTCACCAGCATGACAAACGACATCATCCGCGAGAACTCTTCTTTCGGGTAAATATCGCGCATCAGGGCGTTAATGACCACGCTGGCCGCAGCCGCAGCCAGCCCGTGGAAGAAACGCATCACAATCAGCTGATCGATg                                                                                                   <  1:1164421/150‑1 (MQ=255)
                                              gTGCAATGGTTGTCACCAGCATGACAAACGACATCATCCGCGAGAACTCTTCTTTCGGGTAAATATCGCGCATCAGGGCGTTAATGACCACGCTGGCCGCAGCCGCAGCCAGCCCGTGGAAGAAACGCATCACAATCAGCTGATCGATgg                                                                                                  >  1:1651216/1‑150 (MQ=255)
                                              gTGCAATGGTTGTCACCAGCATGACAAACGACATCATCCGCGAGAACTCTTCTTTCGGGTAAATATCGCGCATCAGGGCGTTAATGACCACGCTGGCCGCAGCCGCAGCCAGCCCGTGGAAGAAACGCATCACAATCAGCTGATCGATgg                                                                                                  >  1:947987/1‑150 (MQ=255)
                                                gCAATGGTTGTCACCAGCATGACAAACGACATCATCCGCGAGAACTCTTCTTTCGGGTAAATATCGCGCATCAGGGCGTTAATGACCACGCTGGCCGCAGCCGCAGCCAGCCCGTGGAAGAAACGCATCACAATCAGCTGATCGATGgtg                                                                                                >  2:928626/1‑150 (MQ=255)
                                                  aaTGGTTGTCACCAGCATGACAAACGACATCATCCGCGAGAACTCTTCTTTCGGGTAAATATCGCGCATCAGGGCGTTAATGACCACGCTGGCCGCAGCCGCAGCCAGCCCGTGGAAGAAACGCATCACAATCAGCTGATCGATGGTGtt                                                                                              >  1:1279665/1‑150 (MQ=255)
                                                   aTGGTTGTCACCAGCATGACAAACGACATCATCCGCGAGAACTCTTCTTTCGGGTAAATATCGCGCATCAGGGCGTTAATGACCACGCTGGCCGCAGCCGCAGCCAGCCCGTGGAAGAAACGCATCACAATCAGCTGATCGATGGTGttt                                                                                             >  1:1505066/1‑150 (MQ=255)
                                                     ggTTGTCACCAGCATGACAAACGACATCATCCGCGAGAACTCTTCTTTCGGGTAAATATCGCGCATCAGGGCGTTAATGACCACGCTGGCCGCAGCCTCAGCCAGCCCGTGGAAGAAACGCATCACAATCAGCTGATCGATGGTGTTTGc                                                                                           <  2:165636/150‑1 (MQ=255)
                                                      gTTGTCACCAGCATGACAAACGACATCATCCGCGAGAACTCTTCTTTCGGGTAAATATCGCGCATCAGGGCGTTAATGACCACGCTGGCCGCAGCCGCAGCCAGCCCGTGGAAGAAACGCATCACAATCAGCTGATCGATGGTGTTTGcc                                                                                          <  2:1703818/150‑1 (MQ=255)
                                                      gTTGTCACCAGCATGACAAACGACATCATCCGCGAGAACTCTTCTTTCGGGTAAATATCGCGCATCAGGGCGTTAATGACCACGCTGGCCGCAGCCGCAGCCAGCCCGTGGAAGAAACGCATCACAATCAGCTGATCGATGGTGTTTGcc                                                                                          <  1:1254250/150‑1 (MQ=255)
                                                      gTTGTCACCAGCATGACAAACGACATCATCCGCGAGAACTCTTCTTTCGGGGAAATATCGCGCATCAGGGCGTTAATGACCACGCTGGCCGCAGCCGCAGCCAGCCCGTGGAAGAAACGCATCACAATCAGCTGATCGATGGTGTTTGcc                                                                                          <  2:1228271/150‑1 (MQ=255)
                                                       ttGTCACCAGCATGACAAACGACATCATCCGCGAGAACTCTTCTTTCGGGTAAATATCGCGCATCAGGGCGTTAATGACCACGCTGGCCGCAGCCGCAGCCAGCCCGTGGAAGAAACGCATCACAATCAGCTGATCGATGGTGTTTGCCa                                                                                         <  2:122473/150‑1 (MQ=255)
                                                         gTCACCAGCATGACAAACGACATCATCCGCGAGAACTCTTCTTTCGGGTAAATATCGCGCATCAGGGCGTTAATGACCACGCTGGCCGCAGCCGCAGCCAGCCCGTGGAAGAAACGCATCACAATCAGCTGATCGATGGTGTTTGCCAAc                                                                                       >  2:2192/1‑150 (MQ=255)
                                                           cACCAGCATGACAAACGACATCATCCGCGAGAACTCTTCTTTCGGGTAAATATCGCGCATCAGGGCGTTAATGACCACGCTGGCCGCAGCCGCAGCCAGCCCGTGGAAGAAACGCATCCCAATCAGCTGATCGATGGTGTTTGCCACCGc                                                                                     >  1:530294/1‑150 (MQ=255)
                                                            aCCAGCATGACAAACGACATCATACGCGAGAACTCTTCTTTCTGTTAAATATCGCGCAACAGTGCGATAATGACCACGCAGGCAGCAGCAGCAGCCAGCCAGTGGAAGAAACGCATCACAATCAGCTGCTCGCTCGTGTATGCCAACGcc                                                                                    >  1:294249/1‑149 (MQ=255)
                                                              cAGCATGACAAACGACATCATCCGCGAGAACTCTTCTTTCGGGTAAATATCGCGCATCAGGGCATTAATGACCACGCTGGCCGCAGCCGCAGCCAGCCCGTGGAAGAAACGCATCACAATCAGCTGATCGATGGTGTTTGCCAACGcacc                                                                                  >  1:1706389/1‑149 (MQ=255)
                                                                 cATGACAAACGACATCATCCGCGAGAACTCTTCTTTCGGGTAAATATCGCGCATCAGGGCGTTAAAGACCACGCTGGCCGCAGCCGCATCCAGCCCGTGGAAGCACCGCATCACAATCCGCTGATCGATGGTGATTGCCAACGCACACGc                                                                               >  1:270551/1‑150 (MQ=255)
                                                                      cAAACGACATCATCCGCGAGAACTCTTCTTTCGGGTAAATATCGCGCATCAGGGCGTTAATGACCACGCTGGCCGCAGCCGCAGCCAGCCCGTGGAAGAAACGCATCACAATCAGCTGATCGATGGTGTTTGCCAACGCACACGCCAccc                                                                          >  1:682374/1‑149 (MQ=255)
                                                                      cAAACGACATCATCCGCGAGAACTCTTCTTTCGGGTAAATATCGCGCATCAGGGCGTGAATGACCACGCTGGCCGCAGCCGCAGCCAGCCCGTGGAAGAAACGCCTCACAAACAGCTGATCGCTGGTGTTTGCCAACGCCCACGCCCCcg                                                                          >  1:330258/1‑150 (MQ=255)
                                                                          cGACATCATCCGCGAGAACTCTTCTTTCGGGTAAATATCGCGCATCAGGGCGTTAATGACCACGCTGGCCGCAGCCGCAGCCAGCCCGTGGAAGAAACGCATCACAATCAGCTGATCGATGGTGTTTGCCAACGCACACGCCCCCgcggc                                                                      >  1:1409147/1‑150 (MQ=255)
                                                                           gACATCATCCGCGAGAACTCTTCTTTCGGGTAAATATCGCGCATCAGGGCGTTAATGACCACGCTGGCCGCAGCCGCAGCCAGCCCGTGGAAGAAACGCATCACAATCAGCTGATCGATGGTGGTTGCCAACGCACACGCCACCgcggcg                                                                     <  1:687420/150‑1 (MQ=255)
                                                                            aCATCATCCGCGAGAACTCGTCTTTCGGGTAAATATCGCGCATCAGGGCGTTAATGCCCACGCTGGCCGCCGCCGCAGCCAGCCCGTGGAAGAAACGCATCCCCATCAGCGGATCGATGGTGTTAGACCACGCACACGCCCCCgcggcgg                                                                    >  1:324692/1‑150 (MQ=255)
                                                                                 atcCGCGAGAACTCTTCTTTCGGGTAAATATCGCGCATCAGGGCGTTAATGACCACGCTGGCCGCAGCCGCAGCCAGCCCGTGGAAGAAACGCATCACAATCAGCTGATCGATGGTGTTTGCCAACGCACACGCCACCGCGGCGGCGGCa                                                               >  1:26313/1‑150 (MQ=255)
                                                                                       gagaACTCTTCTTTCGGGTAAATATCGCGCATCAGGGCGTTAATGACCACGCTGGCCGCAGCCGCAGCCAGCCCGTGGAAGAAACGCATCACAATCAGCTGATCGATGGTGTTTGCCAACGCACACGCCACCGCGGCGGCGGCAAACAcc                                                         <  1:550891/150‑1 (MQ=255)
                                                                                        agaACTCTTCTTTCGGGTAAATATCGCGCATCAGGGCGTTAATGACCACGCTGGCCGCAGCCGCAGCCAGCCCGTGGAAGAAACGCATCACAATTAGCTGATCGATGGTGTTTGCCAACGCACACGCCCCCGCGGCGGCGGCAAACACCa                                                        >  2:439556/1‑150 (MQ=255)
                                                                                         gaACTCTTCTTTCGGGTAAATATCGCGCATCAGGGCGTTAATGTCCAGGCTGGCCGCAGCCGCAGGCAGCCCGTGGAAGAAACGCATCACAATCAGCTGATCGATGGTGTTTGCCAACGCACACGCCACCGCGGCGGCGGCAAACACCAg                                                       <  2:282638/150‑1 (MQ=255)
                                                                                            ctctTCTTTCGGGTAAATATCGCGCATCAGGGCGTTAATGACCACGCTGGCCGCAGCCGCAGCCAGCCCGTGGAAGAAACGCATCCCAATCAGCTGATCGATGGTGTTTGCCAACGCACACGCCACCGCGGCGGCGGCAAACACCAGCGt                                                    >  2:304998/1‑150 (MQ=255)
                                                                                            ctctTCTTTCGGGTAAATATCGCGAATCAGGGCGTTAATGACCACCCTGGAAGCCGCCGCAGACAGCCCGTGGAAGCAACGAATAACAAACAGCTGATCGATGGTGTTTGACAACGCACACGCCACCGCCGCGGCGGCAAACAACAGCGt                                                    >  1:955241/1‑150 (MQ=255)
                                                                                                 cttTCGGGTAAATATCGCGCATCAGGGCGTTAATGACCACGCTGGCCGCAGCCGCAGCCAGCCCGTGGAAGAAACGCATCACAATCAGCTGATCGATGGTGTTTGCCAACGCACACGCCACCGCGGCGGCGGCAAACACCAGCGTAccgc                                               >  2:669742/1‑150 (MQ=255)
                                                                                                 cttTCGGGTAAATATCGCGCATCAGGGCGTTAATGACCACGCTGGCCGCAGCCGCAGCCAGCCCGTGGAAGAAACGCATCACAATCAGCTGATCGATGGTGTTTGCCAACGCACACGCCACCGCGGCGGCGGCAAACACCAGCGTAccgc                                               >  2:1194556/1‑150 (MQ=255)
                                                                                                     cGGGTAAATATCGCGCATCAGGGCGTTAATGACCACGCTGGCCGCAGCCGCAGCCAGCCCGTGGAAGAAACGCATCACAATCAGCTGATCGATGGTGTGTGCCAACGCACACGCCACCGCGGCGGCGGCAAACACCAGCGTACCGCCGAg                                           <  1:694387/150‑1 (MQ=255)
                                                                                                       ggTAAATATCGCGCATCAGGGCGTTAATGACCACGCTGGCCGCAGCCGCAGCCAGCCCGTGGAAGAAACGCATCACAATCAGCTGATCGATGGTGTTTGCCAACGCACACGCCACCGCGGCGGCGGCAAACACCAGCGTACCGCCGAGcc                                         >  2:14485/1‑149 (MQ=255)
                                                                                                         tAAATATCGCGCATCAGGGCGTTAATGACCACGCTGGCCGCAGCCGCAGCCAGCCCGTGGAAGAAACGCATCACAATCAGCTGATCGATGGTGTTTGCCAACGCACACGCCACCGCGGCGGCGGCAAACACCAGCGTACCGCCGAGcacc                                       <  1:1187498/150‑1 (MQ=255)
                                                                                                          aagtatCGCGCATCAGGGCGTTAATGACCACGCTGGCCGCAGCCGCAGCCAGCCCGTGGAAGAAACGCATCACAATCAGCTGATCGATGGTGTTTGCCAACGCACACGCCACCGCGGCGGCGGCAAACACCAGCGTACCGCCGAGcacca                                      <  1:440460/147‑1 (MQ=255)
                                                                                                              atCGCGCATCAGGGCGTTAATGACCACGCTGGCCGCAGCCGCAGCCAGCCCGTGGAAGAAACGCATCACAATCAGCTGATCGATGGTGTTTGCCAACGCACACGCCACCgcggcggcg                                                                  <  1:496081/118‑1 (MQ=255)
                                                                                                              atCGCGCATCAGGGCGTTAATGACCACGCTGGCCGCAGCCGCAGCCAGCCCGTGGAAGAAACGCATCACAATCAGCTGATCGATGGTGTTTGCCAACGCACACGCCACCgcggcggcg                                                                  >  2:496081/1‑118 (MQ=255)
                                                                                                              atCGCGCATCAGGGCGTTAATGACCACGCTGGCCGCAGCCGCAGCCAGCCCGTGGAAGAAACGCATCACAATCAGCTGATCGATGGTGTTTGCCAACGCACACGCCACCGCGGCGGCGGCAAACACCAGCGTACCGCCGAGCACCACCgg                                  <  2:369557/150‑1 (MQ=255)
                                                                                                               tCGCGCATCAGGGCGGTAATGACCACGCTGGCCGCAGCCGCAGCCAGCCCGTGGAAGAAACGCATCACAATCAGCTGCTCGATGGTGTTTGCCAACGCACACGCCACCGCGGCGGCGGCAAACACCAGCGTATCGCCGAGCACCACCGGc                                 <  1:843760/150‑1 (MQ=255)
                                                                                                                cgcgCATCAGGGCGTTAATGACCACGCTGGCCGCAGCCGCAGCCAGCCCGGGGAAGAAACGCCTCACAATCAGCTGATCGATGGTGTTTGCCAACGCACACGCCACCGCGGCGGCGGCACACCCCAGCGTACCGCCGAGCACCACCGGCt                                >  1:313814/1‑150 (MQ=255)
                                                                                                                  cgcATCAGGGCGTTAATGACCACGCTGGCCGCAGCCGCAGCCAGCCCGTGGAAGATACGCATCACAATCAGCAGATCGATGGTGTTTGTCAACGCACACGCCACCGCGGCGGCGGCAAACACCAGCGTTCCGCCGAGCACCACCGGCTTa                              <  1:57658/150‑1 (MQ=255)
                                                                                                                   gcATCAGGGCGTTAATGACCACGCTGGCCGCAGCCGCAGCCAGCCCGTGGAAGAAACGCATCACAATCAGCTGATCGATGGTGTTTGCCAACGCACACGCCACCGCGGCGGCGGCAAACACCAGCGTACCGCCGAGCACCACCGGCTTAc                             >  2:637436/1‑150 (MQ=255)
                                                                                                                     aTCAGGGCGTTAATGACCACGCTGGCCGCAGCCGCAGCCAGCCCGTGGAAGAAACGCATCACAATCAGCTGATCGATGGTGTTTGCCAACGCACACGCCACCGCGGCGGCGGCAAACACCAGCGTACCGCCGAGCACCACCGGCTTACGc                           >  1:1425183/1‑150 (MQ=255)
                                                                                                                            cGTTAATGACCACGCTGGCCGCAGCCGCAGCCAGCCCGTGGAAGAAACGCATCACAATCAGCTGATCGATGGGGTTTGCCAACGCACACGCCACCGCGGCGGCGGCAAACACCAGCGTACCGCCGAGCACCACCGGCTTACGCCCGAAGc                    <  2:1515475/150‑1 (MQ=255)
                                                                                                                              ttAATGACCACGCTGGCCGCAGCCGCAGCCAGCCCGTGGAAGAAACGCATCACAATCAGCTGATCGATGGTGTTTGCCAACGCACACGCCACCGCGGCGGCGGCAAACACCAGCGTACCGCCGAGCACCACCGGCTTACGCCCGAAGCTg                  >  1:659539/1‑150 (MQ=255)
                                                                                                                                   gACCACGCTGGCCGCAGCCGCAGCCAGCCCGTGGAAGAAACGCATCACAATCAGCTGATCGATGGTGTTTGCCAACGCACACGCCACCGCGGCGGCGGCAAACACCCGCGTCCCGCCGACCCCCACCGGCTTACGCCCGAAGCTGTCTGc             >  2:1555859/1‑150 (MQ=255)
                                                                                                                                   gACCACGCTGGCCGCAGCCGCAGCCAGCCCGTGGAAGAAACGCATCACAATCAGCTGATCGATGGTGTTTGCCAACGCACACGCCACCGCGGCGGCGGCAAACACCAGCGTACCGCCGAGCACCCCCGGCTTACGCCCGACGCTGTCTGc             >  1:178512/1‑150 (MQ=255)
                                                                                                                                    aCCACGCTGGCCGCAGCCGCAGCCAGCCCGTGGAAGAAACGCATCACAATCAGCTGATCGATGGTGTTTGCCAACGCACACGCCACCGCGGCGGCGGCAAACACCAGCGTACCGCCGAGCACCACCGGCTTACGCCCGAAGCTGTCTGcc            >  1:452012/1‑150 (MQ=255)
                                                                                                                                     acgggctggctgCGGCTGCGGCCAGCCCGTGGAAGAAACGCATCACAATCAGCTGATCGATGGTGTTTGCCACCGCACACGCCCCCGCGGCGGCGGCAACCCCCCGCGTCCCGCCGGGCCCCCCCGGCTTACGCCCGAAGCGGTCTGccg           >  1:660700/5‑149 (MQ=255)
                                                                                                                                       aCGCTGGCCGCAGCCGCAGCCAGCCCGTGGAAGAAACGCATCACAATCAGCTGATCGATGGTGTTTGCCAACGCACACGCCACCGCGGCGGCGGCAAACACCAGCGTACCGCCGAGCACCACCGGCTTACGCCCGAAGCTGTCTGCCATc         <  2:962326/150‑1 (MQ=255)
                                                                                                                                        cGCTGGCTGCGGCTGCGGCCAGCCCGTGGAAGAAACGCATCACAATCAGCTGATCGATGGTGTTTGCCAACGCACACTCCACCGCGGCGGCGGCATCCCCCAGCGTACCGCCGAGCACCACCGGCTTACGCCCGAAGCTGTCAGCCATcc        >  1:1127291/1‑149 (MQ=255)
                                                                                                                                         gCTGGCCGCAGCCGCAGCCAGCCCGTGGAAGAAACGCATCACAATCAGCTGATCGATGGTGTTTGCCAACGCACACGCCACCGCGGCGGCGGCAAACACCAGCGTACCGCCGAGCACCACCGGCTTACGCCCGAAGCTGTCTGCCATCgg       >  2:1675013/1‑150 (MQ=255)
                                                                                                                                         gCTGGCCGCAGCCGCAGCCAGCCCGTGGAAGAAACGCATCACAATCAGCTGATCGATGGTGTTTGCCAACGCACACGCCACCGCGGCGGCGGCAAACACCAGCGTACCGCCGAGCACCACCGGCTTACGCCCGAAGCTGTCTGCCATCgg       >  2:224220/1‑150 (MQ=255)
                                                                                                                                         gCTGGCAGCGGCTGCGGCCAGCCCGTGGAAGAAACGCATCACAATCAGCTGATCGATGGTGTTTGCCCACGCACCCGCCACCGCGGCGGCGGCCACCACCCGCGTACCGCCGCGCCCCCCCGGCTTACGCCCGCAGCTGTCGGCCCTccg       >  1:599052/1‑148 (MQ=255)
                                                                                                                                          cTGGCCGCAGCCGCAGCCAGCCCGTGGAAGAAACGCATCACAATCAGCTGATCGATGGTGTTTGCCAACGCACACGCCACCGCGGCGGCGGCAAACACCAGCGTACCGCCGAGCACCACCGGCTTACGCCCGAAGCTGTCTGCCATCGGc      <  2:1505066/150‑1 (MQ=255)
                                                                                                                                          cTGGCCGCAGCCGCAGCCAGCCCGTGGAAGAAACGCATCACAATCAGCTGATCGATGGTGTTTGCCAACGCACACGCCACCGCGGCGGCGGCAAACACCAGCGTACCGCCGAGCACCACCGCCTTACGCCCGAAGCAGTCTGCCATCGGc      >  2:210987/1‑150 (MQ=255)
                                                                                                                                             gccgcagccgcaGCCAGCCCGTGGAAGAAACGCATCACAATCAGCTGATCGATGGTGTTTGCCAACGCACACGCCACCGCGGCGGCGGCAAACACCAGCGTACCGCCGAGCACCACCGGCTTACGCCCGAAGCTGTCTGCCATCGGCCCg   >  1:1129120/1‑150 (MQ=255)
                                                                                                                                              ccgcagccgcagCCAGCCCGTGGAAGAAACGCATCACAATCAGCTGATCGATGGTGTTTGCCAACGCACACGCCACCGCGGCGGCGGCAAACACCAGCGTACCGCCGAGCACCACCGGCTTACGCCCGAAGCTGTCTGCCATCGGCCCGt  >  1:595328/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CTCAGCCACACCAGCACCCAGCCGCCAACTATCGGTGCCATCAGCGGTGCAATGGTTGTCACCAGCATGACAAACGACATCATCCGCGAGAACTCTTCTTTCGGGTAAATATCGCGCATCAGGGCGTTAATGACCACGCTGGCCGCAGCCGCAGCCAGCCCGTGGAAGAAACGCATCACAATCAGCTGATCGATGGTGTTTGCCAACGCACACGCCACCGCGGCGGCGGCAAACACCAGCGTACCGCCGAGCACCACCGGCTTACGCCCGAAGCTGTCTGCCATCGGCCCGT  >  NC_000913/2279276‑2279567

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