Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,639,231 T→C 7.3% intergenic (+90/+154) pitA → / ← uspB phosphate transporter, low‑affinity; tellurite importer/universal stress (ethanol tolerance) protein B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,639,2310TC 7.3% 229.1 / 6.3 68intergenic (+90/+154)pitA/uspBphosphate transporter, low‑affinity; tellurite importer/universal stress (ethanol tolerance) protein B
Reads supporting (aligned to +/- strand):  ref base T (29/34);  new base C (2/3);  total (31/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

TTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATT  >  NC_000913/3639083‑3639374
                                                                                                                                                    |                                                                                                                                               
ttCGGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCtt                                                                                                                                                 >  4:546459/1‑149 (MQ=255)
ttCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCtt                                                                                                                                                 >  2:161921/1‑149 (MQ=255)
ttCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCtt                                                                                                                                                 >  3:337659/1‑149 (MQ=255)
   cGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATc                                                                                                                                              <  4:337659/149‑1 (MQ=255)
      cTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                                                           <  1:161921/149‑1 (MQ=255)
       tGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                          <  3:546459/149‑1 (MQ=255)
         cGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCt                                                                                                                                        >  2:1041527/1‑149 (MQ=255)
               tGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACaaaa                                                                                                                                  >  3:665371/1‑149 (MQ=255)
                             cTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCaa                                                                                                                    <  4:665371/149‑1 (MQ=255)
                             cTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCaa                                                                                                                    <  1:448979/149‑1 (MQ=255)
                              tGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAAt                                                                                                                   <  1:740486/149‑1 (MQ=255)
                                      ctctCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAATAACATGCAAATTCaa                                                                                                                    <  1:1057575/140‑1 (MQ=255)
                                      ctctCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGAGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGTAAATTCaa                                                                                                                    >  2:1057575/1‑140 (MQ=255)
                                        ctcCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCaga                                                                                                         >  2:890792/1‑149 (MQ=255)
                                        ctcCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCaga                                                                                                         <  2:746616/149‑1 (MQ=255)
                                                aGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGt                                                                                                 <  1:1041527/149‑1 (MQ=255)
                                                  ttCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAg                                                                                               >  2:1337413/1‑149 (MQ=255)
                                                  ttCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAg                                                                                               <  1:1337413/149‑1 (MQ=255)
                                                              tACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAgaga                                                                                                       >  3:1319030/1‑129 (MQ=255)
                                                              tACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAgaga                                                                                                       <  4:1319030/129‑1 (MQ=255)
                                                                cGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTa                                                                                 <  1:890792/149‑1 (MQ=255)
                                                                   aCCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTAt                                                                                                   >  1:756767/1‑128 (MQ=255)
                                                                   aCCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTAt                                                                                                   <  2:756767/128‑1 (MQ=255)
                                                                    ccAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGAt                                                                                                               >  4:41963/1‑115 (MQ=255)
                                                                    ccAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGAt                                                                                                               <  3:41963/115‑1 (MQ=255)
                                                                     cAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGa                                                                                         >  3:231021/1‑136 (MQ=255)
                                                                     cAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGa                                                                                         <  4:231021/136‑1 (MQ=255)
                                                                     cAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGc                                                                                    <  2:556803/141‑1 (MQ=255)
                                                                     cAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGc                                                                                    >  1:556803/1‑141 (MQ=255)
                                                                      aaaaCGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTa                                                                                                <  3:197979/128‑1 (MQ=255)
                                                                      aaaaCGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTa                                                                                                >  4:197979/1‑128 (MQ=255)
                                                                           gAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAg                                                                                               <  2:234709/124‑1 (MQ=255)
                                                                           gAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAg                                                                                               >  1:234709/1‑124 (MQ=255)
                                                                                 gTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCa                                                                                                           <  2:228682/106‑1 (MQ=255)
                                                                                 gTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCa                                                                                                           >  1:228682/1‑106 (MQ=255)
                                                                                        ggCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATaa                                                                                      <  2:341781/120‑1 (MQ=255)
                                                                                        ggCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATaa                                                                                      >  1:341781/1‑120 (MQ=255)
                                                                                               ttCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGttt                                                          <  4:540420/141‑1 (MQ=255)
                                                                                               ttCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGttt                                                          >  3:540420/1‑141 (MQ=255)
                                                                                                 cAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGtc                                                >  3:754823/1‑149 (MQ=255)
                                                                                                     ttGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCaa                                                                                                                    <  4:49570/77‑1 (MQ=255)
                                                                                                     ttGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCaa                                                                                                                    >  3:49570/1‑77 (MQ=255)
                                                                                                     ttGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGa                                            <  4:754823/149‑1 (MQ=255)
                                                                                                          aCATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAAtt                                                                  <  4:154137/122‑1 (MQ=255)
                                                                                                          aCATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAAtt                                                                  >  3:154137/1‑122 (MQ=255)
                                                                                                              aGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAg                                                                                               >  4:1294228/1‑89 (MQ=255)
                                                                                                              aGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAg                                                                                               <  3:1294228/89‑1 (MQ=255)
                                                                                                                  cgcTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCg                                                             >  1:456800/1‑119 (MQ=255)
                                                                                                                  cgcTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCg                                                             <  2:456800/119‑1 (MQ=255)
                                                                                                                        gTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAg                                                                                               >  3:686262/1‑79 (MQ=255)
                                                                                                                        gTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAg                                                                                               <  4:686262/79‑1 (MQ=255)
                                                                                                                             tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGAtt                                                                                                              >  4:745558/1‑59 (MQ=255)
                                                                                                                             tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGAtt                                                                                                              <  3:745558/59‑1 (MQ=255)
                                                                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCa                                                                            <  1:1238698/92‑1 (MQ=17)
                                                                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCa                                                   <  4:1292289/117‑1 (MQ=17)
                                                                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCa                                                   >  3:1292289/1‑117 (MQ=17)
                                                                                                                                  gATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATg                                                                                                  >  3:1121279/1‑66 (MQ=255)
                                                                                                                                  gATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATg                                                                                                  <  4:1121279/66‑1 (MQ=255)
                                                                                                                                   ataaggcgTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAAt                                                     <  2:1403201/106‑1 (MQ=17)
                                                                                                                                   ataaggcgTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAAt                                                     >  1:1403201/5‑110 (MQ=17)
                                                                                                                                          gTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGtta       >  3:915621/1‑149 (MQ=255)
                                                                                                                                             aaCGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAg                                                                        <  3:87125/81‑1 (MQ=255)
                                                                                                                                             aaCGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAg                                                                        >  4:87125/1‑81 (MQ=255)
                                                                                                                                              aCGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGttt                                                          <  3:856808/94‑1 (MQ=255)
                                                                                                                                              aCGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGttt                                                          >  4:856808/1‑94 (MQ=255)
                                                                                                                                                gCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGCCAGtaaact                             >  4:148434/1‑116 (MQ=255)
                                                                                                                                                gCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGCCAGtaaact                             <  3:148434/121‑6 (MQ=255)
                                                                                                                                                 ccTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTattatt  >  2:163124/1‑147 (MQ=255)
                                                                                                                                                 ccTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTattatt  <  1:163124/147‑1 (MQ=255)
                                                                                                                                                  cTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTAc                              <  3:757780/118‑1 (MQ=255)
                                                                                                                                                  cTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTAc                              >  4:757780/1‑118 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
TTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATT  >  NC_000913/3639083‑3639374

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: