Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,565,730 T→C T41A (ACC→GCC)  yjiM ← putative dehydratase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,565,7300TC98.4% 207.1 / NA 61T41A (ACC→GCC) yjiMputative dehydratase subunit
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (31/29);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTAGCTGCTTTTAATCAGCGGGCAGAGGTTGCGCGGCAGATCTTTCTCCGCTTCTTCAATGGTTTCATCAGAGGTGGAACAGAGCGAAACCACAACCGCACCGGCTGCCATCGGGATCTCTTGCGGCATAAAGGTGCAGTAAGTGCCAACCAGCGGAATGCCGCGCTCCTTGAGATCCATGACGGTGAGAAAGCCTGTCTGGCGAGCTTCAGAGAACTGATCGAAAATGGCGGGTAGATCGGTGACAAGTGACATGATTTTCCTTCCCCG  >  NC_000913/4565595‑4565865
                                                                                                                                       |                                                                                                                                       
cGTAGCTGCTTTTAATCAGCGGGCAGAGGTTGCGCGGCAGATCTTTCTCCGCTTCTTCAATGGTTTCATCAGAGGTGGAACAGAGCGAAACCACAACCGCACAGGCTGCCTTCGAGATCTCTTGCGGCATAAAGGCGCAg                                                                                                                                     >  1:300033/1‑140 (MQ=255)
 gTAGCTGCTTTTAATCAGCGGGCAGAGGTTGCGCGGCAGATCTTTCTCCGCTTCTTCAATGGTTTCATCAGAGGTGGAACAGAGCGAAACCACAACCGCACCGGCTGCCATCGGGATCTCTTGCGGCATAAAGGCGCAGt                                                                                                                                    <  1:76729/140‑1 (MQ=255)
        ctTTTAATCAGCGGGCAGAGGTTGCGCGGCAGATCTTTCTCCGCTTCTTCAATGGTTTCATCAGAGGTGGAACAGAGCGAAACCACAACCGCACCGGCTGCCATCGGGATCTCTTGCGGCATAAAGGCGCAGTAAGTGcc                                                                                                                             >  1:2412/1‑140 (MQ=255)
         ttttAATCAGCGGGCAGAGGTTGCGCGGCAGATCTTTCTCCGCTTCTTCAATGGTTTCATCAGAGGCGGAACAGAGCGAAACCACAACCGCACCGGCTGCCATCCGGAGCTCTTGCGGCATAAAGGCGCAGTAAGTGCCa                                                                                                                            >  1:34230/1‑140 (MQ=255)
             aaTCAGCGGGCAGAGGTTGCGCGGCAGATCTTTCTCCGCTTCTTCAATGGTTTCATCAGAGGTGGAACAGAGCGAAACCACAACCGCACCGGCTGCCATCGGGATCTCTTGCGGCATAAAGGCGCAGTAAGTGCCAACCa                                                                                                                        <  1:2338/140‑1 (MQ=255)
                cAGCGGGCAGAGGTTGCGCGGCAGATCTTTCTCCGCTTCTTCAATGGTTTCATCAGAGGTGGAACAGAGCGAAACCACAACCGCACCGGCTGCCATCGGGATCTCTTGCGGCATAAAGGCGCAGTAAGTGCCAACCAGCg                                                                                                                     <  2:102488/140‑1 (MQ=255)
                         gagGTTGCGCGGCAGATCTTTCTCCGCTTCTTCAATGGTTTCATCAGAGGTGGAACAGAGCGAAACCACAACCGCACCGGCTGCCATCGGGATCTCTTGCGGCATAAAGGCGCAGTAAGTGCCAACCAGCGGAATGCcgc                                                                                                            <  2:8079/140‑1 (MQ=255)
                          agGTTGCGCGGCAGATCTTTCTCCGCTTCTTCAATGGTTTCATCAGAGGTGGAACAGAGCGAAACCACAACCGCACCGGCTGCCATCGGGATCTCTTGCGGCATAAAGGCGCAGTAAGTGCCAACCAGCGGAATGCcgcg                                                                                                           <  1:563539/140‑1 (MQ=255)
                          agGTTGCGCGGCAGATCTTTCTCCGCTTCTTCAATGGTTTCATCAGAGGTGGAACAGAGCGAAACCACAACCGCACCGGCTGCCATCGGGATCTCTTGCGGCATAAAGGCGCAGTAAGTGCCAACCAGCGGAATGCcgcg                                                                                                           <  2:131748/140‑1 (MQ=255)
                           ggTTGCGCGGCAGATCTTTCTCCGCTTCTTCAATGGTTTCATCAGAGGTGGAACAGAGCGAAACCACAACCGCACCGGCTGCCATCGGGATCTCTTGCGGCATAAAGGCGCAGTAAGTGCCAACCAGCGGAATGCCgcgc                                                                                                          <  1:126685/140‑1 (MQ=255)
                               gcgcGGCAGATCTTTCTCCGCTTCTTCAATGGTTTCATCAGAGGTGGAACAGAGCGAAACCACAACCGCACCGGCTGCCATCGGGATCTCTTGCGGCATAAAGGCGCAGTAAGTGCCAACCAGCGGAATGCCGCGCTCCt                                                                                                      >  1:129592/1‑140 (MQ=255)
                                cgcgGCAGATCTTTCTCCGCTTCTTCAATGGTTTCATCAGAGGTGGAACAGAGCGAAACCACAACCGCACCGGCTGCCATCGGGATCTCTTGCGGCATAAAGGCGCAGTACGTGACAGCCAGCGAAATGCCGCGCTCCtt                                                                                                     >  1:247976/1‑140 (MQ=255)
                                cgcgGCAGATCTTTCTCCGCTTCTTCAATGGTTTCATCAGAGGTGGAACAGAGCGAAACCACAACCGCACCGGCTGCCATCGGGATCTCTTGCGGCATAAAGGCGCAGTAAGTGCCAACCAGCGGAATGCCGCGCTCCtt                                                                                                     >  1:126392/1‑140 (MQ=255)
                                cgcgGCAGATCTTTCTCCGCTTCTTCAATGGTTTCATCAGAGGTGGAACAGAGCGAAACCACAACCGCACCGGCTGCCATCGGGATCTCTTGCGGCATAAAGGCGCAGTAAGTGCCAACCAGCGGAATGCCGCGCTCCtt                                                                                                     >  2:17563/1‑140 (MQ=255)
                                 gcgGCAGATCTTTCTCCGCTTCTTCAATGGTTTCATCAGAGGTGGAACAGAGCGAAACCACAACCGCACCGGCTGCCATCGGGATCTCTTGCGGCATAAAGGCGCAGTAAGTGCCAACCAGCGGAATGCCGCGCTCCTTg                                                                                                    <  1:529527/140‑1 (MQ=255)
                                 gcgGCAGATCTTTCTCCCCTTCGTCAATGGTGTCAGCAGTGGTGGAACTGAGCGAAACCACAACCGCACCGGCTGCCATCGGGATCTCTTGCGGCATAAAGGCGCAGTAAGTGCCAACCAGCGGAATGCCGCGCTCCTTg                                                                                                    <  2:235319/140‑1 (MQ=255)
                                 gcgGCAGACCTTTCTCTACTTCTTCAATGGATTTATCAGAGGTGGAACCGAGCGAAACCACAACCGCACCGGCTGCCATCGGGATCTCTTGCGGCATAAAGGCGCAGTAAGTGCCAACCAGCGGAATGCCGCGCTCCTTg                                                                                                    <  2:325081/140‑1 (MQ=255)
                                   ggCAGATCTTTCTCCGCTTCTTCAATGGTTTCATCAGAGGTGGAACAGAGCGAAACCACAACCGCACCGGCTGCCATCGGGATCTCTTGCGGCATAAAGGCGCAGTAAGTGCCAACCAGCGGAATGCCGCGCTCCTTgag                                                                                                  <  2:109953/140‑1 (MQ=255)
                                   ggCAGATCTTTCTCCGCATCTTCAATGGTTTTATCCGAGATGGAACAGGGCGAAACCACAACCGCACCGGCTGCCATCGGGATCTCTTGCGGCATAAAGGCGCAGTAAGTGCCAACCAGCGGAATGCCGCGCTCCTTgag                                                                                                  <  2:247976/140‑1 (MQ=255)
                                    gCAGATCTTTCTCCGCTTCTTCAATGGTTTCATCAGAGGTGGAACAGAGCGAAACCACAACCGCACCGGCTGCCATCGGGATCTCTTGCGGCATAAAGGCGCAGTAAGTGCCAACCAGCGGAATGCCGCGCTCCTTgaga                                                                                                 >  1:139865/1‑140 (MQ=255)
                                               tcCGCTTCTTCAATGGTTTCATCAGAGGTGGAACAGAGCGAAACCACAACCGCACCGGCTGCCATCGGGATCTCTTGCGGCATAAAGGCGCAGTAAGTGCCAACCAGCGGAATGCCGCGCTCCTTGAGATCCAt                                                                                            <  2:120079/134‑1 (MQ=255)
                                               tcCGCTTCTTCAATGGTTTCATCAGAGGTGGAACAGAGCGAAACCACAACCGCACCGGCTGCCATCGGGATCTCTTGCGGCATAAAGGCGCAGTAAGTGCCAACCAGCGGAATGCCGCGCTCCTTGAGATCCAt                                                                                            >  1:120079/1‑134 (MQ=255)
                                                  gCTTCTTCAATGGTTTCATCAGAGGTGGAACAGAGCGAAACCACAACCGCACCGGCTGCCATCGGGATCTCTTGCGGCATAAAGGCGCAGTAAGTGCCAACCAGCGGAATGCCGCGCTCCTTGAGATCCATGACGGTgag                                                                                   >  1:142894/1‑140 (MQ=255)
                                                   cttcttCAATGGTTTCATCAGAGGTGGAACAGAGCGAAACCACAACCGCACCGGCTGCCATCGGGATCTCTTGCGGCATAAAGGCGCAGTAAGTGCCAACCAGCGGAATGCCGCGCTCCTTGAGATCCATGACGGTgaga                                                                                  >  2:392981/1‑140 (MQ=255)
                                                   cttcttCAATGGTTTCATCAGAGGTGGAACAGAGCGAAACCACAACCGCACCGGCTGCCATCGGGATCTCTTGCGGCATAAAGGCGCAGTAAGTGCCAACCAACGGAATGGCGCTCTTCTTGAGATCCATGACGGTgaga                                                                                  >  1:245476/1‑140 (MQ=255)
                                                      cttcAATGGTTTCATCAGAGGTGGAACAGAGCGAAACCACAACCGCACCGGCTGCCATCGGGATCTCTTGCGGCATAAAGGCGCAGTAAGTGCCAACCAGCGGAATGCCGCGCTCCTTGAGATCCATGACGGTGAGAAAg                                                                               >  1:120957/1‑140 (MQ=255)
                                                       ttcAATGGTTTCATCAGAGGTGGAACAGAGCGAAACCACAACCGCACCGGCTGCCATCGGGATCTCTTGCGGCATAAAGGCGCAGTAAGTGCCAACCAGCGGAATGCCGCGCTCCTTGAGATCCATGACGGTGAGAAAGc                                                                              >  2:565661/1‑140 (MQ=255)
                                                                ttCATCAGAGGTGGAACAGAGCGAAACCACAACCGCACCGGCTGCCATCGGGATCTCTTGCGGCATAAAGGCGCAGTAAGTGCCAACCAGCGGAATGCCGCGCTCCTTGAGATCCATGACGGTGAGAAAGCCTGTCTGGc                                                                     <  2:129592/140‑1 (MQ=255)
                                                                ttCATCAGAGGTGGAACAGAGCGAAACCACAACCGCACCGGCTGCCATCGGGATCTCTTGCGGCATAAAGGCGCAGTAAGTGCCAACCAGCGGAATGCCGCGCTCCTTGAGATCCATGACGGTGAGAAAGCCTGTCTGGc                                                                     >  1:150066/1‑140 (MQ=255)
                                                                 tcatcaGAGGTGGAACATAGCGCAACCACAATCGCGCCGACTGCCATCCGGATCTCTTACGGCATAAAGGCGCAGTAAGTGCCAACCAGCGGAATGCCGCGCTCCTTGAGATCCATGACGGTGAGAAAGCCTGTCTGGCg                                                                    <  2:225249/140‑1 (MQ=255)
                                                                   atcaGAGGTGGAACAGAGCGAAACCACAACCGCACCGGCTGCCATCGGGATCTCTTGCGGCATAAAGGCGCAGTAAGTGCCAACCAGCGGAATGCCGCGCTCCTTGAGATCCATGACGGTGAGAAAGCCTGTCTGGCGAg                                                                  <  2:150066/140‑1 (MQ=255)
                                                                   atcaGAGGTGGAACAGAGCGAAACCACAACCGCACCGGCTGCCATCGGGATCTCTTGCGGCATAAAGGCGCAGTAAGTGCCAACCAGCGGAATGCCGCGCTCCTTGAGATCCATGACGGTGAGAAAGCCTGTCTGGCGAg                                                                  <  2:407298/140‑1 (MQ=255)
                                                                   atcaGAGGTGGAACAGAGCGAAACCACAACCGCACCGGCTGCCATCGGGATCTCTTGCGGCATAAAGGCGCAGTAAGTGCCAACCAGCGGAATGCCGCGCTCCTTGAGATCCATGACGGTGAGAAAGCCTGTCTGGCGAg                                                                  <  2:139865/140‑1 (MQ=255)
                                                                     caGAGGTGGAACAGAGCGAAACCACAACCGCACCGGCTGCCATCGGGATCTCTTGCGGCATAAAGGCGCAGTAAGTGCCAACCAGCGGAATGCCGCGCTCCTTGAGATCCATGACGGTGAGAAAGCCTGTCTGGCGAGCt                                                                <  2:667586/140‑1 (MQ=255)
                                                                        agGTGGAACAGAGCGAAACCACAACCGCACCGGCTGCCATCGGGATCTCTTGCGGCATAAAGGCGCAGTAAGTGCCAACCAGCGGAATGCCGCGCTCCTTGAGATCCATGACGGTGAGAAAGCCTGTCTGGCGAGCTTCa                                                             >  2:543521/1‑140 (MQ=255)
                                                                          gTGGAACAGAGCGAAACCACAACCGCACCGGCTGCCATCGGGATCTCTTGCGGCATAAAGGCGCAGTAAGTGCCAACCAGCGGAATGCCGCGCTCCTTGAGATCCATGACGGTGAGAAAGCCTGTCTGGCGAGCTTCaga                                                           >  1:642818/1‑140 (MQ=255)
                                                                            ggAACAGAGCGAAACCACAACCGCACCGGCTGCCATCGGGATCTCTTGCGGCATAAAGGCGCAGTAAGTGCCAACCAGCGGAATGCCGCGCTCCTTGAGATCCATGACGGTGAGAAAGCCTGTCTGGCGAGCTTCAgaga                                                         <  1:515086/140‑1 (MQ=255)
                                                                             gAACAGAGCGAAACCACAACCGCACCGGCTGCCATCGGGATCTCTTGCGGCATAAAGGCGCAGTAAGTGCCAACCAGCGGAATGCCGCGCTCCTTGAGATCCATGACGGTGAGAAAGCCTGTCTGGCGAGCTTCAGAGaa                                                        >  2:138583/1‑140 (MQ=255)
                                                                             gAACAGAGCGAAACCACAACCGCACCGGCTGCCATCGGGATCTCTTGCGGCATAAAGGCGCAGTAAGTGCCAACCAGCGGAATGCCGCGCTCCTTGAGATCCATGACGGTGAGAAAGCCTGTCTGGCGAGCTTCAGAGaa                                                        >  2:10496/1‑140 (MQ=255)
                                                                                  gagCGAAACCACAACCGCACCGGCTGCCATCGGGATCTCTTGCGGCATAAAGGCGCAGTAAGTGCCAACCAGCGGAATGCCGCGCTCCTTGAGATCCATGACGGTGAGAAAGCCTGTCTGGCGAGCTTCAGAGAACTGAt                                                   >  1:478374/1‑140 (MQ=255)
                                                                                   agCGAAACCACAACCGCACCGGCTGCCATCGGGATCTCTTGCGGCATAAAGGCGCAGTAAGTGCCAACCAGCGGAATGCCGCGCTCCTTGAGATCCATGACGGTGAGAAAGCCTGTCTGGCGAGCTTCAGAGAACTGATc                                                  >  1:44441/1‑140 (MQ=255)
                                                                                   agCGAAACCACAACCGCACCGGCTGCCATCGGGATCTCTTGCGGCATAAAGGCGCAGTAAGTGCCAACCAGCGGAATGCCGCGCTCCTTGAGATCCATGACGGTGAGAAAGCCTGTCTGGCGAGCTTCAGAGAACTGATc                                                  >  1:143085/1‑140 (MQ=255)
                                                                                    gCGAAACCACAACCGCACCGGCTGCCATCGGGATCTCTTGCGGCATAAAGGCGCAGTAAGTGCCAACCAGCGGAATGCCGCGCTCCTTGAGGTCCATGACAGTGAGAATGCCTGTCTCGCGAACTTCAGAGAACTGATCg                                                 >  2:195069/1‑140 (MQ=255)
                                                                                    gCGAAACCACAACCGCACCGGCTGCCATCGGGATCTCTTGCGGCATAAAGGCGCAGTAAGTGCCAACCAGCGGAATGCCGCGCTCCTTGAGATCCATGACGGTGAGAAAGCCTGTCTGGCGAGCTTCAGAGAACTGATCg                                                 >  1:333358/1‑140 (MQ=255)
                                                                                    gCGAAACCACAACCGCACCGGCTGCCATCGGGATCTCTTGCGGCATAAAGGCGCAGTAAGTGCCAACCAGCGGAATGCCGCGCTCCTTGAGATCCATGACGGTGAGAAAGCCTGTCTGGCGAGCTTCAGAGAACTGATCg                                                 >  1:34774/1‑140 (MQ=255)
                                                                                    gCGAAACCACAACCGCACCGGCTGCCATCGGGATCTCTTGCGGCATAAAGGCGCAGTAAGTGCCAACCAGCGGAATGCCGCGCTCCTTGAGATCCATGACGGTGAGAAAGCCTGTCTGGCGAGCTTCAGAGAACTGATCg                                                 >  2:148941/1‑140 (MQ=255)
                                                                                            acaACCGCACCGGCTGCCATCGGGATCTCTTGCGGCATAAAGGCGCAGTAAGTGCCAACCAGCGGAATGCCGCGCTCCTTGAGATCCATGACGGTGCGAAAGCCTGTCTGGCGAGCTTCAGAGAACTGATCGAAAATGGc                                         <  2:359928/140‑1 (MQ=255)
                                                                                                ccGCACCGGCTGCCATCGGGATCTCTTGCGGCATAAAGGCGCAGTAAGTGCCAACCAGCGGAATGCCGCGCTCCTTGAGATCCATGACGGTGAGAAAGCCTGTCTGGCGAGCTTCAGAGAACTGATCGAAAATGGCGGGt                                     <  1:74092/140‑1 (MQ=255)
                                                                                                   cACCGGCTGCCATCGGGATCTCTTGCGGCATAAAGGCGCAGTAAGTGCCAACCAGCGGAATGCCGCGCTCCTTGAGATCCATGACGGTGAGAAAGCCTGTCTGGCGAGCTTCAGAGAACTGATCGAGAAGGGCGGGTTGa                                  >  2:615145/1‑140 (MQ=255)
                                                                                                   cACCGGCTGCCATCGGGATCTCTTGCGGCATAAAGGCGCAGTAAGTGCCAACCAGCGGAATGCCGCGCTCCTTGAGATCCATGACGGTGAGAAAGCCTGTCTGGCGAGCTTCAGAGAACTGATCGAAAATGGCGGGTAGa                                  >  1:518306/1‑140 (MQ=255)
                                                                                                   cACCGGCTGCCATCGGGATCTCTTGCGGCATAAAGGCGCAGTAAGTGCCAACCAGCGGAATGCCGCGCTCCTTGAGATCCATGACGGTGAGAAAGCCTGTCTGGCGAGCTTCAGAGAACTGATCGAAAATGGCGGGTAGa                                  <  2:333358/140‑1 (MQ=255)
                                                                                                   cACCGGCTGCCATCGGGATCTCTTGCGGCATAAAGGCGCAGTAAGTGCCAACCAGCGGAATGCCGCGCTCCTTGAGATCCATGACGGTGAGAAAGCCTGTCTGGCGAGCTTCAGAGAACTGATCGAAAATGGCGGGTAGa                                  >  2:359348/1‑140 (MQ=255)
                                                                                                        gCTGCCATCGGGATCTGTTGCGTCATAAAGGTGTATTAACTGCCAACCCGCGGAATGCCGCGCTCCTTGAGATCCATGACTGTGAGAAAGCCTGTCTGGCGAGCTTCAGAGAACTGATCGAAAATGGCGGGTAGATCGGt                             <  2:305696/140‑1 (MQ=255)
                                                                                                        gCTGCCATCGGGATCTCTTGCGGCATAAAGGCGCAGTAAGTGCCAACCAGCGGAATGCCGCGCTCCTTGAGATCCATGACGGTGAGAAAGCCTGTCTGGCGAGCTTCAGAGAACTGATCGAAAATGGCGGGTAGATCGGt                             <  1:615145/140‑1 (MQ=255)
                                                                                                        gCTGCCATCGGGATCTCTTGCGGCATAAAGGCGCAGTAAGTGCCAACCAGCGGAATGCCGCGCTCCTTGAGATCCATGACGGTGAGAAAGCCTGTCTGGCGAGCTTCAGAGAACTGAGCGAAAATGGCGGGTAGATCggg                             >  2:241617/1‑139 (MQ=255)
                                                                                                        gCTGCCATCGGGATCTCTTGCGGCATAAAGGCGCAGTAAGTGCCAACCAGCGGAATGCCGCGCTCCTTGAGATCCATGACGGTGAGAAAGCCTGTCTGGCCAGCATCATAGAACTGAGCGAATATGGCGGGTAgtgcggg                             >  2:303935/1‑134 (MQ=255)
                                                                                                               tCGGGATCTCTTGCGGCATAAAAGCGCAGTAAGTGCCCACCAGCGGAATGCCGCGCTCCTTGAGATCCATGACGGTGAGAAAGCCTGTCTGGCGAGCTTCAGAGAACTGATCGAAAATGGCGGGTAGATCGGTGACAAGt                      <  1:303935/140‑1 (MQ=255)
                                                                                                                 gggATCTCTTGCGGCATAAAGGCGCAGTAAGTGCCAACCAGCGGAATGCCGCGCTCCTTGAGATCCATGACGGTGAGAAAGCCTGTCTGGCGAGCTTCAGAGAACTGATCGAAAATGGCGGGTAGATCGGTGACAAGTGa                    <  2:555475/140‑1 (MQ=255)
                                                                                                                  ggATCTCTTGCGGCATAAAGGCGCAGTAAGTGCCAACCAGCGGAATGCCGCGCTCCTTGAGATCCATGACGGTGAGAAAGCCTGTCTGGCGAGCTTCAGAGAACTGATCGAAAATGGCGGGTAGATCGGTGACAAGTGAc                   <  2:518306/140‑1 (MQ=255)
                                                                                                                      ctctTGCGGCATAAAGGCGCAGTAAGTGCCAACCAGCGGAATGCCGCGCTCCTTGAGATCCATGACGGTGAGAAAGCCTGTCTGGCGAGCTTCAGAGAACTGATCGAAAATGGCGGGTAGATCGGTGACAAGTGACATGa               <  1:392981/140‑1 (MQ=255)
                                                                                                                        ctTGCGGCATAAAGGCGCAGTACGTCCCATCCAGCGGTATGCCGCGCTCCTTGAGATCCATGACGGTGAGAAAGCCTGTCTGGCGAGCTTCAGAGAACTGATCGAAAATGGCGGGTAGATCGGTGACAAGTGACATGAtt             <  1:195069/140‑1 (MQ=255)
                                                                                                                         ttGCGGCATAAAGGCGCAGTAAGTGCCAACCAGCGGAATGCCGCGCTCCTTGAGATCCATGACGGTGAGAAAGCCTGTCTGGCGAGCTTCAGAGAACTGATCTAAAATGGGGGGGAGGTCGGTGAGAAGGGACATGAttt            >  1:238795/1‑140 (MQ=255)
                                                                                                                                  aaaGGCGCAGTAAGTGCCAACCTGCGGAATGCCGCTCTCTTTGAGATCGTGACGGTGAGAAAGCCTGTCTGGCGAGCTTCAGAGAACTGATCGAAAATGGCGGGTAGATCGGTGACAAGTGACATGATTTTCCTTCCCCg  <  2:227856/140‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                       
CGTAGCTGCTTTTAATCAGCGGGCAGAGGTTGCGCGGCAGATCTTTCTCCGCTTCTTCAATGGTTTCATCAGAGGTGGAACAGAGCGAAACCACAACCGCACCGGCTGCCATCGGGATCTCTTGCGGCATAAAGGTGCAGTAAGTGCCAACCAGCGGAATGCCGCGCTCCTTGAGATCCATGACGGTGAGAAAGCCTGTCTGGCGAGCTTCAGAGAACTGATCGAAAATGGCGGGTAGATCGGTGACAAGTGACATGATTTTCCTTCCCCG  >  NC_000913/4565595‑4565865

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