Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,048,429 C→T 16.7% W88* (TGG→TAG) ‡ gfcC ← putative O‑antigen capsule production periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,048,4290CT16.7% 219.9 / 33.4 96W88* (TGG→TAG) gfcCputative O‑antigen capsule production periplasmic protein
Reads supporting (aligned to +/- strand):  ref base C (40/40);  new base T (8/8);  total (48/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.89e-01

Only 100 of 101 total aligned reads displayed.

TCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGC  >  NC_000913/1048330‑1048528
                                                                                                   |                                                                                                     
tCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTccc                                                                                                       <  2:442211/100‑1 (MQ=255)
 ccAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCa                                                                                                      <  1:2448948/100‑1 (MQ=255)
  cAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAg                                                                                                     <  2:901090/100‑1 (MQ=255)
  cAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCCTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAg                                                                                                     <  1:706486/100‑1 (MQ=255)
   aGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAgg                                                                                                    <  2:802645/100‑1 (MQ=255)
    gTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAgg                                                                                                    <  2:2347179/99‑1 (MQ=255)
      ttGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGaag                                                                                                 <  2:2132766/100‑1 (MQ=255)
      ttGACCGGTAAGCGACCAGTGATATTGAGAGTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACAGCGTCGTCAGCTTCTGCTTCCCAGGgag                                                                                                 <  1:1467973/100‑1 (MQ=255)
       tGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGaagc                                                                                                <  2:165667/100‑1 (MQ=255)
       tGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGaagc                                                                                                >  1:1909072/1‑100 (MQ=255)
        gACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCTGGGaagc‑‑a                                                                                             <  2:1506829/100‑1 (MQ=255)
          ccGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGaagc‑‑aag                                                                                           >  2:1526006/1‑100 (MQ=255)
          ccGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGaagc‑‑aag                                                                                           <  1:2061891/100‑1 (MQ=255)
           cGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCTGGGaagc‑‑agaa                                                                                          <  1:878692/100‑4 (MQ=255)
           cGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGaagc‑‑aagc                                                                                          >  2:183500/1‑100 (MQ=255)
            ggTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCTGGGaagc‑‑agaag                                                                                         <  1:2208216/100‑5 (MQ=255)
            ggTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCTGGGaagc‑‑agaag                                                                                         <  1:1476801/100‑5 (MQ=255)
            ggTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCTGGGaagc‑‑agaag                                                                                         <  2:1552893/100‑5 (MQ=255)
            ggTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCt                                                                                         >  2:706307/1‑100 (MQ=255)
             gTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCTGGGaagc‑‑agaagc                                                                                        <  2:1261592/100‑6 (MQ=255)
             gTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCTGGGaagc‑‑agaagc                                                                                        <  2:1783059/100‑6 (MQ=255)
                aGCGACCAGTGATATTGAGATTGAGTAACTGCTGTCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGc                                                                                     <  1:1603750/100‑1 (MQ=255)
                aGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCTGGGaagc‑‑agaagctga                                                                                     <  1:533783/100‑9 (MQ=255)
                aGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGaagc‑‑aagc                                                                                          >  1:393159/1‑95 (MQ=255)
                   gACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCAt                                                                                  <  2:1493307/100‑1 (MQ=255)
                   gACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCAt                                                                                  <  1:356862/100‑1 (MQ=255)
                   gACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCAt                                                                                  >  1:281381/1‑100 (MQ=255)
                    aCCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATg                                                                                 <  2:550551/100‑1 (MQ=255)
                    aCCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATg                                                                                 <  2:1840439/100‑1 (MQ=255)
                          gATATTGAGATTGAGTAACTGCTGGCGCATGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCCGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGt                                                                           <  1:1552072/100‑1 (MQ=255)
                           atatTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGtt                                                                          <  1:1887755/100‑1 (MQ=255)
                            tatTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTg                                                                         >  2:1121820/1‑100 (MQ=255)
                             attgagattgagTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTgg                                                                        >  2:165084/1‑100 (MQ=255)
                             attgagattgagTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTgg                                                                        <  1:1523858/100‑1 (MQ=255)
                              ttgagattgagTAACGGCTGGCGCACGGATTTAATCGTTGCCGCTACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTgg                                                                        >  1:381293/1‑99 (MQ=255)
                                gagattgagTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAg                                                                     <  1:694432/100‑1 (MQ=255)
                                 agattgagTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGt                                                                    <  1:2023511/100‑1 (MQ=255)
                                  gattgagTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTc                                                                   <  1:2435230/100‑1 (MQ=255)
                                  gattgagTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTc                                                                   <  2:1496568/100‑1 (MQ=255)
                                  gattgagTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTc                                                                   >  2:495907/1‑100 (MQ=255)
                                   attgagTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCt                                                                  >  2:2055557/1‑100 (MQ=255)
                                     tgagTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCttt                                                                >  1:798395/1‑100 (MQ=255)
                                     tgagTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCtt                                                                 <  2:2117488/99‑1 (MQ=255)
                                      gagTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCTTTa                                                               >  1:1854027/1‑100 (MQ=255)
                                       agTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCTTTa                                                               >  1:1526617/1‑99 (MQ=255)
                                       agTAACGGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCTTTaa                                                              <  1:844074/100‑1 (MQ=255)
                                        gTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCTTTAAg                                                             <  1:1245387/100‑1 (MQ=255)
                                        gTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCTTTAAg                                                             <  1:1183295/100‑1 (MQ=255)
                                              gctgGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCttt                                                       >  1:2001006/1‑100 (MQ=255)
                                               ctgGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTa                                                      <  2:651119/100‑1 (MQ=255)
                                                tgGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAg                                                     >  1:1979920/1‑100 (MQ=255)
                                                   cgcACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCtt                                                  >  2:1236411/1‑100 (MQ=255)
                                                   agcACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTACCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCtt                                                  >  2:266007/2‑100 (MQ=255)
                                                   agcACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTACCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCtt                                                  >  2:266088/2‑100 (MQ=255)
                                                    gcACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCttt                                                 >  1:1127212/1‑100 (MQ=255)
                                                    gcACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCttt                                                 >  2:425370/1‑100 (MQ=255)
                                                      aCGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTcg                                               >  1:1759755/1‑100 (MQ=255)
                                                       cGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTcgc                                              >  2:783584/1‑100 (MQ=255)
                                                           tttAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCa                                          >  2:916160/1‑100 (MQ=255)
                                                           tttAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTCCCATGACATGTTGGTAGTCTTTAAGCGCTATAGCTTTCGCCGCa                                          >  1:1858728/1‑100 (MQ=255)
                                                             tAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCaga                                        >  1:1834200/1‑100 (MQ=255)
                                                             tAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCaga                                        >  2:1815717/1‑100 (MQ=255)
                                                                tCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTc                                     <  1:1150906/100‑1 (MQ=255)
                                                                tCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTc                                     <  1:1689708/100‑1 (MQ=255)
                                                                 cGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCg                                    <  1:998994/100‑1 (MQ=255)
                                                                 cGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTGCGCCGCAGAGTCg                                    <  1:1603434/100‑1 (MQ=255)
                                                                  gTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCgg                                   <  1:219117/100‑1 (MQ=255)
                                                                  gTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCgg                                   >  2:2256820/1‑100 (MQ=255)
                                                                   ttGCCGCCACATCGTCGTCAGCTTCTGCTGCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGGAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGt                                  <  1:1733255/100‑1 (MQ=255)
                                                                    tgtcgccACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTc                                 <  1:554490/97‑1 (MQ=255)
                                                                    tGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTc                                 <  2:2495231/100‑1 (MQ=255)
                                                                      ccgccACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAg                               >  1:1717218/1‑100 (MQ=255)
                                                                      ccgccACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAg                               >  2:2102841/1‑100 (MQ=255)
                                                                       cgccACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGc                              >  2:1928183/1‑100 (MQ=255)
                                                                         ccACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCaa                            >  2:407311/1‑100 (MQ=255)
                                                                            caTCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGc                         >  1:1985933/1‑100 (MQ=255)
                                                                                  cgtcAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAgg                   <  2:1253848/100‑1 (MQ=255)
                                                                                  cgtcAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAgg                   >  2:1669733/1‑100 (MQ=255)
                                                                                   gtcAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGGTGGTAGTCTTTAAGCACTGTAGCATTCGCTGCAGAGACGGGAAGCAAGGCACCAgga                  >  1:1985040/1‑99 (MQ=255)
                                                                                    tcAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGcc                 >  2:1530755/1‑100 (MQ=255)
                                                                                    tcAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGcc                 >  1:2222443/1‑100 (MQ=255)
                                                                                     cacagcttGCTTCCTGGGAAGCAGAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGc                  >  2:1065376/8‑100 (MQ=37)
                                                                                      acagcttGCTTCCTGGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGccac               >  2:1907344/7‑100 (MQ=255)
                                                                                       gCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGccacc              <  2:1274349/100‑1 (MQ=255)
                                                                                       cagcttGCTTCCTGGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGccacc              >  2:82593/6‑100 (MQ=255)
                                                                                       cagcttGCTTCCTGGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGccacc              >  1:297447/6‑100 (MQ=255)
                                                                                        agcttGCTTCCTGGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGccacca             >  1:2008820/5‑100 (MQ=255)
                                                                                         ttCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCag            >  1:1919749/1‑100 (MQ=255)
                                                                                         ttCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCag            >  1:1098026/1‑100 (MQ=255)
                                                                                         ttCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCag            <  1:2343903/100‑1 (MQ=255)
                                                                                         gcttGCTTCCTGGGAAGC‑‑AGGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCag            >  1:1440758/4‑100 (MQ=255)
                                                                                         gcttGCTTCCTGGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTTGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCag            >  1:786324/4‑100 (MQ=255)
                                                                                         gcttGCTTCCTGGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGACACCAGGCCACCag            >  2:842119/4‑100 (MQ=255)
                                                                                             gCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGCGGc        >  1:1246114/1‑100 (MQ=255)
                                                                                             gCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGc        <  1:1525997/100‑1 (MQ=255)
                                                                                              cTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCg       <  1:89105/100‑1 (MQ=255)
                                                                                                  ccAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATcg   >  1:360111/1‑100 (MQ=255)
                                                                                                  ccAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATcg   >  2:152074/1‑100 (MQ=255)
                                                                                                  ccAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATcg   >  1:1873510/1‑100 (MQ=255)
                                                                                                   cAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATcgc  >  2:816377/1‑100 (MQ=255)
                                                                                                   |                                                                                                     
TCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGC‑‑AAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGC  >  NC_000913/1048330‑1048528

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