Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 295,247 G→A 7.2% G217E (GGG→GAG)  ssuB → aliphatic sulfonates ABC transporter ATP‑binding protein SsuB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947295,2470GA 7.2% 204.8 / 4.9 69G217E (GGG→GAG) ssuBaliphatic sulfonates ABC transporter ATP‑binding protein SsuB
Reads supporting (aligned to +/- strand):  ref base G (34/30);  new base A (2/3);  total (36/33)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CTGGATGCGTTGACCCGTATCGAGATGCAGCAACTGATCGAACGCCTGTGGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTGTCCACCCCGGGCACTGCGCCCGAGCCGGATCCTGTAGCCCCTCTGCCCACGCAACTGC  >  NC_002947/295099‑295396
                                                                                                                                                    |                                                                                                                                                     
cTGGTTGCGTTGACCGGTTTCGAGATGCAGCAACTGATCGAACGCCTGTGGCGTCAGCACGGCTTCACTGTACTGCTCGTCACCCGCGACGTCAGCGAGGCGGTTGCCGTGGCTGTCCGGGTGATCCTGATCGATGACGGCGAGGTCggg                                                                                                                                                      <  2:28446/150‑1 (MQ=255)
 tGGATGCGTTGACCCGTATCGAGATGCAGCAACTGATCGAACGCCTGTGGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGAGc                                                                                                                                                     <  2:709197/150‑1 (MQ=255)
  ggATGCGTTGACCCGTATCGAGATGCAGCAACTGATCGAACGCCTGTGGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCt                                                                                                                                                    >  1:869534/1‑150 (MQ=255)
             cccGTATCGAGATGCAGCAACTGATCGAACGCCTGTGGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTg                                                                                                                                         <  2:1106964/150‑1 (MQ=255)
              ccGTATCGAGATGCAGCAACTGATCGAACGCCTGTGGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGt                                                                                                                                        >  1:435890/1‑150 (MQ=255)
              ccGTATCGAGATGCAGCAACTGATCGAACGCCTGTGGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGt                                                                                                                                        >  2:415589/1‑150 (MQ=255)
              ccGTATCGAGATGCAGCAACTGATCGAACGCCTGTGGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGAAGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGt                                                                                                                                        >  1:1051058/1‑150 (MQ=255)
                     gagaTGCAGCAACTGATCGAACGCCTGTGGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCAATAATGACCTg                                                                                                                                 >  1:653923/1‑150 (MQ=255)
                      agaTGCAGCAACTGATCGAACGCCTGTGGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCCGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTgg                                                                                                                                >  2:41137/1‑150 (MQ=255)
                      agaTGCAGCAACTGATCGAACGCCTGTGGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGAGCCCGACCTcgc                                                                                                                                           <  1:649939/139‑3 (MQ=255)
                      agaTGCAGCAACTGATCGAACGCCTGTGGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGAGCCCGACCTcgc                                                                                                                                           >  2:649939/1‑137 (MQ=255)
                          gcagcaACTGATCGAACGCCTGTGGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACg                                                                                                                            >  2:1050059/1‑150 (MQ=255)
                             gcaACTGATCGAACGCCTGTGGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGcc                                                                                                                         >  1:599144/1‑150 (MQ=255)
                              caACTGATCGAACGCCTGTGGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGTTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGccc                                                                                                                        >  2:428932/1‑150 (MQ=255)
                              caACTGATCGAACGCCTGTGGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCAGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACATGGCACGGccc                                                                                                                        <  2:273511/150‑1 (MQ=255)
                                aCTGATCGAACGCCTGTGGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCgg                                                                                                                                                       <  2:628281/117‑1 (MQ=255)
                                aCTGATCGAACGCCTGTGGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCgg                                                                                                                                                       >  1:628281/1‑117 (MQ=255)
                                     tCGAACGCCTGTGGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGtcgc                                                                                                                 >  2:440875/1‑150 (MQ=255)
                                       gAACGCCTGTGGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGgcgc                                                                                                                 >  2:930814/1‑148 (MQ=255)
                                       gAACGCCTGTGGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGgcgc                                                                                                                 <  1:930814/148‑1 (MQ=255)
                                        aaCGCCTGTGGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTg                                                                                                              <  2:599144/150‑1 (MQ=255)
                                              tgtgGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTc                                                                                                                                                   >  1:357575/1‑107 (MQ=255)
                                              tgtgGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTc                                                                                                                                                   <  2:357575/107‑1 (MQ=255)
                                               gtgGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCa                                                                                                                              <  1:414819/127‑1 (MQ=255)
                                               gtgGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCa                                                                                                                              >  2:414819/1‑127 (MQ=255)
                                                    gTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCCAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTc                                                                                                  <  2:869534/150‑1 (MQ=255)
                                                       aGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTgg                                                                                               <  1:1050059/150‑1 (MQ=255)
                                                        gCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGc                                                                                              >  1:1000563/1‑150 (MQ=255)
                                                         cACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGt                                                                                                            <  1:1047872/135‑1 (MQ=255)
                                                         cACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGt                                                                                                            >  2:1047872/1‑135 (MQ=255)
                                                            ggCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCcgcg                                                                                          >  2:1171637/1‑150 (MQ=255)
                                                               ttCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCAcc                                                                                                     <  1:121304/136‑1 (MQ=255)
                                                               ttCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCAcc                                                                                                     >  2:121304/1‑136 (MQ=255)
                                                                      tACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGAc                                                                                                                                                >  1:1107301/1‑86 (MQ=255)
                                                                      tACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGAc                                                                                                                                                <  2:1107301/86‑1 (MQ=255)
                                                                      tACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCg                                                                                >  1:583568/1‑150 (MQ=255)
                                                                       aCTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGa                                                                               >  1:1076188/1‑150 (MQ=255)
                                                                         tgctgGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTc                                                                                                                                                   >  2:179133/1‑80 (MQ=255)
                                                                         tgctgGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTc                                                                                                                                                   <  1:179133/80‑1 (MQ=255)
                                                                         tgctgGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAg                                                                             >  2:990596/1‑150 (MQ=255)
                                                                          gctgGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGt                                                                            <  2:1051058/150‑1 (MQ=255)
                                                                            tgGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTgg                                                                                                             <  1:832551/115‑1 (MQ=255)
                                                                            tgGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTgg                                                                                                             >  2:832551/1‑115 (MQ=255)
                                                                               tCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCt                                                                                        >  2:139100/1‑133 (MQ=255)
                                                                               tCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCt                                                                                        <  1:139100/133‑1 (MQ=255)
                                                                                   ccACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTg                                                                                                              >  2:472819/1‑107 (MQ=255)
                                                                                   ccACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTg                                                                                                              <  1:472819/107‑1 (MQ=255)
                                                                                   ccACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCa                                                                                                       <  2:962217/114‑1 (MQ=255)
                                                                                   ccACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCa                                                                                                       >  1:962217/1‑114 (MQ=255)
                                                                                   ccACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTg                                                                                                <  1:1010595/121‑1 (MQ=255)
                                                                                   ccACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTg                                                                                                >  2:1010595/1‑121 (MQ=255)
                                                                                    cACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCgt                                                                  >  2:248655/1‑150 (MQ=255)
                                                                                    cACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCgt                                                                  >  2:738031/1‑150 (MQ=255)
                                                                                    cACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCgt                                                                  <  2:243070/150‑1 (MQ=255)
                                                                                     acgacgTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCgtg                                                                 >  1:1117358/1‑150 (MQ=255)
                                                                                          gTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTg                                                            <  1:1191142/150‑1 (MQ=255)
                                                                                                  ggCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTGTCCAcccc                                                    <  1:415589/150‑1 (MQ=255)
                                                                                                    cGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTGTCCACCCCgg                                                  <  1:41137/150‑1 (MQ=255)
                                                                                                    cGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTGTCCACCCCgg                                                  <  2:583568/150‑1 (MQ=255)
                                                                                                                 tGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTGTCCACCCCGGGCACTGCGCCCGa                                     >  1:1147750/1‑150 (MQ=255)
                                                                                                                       ggTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGgc                                                                                                                   <  1:548976/66‑1 (MQ=255)
                                                                                                                       ggTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGgc                                                                                                                   >  2:548976/1‑66 (MQ=255)
                                                                                                                            tCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTGTCCACCCCGGGCACTGCGCCCGAGCCGGATCCTg                          <  2:480031/150‑1 (MQ=255)
                                                                                                                                 aTCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGaa                                                                              <  1:720710/93‑1 (MQ=255)
                                                                                                                                 aTCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGaa                                                                              >  2:720710/1‑93 (MQ=255)
                                                                                                                                     aGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTGTCCACCCCGGGCACTGCGCCCGAGCCGGATCCTGTAGCCCctc                 <  2:1076188/150‑1 (MQ=255)
                                                                                                                                        aCGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTGTCCACCCCGGGCACTGCGCCCGAGCCGGATCCTGTAGCCCCTCTGc              <  1:428932/150‑1 (MQ=255)
                                                                                                                                         cGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTGTCCACCCCGGGCAc                                              <  2:788778/117‑1 (MQ=255)
                                                                                                                                         cGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTGTCCACCCCGGGCAc                                              >  1:788778/1‑117 (MQ=255)
                                                                                                                                            tgAGGTCGAGCCCGACCTCGCTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGaaa                                                                                   <  1:181213/76‑1 (MQ=38)
                                                                                                                                            tgAGGTCGAGCCCGACCTCGCTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGaaa                                                                                   >  2:181213/2‑77 (MQ=38)
                                                                                                                                            ggAGGTCGGGCTCGACCTCACTGTTGACCTGGCCCGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTGTCCACCCCGGGCACTGCGCCCGAGCCGGATCCTGTAGCCCCTCTGCCCAc          <  1:440875/149‑1 (MQ=255)
                                                                                                                                                    ggCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTGTCCACCCCGGGCACTGCGCCCGAGCCGGATCCTGTAGCCCCTCTGCCCACGCAACTgc  >  1:1040354/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTGGATGCGTTGACCCGTATCGAGATGCAGCAACTGATCGAACGCCTGTGGCGTCAGCACGGCTTCACTGTACTGCTGGTCACCCACGACGTCAGCGAGGCGGTTGCCGTGGCTGACCGGGTGATCCTGATCGAGGACGGCGAGGTCGGGCTCGACCTCACTGTTGACCTGGCACGGCCCCGGGCGCGTGGTTCGCACCGTCTGGCCGCGCTGGAAAGCGAAGTGCTCAACCGTGTTCTGTCCACCCCGGGCACTGCGCCCGAGCCGGATCCTGTAGCCCCTCTGCCCACGCAACTGC  >  NC_002947/295099‑295396

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