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,0540TC97.6% 287.1 / ‑3.0 84S83P (TCG→CCG) fimHminor component of type 1 fimbriae
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (38/44);  total (39/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.48e-01

GTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGG  >  NC_000913/4548928‑4549186
                                                                                                                              |                                                                                                                                    
gTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATggcg                                                                                                                          <  1:56646/139‑1 (MQ=255)
  ttATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATggcggc                                                                                                                        <  2:623121/139‑1 (MQ=255)
  ttATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATggcggc                                                                                                                        <  1:144120/139‑1 (MQ=255)
     tGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCgtg                                                                                                                     <  2:572433/139‑1 (MQ=255)
         aaCCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTAt                                                                                                                 >  2:1492944/1‑139 (MQ=255)
           ccTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCt                                                                                                               >  2:1695094/1‑139 (MQ=255)
              tGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAAt                                                                                                            >  2:464004/1‑139 (MQ=255)
              tGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAAt                                                                                                            >  2:1421177/1‑139 (MQ=255)
               gcgcCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAAtt                                                                                                           >  2:1059946/1‑139 (MQ=255)
                cgcCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAAttt                                                                                                          >  1:1639790/1‑139 (MQ=255)
                      tcgtGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCgg                                                                                                    <  1:1791068/139‑1 (MQ=255)
                      tcgtGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCgg                                                                                                    <  2:1183133/139‑1 (MQ=255)
                      tcgtGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCgg                                                                                                    <  2:1345187/139‑1 (MQ=255)
                        gtGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                                  <  2:1006874/139‑1 (MQ=255)
                        gtGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                                  <  1:1607369/139‑1 (MQ=255)
                        gtGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                                  <  2:214931/139‑1 (MQ=255)
                           aaTGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCg                                                                                               <  2:620436/139‑1 (MQ=255)
                             tgtgGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTa                                                                                             <  1:255608/139‑1 (MQ=255)
                                ggggCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTaaaa                                                                                          >  1:213460/1‑139 (MQ=255)
                                 gggCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTa                                                                                                                               >  1:1535861/1‑101 (MQ=255)
                                 gggCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTa                                                                                                                               <  2:1535861/101‑1 (MQ=255)
                                 gggCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAat                                                                                         >  1:758149/1‑139 (MQ=255)
                                 gggCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAat                                                                                         >  1:466880/1‑139 (MQ=255)
                                 gggCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAat                                                                                         >  2:797434/1‑139 (MQ=255)
                                 gggCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAat                                                                                         >  1:1375331/1‑139 (MQ=255)
                                 gggCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAat                                                                                         >  2:652857/1‑139 (MQ=255)
                                    cAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAAtata                                                                                      <  1:464004/139‑1 (MQ=255)
                                         ccTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGc                                                                                 <  2:758149/139‑1 (MQ=255)
                                               cGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGc                                                                           <  1:444963/139‑1 (MQ=255)
                                                   gATTTTTCGACCCAAATTTTTTTGCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATc                                                                       <  1:105531/139‑1 (MQ=255)
                                                        ttCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                                  <  2:42306/107‑1 (MQ=255)
                                                        ttCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                                  >  1:42306/1‑107 (MQ=255)
                                                             cGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTac                                                             >  1:235389/1‑139 (MQ=255)
                                                              gCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTTGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTacc                                                            >  1:35389/1‑139 (MQ=255)
                                                                aaaTCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTaccac                                                          <  1:797434/139‑1 (MQ=255)
                                                                aaaTCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTaccac                                                          <  1:1695094/139‑1 (MQ=255)
                                                                 aaTCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                                  >  1:491405/1‑98 (MQ=255)
                                                                 aaTCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                                  <  2:491405/98‑1 (MQ=255)
                                                                       ttGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGaa                                                   >  2:1600277/1‑139 (MQ=255)
                                                                              aaCGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGTTTATGGCGGCGTGTTATCTATTTTTTCCGGGACCGTAAAATATCGTGTCAGTAGCTATCCATTTCCTACCACCAGCGAAAcgccgc                                            <  1:1708120/139‑1 (MQ=255)
                                                                              aaCGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                                  <  1:486623/85‑1 (MQ=255)
                                                                              aaCGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                                  >  2:486623/1‑85 (MQ=255)
                                                                                cGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTTGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCcgcgc                                          <  2:35389/139‑1 (MQ=255)
                                                                                cGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCcgcgc                                          <  1:1492944/139‑1 (MQ=255)
                                                                                 gATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCcgcgcg                                         <  1:1059946/139‑1 (MQ=255)
                                                                                 gATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCcgcgcg                                         <  2:466880/139‑1 (MQ=255)
                                                                                 gATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCcgcgcg                                         <  2:418805/139‑1 (MQ=255)
                                                                                 gATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCcgcgcg                                         <  2:1375331/139‑1 (MQ=255)
                                                                                           aaaCCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                                  <  1:1649524/72‑1 (MQ=255)
                                                                                           aaaCCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                                  >  2:1649524/1‑72 (MQ=255)
                                                                                           aaaCCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAAtata                                                                                      >  1:1391243/1‑84 (MQ=255)
                                                                                           aaaCCATTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAAtata                                                                                      <  2:1391243/84‑1 (MQ=255)
                                                                                                gttACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTacc                                                            <  2:1151360/104‑1 (MQ=255)
                                                                                                aTTACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTacc                                                            >  1:1151360/1‑105 (MQ=255)
                                                                                                  tACAGACTATGTCACACTGCAACGATGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCgttgtt                                    >  2:113510/1‑127 (MQ=255)
                                                                                                  tACAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCgttgtt                                    <  1:113510/127‑1 (MQ=255)
                                                                                                   aCAGACTATGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGaa                       <  1:1421177/139‑1 (MQ=255)
                                                                                                           tGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCgg                                                                                                    <  1:1365043/54‑1 (MQ=255)
                                                                                                           tGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCgg                                                                                                    >  2:1365043/1‑54 (MQ=255)
                                                                                                           tGTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAg               >  2:257063/1‑139 (MQ=255)
                                                                                                            gTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                                  <  1:1025207/55‑1 (MQ=255)
                                                                                                            gTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                                  <  1:1333287/55‑1 (MQ=255)
                                                                                                            gTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                                  >  2:1333287/1‑55 (MQ=255)
                                                                                                            gTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGa                                                                                                  >  2:1025207/1‑55 (MQ=255)
                                                                                                            gTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAg                                                                            >  2:341770/1‑77 (MQ=255)
                                                                                                            gTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAg                                                                            <  1:341770/77‑1 (MQ=255)
                                                                                                            gTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTacc                                                            <  1:617473/93‑1 (MQ=255)
                                                                                                            gTCACACTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTacc                                                            >  2:617473/1‑93 (MQ=255)
                                                                                                                 acTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCa                                                                     >  1:1586351/1‑79 (MQ=255)
                                                                                                                 acTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCa                                                                     <  2:1586351/79‑1 (MQ=255)
                                                                                                                 acTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCcgcgcg                                         <  1:660659/107‑1 (MQ=255)
                                                                                                                 acTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCcgcgcg                                         >  2:660659/1‑107 (MQ=255)
                                                                                                                 acTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGaa                       <  2:1294874/125‑1 (MQ=255)
                                                                                                                 acTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGaa                       >  1:1294874/1‑125 (MQ=255)
                                                                                                                 acTGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTgg         >  2:1526578/1‑139 (MQ=255)
                                                                                                                   tGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGaa                       >  2:793870/1‑123 (MQ=255)
                                                                                                                   tGCAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGaa                       <  1:793870/123‑1 (MQ=255)
                                                                                                                     cAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCa                                                                     >  1:1337704/1‑75 (MQ=255)
                                                                                                                     cAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCa                                                                     <  2:1337704/75‑1 (MQ=255)
                                                                                                                     cAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCgg     >  1:61730/1‑139 (MQ=255)
                                                                                                                     cAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCgg     >  2:1262586/1‑139 (MQ=255)
                                                                                                                     cAACGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCgg     >  2:1329037/1‑139 (MQ=255)
                                                                                                                        cGAGGCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTgg  >  2:403555/1‑139 (MQ=255)
                                                                                                                           ggCCCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAg                                                                            <  1:855552/62‑1 (MQ=255)
                                                                                                                              |                                                                                                                                    
GTTTATGTAAACCTTGCGCCCGTCGTGAATGTGGGGCAAAACCTGGTCGTGGATCTTTCGACGCAAATCTTTTGCCATAACGATTATCCGGAAACCATTACAGACTATGTCACACTGCAACGAGGCTCGGCTTATGGCGGCGTGTTATCTAATTTTTCCGGGACCGTAAAATATAGTGGCAGTAGCTATCCATTTCCTACCACCAGCGAAACGCCGCGCGTTGTTTATAATTCGAGAACGGATAAGCCGTGGCCGGTGG  >  NC_000913/4548928‑4549186

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: