Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,549,054 T→C S83P (TCG→CCG)  fimH → minor component of type 1 fimbriae

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,549,0540TC100.0% 335.8 / NA 95S83P (TCG→CCG) fimHminor component of type 1 fimbriae
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (44/51);  total (44/51)

CGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGG  >  NC_000913/4548921‑4549186
                                                                                                                                     |                                                                                                                                    
cgcCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGATTATGTCACACTGCAACGAGGCCCGGCt                                                                                                                                 >  1:71180/1‑139 (MQ=255)
cgcCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCt                                                                                                                                 >  1:154797/1‑139 (MQ=255)
       gTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATggcg                                                                                                                          <  2:490253/139‑1 (MQ=255)
       gTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATggcg                                                                                                                          <  2:154797/139‑1 (MQ=255)
          tATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGgcggcg                                                                                                                       >  1:536294/1‑139 (MQ=255)
                    ttGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTaa                                                                                                             >  2:1486624/1‑139 (MQ=255)
                     tGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAAt                                                                                                            >  1:1526161/1‑139 (MQ=255)
                     tGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAAt                                                                                                            >  1:923585/1‑139 (MQ=255)
                       cgcCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAAttt                                                                                                          >  1:349372/1‑139 (MQ=255)
                       cgcCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAAttt                                                                                                          >  1:1500749/1‑139 (MQ=255)
                       cgcCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAAttt                                                                                                          >  1:921530/1‑139 (MQ=255)
                       cgcCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAAt                                                                                                            >  1:1184336/1‑137 (MQ=255)
                       cgcCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAAt                                                                                                            <  2:1184336/137‑1 (MQ=255)
                          ccGTCGTGAATCTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTc                                                                                                       <  1:1082145/139‑1 (MQ=255)
                               gtGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                                  <  2:1374433/139‑1 (MQ=255)
                               gtGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                                  <  1:668692/139‑1 (MQ=255)
                               gtGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                                  >  1:1374433/1‑139 (MQ=255)
                               gtGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                                  <  2:1427445/139‑1 (MQ=255)
                               gtGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                                  <  2:1500749/139‑1 (MQ=255)
                               gtGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                                  <  2:921530/139‑1 (MQ=255)
                               gtGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                                  <  1:1112980/139‑1 (MQ=255)
                                   aTGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCGACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGt                                                                                              <  2:536768/139‑1 (MQ=255)
                                     gtgGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTaa                                                                                            <  1:707914/139‑1 (MQ=255)
                                       ggggCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCCTTACAGACTATGTAACACTGCAACGAGGCCCGGCCTATGGCGGCGTGTGATCTAATTTTTCCGGAACCGTaaaa                                                                                          >  2:995947/1‑139 (MQ=255)
                                         ggCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTa                                                                                             <  2:1206865/134‑1 (MQ=255)
                                         ggCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTa                                                                                             >  1:1206865/1‑134 (MQ=255)
                                          gCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAatat                                                                                       >  1:954887/1‑139 (MQ=255)
                                           cAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAAtata                                                                                      <  1:106496/139‑1 (MQ=255)
                                           cAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAAtata                                                                                      <  2:536294/139‑1 (MQ=255)
                                              aaCCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                                  >  1:1130003/1‑124 (MQ=255)
                                              aaCCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCAGTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                                  <  2:1130003/124‑1 (MQ=255)
                                                ccTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGc                                                                                 <  2:263379/139‑1 (MQ=255)
                                                 cTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCa                                                                                <  2:1161890/139‑1 (MQ=255)
                                                     tCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAg                                                                            <  2:71180/139‑1 (MQ=255)
                                                      cGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGc                                                                           <  2:1526161/139‑1 (MQ=255)
                                                        tGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTaa                                                                                            <  1:787210/120‑1 (MQ=255)
                                                        tGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTaa                                                                                            >  2:787210/1‑120 (MQ=255)
                                                          gATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATc                                                                       >  2:956952/1‑139 (MQ=255)
                                                              tttCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATggcg                                                                                                                          <  1:286179/84‑1 (MQ=255)
                                                              tttCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATggcg                                                                                                                          >  2:286179/1‑84 (MQ=255)
                                                              tttCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                                  <  1:557045/108‑1 (MQ=255)
                                                              tttCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                                  >  2:557045/1‑108 (MQ=255)
                                                               ttCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCAttt                                                                  >  1:1250329/1‑139 (MQ=255)
                                                                       aaaTCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                                  <  2:1349846/99‑1 (MQ=255)
                                                                       aaaTCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                                  >  1:1349846/1‑99 (MQ=255)
                                                                       aaaTCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAAtata                                                                                      >  1:947311/1‑111 (MQ=255)
                                                                       aaaTCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAAtata                                                                                      <  2:947311/111‑1 (MQ=255)
                                                                              ttGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTAACCATTTCCTACCACCAGCGaa                                                   >  2:414329/1‑139 (MQ=255)
                                                                               tGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAAt                                                                                                            <  2:1446781/81‑1 (MQ=255)
                                                                               tGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAAt                                                                                                            >  1:1446781/1‑81 (MQ=255)
                                                                                     aaCGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATggcg                                                                                                                          <  2:1602405/61‑1 (MQ=255)
                                                                                     aaCGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATggcg                                                                                                                          >  1:1602405/1‑61 (MQ=255)
                                                                                     aaCGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAAcgccgc                                            >  2:73735/1‑139 (MQ=255)
                                                                                     aaCGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAAcg                                                <  2:733625/135‑1 (MQ=255)
                                                                                     aaCGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAAcg                                                >  1:733625/1‑135 (MQ=255)
                                                                                     aaCGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAAcg                                                <  2:497520/135‑1 (MQ=255)
                                                                                     aaCGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAAcg                                                >  1:497520/1‑135 (MQ=255)
                                                                                      aCGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTAACCATTTCCTACCACCAGCGAAACGCcgcg                                           <  1:414329/139‑1 (MQ=255)
                                                                                        gATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCcgcgcg                                         <  1:956952/139‑1 (MQ=255)
                                                                                         aTTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCgt                                        <  1:536750/139‑1 (MQ=255)
                                                                                             tCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCgttgtt                                    <  2:954887/139‑1 (MQ=255)
                                                                                               cGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTta                                  <  1:1486624/139‑1 (MQ=255)
                                                                                                    aCCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGAcc                                                                                                >  2:918961/1‑72 (MQ=255)
                                                                                                    aCCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGAcc                                                                                                <  1:918961/72‑1 (MQ=255)
                                                                                                          aCAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                                  <  2:1050202/64‑1 (MQ=255)
                                                                                                          aCAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                                  >  1:1050202/1‑64 (MQ=255)
                                                                                                              aCTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCcgcgcg                                         <  1:1333081/117‑1 (MQ=255)
                                                                                                              aCTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCcgcgcg                                         >  2:1333081/1‑117 (MQ=255)
                                                                                                              aCTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGa                   <  2:1250329/139‑1 (MQ=255)
                                                                                                               cTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGAt                  <  1:73735/139‑1 (MQ=255)
                                                                                                                tATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGc                                                                                 <  1:1234801/75‑1 (MQ=255)
                                                                                                                tATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGc                                                                                 >  2:1234801/1‑75 (MQ=255)
                                                                                                                  tGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAAtata                                                                                      >  1:1232450/1‑68 (MQ=255)
                                                                                                                  tGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTAt                                                                        >  1:1401528/1‑82 (MQ=255)
                                                                                                                  tGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTAt                                                                        <  2:1401528/82‑1 (MQ=255)
                                                                                                                  tGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCa                                                                     <  2:420131/85‑1 (MQ=255)
                                                                                                                  tGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCa                                                                     >  1:420131/1‑85 (MQ=255)
                                                                                                                  tGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAg               >  2:425044/1‑139 (MQ=255)
                                                                                                                  tGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAg               >  2:1291966/1‑139 (MQ=255)
                                                                                                                  tGTCACACGGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAAtata                                                                                      <  2:1232450/68‑1 (MQ=255)
                                                                                                                   gTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCgt                                        <  1:1476423/113‑1 (MQ=255)
                                                                                                                   gTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCgt                                        >  2:1476423/1‑113 (MQ=255)
                                                                                                                        acTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTgg         >  2:691791/1‑139 (MQ=255)
                                                                                                                          tGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCcgcgcg                                         >  2:1263953/1‑105 (MQ=255)
                                                                                                                          tGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCcgcgcg                                         <  1:1263953/105‑1 (MQ=255)
                                                                                                                          tGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGcc       >  1:1169700/1‑139 (MQ=255)
                                                                                                                            cAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTc                                                                                                       <  1:80213/41‑1 (MQ=255)
                                                                                                                            cAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTc                                                                                                       >  2:80213/1‑41 (MQ=255)
                                                                                                                            cAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTtat                                 <  2:1615077/111‑1 (MQ=255)
                                                                                                                            cAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTtat                                 >  1:1615077/1‑111 (MQ=255)
                                                                                                                             aaCGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGt    <  1:1291966/139‑1 (MQ=255)
                                                                                                                               cGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCcgcgcg                                         <  1:355893/100‑1 (MQ=255)
                                                                                                                               cGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTgg  <  1:425044/139‑1 (MQ=255)
                                                                                                                               cGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGATATCCATTTCCTACCACCAGCGAAACGCcgcgcg                                         >  2:355893/1‑100 (MQ=255)
                                                                                                                                  ggCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCcgcgcg                                         <  1:1646564/97‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                    
CGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGG  >  NC_000913/4548921‑4549186

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