Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 155,935 T→C K89K (AAA→AAG)  yadV ← putative periplasmic pilin chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913155,9350TC100.0% 295.1 / NA 80K89K (AAA→AAG) yadVputative periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (42/38);  total (42/38)

GCAGGCTTTGATTCGCGACCTTTTCTGCATCTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCA  >  NC_000913/155802‑156062
                                                                                                                                     |                                                                                                                               
gcagGCTTTGATTCGCGACCTTTTCTGCATCTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTACTTTGATTGTTTGCCCACGCTTGGc                                                                                                                            <  1:1217783/139‑1 (MQ=255)
   ggCTTTGATTCGCGACCTTTTCTGCATCTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCAt                                                                                                                          >  2:1508191/1‑138 (MQ=255)
   ggCTTTGATTCGCGACCTTTTCTGCATCTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCAt                                                                                                                          <  1:1508191/138‑1 (MQ=255)
              gcgACCTTTTCTGCATCTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATa                                                                                                              >  2:731679/1‑139 (MQ=255)
                          gCATCTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAg                                                                                                  <  2:228865/139‑1 (MQ=255)
                          gCATCTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATGAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGCAg                                                                                                  <  1:1040796/139‑1 (MQ=255)
                                    ttGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAgg                                                                                        <  1:867606/139‑1 (MQ=255)
                                    ttGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAgg                                                                                        <  1:804922/139‑1 (MQ=255)
                                      ggtggAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGa                                                                                      >  2:712738/1‑139 (MQ=255)
                                           aaCTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATATGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTa                                                                                 >  1:506246/1‑139 (MQ=255)
                                                ccAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACt                                                                            >  2:188456/1‑139 (MQ=255)
                                                      cGTTAAGCCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACGGGCGGCGTAGCAGTAAAAGGGACTTTAATACTGCCAgg                                                                      <  1:192594/139‑1 (MQ=255)
                                                       gTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATATGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGc                                                                     <  2:506246/139‑1 (MQ=255)
                                                       gTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTggcgg                                                                                                       >  1:1702207/1‑105 (MQ=255)
                                                       gTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTggcgg                                                                                                       <  2:1702207/105‑1 (MQ=255)
                                                       gTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAgg                                                                                        <  2:1037367/120‑1 (MQ=255)
                                                       gTTAAACCAGAACACGCTCTCTCTGTCTTAAGGCAGTGAGGTGCTTGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAgg                                                                                        >  1:1037367/1‑120 (MQ=255)
                                                         tAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAg                                                                                                  <  2:1492683/108‑1 (MQ=255)
                                                         tAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAg                                                                                                  >  1:1492683/1‑108 (MQ=255)
                                                         tAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTc                                                                   >  2:191151/1‑139 (MQ=255)
                                                         tAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTc                                                                   >  1:584379/1‑139 (MQ=255)
                                                           aaCCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAg                                                                 >  1:1278688/1‑139 (MQ=255)
                                                              cAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAg                                                                                                  <  2:62206/103‑1 (MQ=255)
                                                              cAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAg                                                                                                  >  1:62206/1‑103 (MQ=255)
                                                                gAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTg                                                            >  1:837564/1‑139 (MQ=255)
                                                                 aaCACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGt                                                           >  1:323798/1‑139 (MQ=255)
                                                                 aaCACGCGCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGt                                                           >  1:1174903/1‑139 (MQ=255)
                                                                  acacGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTc                                                          <  1:1431327/139‑1 (MQ=255)
                                                                      gCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGt                                                                                                    <  1:1787761/93‑1 (MQ=255)
                                                                      gCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGt                                                                                                    >  2:1787761/1‑93 (MQ=255)
                                                                        tctctctGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGcc                                                    >  2:1577057/1‑139 (MQ=255)
                                                                        tctctctGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGcc                                                    <  2:323798/139‑1 (MQ=255)
                                                                           ctctGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGctt                                                                                                                               >  1:422416/1‑61 (MQ=255)
                                                                             ctGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTAt                                               <  1:731679/139‑1 (MQ=255)
                                                                                tCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTa                                            >  1:294429/1‑139 (MQ=255)
                                                                                  tttAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAAc                                          <  1:712738/139‑1 (MQ=255)
                                                                                    tAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGt                                                 >  1:1733087/1‑130 (MQ=255)
                                                                                    tAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGt                                                 <  2:1733087/130‑1 (MQ=255)
                                                                                          gTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTaaaa                                                                                          >  2:428165/1‑83 (MQ=255)
                                                                                          gTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTaaaa                                                                                          <  1:428165/83‑1 (MQ=255)
                                                                                          gTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAgg                                                                                        <  1:32255/85‑1 (MQ=255)
                                                                                          gTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAgg                                                                                        >  2:32255/1‑85 (MQ=255)
                                                                                          gTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCAt                                                        >  1:759761/1‑117 (MQ=255)
                                                                                          gTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCAt                                                        <  2:759761/117‑1 (MQ=255)
                                                                                           tGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGc                                                                     <  2:1011765/103‑1 (MQ=255)
                                                                                           tGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGc                                                                     >  1:1011765/1‑103 (MQ=255)
                                                                                           tGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTgg                                 >  1:1735957/1‑139 (MQ=255)
                                                                                            gAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGa                                <  2:1278688/139‑1 (MQ=255)
                                                                                            gAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGa                                >  2:180026/1‑139 (MQ=255)
                                                                                            gAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACACGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCg                                                               <  2:1547376/108‑1 (MQ=255)
                                                                                            gAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACACGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCg                                                               >  1:1547376/1‑108 (MQ=255)
                                                                                             aGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGAc                               <  1:191151/139‑1 (MQ=255)
                                                                                                 gctggctgTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGc                           >  1:418642/1‑139 (MQ=255)
                                                                                                  ctggctgTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCa                                                         <  2:568306/108‑1 (MQ=255)
                                                                                                  ctggctgTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCa                                                         >  1:568306/1‑108 (MQ=255)
                                                                                                   tggctgTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCaa                         <  1:1577057/139‑1 (MQ=255)
                                                                                                     gctgTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACt                                                                                    >  2:800945/1‑78 (MQ=255)
                                                                                                     gctgTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACt                                                                                    <  1:800945/78‑1 (MQ=255)
                                                                                                         tgtACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCa                          >  1:1432182/1‑132 (MQ=255)
                                                                                                         tgtACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCa                          <  2:1432182/132‑1 (MQ=255)
                                                                                                         tgtACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGtt                   >  2:115222/1‑139 (MQ=255)
                                                                                                          gtACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTa                  >  2:1367120/1‑139 (MQ=255)
                                                                                                          gtACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTCTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTa                  <  2:1174903/139‑1 (MQ=255)
                                                                                                               ttAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTccc             >  2:1757700/1‑139 (MQ=255)
                                                                                                                tAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTaa                                                                                >  2:1606323/1‑71 (MQ=255)
                                                                                                                tAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTaa                                                                                <  1:1606323/71‑1 (MQ=255)
                                                                                                                     tGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAg                                                                                                  <  2:534815/48‑1 (MQ=255)
                                                                                                                     tGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAg                                                                                                  >  1:534815/1‑48 (MQ=255)
                                                                                                                       aTTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTg                                                                           <  2:1708762/69‑1 (MQ=255)
                                                                                                                       aTTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTg                                                                           >  1:1708762/1‑69 (MQ=255)
                                                                                                                       aTTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAAc                                          <  2:502985/102‑1 (MQ=255)
                                                                                                                       aTTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAAc                                          >  1:502985/1‑102 (MQ=255)
                                                                                                                       aTTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTAtttt     >  2:297548/1‑139 (MQ=255)
                                                                                                                        ttGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTc    >  2:1627106/1‑139 (MQ=255)
                                                                                                                          gTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCa  <  2:418642/139‑1 (MQ=255)
                                                                                                                            ttGCCCACGCTTGGCATCAATACGCGATACTGgcggcg                                                                                                     >  2:1516309/1‑38 (MQ=38)
                                                                                                                            ttGCCCACGCTTGGCATCAATACGCGATACTGgcggcg                                                                                                     <  1:1516309/38‑1 (MQ=38)
                                                                                                                             tGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGtt                   >  2:40416/1‑119 (MQ=255)
                                                                                                                             tGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGtt                   <  1:40416/119‑1 (MQ=255)
                                                                                                                             tGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTAt                 <  1:1709267/121‑1 (MQ=255)
                                                                                                                             tGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTAt                 >  2:1709267/1‑121 (MQ=255)
                                                                                                                                     |                                                                                                                               
GCAGGCTTTGATTCGCGACCTTTTCTGCATCTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCA  >  NC_000913/155802‑156062

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