Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,588,476 T→A 6.4% M122L (ATG→TTG)  ydeS ← putative fimbrial‑like adhesin protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,588,4760TA 6.4% 587.7 / 19.5 169M122L (ATG→TTG) ydeSputative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base T (74/84);  new base A (7/4);  total (81/88)
Fisher's exact test for biased strand distribution p-value = 3.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 176 total aligned reads displayed.

TTGCCGAAGCCCTCACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATC  >  NC_000913/1588332‑1588624
                                                                                                                                                |                                                                                                                                                    
ttGCCGAAGCCCTCACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCaaatat                                                                                                                                                  <  3:707224/149‑3 (MQ=255)
  gCCGAAGCCCTCACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGc                                                                                                                                                <  4:819345/149‑1 (MQ=255)
    cGAAGCCCTCACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCAAATATGCGt                                                                                                                                              <  2:1519004/149‑1 (MQ=255)
     gAAGCCCTCACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTc                                                                                                                                             >  2:1047846/1‑149 (MQ=255)
     gAAGCCCTCACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTc                                                                                                                                             >  2:821170/1‑149 (MQ=255)
      aaGCCCTCACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTcc                                                                                                                                            >  1:379798/1‑149 (MQ=255)
       aGCCCTCACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCAAATATGCGTccg                                                                                                                                           <  2:1120413/149‑2 (MQ=255)
           cTCACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAatat                                                                                                                                       >  3:826739/1‑149 (MQ=255)
           cTCACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAatat                                                                                                                                       >  2:1638708/1‑149 (MQ=255)
             cACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATAttt                                                                                                                                     <  3:1703101/149‑1 (MQ=255)
              aCCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTc                                                                                                                                    >  2:1430524/1‑149 (MQ=255)
                cAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCta                                                                                                                                  >  4:1234912/1‑149 (MQ=255)
                  gTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCtata                                                                                                                                <  1:1575416/149‑1 (MQ=255)
                     ccGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATAccc                                                                                                                             <  1:353453/149‑1 (MQ=255)
                     ccGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATAccc                                                                                                                             >  4:923187/1‑149 (MQ=255)
                            gACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCgg                                                                                                                      <  4:1616073/149‑1 (MQ=255)
                             aCGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAATCGGa                                                                                                                     <  4:883884/149‑1 (MQ=255)
                             aCGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGa                                                                                                                     >  2:588010/1‑149 (MQ=255)
                              cGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATGTTGTTCTGTTCTGGTACCAGTGGGAGCCACTGCCTCCCGGCATGAAGATGATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAg                                                                                                                    <  4:1447502/149‑1 (MQ=255)
                               ggACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAgg                                                                                                                   >  3:311817/1‑149 (MQ=255)
                                gACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATAttt                                                                                                                                     <  2:1096245/130‑1 (MQ=255)
                                 aCTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGcc                                                                                                                 <  4:1279294/149‑1 (MQ=255)
                                 aCTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGcc                                                                                                                 <  4:689788/149‑1 (MQ=255)
                                   ttCTGAGTTGACTTCAGTCGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACAGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGt                                                                                                               >  4:1067744/1‑149 (MQ=255)
                                   ttCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGt                                                                                                               >  1:891777/1‑149 (MQ=255)
                                      tGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTAtt                                                                                                            <  3:1276101/149‑1 (MQ=255)
                                       gAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTg                                                                                                           >  3:71681/1‑149 (MQ=255)
                                        aGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGc                                                                                                          <  1:108898/149‑1 (MQ=255)
                                          ttGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGctt                                                                                                        >  1:92876/1‑149 (MQ=255)
                                          ttGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGctt                                                                                                        <  1:991796/149‑1 (MQ=255)
                                            gACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGcttcc                                                                                                      <  3:1435082/149‑1 (MQ=255)
                                             aCTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGcttcct                                                                                                     <  4:642674/149‑1 (MQ=255)
                                             aCTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGcttcct                                                                                                     <  1:1027052/149‑1 (MQ=255)
                                              cTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGcttcctt                                                                                                    <  4:390107/149‑1 (MQ=255)
                                               ttCAGACGAGCGGAGTACGGCAATATTTTTTTCTGTTCTGGTACCAGTGGGAGCGACTGCATGCCGGCAGGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACTGGAGGCCGTATTGcttccttc                                                                                                   >  2:644434/1‑149 (MQ=255)
                                               ttCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCgg                                                                                                                      >  3:215025/1‑130 (MQ=255)
                                               ttCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGcttccttc                                                                                                   >  4:28806/1‑149 (MQ=255)
                                               ttCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGcttccttc                                                                                                   >  3:634152/1‑149 (MQ=255)
                                               ttCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGcttccttc                                                                                                   >  2:495093/1‑149 (MQ=255)
                                                       aGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGtt                                                                                           >  2:336379/1‑149 (MQ=255)
                                                          ggAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTca                                                                                        <  1:1212670/149‑1 (MQ=255)
                                                           gAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGcttcct                                                                                                     >  1:521546/1‑135 (MQ=255)
                                                           gAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTcaa                                                                                       <  4:1648011/149‑1 (MQ=255)
                                                            aGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTcaac                                                                                      <  1:18184/149‑1 (MQ=255)
                                                              tAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTcaacaa                                                                                    >  4:213892/1‑149 (MQ=255)
                                                               aaGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACaaa                                                                                   <  2:220314/149‑1 (MQ=255)
                                                                aGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAg                                                                                  >  1:1224545/1‑149 (MQ=255)
                                                                 ggCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGt                                                                                 >  4:414431/1‑149 (MQ=255)
                                                                     aaaTATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGtatta                                                                             >  1:740476/1‑149 (MQ=255)
                                                                       atatTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTattatt                                                                           <  4:562773/149‑1 (MQ=255)
                                                                             ttctgttctgGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGc                                                                     >  4:1554632/1‑149 (MQ=255)
                                                                                 gttctgGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATcc                                                                 <  4:1699020/149‑1 (MQ=255)
                                                                                     tgGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCaa                                                             >  1:1471121/1‑149 (MQ=255)
                                                                                     tgGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTACACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCaa                                                             <  1:819940/149‑1 (MQ=255)
                                                                                        tACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACaaa                                                                                   >  1:1466560/1‑124 (MQ=255)
                                                                                         aCCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATAcc                                                         <  3:923187/149‑1 (MQ=255)
                                                                                          ccAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCg                                                        >  1:283973/1‑149 (MQ=255)
                                                                                            aGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTcaa                                                                                       >  2:1319747/1‑116 (MQ=255)
                                                                                             gTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTa                                                     <  3:414431/149‑1 (MQ=255)
                                                                                               gggATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTaaa                                                   <  2:567883/149‑1 (MQ=255)
                                                                                               gggATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTaaa                                                   <  4:1334761/149‑1 (MQ=255)
                                                                                               gggATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTaaa                                                   >  4:19737/1‑149 (MQ=255)
                                                                                                ggATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAAt                                                  <  2:117803/149‑1 (MQ=255)
                                                                                                ggATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAAt                                                  <  3:836145/149‑1 (MQ=255)
                                                                                                ggATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAAt                                                  <  2:359495/149‑1 (MQ=255)
                                                                                                   tCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCa                                               >  1:340344/1‑149 (MQ=255)
                                                                                                    ccACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCaa                                              <  3:1336868/149‑1 (MQ=255)
                                                                                                    ccACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCaa                                              <  2:318828/149‑1 (MQ=255)
                                                                                                     cACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAAc                                             <  4:757153/149‑1 (MQ=255)
                                                                                                      aCTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAAcc                                            >  3:595101/1‑149 (MQ=255)
                                                                                                      aCTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAAcc                                            >  3:879149/1‑149 (MQ=255)
                                                                                                              ccGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAAcccc                                    <  3:500484/149‑1 (MQ=255)
                                                                                                              ccGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAAcccc                                    >  4:804106/1‑149 (MQ=255)
                                                                                                                    tGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGt                              <  3:19737/149‑1 (MQ=255)
                                                                                                                        gATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTcaa                                                                                       >  2:1759125/1‑88 (MQ=255)
                                                                                                                        gATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCt                                                                       >  2:905845/1‑104 (MQ=255)
                                                                                                                           cATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTtg                       <  1:1100065/149‑1 (MQ=255)
                                                                                                                            aTTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTg                               <  4:360419/140‑1 (MQ=255)
                                                                                                                            aTTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTg                               >  3:360419/1‑140 (MQ=255)
                                                                                                                             ttCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTtgct                     >  1:1696283/1‑149 (MQ=255)
                                                                                                                                 gTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCAt                 <  4:231921/149‑1 (MQ=255)
                                                                                                                                  tttCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCGGTTGTCCCTTTGCTGCAtt                >  4:1176651/1‑149 (MQ=255)
                                                                                                                                   ttCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTc               <  3:488980/149‑1 (MQ=255)
                                                                                                                                     cACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCAc             <  2:1224545/149‑1 (MQ=255)
                                                                                                                                     cACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCAc             <  1:145697/149‑1 (MQ=255)
                                                                                                                                       ttggACGCAAATATGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACtt           >  1:1239597/3‑149 (MQ=255)
                                                                                                                                       ccGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACtt           <  3:392283/149‑1 (MQ=255)
                                                                                                                                         ggACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTaa         <  4:1358664/149‑1 (MQ=255)
                                                                                                                                         ggACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTaa         <  1:1689943/149‑1 (MQ=255)
                                                                                                                                         ggACGCAAATATGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTaa         >  1:491841/1‑149 (MQ=255)
                                                                                                                                          ggcGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAAc        <  4:827553/147‑1 (MQ=255)
                                                                                                                                          gACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGCACCCCCGTTGTCCCTTTGCTGCATTCACTTAAc        <  2:1471121/149‑1 (MQ=255)
                                                                                                                                          gACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAAc        <  1:1047846/149‑1 (MQ=255)
                                                                                                                                             gCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTa     >  1:343798/1‑149 (MQ=255)
                                                                                                                                             gCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTa     >  2:1671363/1‑149 (MQ=255)
                                                                                                                                             gCATATTTGCGTCCAATATTTCTATACCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTaa    <  4:1303409/149‑1 (MQ=255)
                                                                                                                                              cATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTaa         <  1:1332469/144‑1 (MQ=255)
                                                                                                                                              cATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTaa         >  2:1332469/1‑144 (MQ=255)
                                                                                                                                              cATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTaa    <  2:740476/149‑1 (MQ=255)
                                                                                                                                                tatTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATc  >  3:1564166/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
TTGCCGAAGCCCTCACCAGTCCCGGATTGACGGACTTCTGAGTTGACTTCAGACGAGCGGAGTAAGGCAAAATATTGTTCTGTTCTGGTACCAGTGGGATCCACTGCATCCCGGCATGAAGATCATTCAGTTTCACCGGACGCATATTTGCGTCCAATATTTCTATACCCAAACCGGAGGCCGTATTGCTTCCTTCATCCAGTTTCAACAAAGTATTATTTTCTGCATCCTCAATACCGTTAAATGCAACCCGAACCCCCGTTGTCCCTTTGCTGCATTCACTTAACGTAATC  >  NC_000913/1588332‑1588624

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