Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 560,553 C→T A286V (GCC→GTC)  sfmD → putative outer membrane export usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913560,5530CT96.6% 388.1 / 0.8 118A286V (GCC→GTC) sfmDputative outer membrane export usher protein
Reads supporting (aligned to +/- strand):  ref base C (3/1);  new base T (61/53);  total (64/54)
Fisher's exact test for biased strand distribution p-value = 6.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 121 total aligned reads displayed.

GGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTT  >  NC_000913/560456‑560648
                                                                                                 |                                                                                               
ggCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGtcc                                                                                               <  2:1378859/100‑1 (MQ=255)
 gCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGtccc                                                                                              >  1:2586591/1‑100 (MQ=255)
 gCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGtccc                                                                                              <  1:2161744/100‑1 (MQ=255)
  cAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCa                                                                                             >  2:2273075/1‑100 (MQ=255)
  cAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCa                                                                                             >  1:386757/1‑100 (MQ=255)
  cAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCa                                                                                             <  1:2329666/100‑1 (MQ=255)
   aaTACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCaa                                                                                            >  1:914784/1‑100 (MQ=255)
   aaTACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCaa                                                                                            <  1:36196/100‑1 (MQ=255)
    aTACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCAAc                                                                                           <  2:1139929/100‑1 (MQ=255)
    aTACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGCCCCAAc                                                                                           <  1:2078038/100‑1 (MQ=255)
      aCAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACgg                                                                                         >  1:2378432/1‑100 (MQ=255)
       cAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGt                                                                                        >  2:404440/1‑100 (MQ=255)
        aGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTa                                                                                       <  2:53344/100‑1 (MQ=255)
           aaGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCTTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGt                                                                                    >  1:1576756/1‑100 (MQ=255)
             gTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTgg                                                                                  <  2:1619065/100‑1 (MQ=255)
              tGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTggg                                                                                 >  2:2798922/1‑100 (MQ=255)
               gATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGa                                                                                >  2:1981669/1‑100 (MQ=255)
               gATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGa                                                                                >  1:440841/1‑100 (MQ=255)
                  atTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTg                                                                             >  2:2835067/1‑100 (MQ=255)
                   ttttCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGc                                                                            <  2:2814531/100‑1 (MQ=255)
                    tttCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGcc                                                                           <  2:1531245/100‑1 (MQ=255)
                    tttCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGcc                                                                           <  1:2363232/100‑1 (MQ=255)
                    tttCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGcc                                                                           <  2:90378/100‑1 (MQ=255)
                     ttCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGcc                                                                           <  1:2387635/99‑1 (MQ=255)
                       cGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGt                                                                        >  2:2546044/1‑100 (MQ=255)
                        gATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCACAACGGTACGTGGGATTGCCCGTa                                                                       >  1:1255134/1‑100 (MQ=255)
                            gCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTAcggc                                                                   <  2:2503232/100‑1 (MQ=255)
                            gCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTAcggc                                                                   >  1:1088547/1‑100 (MQ=255)
                             cGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTAcggcg                                                                  <  2:2020217/100‑1 (MQ=255)
                             cGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTAcggcg                                                                  <  1:2688713/100‑1 (MQ=255)
                               ttGGATTTCTTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACggcggc                                                                <  2:2470231/100‑1 (MQ=255)
                               ttGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACggcggc                                                                <  1:1699810/100‑1 (MQ=255)
                                  gATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCa                                                             <  2:2296227/100‑1 (MQ=255)
                                   aTTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAg                                                            >  1:536540/1‑100 (MQ=255)
                                    tttCGTGGTGTACGGCTTTATTCTTCTGTTAATATGTATCCTGATAGCGATCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCa                                                             >  2:554819/1‑98 (MQ=255)
                                    tttCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTGCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGc                                                           <  2:1551232/100‑1 (MQ=255)
                                    tttCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGc                                                           <  2:2058088/100‑1 (MQ=255)
                                      tCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTa                                                         <  2:1658787/100‑1 (MQ=255)
                                      tCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTa                                                         >  1:1313473/1‑100 (MQ=255)
                                        gtggtgTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAAc                                                       >  1:1681764/1‑100 (MQ=255)
                                        gtggtgTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAAc                                                       <  1:2205134/100‑1 (MQ=255)
                                        gtggtgTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAAc                                                       <  1:24506/100‑1 (MQ=255)
                                        gtggtgTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAAc                                                       <  1:24561/100‑1 (MQ=255)
                                             gtACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTc                                                  >  2:2482539/1‑100 (MQ=255)
                                                cGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGc                                               >  2:2808897/1‑100 (MQ=255)
                                                 ggCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCa                                              >  1:2273321/1‑100 (MQ=255)
                                                  gCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCaa                                             >  1:1469525/1‑100 (MQ=255)
                                                  gCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCaa                                             >  1:1114047/1‑100 (MQ=255)
                                                   cTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCCGCTAACGATTCGGCaaa                                            >  2:1166134/1‑100 (MQ=255)
                                                          cttctGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGttt                                     <  2:94105/100‑1 (MQ=255)
                                                          cttctGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGttt                                     >  1:416559/1‑100 (MQ=255)
                                                          cttctGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGttt                                     <  1:1809310/100‑1 (MQ=255)
                                                          cttctGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGttt                                     <  1:369154/100‑1 (MQ=255)
                                                          cttctGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGtt                                      >  1:1555014/1‑99 (MQ=255)
                                                           ttctGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGtttt                                    >  1:2654820/1‑100 (MQ=255)
                                                             ctGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTttat                                  >  2:2678194/1‑100 (MQ=255)
                                                             ctGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTttat                                  <  2:914782/100‑1 (MQ=255)
                                                              tGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTttatt                                 <  1:1400535/100‑1 (MQ=255)
                                                              tGATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTttatt                                 >  1:188036/1‑100 (MQ=255)
                                                               gATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTttatta                                >  1:812179/1‑100 (MQ=255)
                                                               gATAATATGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTttatta                                >  1:1821048/1‑100 (MQ=255)
                                                                ataataTGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTttatta                                <  2:636818/99‑1 (MQ=255)
                                                                ataataTGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTtattat                               >  1:2010438/1‑100 (MQ=255)
                                                                 taataTGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATtatc                              <  2:79365/100‑1 (MQ=255)
                                                                 taataTGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATtatc                              <  1:1600944/100‑1 (MQ=255)
                                                                     atGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATtatctatc                          <  1:2798943/100‑1 (MQ=255)
                                                                      tGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCa                         >  1:11651/1‑100 (MQ=255)
                                                                      tGTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCa                         >  1:1198312/1‑100 (MQ=255)
                                                                       gTATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCaa                        <  1:385065/100‑1 (MQ=255)
                                                                        tATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTttattaagatcggaa                       >  2:2111815/1‑91 (MQ=255)
                                                                        tATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCaaa                       <  1:1981683/100‑1 (MQ=255)
                                                                        tATCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCaaa                       >  2:731797/1‑100 (MQ=255)
                                                                         aTCCTGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAg                      >  2:2335185/1‑100 (MQ=255)
                                                                         aTCCTGATAGCCAGCAAGGATTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAg                      >  2:58752/1‑100 (MQ=255)
                                                                          tCCTGATAGCCAGCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGc                     >  1:2010730/1‑100 (MQ=255)
                                                                             tGATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTAt                  >  1:121304/1‑100 (MQ=255)
                                                                              gATAGCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATg                 <  1:2188513/100‑1 (MQ=255)
                                                                                  gCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTc             >  1:467533/1‑100 (MQ=255)
                                                                                  gCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTc             <  2:1302062/100‑1 (MQ=255)
                                                                                  gCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTc             <  1:2622610/100‑1 (MQ=255)
                                                                                  gCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTc             >  2:444078/1‑100 (MQ=255)
                                                                                  gCCAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTc             >  1:2086070/1‑100 (MQ=255)
                                                                                   ccAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTcc            <  2:440840/100‑1 (MQ=255)
                                                                                   ccAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTcc            >  2:570275/1‑100 (MQ=255)
                                                                                   ccAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATNATCTATCAAAGCTATGTTTcc            >  2:1661386/1‑100 (MQ=255)
                                                                                    cAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTccc           >  1:1477570/1‑100 (MQ=255)
                                                                                    cAGCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTccc           <  2:714071/100‑1 (MQ=255)
                                                                                      gCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTccccc         >  2:2265760/1‑100 (MQ=255)
                                                                                      gCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTccccc         >  1:1556678/1‑100 (MQ=255)
                                                                                      gCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTccccc         <  1:1718447/100‑1 (MQ=255)
                                                                                      gCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTccccc         >  1:1860133/1‑100 (MQ=255)
                                                                                      gCAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAACGGTTTTATTATCTATCAAAGCTATGTTTccccc         <  2:1903472/100‑1 (MQ=255)
                                                                                       cAAGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCg        >  1:1576240/1‑100 (MQ=255)
                                                                                        aaGGGTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCgg       <  2:188035/100‑1 (MQ=255)
                                                                                         aGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGgc      >  1:809113/1‑100 (MQ=255)
                                                                                          gggTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCTCCCGgc      >  2:1624180/1‑99 (MQ=255)
                                                                                          gggTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGgcg     >  1:202057/1‑100 (MQ=255)
                                                                                           ggTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGgcgc    >  1:554371/1‑100 (MQ=255)
                                                                                            gTTTGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCt   <  1:1865431/100‑1 (MQ=255)
                                                                                             tttGTCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCtt  >  2:177766/1‑100 (MQ=255)
                                                                                                 |                                                                                               
GGCAATACAGGAAGTGATATTTTCGATGGCGTTGGATTTCGTGGTGTACGGCTTTATTCTTCTGATAATATGTATCCTGATAGCCAGCAAGGGTTTGCCCCAACGGTACGTGGGATTGCCCGTACGGCGGCCCAGCTAACGATTCGGCAAAATGGTTTTATTATCTATCAAAGCTATGTTTCCCCCGGCGCTT  >  NC_000913/560456‑560648

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