Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,549,001 G→A C65Y (TGC→TAC)  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,0010GA100.0% 249.0 / NA 76C65Y (TGC→TAC) fimHminor component of type 1 fimbriae
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (37/39);  total (37/39)

GTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCA  >  NC_000913/4548867‑4549132
                                                                                                                                      |                                                                                                                                   
gTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCAt                                                                                                                                 >  1:294309/1‑139 (MQ=255)
gTCATTCGCCTGTAAAACCGCCAACGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCAt                                                                                                                                 >  2:304316/1‑139 (MQ=255)
     tCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGa                                                                                                                            >  2:139459/1‑139 (MQ=255)
     tCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGa                                                                                                                            >  2:337703/1‑139 (MQ=255)
      cGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGAt                                                                                                                           >  2:24560/1‑139 (MQ=255)
           gTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATcc                                                                                                                      >  1:464648/1‑139 (MQ=255)
           gTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTATATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCAACGCAAATCTTTTACCATAACGATTATcc                                                                                                                      >  2:926222/1‑139 (MQ=255)
            tAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCACATCTTTTACCATAACGATTATCCg                                                                                                                     <  1:926222/139‑1 (MQ=255)
              aaaCCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGa                                                                                                                   <  2:597478/139‑1 (MQ=255)
                aCCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGaaa                                                                                                                 <  1:304316/139‑1 (MQ=255)
                 ccGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAAc                                                                                                                <  2:175590/139‑1 (MQ=255)
                        tGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACa                                                                                                         <  2:464648/139‑1 (MQ=255)
                        tGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACa                                                                                                         <  1:764680/139‑1 (MQ=255)
                        tGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACa                                                                                                         <  2:294309/139‑1 (MQ=255)
                          gTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGa                                                                                                       <  2:791519/139‑1 (MQ=255)
                            aCCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACt                                                                                                     <  1:337703/139‑1 (MQ=255)
                               gCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGaaa                                                                                                                 >  1:949005/1‑124 (MQ=255)
                               gCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGaaa                                                                                                                 <  2:949005/124‑1 (MQ=255)
                                    cccTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTcaca                                                                                             <  1:411193/139‑1 (MQ=255)
                                     ccTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCacac                                                                                            >  2:571872/1‑139 (MQ=255)
                                        aTTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGc                                                                                         >  1:941130/1‑139 (MQ=255)
                                             cGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGa                                                                                    >  2:510802/1‑139 (MQ=255)
                                                 ggCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCt                                                                                <  1:139459/139‑1 (MQ=255)
                                                            tGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCt                                                                                >  1:487935/1‑128 (MQ=255)
                                                            tGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCt                                                                                <  2:487935/128‑1 (MQ=255)
                                                            tGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATggc                                                                     <  2:941130/139‑1 (MQ=255)
                                                             gTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATggcg                                                                    <  2:1189968/139‑1 (MQ=255)
                                                                     aaaCCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTa                                                            >  2:189931/1‑139 (MQ=255)
                                                                      aaCCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTAt                                                           >  2:616663/1‑139 (MQ=255)
                                                                       aCCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATc                                                          >  2:1206867/1‑139 (MQ=255)
                                                                        ccTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCt                                                         >  2:235103/1‑139 (MQ=255)
                                                                          ttGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTa                                                                                                           >  1:860766/1‑87 (MQ=255)
                                                                          ttGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTa                                                                                                           <  2:860766/87‑1 (MQ=255)
                                                                          ttGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGa                                                                                                       <  1:220250/91‑1 (MQ=255)
                                                                          ttGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGa                                                                                                       >  2:220250/1‑91 (MQ=255)
                                                                           tGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCa                                                                                                                                  <  1:875024/63‑1 (MQ=255)
                                                                           tGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCa                                                                                                                                  >  2:875024/1‑63 (MQ=255)
                                                                           tGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGaaa                                                                                                                 >  1:1206765/1‑80 (MQ=255)
                                                                           tGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGaaa                                                                                                                 <  2:1206765/80‑1 (MQ=255)
                                                                            gcgcCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGaaa                                                                                                                 >  1:165896/1‑79 (MQ=255)
                                                                            gcgcCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGaaa                                                                                                                 <  2:165896/79‑1 (MQ=255)
                                                                            gcgcCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAAtt                                                     >  1:16075/1‑139 (MQ=255)
                                                                             cgctcgtcgtGAATGTGGGGCAAAACCTGGTCTTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTa                                                                         >  1:1046746/5‑118 (MQ=255)
                                                                             cgcTCGTCGTGAATGTGGGGCAAAACCTGGTCTTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTa                                                                         <  2:1046746/118‑1 (MQ=255)
                                                                                           gtgGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTaa                                      <  2:16075/139‑1 (MQ=255)
                                                                                              gggCAAAACCTGGTCGTGGATCTTTCGACGCAGATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAat                                   >  1:244161/1‑139 (MQ=255)
                                                                                              gggCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATaa                                                                                                                               >  2:851049/1‑47 (MQ=255)
                                                                                              gggCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATaa                                                                                                                               <  1:851049/47‑1 (MQ=255)
                                                                                              gggCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                            >  1:351507/1‑130 (MQ=255)
                                                                                              gggCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                            <  2:351507/130‑1 (MQ=255)
                                                                                              gggCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAat                                   >  1:928455/1‑139 (MQ=255)
                                                                                                 cAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAAtata                                <  2:928455/139‑1 (MQ=255)
                                                                                                 cAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAAtata                                <  1:24560/139‑1 (MQ=255)
                                                                                                      ccTGGTCGTGGATCTTTCGACGCAAATCTTTTACCa                                                                                                                                  <  1:761722/36‑1 (MQ=255)
                                                                                                      ccTGGTCGTGGATCTTTCGACGCAAATCTTTTACCa                                                                                                                                  >  2:761722/1‑36 (MQ=255)
                                                                                                      ccTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGc                           <  1:571872/139‑1 (MQ=255)
                                                                                                      ccTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGc                           <  1:616663/139‑1 (MQ=255)
                                                                                                      ccTGGTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGc                           <  1:1206867/139‑1 (MQ=255)
                                                                                                       cTGGTCGTGGATCTTTCGACGCAGATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCa                          <  2:244161/139‑1 (MQ=255)
                                                                                                         ggTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTca                                                                                               >  1:908654/1‑68 (MQ=255)
                                                                                                         ggTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTca                                                                                               <  2:908654/68‑1 (MQ=255)
                                                                                                         ggTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCt                                                                                <  1:1060344/83‑1 (MQ=255)
                                                                                                         ggTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCt                                                                                >  2:1060344/1‑83 (MQ=255)
                                                                                                          gTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGaaa                                                                                                                 >  2:1307449/1‑49 (MQ=255)
                                                                                                          gTCGTGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGaaa                                                                                                                 <  1:1307449/49‑1 (MQ=255)
                                                                                                              tGGATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTa                   <  2:1206456/139‑1 (MQ=255)
                                                                                                                gATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCgg                                              <  2:621548/110‑1 (MQ=255)
                                                                                                                gATCTTTCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCgg                                              >  1:621548/1‑110 (MQ=255)
                                                                                                                     ttCGACGCAAATCTTTTACCATAACGATTATCCGGaaa                                                                                                                 <  2:261922/38‑1 (MQ=255)
                                                                                                                     ttCGACGCAAATCTTTTACCATAACGATTATCCGGaaa                                                                                                                 >  1:261922/1‑38 (MQ=255)
                                                                                                                      tCGACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTc           >  2:682216/1‑139 (MQ=255)
                                                                                                                        gACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCt         >  2:638474/1‑139 (MQ=255)
                                                                                                                        gACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCt         <  1:510802/139‑1 (MQ=255)
                                                                                                                        gACGCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCt         <  1:235103/139‑1 (MQ=255)
                                                                                                                           gCAAATCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTacc      >  1:25900/1‑139 (MQ=255)
                                                                                                                               aTCTTTTACCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTAccacca  >  2:649558/1‑139 (MQ=255)
                                                                                                                                      |                                                                                                                                   
GTCATTCGCCTGTAAAACCGCCAATGGTACCGCTATCCCTATTGGCGGTGGCAGCGCCAATGTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCA  >  NC_000913/4548867‑4549132

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: