Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,188,968 G→A A648V (GCC→GTC)  yehB ← putative fimbrial usher protein YehB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,188,9680GA98.6% 229.9 / ‑4.4 69A648V (GCC→GTC) yehBputative fimbrial usher protein YehB
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (38/30);  total (39/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

TTTGTTGTGCGATATTTTTGCCCATTGACATAAGCATCTTTAATTCCTGGCGCATTCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGGCTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTCCAGGTCAAATTCGCCCCAGCTGTCGTTTCATTTCCCTGATGTTGATGACTCAGGTTGACACCATAATTGAAC  >  NC_000913/2188833‑2189105
                                                                                                                                        |                                                                                                                                         
tttGTTGTGCGATATTTTTGCCCATTGACATAAGCATCTTTAATTCCTGGCGCATTCATCACAGCAAACGTTTCGGAAATACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGACTcc                                                                                                                                       >  1:753833/1‑140 (MQ=255)
tttGTTGTGCGATATTTTTGCCCATTGACATAAGCATCTTTAATTCCTGGCGCATTCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGACTcc                                                                                                                                       >  2:548871/1‑140 (MQ=255)
tttGTTGTGCGATATTTTTGCCCATTGACATAAGCATCTTTAATTCCTGGCGCATTCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGACTcc                                                                                                                                       >  1:340343/1‑140 (MQ=255)
tttGTTGTGCGATATTTTTGCCCATTGACATAAGCATCTTTAATTCCTGGCGCATTCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGACTcc                                                                                                                                       >  2:105806/1‑140 (MQ=255)
 ttgttgTGCGATATTTTTGCCCATTGACATAAGCATCTTTAATTCCTGGCGCATTCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGACTCCa                                                                                                                                      <  2:524105/140‑1 (MQ=255)
  tgttgtGCGATATTTTTGCCCATTGACATAAGCATCTTTAATTCCTGGCGCATTCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGATAATGCCCCCTGAAACACTGACTCCAg                                                                                                                                     >  2:861586/1‑140 (MQ=255)
  tgttgtGCGATATTTTTGCCCATTGACATAAGCATCTTTAATTCCTGGCGCATTCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGACTCCAg                                                                                                                                     >  1:663061/1‑140 (MQ=255)
   gttgtGCGATATTTTTGCCCATTGACATAAGCATCTTTAATTCCTGGCGCATTCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGACTCCAGc                                                                                                                                    <  2:940824/140‑1 (MQ=255)
    ttgtGCGATATTTTTGCCCATTGACATAAGCATCTTTAATTCCTGGCGCATTCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGACTCCAGcc                                                                                                                                   >  2:583653/1‑140 (MQ=255)
      gtgCGATATTTTTGCCCATTGACATAAGCATCTTTAATTCCTGGCGCATTCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGACTCCAGCCTg                                                                                                                                 <  1:351241/140‑1 (MQ=255)
         cGATATTTTTGCCCATTGACATAAGCATCTTTAATTCCTGGCGCATTCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGACTCCAGCCTGTCg                                                                                                                              >  2:756368/1‑140 (MQ=255)
          aatatTTTTGCCATTTGAATAAGGCTTCTTTAATCCCTGGCGCATTCATCTCACCAACCGTTTCGCAAAACCGGTCCGCCGGTTTACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGACTCCAGCCTGTCGAt                                                                                                                            <  2:193715/139‑1 (MQ=255)
           atatTTTTGCCCATTGACATAAGCATCTTTAATTCCTGGCGCATTCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCCCCCGACCAGGCGACAATGCCCCCTGAAACACTGACTCCAGCCTGTCGAt                                                                                                                            <  2:979406/140‑1 (MQ=255)
           atatTTTTGCCCATTGACATAAGCATCTTTAATTCCTGGCGCATTCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGACTCCAGCCTGTCGAt                                                                                                                            <  1:373525/140‑1 (MQ=255)
           atatTTTTGCCCATTGACATAAGCATCTTTAATTCCTGGCGCATTCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGACTCCAGCCTGTCGAt                                                                                                                            <  1:597/140‑1 (MQ=255)
                  tGCCCATTGACATAAGCATCTTTAATTCCTGGCGCATTCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGACTCCAGCCTGTCGATAAGTACt                                                                                                                     <  1:965187/140‑1 (MQ=255)
                    cccATTGACATAAGCATCTTTAATTCCTGGCGCATTCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGACTCCAGCCTGTCGATAAGTACTCg                                                                                                                   >  2:968733/1‑140 (MQ=255)
                     ccATTGACATAAGCATCTTTAATTCCTGGCGCATTCATCACAGCAAACGTTTCGGAAATACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGACTCCAGCCTGTCGATAAGTACTCga                                                                                                                  <  2:753833/140‑1 (MQ=255)
                     ccATTGACATAAGCATCTTTAATTCCTGGCGCATTCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGACTCCAGCCTGTCGATAAGTACTCga                                                                                                                  <  1:318244/140‑1 (MQ=255)
                      cATTGACATAAGCATCTTTAATTCCTGGCGCATTCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGACTCCAGCCTGTCGATAAGTACTCgac                                                                                                                 <  1:419071/140‑1 (MQ=255)
                      cATTGACATAAGCATCTTTAATTCCTGGCGCATTCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGACTCCAGCCTGTCGATAAGTACTCgac                                                                                                                 >  2:85428/1‑140 (MQ=255)
                        ttGACATAAGCATCTTTAATTCCTGGCGCATTCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGATAATGCCCCCTGAAACACTGACTCCAGCCTGTCGATAAGTACTCgactg                                                                                                               <  1:861586/140‑1 (MQ=255)
                        ttGACATAAGCATCTTTAATTCCTGGCGCATTCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGACTCCAGCCTGTCGATAAGTACTCgactg                                                                                                               <  1:27340/140‑1 (MQ=255)
                             aTAAGCATCTTTAATTCCTGGCGCATTCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGACTCCAGCCTGTCGATAAGTACTCGACTGACtat                                                                                                          <  2:150977/140‑1 (MQ=255)
                                aGCATCTTTAATTCCTGGCGCATTCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGACTCCAGCCTGTCGATAAGTACTCGACTGACTATAAc                                                                                                       <  2:543951/140‑1 (MQ=255)
                                   cATCTTTAATTCCGGGCGCATTCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGACTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTg                                                                                                     <  2:609598/140‑1 (MQ=255)
                                    aTCTTTAATTCCTGGCGCATTCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGACTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGc                                                                                                    >  2:266942/1‑140 (MQ=255)
                                     tCTTTAATTCCTGGCGCATTCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGACTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGcc                                                                                                   >  1:352783/1‑140 (MQ=255)
                                      cTTTAATTCCTGGCGCATTCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGACTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCa                                                                                                  >  1:165297/1‑140 (MQ=255)
                                      cTTTAATTCCTGGCGCATTCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGACTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCa                                                                                                  >  2:314064/1‑140 (MQ=255)
                                          aaTTCCTGGCGCATTCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACTCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGACTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCa                                                                                              <  1:577320/140‑1 (MQ=255)
                                             tCCTGGCGCATTCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGACTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTg                                                                                           >  1:1021271/1‑140 (MQ=255)
                                             tCCTGGCGCATTCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGACTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTg                                                                                           >  1:985850/1‑140 (MQ=255)
                                              ccTGGCGCATTCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGACTCCAGCCTGTCGATAAGTACTCGACTGACtat                                                                                                          >  2:407081/1‑124 (MQ=255)
                                              ccTGGCGCATTCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGACTCCAGCCTGTCGATAAGTACTCGACTGACtat                                                                                                          <  1:407081/124‑1 (MQ=255)
                                                 ggCGCATTCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGACTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGc                                                                                       >  1:705517/1‑140 (MQ=255)
                                                 ggCGCATTCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGACTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGc                                                                                       >  2:660332/1‑140 (MQ=255)
                                                         catcaCAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGACTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGgcg                                                                               <  1:266942/140‑1 (MQ=255)
                                                                    aCGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGACTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTCCAGGTc                                                                    >  2:717271/1‑140 (MQ=255)
                                                                      gTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGACTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTCCAGGTCaa                                                                  >  2:248613/1‑140 (MQ=255)
                                                                         tCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGACTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTCCAGGTCAAAtt                                                               >  1:326385/1‑140 (MQ=255)
                                                                         tCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGACTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTCCAGGTCAAAtt                                                               >  2:835401/1‑140 (MQ=255)
                                                                         tCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGACTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTCCAGGTCAAAtt                                                               >  1:842478/1‑140 (MQ=255)
                                                                              aaGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGACTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTCCAGGTCAAATTCGccc                                                          <  2:111093/140‑1 (MQ=255)
                                                                                gACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTTACTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTCCAGGTCAAATTCGCCCCa                                                        >  1:886944/1‑140 (MQ=255)
                                                                                gACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGACTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTCCAGGTCAAATTCGCCCCa                                                        >  2:840059/1‑140 (MQ=255)
                                                                                gACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGACTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTCCAGGTCAAATTCGCCCCa                                                        >  1:522495/1‑140 (MQ=255)
                                                                                gACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGACTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTCCAGGTCAAATTCGCCCCa                                                        >  1:984797/1‑140 (MQ=255)
                                                                                gACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGACTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTCCAGGTCAAATTCGCCCCa                                                        >  2:167166/1‑140 (MQ=255)
                                                                                     ttCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGACTCCAGCCTGACGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTCCAGGTCAAATTCGCCCCAGCTGt                                                   >  2:452434/1‑140 (MQ=255)
                                                                                       cGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGACTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTCCAGGTCAAATTCGCCCCAGCTGTCg                                                 <  2:541694/140‑1 (MQ=255)
                                                                                       cGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGACTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTCCAGGTCAAATTCGCCCCAGCTGTCg                                                 >  1:588726/1‑140 (MQ=255)
                                                                                               tAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGACTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTCCAGGTCAAATTCGCCCCAGCTGTCGTTTCAttt                                         >  2:1021752/1‑140 (MQ=255)
                                                                                                 aCGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGACTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTCCAGGTCAAATTCGCCCCAGCTGTCGTTTCATTTcc                                       >  1:821323/1‑140 (MQ=255)
                                                                                                      aCCCGACCAGGCGACAATGCCCCCTGAAACACTGACTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTCCAGGTCAAATTCGCCCCAGCTGTCGTTTCATTTCCCTGAt                                  >  1:369086/1‑140 (MQ=255)
                                                                                                       cccGACCAGGCGACAATGCCCCCTGAAACACTGACTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTCCAGGTCAAATTCGCCCCAGCTGTCGTTTCATTTccc                                      <  1:706142/135‑1 (MQ=255)
                                                                                                       cccGACCAGGCGACAATGCCCCCTGAAACACTGACTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTCCAGGTCAAATTCGCCCCAGCTGTCGTTTCATTTccc                                      >  2:706142/1‑135 (MQ=255)
                                                                                                            ccAGGCGACAATGCCCCCTGAAACACTGACTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTCCAGGTCAAATTCGCCCCAGCTGTCGTTTCATTTCCCTGATGTtgat                            <  2:690871/140‑1 (MQ=255)
                                                                                                            ccAGGCGACAATGCCCCCTGAAACACTGACTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTCCAGGTCAAATTCGCCCCAGCTGTCGTTTCATTTCCCTGATGTtgat                            >  2:688215/1‑140 (MQ=255)
                                                                                                              aGGCGACAATACCCCCTGAAACACTGACTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTCCAGGTCAAATTCGCCCCAGCTGTCGTTTCATTTCCCTGATGTtgatga                          <  1:915063/140‑1 (MQ=255)
                                                                                                                        gCCCCCTGAAACGCTGACTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTCCAGGTCAAATTCGCCCCAGCTGTCGTTTCATTTCCCTGATGTTGATGACTCAGGTTGa                <  2:622921/140‑1 (MQ=255)
                                                                                                                        gCCCCCTGAAACACTGACTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTCCAGGTCAAATTCGCCCCAGCTGTCGTTTCATTTCCCTGATGTTGATGACTCAGGTTGa                <  2:524122/140‑1 (MQ=255)
                                                                                                                         cccccTGAAACACTGACTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTCCAGGTCAAATTCGCCCCAGCTGTCGTTTCATTTCCCTGATGTTGATGACTCAGGTTGac               >  1:773701/1‑140 (MQ=255)
                                                                                                                           cccTGAAACACTGACTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTCCAGGTCAAATTCGCCCCAGCTGTCGTTTCATTTCCCTGATGTTGATGACTCAGGTTGacac             <  2:842478/140‑1 (MQ=255)
                                                                                                                                 aaCACTGACTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTCCAGGTCAAATTCGCCCCAGCTGTCGTTTCATTTCCCTGATGTTGATGACTCAGGTTGACACCATAAt       <  2:326385/140‑1 (MQ=255)
                                                                                                                                 aaCACTGACTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTCCAGGTCAAATTCGCCCCAGCTGTCGTTTCATTTCCCTGATGTTGATGACTCAGGTTGACACCATAAt       >  1:152871/1‑140 (MQ=255)
                                                                                                                                  acacTGACTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTCCAGGTCAAATTCGCCCCAGCTGTCGTTTCATTTCCCTGATGTTGATGACTCAGGTTGACACCATAAtt      <  1:1025102/140‑1 (MQ=255)
                                                                                                                                     cTGACTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTCCAGGTCAAATTCGCCCCAGCTGTCGTTTCATTTCCCTGATGTTGATGACTCAGGTTGACACCATAATTGaa   <  2:877297/140‑1 (MQ=255)
                                                                                                                                      tGGCTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTCCAGGTCAAATTCGCCCCAGCTGTCGTTTCATTTCCCTGATGTTGATGACTCAGGTTGACACCATAATTGAAc  >  1:542585/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                         
TTTGTTGTGCGATATTTTTGCCCATTGACATAAGCATCTTTAATTCCTGGCGCATTCATCACAGCAAACGTTTCGGAAAGACGGTTCGCCAGATTAACGCCACCCGACCAGGCGACAATGCCCCCTGAAACACTGGCTCCAGCCTGTCGATAAGTACTCGACTGACTATAACTGCCATTCACTGTCGCAACCGGCGCGTTCCAGGTCAAATTCGCCCCAGCTGTCGTTTCATTTCCCTGATGTTGATGACTCAGGTTGACACCATAATTGAAC  >  NC_000913/2188833‑2189105

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