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% 296.8 / NA 82K89K (AAA→AAGyadVputative periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (45/37);  total (45/37)

TCTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGAC  >  NC_000913/155831‑156065
                                                                                                        |                                                                                                                                  
tCTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAg                                                                                                     <  1:1342491/136‑1 (MQ=255)
tCTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAg                                                                                                     >  2:1342491/1‑136 (MQ=255)
 cTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGt                                                                                                 >  1:246097/1‑139 (MQ=255)
   ggTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACt                                                                                                               <  2:1103449/123‑1 (MQ=255)
   ggTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACt                                                                                                               >  1:1103449/1‑123 (MQ=255)
   ggTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTaa                                                                                               >  2:195116/1‑139 (MQ=255)
       ttGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTTAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAgg                                                                                           <  1:195116/139‑1 (MQ=255)
       ttGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAgg                                                                                           <  2:961062/139‑1 (MQ=255)
          gtggAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGAc                                                                                        >  1:1365122/1‑139 (MQ=255)
                     aGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGTGGCGTAg                                                                                                     >  1:1198680/1‑115 (MQ=255)
                     aGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTATTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAg                                                                                                     <  2:1198680/115‑1 (MQ=255)
                           ttAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAAt                                                                                                                         >  2:215175/1‑89 (MQ=255)
                           ttAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAAt                                                                                                                         <  1:215175/89‑1 (MQ=255)
                            tAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAgg                                                                                           <  2:522230/118‑1 (MQ=255)
                            tAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAgg                                                                                           >  1:522230/1‑118 (MQ=255)
                                 cAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGt                                                                 >  2:1618624/1‑139 (MQ=255)
                                   gAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCAt                                                                                                                             <  1:1647023/77‑1 (MQ=255)
                                   gAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCAt                                                                                                                             >  2:1647023/1‑77 (MQ=255)
                                    aaCACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGgcggcg                                                                                                        <  1:1794576/97‑1 (MQ=255)
                                    aaCACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGgcggcg                                                                                                        >  2:1794576/1‑97 (MQ=255)
                                     acacGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCAGGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTc                                                             <  2:1124901/139‑1 (MQ=255)
                                           tctctctGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTgg                                                                                                             <  1:985601/85‑1 (MQ=255)
                                           tctctctGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTgg                                                                                                             >  2:985601/1‑85 (MQ=255)
                                              ctctGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTc                                                             >  1:1667268/1‑130 (MQ=255)
                                              ctctGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTc                                                             <  2:1667268/130‑1 (MQ=255)
                                                     tttAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCg                                                                  >  2:630379/1‑118 (MQ=255)
                                                     tttAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCg                                                                  <  1:630379/118‑1 (MQ=255)
                                                     tttAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAAc                                             <  2:246097/139‑1 (MQ=255)
                                                       tAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCg                                                                  >  2:762157/1‑116 (MQ=255)
                                                       tAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCg                                                                  <  1:762157/116‑1 (MQ=255)
                                                             gTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGAc                                                                                        >  2:162312/1‑88 (MQ=255)
                                                             gTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGAc                                                                                        <  1:162312/88‑1 (MQ=255)
                                                             gTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCa                                                            >  2:876864/1‑116 (MQ=255)
                                                             gTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCa                                                            <  1:876864/116‑1 (MQ=255)
                                                             gTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCa                                                      <  2:468184/122‑1 (MQ=255)
                                                             gTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCa                                                      >  1:468184/1‑122 (MQ=255)
                                                              tGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTgg                                    >  2:1724912/1‑139 (MQ=255)
                                                              tGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTgg                                    >  2:314279/1‑139 (MQ=255)
                                                                     ctggctgTGTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAAGACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCa                             <  1:314279/139‑1 (MQ=255)
                                                                          tgtgTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAAc                                         <  1:1471337/122‑1 (MQ=255)
                                                                          tgtgTACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAAc                                         >  2:1471337/1‑122 (MQ=255)
                                                                            tgtACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAgg                                                                                           >  1:498342/1‑70 (MQ=255)
                                                                            tgtACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAgg                                                                                           <  2:498342/70‑1 (MQ=255)
                                                                            tgtACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGtt                      <  2:1365122/139‑1 (MQ=255)
                                                                             gtACATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTa                     >  2:1743977/1‑139 (MQ=255)
                                                                                cATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGAc                                  <  2:929554/123‑1 (MQ=255)
                                                                                cATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTc                  >  2:823191/1‑139 (MQ=255)
                                                                                cATTAATTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTCGTCATCGCCAGTATCTAACCAACTCTGGAc                                  >  1:929554/1‑123 (MQ=255)
                                                                                  ttAATTTGATTGTTCGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAAAACTGCCAGGCTCAGCGTTGTCATCGCAAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTccc                >  2:1114160/1‑139 (MQ=255)
                                                                                    aaTTTGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCtt              <  1:463317/139‑1 (MQ=255)
                                                                                      tttGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGtt                                                                >  2:1516633/1‑87 (MQ=255)
                                                                                      tttGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGtt                                                                <  1:1516633/87‑1 (MQ=255)
                                                                                      tttGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAAc                                             <  1:1026081/106‑1 (MQ=255)
                                                                                      tttGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAAc                                             >  2:1026081/1‑106 (MQ=255)
                                                                                      tttGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTa            >  1:472094/1‑139 (MQ=255)
                                                                                        tGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTaaa                                                                                              >  2:448089/1‑55 (MQ=255)
                                                                                        tGATTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTaaa                                                                                              <  1:448089/55‑1 (MQ=255)
                                                                                          aTTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAgg                                                                                           <  2:1598163/56‑1 (MQ=255)
                                                                                          aTTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAgg                                                                                           >  1:1598163/1‑56 (MQ=255)
                                                                                          aTTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGc                                                        <  2:459944/91‑1 (MQ=255)
                                                                                          aTTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGc                                                        >  1:459944/1‑91 (MQ=255)
                                                                                          aTTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTAtttt        >  2:1567693/1‑139 (MQ=255)
                                                                                          aTTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTAtttt        >  2:913521/1‑139 (MQ=255)
                                                                                          aTTGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTAtttt        >  2:614792/1‑139 (MQ=255)
                                                                                           ttGTTTGCCCACGCTTGTCATCAATACGCGATACt                                                                                                               >  1:1050235/1‑35 (MQ=25)
                                                                                           ttGTTTGCCCACGCTTGGCATCAATACGCGATACt                                                                                                               <  2:1050235/35‑1 (MQ=255)
                                                                                           ttGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAgg                                                                                           <  2:129666/55‑1 (MQ=255)
                                                                                           ttGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAgg                                                                                           >  1:129666/1‑55 (MQ=255)
                                                                                           ttGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTaa                                                                                   >  1:1188452/1‑63 (MQ=255)
                                                                                           ttGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTaa                                                                                   <  2:1188452/63‑1 (MQ=255)
                                                                                           ttGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGc                                                        >  1:426604/1‑90 (MQ=255)
                                                                                           ttGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGc                                                        <  2:426604/90‑1 (MQ=255)
                                                                                           ttGTTTGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTc       >  1:1163556/1‑139 (MQ=255)
                                                                                              tttGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTaaaa                                                                                             >  2:1639175/1‑50 (MQ=255)
                                                                                              tttGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTaaaa                                                                                             <  1:1639175/50‑1 (MQ=255)
                                                                                                tGCCCACGCTTGGCATCAATACGTGATACTGGCGGCGTAGCAGTAAAAGGGa                                                                                         <  1:841070/52‑1 (MQ=255)
                                                                                                tGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTATCAGTAAAAGGGa                                                                                         >  2:841070/1‑52 (MQ=255)
                                                                                                tGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGAc                                  >  2:1671817/1‑107 (MQ=255)
                                                                                                tGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGAc                                  <  1:1671817/107‑1 (MQ=255)
                                                                                                tGCCCACGCTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGac  >  1:1455093/1‑139 (MQ=255)
                                                                                                tGCCCACGCTTGGCATCAATACGCGATACCGGCGGCGTAGCAGTaaa                                                                                              >  2:808814/1‑47 (MQ=255)
                                                                                                tGCCCACGCTTGGCATCAATACGCGATACCGGCGGCGTAGCAGTaaa                                                                                              <  1:808814/47‑1 (MQ=255)
                                                                                                        |                                                                                                                                  
TCTGGTTTTGGTGGAACTTCCAGTACGTTAAACCAGAACACGCTCTCTCTGTCTTTAGGCAGTGAGGTGCTGGCTGTGTACATTAATTTGATTGTTTGCCCACGTTTGGCATCAATACGCGATACTGGCGGCGTAGCAGTAAAAGGGACTGTAATACTGCCAGGCTCAGCGTTGTCATCGCCAGTATCTAACCAACTCTGGACAAGCAACGGGTTATTCCCTTTATTTTCCAGAC  >  NC_000913/155831‑156065

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