Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 514,375 G→A G128D (GGC→GAC)  cueR → copper‑responsive regulon transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913514,3750GA98.3% 217.8 / NA 59G128D (GGC→GAC) cueRcopper‑responsive regulon transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (24/34);  minor base T (1/0);  total (25/34)
Fisher's exact test for biased strand distribution p-value = 4.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCG  >  NC_000913/514229‑514517
                                                                                                                                                  |                                                                                                                                              
cGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGactgc                                                                                                                                            >  2:113816/1‑151 (MQ=255)
cGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGactg                                                                                                                                             <  1:678759/150‑1 (MQ=255)
 gTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGactgct                                                                                                                                           <  2:249194/151‑1 (MQ=255)
   cAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGActgctg                                                                                                                                          <  1:499272/150‑1 (MQ=255)
   cAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGACTGCTGt                                                                                                                                         >  1:451048/1‑151 (MQ=255)
      aCGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGACTGCTGtcat                                                                                                                                      <  1:294053/151‑1 (MQ=255)
          cgcACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGACTGCTGTCatcatc                                                                                                                                  >  1:529006/1‑151 (MQ=255)
                  ggAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGACTGCTGTCATCATCGGGCAgg                                                                                                                           <  2:451082/150‑1 (MQ=255)
                   gagaAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGACTGCTGTCATCATCGGGCAggg                                                                                                                          >  1:415055/1‑150 (MQ=255)
                   gagaAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGACTGCTGTCATCATCGGGCAGGGt                                                                                                                         <  2:811407/151‑1 (MQ=255)
                       aGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCTCTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGACTGCTGTCATCATCGGGCAGGGTGAtt                                                                                                                     <  1:260875/151‑1 (MQ=255)
                              ggAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGACTGCTGTCATCATCGGGCAGGGTGATTAagacga                                                                                                              <  2:246893/151‑1 (MQ=255)
                                 gaTCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGACTGCTGTCATCATCGGGCAGGGTGATTAagacgagac                                                                                                           <  2:134506/151‑1 (MQ=255)
                                  aTCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGACTGCTGTCATCATCGGGCAGGGTGATTAagacgagacg                                                                                                          <  1:104446/151‑1 (MQ=255)
                                    cGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGACTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGc                                                                                                        >  1:114824/1‑151 (MQ=255)
                                     gAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGACTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCt                                                                                                       <  2:839912/151‑1 (MQ=255)
                                     gAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGACTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCt                                                                                                       <  1:535739/151‑1 (MQ=255)
                                          acacATTGAGGAGCTGCGATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCTGACTGCCCGATTATCGAAAATCTCTCCGACTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGat                                                                                                  <  1:573754/151‑1 (MQ=255)
                                           cacaTTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGACTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGat                                                                                                  >  2:689810/1‑150 (MQ=255)
                                                   ggagCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGACTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGt                                                                                         <  2:506857/151‑1 (MQ=255)
                                                     agCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGACTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGtt                                                                                        >  2:732310/1‑150 (MQ=255)
                                                           aaTCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGACTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCtt                                                                                 <  1:707684/151‑1 (MQ=255)
                                                            aTCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGACTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTc                                                                                >  1:820675/1‑151 (MQ=255)
                                                              ccATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGACTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTa                                                                              <  2:114715/151‑1 (MQ=255)
                                                                  gcgcGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGACTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCg                                                                           <  1:28474/150‑1 (MQ=255)
                                                                         cAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGACTGCTGTCATCATCGGGCAGGATGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACt                                                                   >  1:433881/1‑151 (MQ=255)
                                                                             tgctgGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGACTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCa                                                                >  1:633107/1‑150 (MQ=255)
                                                                             tgctgGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGACTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCa                                                                >  2:716447/1‑150 (MQ=255)
                                                                              gctgGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGACTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCaa                                                               <  2:725702/150‑1 (MQ=255)
                                                                               ctgGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGACTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAAc                                                              >  1:670158/1‑150 (MQ=255)
                                                                                    aCTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGACTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATGCGt                                                        <  2:117559/151‑1 (MQ=255)
                                                                                      tcgcgaccgcctgccGTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCCCCAACTGCTGTCATCATCAGGAGGGGTGATTAAGACGAGCCGGCTCGGATATGTAGGGTTATCCCTTCCATCGCCATGACTTCAACATGCGCAc                                                      >  1:330220/9‑151 (MQ=255)
                                                                                         cGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGACTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATGCGTACATg                                                   >  1:334185/1‑151 (MQ=255)
                                                                                           aaTGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGACTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATGCGTACCTGCa                                                 <  1:625235/151‑1 (MQ=255)
                                                                                           aaTGACTGCCCTGGCGATGACCGCGCCGACTGCTCGATTACCGCAAATCTCTGCGTCTGCGGTTATCTTTGGGCAGGGTGATTAAGACGAGTGGGCTCGGGTACGTAGGGTTAGCCCTTATATCGCACTGACTTCAACATGCGTAcc                                                     >  1:188339/1‑147 (MQ=255)
                                                                                              gcctgccCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGACTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATGCGTACCTGCAcc                                               <  1:235164/150‑1 (MQ=255)
                                                                                                 tgccCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGACTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCg                                              >  1:729130/1‑148 (MQ=255)
                                                                                                  gccCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGACTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGAt                                          <  1:113793/151‑1 (MQ=255)
                                                                                                       ggCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGACTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTc                                      <  1:28073/150‑1 (MQ=255)
                                                                                                          gATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGACTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTg                                  <  1:500155/151‑1 (MQ=255)
                                                                                                           aTGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGACTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTgg                                 <  1:810500/151‑1 (MQ=255)
                                                                                                                   gcagACTGCCCGATTATCGAAAATCTCTCCGACTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATGCGTACCCGCACCGAGATCCTCGCTGGCGCTGACa                         <  2:630293/148‑1 (MQ=255)
                                                                                                                     cGACTGCCCGATTATCGAAAATCTCTCCGACTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAgg                       >  1:307138/1‑151 (MQ=255)
                                                                                                                     cGACTGCCCGATTATCGAAAATCTCTCCGACTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAgg                       <  2:729173/151‑1 (MQ=255)
                                                                                                                     cGACTGCCCGATTATCGAAAATCTCTCCGACTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAgg                       <  1:672628/151‑1 (MQ=255)
                                                                                                                       aCTGCCCGATTATCGAAAATCTCTCCGACTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGcc                     >  2:799813/1‑151 (MQ=255)
                                                                                                                       aCTGCCCGATTATCGAAAATCTCTCCGACTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGcc                     <  2:777781/151‑1 (MQ=255)
                                                                                                                          gCCCGATTATCGAAAATCTCTCCGACTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATg                  >  1:694943/1‑151 (MQ=255)
                                                                                                                           cccGATTATCGAAAATCTCTCCGACTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGa                 <  2:636551/151‑1 (MQ=255)
                                                                                                                            ccGATTATCGAAAATCTCTCCGACTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGa                 >  1:179296/1‑150 (MQ=255)
                                                                                                                             cGATTATCGAAAATCTCTCCGACTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAAc               >  2:257220/1‑151 (MQ=255)
                                                                                                                              gATTATCGAAAATCTCTCCGACTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACt              <  1:632781/151‑1 (MQ=255)
                                                                                                                               aTTATCGAAAATCTCTCCGACTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTg             >  1:378039/1‑151 (MQ=255)
                                                                                                                                 tATCGAAAATCTCTCCGACTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTc           >  2:438047/1‑151 (MQ=255)
                                                                                                                                  aTCGAAAATCTCTCCGACTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTc           >  1:103227/1‑150 (MQ=255)
                                                                                                                                    cGAAAATCTCTCCGACTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCAc         <  1:723915/150‑1 (MQ=255)
                                                                                                                                    cGAAAATCTCTCCGACTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCAc         <  1:442400/150‑1 (MQ=255)
                                                                                                                                      aaaaTCTCTCCGACTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCg       >  2:565425/1‑150 (MQ=255)
                                                                                                                                       aaaTCTCTCCGACTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGAc     <  1:640078/151‑1 (MQ=255)
                                                                                                                                       aaaTCTCTCCGACTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGAc     <  1:721881/151‑1 (MQ=255)
                                                                                                                                         aTCTCTCCGACTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGAcgc   >  2:217939/1‑151 (MQ=255)
                                                                                                                                             ctcCGACTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGAcgcg  <  1:62734/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
CGTCAAACGGCGCACGCTGGAGAAGGTGGCGGAGATCGAACGACACATTGAGGAGCTGCAATCCATGCGCGACCAGCTGCTGGCACTGGCGAATGCCTGCCCTGGCGATGACAGCGCCGACTGCCCGATTATCGAAAATCTCTCCGGCTGCTGTCATCATCGGGCAGGGTGATTAAGACGAGACGGCTCGGATATGTAGGGTTATCCCTTCTATCGCAATGACTTCAACATGCGTACCTGCACCGAGATCCTCGCTGGCGCTGACAGGCCATGAACTGTCACCGACGCG  >  NC_000913/514229‑514517

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