Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 504,822 G→C 5.3% intergenic (‑126/‑92) fsr ← / → ushA putative fosmidomycin efflux system protein/bifunctional UDP‑sugar hydrolase/5'‑nucleotidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913504,8220GC 5.3% 273.8 / 5.6 76intergenic (‑126/‑92)fsr/ushAputative fosmidomycin efflux system protein/bifunctional UDP‑sugar hydrolase/5'‑nucleotidase
Reads supporting (aligned to +/- strand):  ref base G (36/36);  new base C (2/2);  total (38/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.72e-01

CACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTGTATGCGTCGTCTGAAGGAAGTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTAC  >  NC_000913/504705‑504963
                                                                                                                     |                                                                                                                                             
cACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTGTATGCGTCGTCTGAAGGAAGTCTCAACGCCGAATACAGAATTTCTAATCTGGa                                                                                                                                              <  3:60738/119‑1 (MQ=255)
cACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTGTATGCGTCGTCTGAAGGAAGTCTCAACGCCGAATACAGAATTTCTAATCTGGa                                                                                                                                              >  4:60738/1‑119 (MQ=255)
  cTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTGTATGCGTCGTCTGAAGGAAGTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTAt                                                                                                                                   <  4:804355/128‑1 (MQ=255)
  cTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTGTATGCGTCGTCTGAAGGAAGTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTAt                                                                                                                                   >  3:804355/1‑128 (MQ=255)
   tttttCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTGTATGCGTCGTCTGAAGGAAGTCTCAACGCCGAATACAGAATTTCTAATCTGCATCCAGAttaga                                                                                                                                   >  4:882657/1‑124 (MQ=255)
   tttttCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTGTATGCGTCGTCTGAAGGAAGTCTCAACGCCGAATACAGAATTTCTAATCTGCATCCAGAttaga                                                                                                                                   <  3:882657/127‑4 (MQ=255)
                 gATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTGTATGCGTCGTCTGAAGGAAGTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTc                                                                                                             >  2:28330/1‑135 (MQ=255)
                 gATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTGTATGCGTCGTCTGAAGGAAGTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTc                                                                                                             <  1:28330/135‑1 (MQ=255)
                 gATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTGTATGCGTCGTCTGAAGGAAGTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGt                                                                                                >  1:38158/1‑148 (MQ=255)
                    gTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTGTATGCGTCGTCTGAAGGAAGTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTc                                                                                                                               <  4:827072/114‑1 (MQ=255)
                    gTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTGTATGCGTCGTCTGAAGGAAGTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTc                                                                                                                               >  3:827072/1‑114 (MQ=255)
                        aGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTGTATGCGTCGTCTGAAGGAAGTCTCAACGCCGAATACAGAATTTCTAATCTGGa                                                                                                                                              >  4:889030/1‑95 (MQ=255)
                        aGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTGTATGCGTCGTCTGAAGGAAGTCTCAACGCCGAATACAGAATTTCTAATCTGGa                                                                                                                                              <  3:889030/95‑1 (MQ=255)
                               tCTTATTTGTGATTATTACCAGACTAACATACCTGTATGCGTCGTCTGAAGGAAGTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCa                                                                                  >  2:1113529/1‑148 (MQ=255)
                                cTTATTTGTGATTATTACCAGGCTAACATACCTGTATGCGTCGTCTGAAGGAAGTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCaa                                                                                 >  3:370392/1‑148 (MQ=255)
                                  tATTTGTGATTATTACCAGACTAACATACCTGTATGCGTCGTCTGAAGGAAGTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGa                                                                                                                         >  1:324406/1‑106 (MQ=255)
                                  tATTTGTGATTATTACCAGACTAACATACCTGTATGCGTCGTCTGAAGGAAGTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGa                                                                                                                         <  2:324406/106‑1 (MQ=255)
                                  tATTTGTGATTATTACCAGACTAACATACCTGTATGCGTCGTCTGAAGGAAGTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATa                                                                                            >  4:52351/1‑135 (MQ=255)
                                  tATTTGTGATTATTACCAGACTAACATACCTGTATGCGTCGTCTGAAGGAAGTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATa                                                                                            <  3:52351/135‑1 (MQ=255)
                                         gATTATTACCAGACTAACATACCTGTATGCGTCGTCTGAAGGAAGTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCa                                                                                              >  1:1025569/1‑126 (MQ=255)
                                         gATTATTACCAGACTAACATACCTGTATGCGTCGTCTGAAGGAAGTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCa                                                                                              <  2:1025569/126‑1 (MQ=255)
                                         gATTATTACCAGACTAACATACCTGTATGCGTCGTCTGAAGGAAGTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGaa                                                                        <  2:38158/148‑1 (MQ=255)
                                            tattaCCAGACTAACATACCTGTATGCGTCGTCTGAAGGAAGTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCa                                                                                                                              <  1:10919/91‑1 (MQ=255)
                                            tattaCCAGACTAACATACCTGTATGCGTCGTCTGAAGGAAGTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCa                                                                                                                              >  2:10919/1‑91 (MQ=255)
                                                   aGACTAACATACCTGTATGCGTCGTCTGAAGGAAGTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTcgc                                                                                                  <  3:997940/112‑1 (MQ=255)
                                                   aGACTAACATACCTGTATGCGTCGTCTGAAGGAAGTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTcgc                                                                                                  >  4:997940/1‑112 (MQ=255)
                                                        aaCATACCTGTATGCGTCGTCTGAAGGAAGTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTg                                                                                                               >  1:660664/1‑94 (MQ=255)
                                                        aaCATACCTGTATGCGTCGTCTGAAGGAAGTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTg                                                                                                               <  2:660664/94‑1 (MQ=255)
                                                             aCCTGTATGCGTCGTCTGAAGGAAGTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGt                                                    <  1:1113529/148‑1 (MQ=255)
                                                             aCCTGTATGCGTCGTCTGAAGGAAGTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGt                                                    >  3:491938/1‑148 (MQ=255)
                                                              ccTGTATGCGTCGTCTGAAGGAAGTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCa                                                                                  >  2:267513/1‑117 (MQ=255)
                                                              ccTGTATGCGTCGTCTGAAGGAAGTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCa                                                                                  <  1:267513/117‑1 (MQ=255)
                                                                tggaTGCGTCGTCTGAAGGGAGTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATa                                                                                            >  2:220732/4‑105 (MQ=255)
                                                                tGTATGCGTCGTCTGAAGGAAGTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATa                                                                                            <  1:220732/105‑1 (MQ=255)
                                                                tGTATGCGTCGTCTGAAGGAAGTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACg                                                                             <  3:444170/120‑1 (MQ=255)
                                                                tGTATGCGTCGTCTGAAGGAAGTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACg                                                                             >  4:444170/1‑120 (MQ=255)
                                                                  tATGCGTCGTCTGAAGGAAGTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAGCACGTTGAAATCAGGTCAGGGAGAGAAGTATGaa                                               <  4:370392/148‑1 (MQ=255)
                                                                      cgtcgtCTGAAGGAAGTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGgagaga                                                       >  4:203977/1‑136 (MQ=255)
                                                                      cgtcgtCTGAAGGAAGTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGgagaga                                                       <  3:203977/136‑1 (MQ=255)
                                                                            cTGAAGGAAGTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTatgat                                                                                                       >  3:838293/1‑82 (MQ=255)
                                                                            cTGAAGGAAGTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTatgat                                                                                                       <  4:838293/82‑1 (MQ=255)
                                                                              gaaggaagTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAAtca                                                                    <  1:824627/115‑1 (MQ=255)
                                                                              gaaggaagTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAAtca                                                                    >  2:824627/1‑115 (MQ=255)
                                                                              gaaggaagTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAAt                                                                      <  2:58098/113‑1 (MQ=255)
                                                                              gaaggaagTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAAt                                                                      >  1:58098/1‑113 (MQ=255)
                                                                               aaggaagTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAAtcag                                                                   <  2:37919/115‑1 (MQ=255)
                                                                               aaggaagTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAAtcag                                                                   >  1:37919/1‑115 (MQ=255)
                                                                                aggaagTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCa                                                                                              <  1:138133/87‑1 (MQ=255)
                                                                                aggaagTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCa                                                                                              >  2:138133/1‑87 (MQ=255)
                                                                                   aagTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGt                                                                                                <  2:157671/82‑1 (MQ=255)
                                                                                   aagTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGt                                                                                                >  1:157671/1‑82 (MQ=255)
                                                                                   aagTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGt                                                                            >  1:118868/1‑102 (MQ=255)
                                                                                   aagTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGt                                                                            <  2:118868/102‑1 (MQ=255)
                                                                                      tctcAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGgaga                                                         >  1:1017004/1‑118 (MQ=255)
                                                                                      tctcAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGgaga                                                         <  2:1017004/118‑1 (MQ=255)
                                                                                       ctcAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAatta                                           >  3:19044/1‑131 (MQ=255)
                                                                                       ctcAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAatta                                           <  4:19044/131‑1 (MQ=255)
                                                                                         cAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTatga                                                                                                        >  2:440566/1‑68 (MQ=255)
                                                                                         cAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTatga                                                                                                        <  1:440566/68‑1 (MQ=255)
                                                                                         cAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAAtca                                                                    <  2:699253/104‑1 (MQ=255)
                                                                                         cAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAAtca                                                                    >  1:699253/1‑104 (MQ=255)
                                                                                          aaCGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGtt                       <  4:491938/148‑1 (MQ=255)
                                                                                            cGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCa                                                                                  <  3:877736/87‑1 (MQ=255)
                                                                                            cGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCa                                                                                  >  4:877736/1‑87 (MQ=255)
                                                                                                  aTACAGAATTTCTAATCTGCATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCgg                                  >  2:224959/1‑129 (MQ=255)
                                                                                                  aTACAGAATTTCTAATCTGCATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCgg                                  <  1:224959/129‑1 (MQ=255)
                                                                                                    aCAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGt                                                                            >  1:990539/1‑85 (MQ=255)
                                                                                                    aCAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGt                                                                            <  2:990539/85‑1 (MQ=255)
                                                                                                    aCAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAattat                                          >  1:519086/1‑119 (MQ=255)
                                                                                                    aCAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAattat                                          <  2:519086/119‑1 (MQ=255)
                                                                                                       gAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCgg                                  >  2:43463/1‑124 (MQ=255)
                                                                                                       gAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCgg                                  <  1:43463/124‑1 (MQ=255)
                                                                                                        aaTTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTAt                                                  >  3:937598/1‑107 (MQ=255)
                                                                                                        aaTTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTAt                                                  <  4:937598/107‑1 (MQ=255)
                                                                                                           ttCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAAtcaggtcag                                                              >  4:237340/1‑92 (MQ=255)
                                                                                                           ttCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAAtcaggtcag                                                              <  3:237340/92‑1 (MQ=255)
                                                                                                                     gATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTac  <  4:501384/142‑1 (MQ=255)
                                                                                                                     gATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTac  >  3:501384/1‑142 (MQ=255)
                                                                                                                     |                                                                                                                                             
CACTTTTTCAGGGTTGCGATGTAAAGAATGATCTTATTTGTGATTATTACCAGACTAACATACCTGTATGCGTCGTCTGAAGGAAGTCTCAACGCCGAATACAGAATTTCTAATCTGGATGCAGATTTATCTTCACCGGACGCAGACTTGTCTATGATGTCGCGTCATACTATTTTTCAACACGTTGAAATCAGGTCAGGGAGAGAAGTATGAAATTATTGCAGCGGGGCGTGGCGTTAGCGCTGTTAACCACATTTAC  >  NC_000913/504705‑504963

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