Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,379,266 (C)7→9 coding (732/1128 nt) zapE ← cell division factor ZapE with ATPase activity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,379,2591.C98.0% 180.5 / ‑3.2 51coding (739/1128 nt)zapEcell division factor ZapE with ATPase activity
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (23/27);  total (23/28)
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.38e-01
*NC_0009133,379,2592.C98.0% 181.3 / ‑3.2 51coding (739/1128 nt)zapEcell division factor ZapE with ATPase activity
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (23/27);  total (23/28)
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.88e-01

GTGAGAGCGCAATATAGTCATGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTT  >  NC_000913/3379121‑3379382
                                                                                                                                             ||                                                                                                                            
gTGAGAGCGCAATATAGTCATGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑c                                                                                                                            <  1:877826/140‑1 (MQ=255)
  gagaGCGCAATATAGTCATGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑ccc                                                                                                                          >  2:693758/1‑140 (MQ=255)
       cgcAATATAGTCATGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑cccccccc                                                                                                                     >  2:68012/1‑139 (MQ=255)
         cAATATAGTCATGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑cccccccccg                                                                                                                   >  1:182580/1‑137 (MQ=255)
          aaTATAGTCATGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑cccccccccgc                                                                                                                  >  2:464468/1‑136 (MQ=255)
           atatAGTCATGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑ccccccccccgc                                                                                                                 <  1:1077713/140‑3 (MQ=255)
            tataGTCATGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑cccccccccgcca                                                                                                                <  2:182580/140‑7 (MQ=255)
                gTCATGCTGACTGCGGGCGTCGACGCACAGCGTAGTAGAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCGCCAgcgc                                                                                                              <  1:1051637/140‑1 (MQ=255)
                gTCATGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCGCCAgcgc                                                                                                              <  1:1064198/140‑1 (MQ=255)
                 tcttGCTGACGGCGGGCGTCGACGCAAGGCGTAGAAAAAGAGCCCGCCAGCGTCTGGTTCTCGACGCCCTTTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCGCCAGCGcc                                                                                                             <  2:292924/137‑1 (MQ=255)
                 tCATGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCGCCAGCGcc                                                                                                             <  1:824719/140‑1 (MQ=255)
                    tGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCGCCAGCGCCAAc                                                                                                          >  1:12408/1‑140 (MQ=255)
                    tGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCGCCAGCGCCAAc                                                                                                          >  2:1114219/1‑140 (MQ=255)
                    tGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCCACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCGCCAGCGCCAAc                                                                                                          >  2:1159817/1‑140 (MQ=255)
                      cTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTTGTT‑‑TCCCCCCCCCGCCAGCGCCAACCa                                                                                                        >  2:79716/1‑140 (MQ=255)
                      cTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCCGCCAGCGCCAAcc                                                                                                         >  2:1074137/1‑140 (MQ=255)
                          cTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCGCCAGCGCCAACCATAGt                                                                                                    <  1:901474/140‑1 (MQ=255)
                           tGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCGCCAGCGCCAACCATAGtt                                                                                                   >  1:268677/1‑140 (MQ=255)
                                         cacgGCGTAGAAAAAAGGCCGGCCAGCTTCTGGTTCTCGACGCCCTTTTTTGCAAATGGCCGATGGTTATTTTCTAACGTCGGTGATTTTTCTCTTTCCC‑‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGcc                                                                                     <  2:581823/140‑1 (MQ=255)
                                               gtagtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGtt                                                                               <  2:268677/140‑1 (MQ=255)
                                               gtagtaAAAGAGACCGCCAGCGTCTGGTTCTCCACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGtt                                                                               <  1:1159817/140‑1 (MQ=255)
                                                tagtaAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGtt                                                                               <  1:1074137/140‑1 (MQ=255)
                                                    aaaaGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATc                                                                          <  1:68012/140‑1 (MQ=255)
                                                        gagaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGa                                                                      >  1:426631/1‑140 (MQ=255)
                                                        gagaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGa                                                                      >  2:22913/1‑140 (MQ=255)
                                                          gaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTgcgc                                                                                      <  2:206389/122‑1 (MQ=255)
                                                          gaCCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTgcgc                                                                                      >  1:206389/1‑122 (MQ=255)
                                                             cGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTgg                                                                 >  1:151557/1‑140 (MQ=255)
                                                             cGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTgg                                                                 >  2:606698/1‑140 (MQ=255)
                                                             cGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTgg                                                                 >  2:1004586/1‑140 (MQ=255)
                                                                   cGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAg                                                           >  1:21658/1‑140 (MQ=255)
                                                                   cGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAg                                                           >  1:792284/1‑140 (MQ=255)
                                                                         gTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAg                                                     <  2:418854/140‑1 (MQ=255)
                                                                           tctcGACGCCCTTTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGAt                                                   <  2:557625/140‑1 (MQ=255)
                                                                               gacccgcATTGTTGCAAGTGGCCGGTGGTTGATTTGTAACGGCGGTGAATTTTGTGGTTGCC‑‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCAGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATgcgc                                               <  2:1072236/134‑1 (MQ=255)
                                                                                 cGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCt                                             <  2:385457/140‑1 (MQ=255)
                                                                                   cccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtga                                           <  1:22913/140‑1 (MQ=255)
                                                                                    ccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTTGTT‑‑TCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgag                                          <  1:79716/140‑1 (MQ=255)
                                                                                    ccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgag                                          <  2:1063865/140‑1 (MQ=255)
                                                                                    ccATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgag                                          <  1:693758/140‑1 (MQ=255)
                                                                                              tAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACg                                <  1:1004586/140‑1 (MQ=255)
                                                                                              tAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACg                                >  2:792527/1‑140 (MQ=255)
                                                                                               aaTGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGc                               <  2:426631/140‑1 (MQ=255)
                                                                                                aTGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCa                              <  1:792527/140‑1 (MQ=255)
                                                                                                aTGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCa                              <  1:606698/140‑1 (MQ=255)
                                                                                                aTGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCa                              >  1:108984/1‑140 (MQ=255)
                                                                                                aTGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCa                              <  1:755077/140‑1 (MQ=255)
                                                                                                  ggCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtga                                       <  1:639663/129‑1 (MQ=255)
                                                                                                  ggCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCtgagtga                                       >  2:639663/1‑129 (MQ=255)
                                                                                                  ggCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGa                            <  2:12408/140‑1 (MQ=255)
                                                                                                  ggCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGa                            >  2:608962/1‑140 (MQ=255)
                                                                                                  ggCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGa                            <  2:21658/140‑1 (MQ=255)
                                                                                                    ccGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACg                          >  2:420976/1‑140 (MQ=255)
                                                                                                    ccGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACg                          >  1:809116/1‑140 (MQ=255)
                                                                                                     cGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGa                         <  2:108984/140‑1 (MQ=255)
                                                                                                     cGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGa                         <  2:792284/140‑1 (MQ=255)
                                                                                                             gATTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAac                 >  1:1079110/1‑140 (MQ=255)
                                                                                                              aTTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAaca                >  2:551863/1‑140 (MQ=255)
                                                                                                              aTTTCTAACGTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAaca                >  1:614042/1‑140 (MQ=255)
                                                                                                                       gTCGGTGAATTTTCTCGTT‑‑TCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTcc       <  2:614042/140‑1 (MQ=255)
                                                                                                                            tGAATTTTCTCGTT‑‑TCCCCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGtt  <  1:420976/140‑1 (MQ=255)
                                                                                                                                             ||                                                                                                                            
GTGAGAGCGCAATATAGTCATGCTGACTGCGGGCGTCGACGCACAGCGTAGTAAAAGAGACCGCCAGCGTCTGGTTCTCGACGCCCATTGTTGCTAATGGCCGATGGTTGATTTCTAACGTCGGTGAATTTTCTCGTT‑‑T‑‑‑CCCCCCCGCCAGCGCCAACCATAGTTTATCCATTTGCGCCCGGGTTTCATCGTGAAGTGGCGAAAGCCACAGATGCGCCTGAGTGAGTGTACGCAGACGATAATCAACACCAGCGTCCACGTT  >  NC_000913/3379121‑3379382

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