Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,600,156 C→T 5.1% A12T (GCC→ACC)  yjjA ← putative DUF2501 family periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,600,1560CT 5.1% 329.5 / 5.5 100A12T (GCC→ACC) yjjAputative DUF2501 family periplasmic protein
Reads supporting (aligned to +/- strand):  ref base C (37/57);  new base T (2/3);  total (40/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.39e-01

Only 100 of 110 total aligned reads displayed.

AGCAAGTTAGTTAATGACGCGAGCGACCAACCGCCTTCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTA  >  NC_000913/4600011‑4600304
                                                                                                                                                 |                                                                                                                                                     
aGCAAGTTAGTTAATGACGCGAGCGACCAACCGCCTTCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGgcaca                                                                                                                                                    >  1:772635/1‑149 (MQ=255)
aGCAAGTTAGTTAATGACGCGAGCGACCAACCGCCTTCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGgcaca                                                                                                                                                    <  3:368273/149‑1 (MQ=255)
    agttagttAATGACGCGAGCGACCAACCGCCTTCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGgcacagcac                                                                                                                                                >  3:321804/1‑149 (MQ=255)
      ttagttaATGACGCGAGCGACCAACCGCCTTCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGTGCTGGCGGCAATGGcacagcaca‑g                                                                                                                                             >  1:333270/1‑149 (MQ=255)
      ttagttaATGACGCGAGCGACCAACCGCCTTCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGcacagcaca‑g                                                                                                                                             >  1:640576/1‑149 (MQ=255)
        agttaATGACGCGAGCGACCAACCGCCTTCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACa‑gaa                                                                                                                                           <  2:726702/149‑1 (MQ=255)
        agttaATGACGCGAGCGACCAACCGCCTTCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACa‑gaa                                                                                                                                           <  2:443700/149‑1 (MQ=255)
        agttaATGACGCGAGCGACCAACCGCCTTCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACa‑gaa                                                                                                                                           <  3:61652/149‑1 (MQ=255)
         gttaATGACGCGAGCGACCAACCGCCTTCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACa‑gaag                                                                                                                                          >  2:20678/1‑149 (MQ=255)
          ttaATGACGCGAGCGACCAACCGCCTTCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACa‑gaaga                                                                                                                                         <  2:419492/149‑1 (MQ=255)
          ttaATGACGCGAGCGACCAACCGCCTTCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACa‑gaaga                                                                                                                                         >  2:480910/1‑149 (MQ=255)
           taATGACGCGAGCGACCAACCGCCTTCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGAt                                                                                                                                        <  4:697082/149‑1 (MQ=255)
             aTGACGCGAGCGACCAACCGCCTTCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGt                                                                                                                                      <  3:639896/149‑1 (MQ=255)
              tGACGCGAGCGACCAACCGCCTTCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGtt                                                                                                                                     <  3:525678/149‑1 (MQ=255)
                 cgcgAGCGACCAACCGCCTTCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTaa                                                                                                                                  <  2:763484/149‑1 (MQ=255)
                  gcgagcgaCCAGCCGCCTTCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAAc                                                                                                                                 >  2:246105/1‑149 (MQ=255)
                     agctaCCATCCGCCTGCCTGTGTCGTGCTGTTTTTGTTGCCAAGTTCGCTGGCGGCGCTGGATAGCGCATCTTTCCAGGACGCAGCATGAACCCCGATGTAAATTAATTCGCTGGCGGCAACGGCACCGCACA‑GAAGATGTTTAACAGt                                                                                                                              <  4:235174/149‑1 (MQ=255)
                          ccAACCGCCTTCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTtca                                                                                                                         <  4:1687/149‑1 (MQ=255)
                          ccAACCGCCTTCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTtca                                                                                                                         <  1:87531/149‑1 (MQ=255)
                           cAACCGCCTTCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTtcat                                                                                                                        <  4:222546/149‑1 (MQ=255)
                           cAACCGCCTTCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTtcat                                                                                                                        <  4:835220/149‑1 (MQ=255)
                            aaCCGCCTTCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTtcatc                                                                                                                       >  2:687298/1‑149 (MQ=255)
                               cGCCTTCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTtcatcatc                                                                                                                    >  1:179486/1‑149 (MQ=255)
                               cGCCTTCCTGTGTCGTGCTGTTTTGGTTG‑CAAATTCGGTGGGGGCCCTGGAAAAGAGATTTTTCCCGGGCGCAACATGGACCCCGGGGGGAATTAATGGGCTTGGGGGCATGGGACAGCACAGGAGGATTTTTAAAGGTTttc‑tcatca                                                                                                                   >  1:786476/1‑149 (MQ=255)
                                gCCTTCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTtcatcatca                                                                                                                   >  3:667506/1‑149 (MQ=255)
                                 ccTTCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTaa                                                                                                                                  <  1:142706/133‑1 (MQ=255)
                                 ccTTCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTaa                                                                                                                                  >  2:142706/1‑133 (MQ=255)
                                    tCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTAc                                                                                                               <  3:682768/149‑1 (MQ=255)
                                       tgtgTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCa                                                                                                            <  1:302520/149‑1 (MQ=255)
                                       tgtgTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCa                                                                                                            <  3:347759/149‑1 (MQ=255)
                                         tgtCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTCGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAg                                                                                                          <  2:93041/149‑1 (MQ=255)
                                              tGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAGATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGa                                                                                                     >  2:428348/1‑149 (MQ=255)
                                              tGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGa                                                                                                     >  4:641803/1‑149 (MQ=255)
                                               gCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGaa                                                                                                    <  4:667506/149‑1 (MQ=255)
                                                 tGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGTCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAAtt                                                                                                  <  1:731295/149‑1 (MQ=255)
                                                 tGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAAtt                                                                                                  >  3:674528/1‑149 (MQ=255)
                                                 tGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAAtt                                                                                                  >  3:666096/1‑149 (MQ=255)
                                                   ttttGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATtgt                                                                                                <  4:215040/149‑1 (MQ=255)
                                                   ttttGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATtgt                                                                                                >  1:8976/1‑149 (MQ=255)
                                                     ttGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTgtgt                                                                                              <  1:373004/149‑1 (MQ=255)
                                                         ttGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGTACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCa                                                                                          <  3:730668/149‑1 (MQ=255)
                                                         ttGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCa                                                                                          >  1:139993/1‑149 (MQ=255)
                                                           gCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGt                                                                                        >  1:95819/1‑149 (MQ=255)
                                                           gCCAAGTTCGCTGGCGGCGCTGGAGAGCGCAGCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGt                                                                                        >  4:539913/1‑149 (MQ=255)
                                                             cAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGtat                                                                                      >  2:540588/1‑149 (MQ=255)
                                                             cAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGtat                                                                                      >  2:735944/1‑149 (MQ=255)
                                                              aaGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTgtgt                                                                                              <  4:317390/140‑1 (MQ=255)
                                                              aaGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTgtgt                                                                                              >  3:317390/1‑140 (MQ=255)
                                                                   cGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGcc                                                                                <  2:8976/149‑1 (MQ=255)
                                                                    gCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCa                                                                               <  4:350360/149‑1 (MQ=255)
                                                                    gCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCa                                                                               <  1:873852/149‑1 (MQ=255)
                                                                      tGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGa                                                                             >  2:760286/1‑149 (MQ=255)
                                                                           gcgcTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACgg                                                                        <  4:834421/149‑1 (MQ=255)
                                                                           gcgcTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACgg                                                                        <  3:46631/149‑1 (MQ=255)
                                                                            cgcTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACggg                                                                       <  3:837259/149‑1 (MQ=255)
                                                                              cTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGa                                                                     >  3:233784/1‑149 (MQ=255)
                                                                               tGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGAt                                                                    <  3:627608/149‑1 (MQ=255)
                                                                                  agagCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTaa                                                                 <  1:592911/149‑1 (MQ=255)
                                                                                    agCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATa                                                               <  2:689590/149‑1 (MQ=255)
                                                                                     gcgcATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATAc                                                              >  1:782064/1‑149 (MQ=255)
                                                                                       gcATCTTTCCAGGACGCAGCCTGCACCCCGGTGGCAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACCCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATActc                                                            <  4:685985/149‑1 (MQ=255)
                                                                                         aTCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACgg                                                                        <  3:62791/135‑1 (MQ=255)
                                                                                         aTCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACgg                                                                        >  4:62791/1‑135 (MQ=255)
                                                                                         aTCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCtt                                                          <  3:627704/149‑1 (MQ=255)
                                                                                         aTCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCtt                                                          <  4:666096/149‑1 (MQ=255)
                                                                                         aTCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCtt                                                          >  4:254087/1‑149 (MQ=255)
                                                                                            tttCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTAc                                                       >  1:3819/1‑149 (MQ=255)
                                                                                            tttCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTAc                                                       >  4:31401/1‑149 (MQ=255)
                                                                                            tttCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTAc                                                       >  4:113762/1‑149 (MQ=255)
                                                                                            tttCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTAc                                                       >  4:38727/1‑149 (MQ=255)
                                                                                              tCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCt                                                     <  3:38727/149‑1 (MQ=255)
                                                                                                 aGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGTACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATAc                                                              >  4:565399/1‑137 (MQ=255)
                                                                                                 aGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGTACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATAc                                                              <  1:497968/137‑1 (MQ=255)
                                                                                                 aGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGTACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATAc                                                              >  2:497968/1‑137 (MQ=255)
                                                                                                 aGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGTACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATAc                                                              <  3:565399/137‑1 (MQ=255)
                                                                                                        agcaTGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGc                                           <  1:735944/149‑1 (MQ=255)
                                                                                                          caTGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTa                                         <  1:511238/149‑1 (MQ=255)
                                                                                                           aTGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGgca                                                                                                                                                      >  2:640699/1‑40 (MQ=255)
                                                                                                           aTGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGgca                                                                                                                                                      <  1:640699/40‑1 (MQ=255)
                                                                                                               aCCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCATAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTa                                    <  3:539913/149‑1 (MQ=255)
                                                                                                               aCCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTa                                    >  3:444025/1‑149 (MQ=255)
                                                                                                                      tGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGGTGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACggg                                                                       >  4:509344/1‑107 (MQ=255)
                                                                                                                      tGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACggg                                                                       <  3:509344/107‑1 (MQ=255)
                                                                                                                      tGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATcc                             <  2:516417/149‑1 (MQ=255)
                                                                                                                        gAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCa                           >  4:767238/1‑149 (MQ=255)
                                                                                                                          aaTTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGTAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGt                         >  1:306915/1‑149 (MQ=255)
                                                                                                                           aTTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGTAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGtt                        <  2:306915/149‑1 (MQ=255)
                                                                                                                           aTTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGtt                        >  4:268401/1‑149 (MQ=255)
                                                                                                                              aaTGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGaa                     <  4:799357/149‑1 (MQ=255)
                                                                                                                               aTGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAg                    <  1:333407/149‑1 (MQ=255)
                                                                                                                               aTGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAg                    <  4:674528/149‑1 (MQ=255)
                                                                                                                                 gcgcTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGAt                  >  2:20379/1‑149 (MQ=255)
                                                                                                                                   gcTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCa                <  2:870046/149‑1 (MQ=255)
                                                                                                                                    cTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCAc               <  2:744155/149‑1 (MQ=255)
                                                                                                                                    cTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCAc               >  3:219833/1‑149 (MQ=255)
                                                                                                                                            aaTGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAAc       <  2:179486/149‑1 (MQ=255)
                                                                                                                                               ggCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTGCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGt    <  3:90847/149‑1 (MQ=255)
                                                                                                                                                gcacagcaca‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGtt   <  1:246105/149‑1 (MQ=255)
                                                                                                                                                 cacagcaca‑gAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTa  <  4:554209/149‑1 (MQ=255)
                                                                                                                                                 cacagcaca‑gAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTa  <  4:852030/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
AGCAAGTTAGTTAATGACGCGAGCGACCAACCGCCTTCCTGTGTCGTGCTGTTTTGGTTGCCAAGTTCGCTGGCGGCGCTGGAGAGCGCATCTTTCCAGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACA‑GAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTA  >  NC_000913/4600011‑4600304

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