Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,191,432 C→T 100% G72S (GGC→AGC)  yehC ← putative periplasmic pilin chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,191,4320CT97.6% 314.2 / ‑3.2 103G72S (GGC→AGC) yehCputative periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (50/50);  total (52/51)
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.00e+00

Only 100 of 105 total aligned reads displayed.

TTTTTACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCCATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGCGCGT  >  NC_000913/2191336‑2191528
                                                                                                |                                                                                                
tttttACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTTGTTGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCtatc                                                                                               <  2:2855444/100‑1 (MQ=255)
tttttACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCtatc                                                                                               >  2:1588160/1‑100 (MQ=255)
   ttACTTGCTGCCCGGAATTTGCCCCTATTTTTTTTTTTTTTTTCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTatcatc                                                                                            >  2:2071355/1‑100 (MQ=255)
   ttACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTatcatc                                                                                            <  1:433570/100‑1 (MQ=255)
   gtACTTGCTGCCCGGAATTTACCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTatcatc                                                                                            >  2:574781/2‑100 (MQ=255)
      cTTGCTGCCCGGAATTTGCCCCTAATTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGATGATAATGACGTATCGCTATCATCAAt                                                                                         >  2:2863626/1‑100 (MQ=255)
           tgcCGGGAATTTGCCCCTATTTTTGCTGCTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATTACGTATAGCTATCATCAATCCAc                                                                                     <  2:549237/99‑1 (MQ=255)
           tgcCCGGAATTTGCCCGTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGATGTATCGCTATCATCAATCCACg                                                                                    <  2:389647/100‑1 (MQ=255)
            gcCCGGAATTTGCCCCTATTTTTTCTATTGTTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCCATCCACg                                                                                    >  2:1544377/1‑99 (MQ=255)
            gcCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGc                                                                                   >  1:1514442/1‑100 (MQ=255)
            gcCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGc                                                                                   >  1:391512/1‑100 (MQ=255)
               cGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTATGATGGTAATGACGTATCGCTATCATCAATCCACGCCTg                                                                                >  1:1304424/1‑100 (MQ=255)
                  aaTTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTgcag                                                                             <  2:862841/100‑1 (MQ=255)
                   aTTTGCCCGTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTgcagc                                                                            <  2:849732/100‑1 (MQ=255)
                   aTTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTgcagc                                                                            >  1:2851110/1‑100 (MQ=255)
                    tttGCCCCTATTTTTGCCACTGTTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTgcagc                                                                            <  2:340404/99‑1 (MQ=255)
                     ttGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGcagcag                                                                          <  2:877329/100‑1 (MQ=255)
                      tGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCagcagc                                                                         >  2:2017383/1‑100 (MQ=255)
                        ccccTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCga                                                                       >  1:3033680/1‑100 (MQ=255)
                        ccccTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCga                                                                       >  1:816431/1‑100 (MQ=255)
                        ccccTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCga                                                                       <  2:2859228/100‑1 (MQ=255)
                         cccTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCgaa                                                                      >  1:2432518/1‑100 (MQ=255)
                         cccTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCgaa                                                                      >  1:848828/1‑100 (MQ=255)
                          ccTATTTTTGCTATTTGTTGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCATGCATCAGCgaa                                                                      >  2:825889/1‑99 (MQ=255)
                              tttttGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACg                                                                 <  1:17083/100‑1 (MQ=255)
                                   gCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTc                                                            >  1:2288555/1‑100 (MQ=255)
                                   gCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTc                                                            >  1:214379/1‑100 (MQ=255)
                                    ncACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTcc                                                           >  2:126/2‑100 (MQ=255)
                                    ccACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTcc                                                           <  2:1446642/100‑1 (MQ=255)
                                    ccACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTcc                                                           >  1:451279/1‑100 (MQ=255)
                                    ccACTGGTGGAGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATTACGTATCGCTATCATCAATCCAGGCCTGCAGCAGCGAAAAGCGGctttc                                                           >  2:297667/1‑98 (MQ=255)
                                     cACTGGTGGCGTTAACATGAAAGGTACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTccc                                                          <  2:112789/100‑1 (MQ=255)
                                      aCTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTTGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCt                                                         <  1:72829/100‑1 (MQ=255)
                                      aCTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCt                                                         >  2:2601696/1‑100 (MQ=255)
                                        tggtggCGTTAACATGAAAGGAACCTGATTTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTgg                                                       <  1:1051569/100‑1 (MQ=255)
                                         ggtggCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGt                                                      <  1:1012530/100‑1 (MQ=255)
                                         ggtggCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGt                                                      >  2:2626798/1‑100 (MQ=255)
                                         ggtggCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGt                                                      <  1:914176/100‑1 (MQ=255)
                                          gttgtgTTAACATGATAGGAGCCTGATTTTTTTCTGGTGGTAATGACGTATGGCTATCATCAGTCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGtt                                                     <  2:2161376/95‑1 (MQ=255)
                                           tggCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTtc                                                    >  1:1390144/1‑100 (MQ=255)
                                            ggCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTtca                                                   >  2:3023091/1‑100 (MQ=255)
                                            ggCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTtca                                                   <  1:2194131/100‑1 (MQ=255)
                                            ggCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACAGTTTCCCTGGTtca                                                   >  1:1323755/1‑100 (MQ=255)
                                             gCGTTAACATGTAAGTAACCTGATTTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTtca                                                   <  1:1476657/99‑1 (MQ=255)
                                              cGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTtcatc                                                 >  1:362275/1‑100 (MQ=255)
                                              cGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTtcatc                                                 >  1:2186984/1‑100 (MQ=255)
                                              cGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTtcatc                                                 <  1:1031848/100‑1 (MQ=255)
                                               gTTAACCTATAATGAACCTAAATTTTTTCTGGTGGTAATTACGTATCGATATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTtcatca                                                <  1:2655071/100‑1 (MQ=255)
                                               gTTAACATGAAAGGAACCTGAATTTTTTCTTGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTtcatca                                                <  1:1380902/100‑1 (MQ=255)
                                               gTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATAAATACACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTtcatca                                                >  1:939610/1‑100 (MQ=255)
                                                ttAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCaa                                               <  1:894007/100‑1 (MQ=255)
                                                ttAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCaa                                               <  1:2510714/100‑1 (MQ=255)
                                                 tAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCATGGTTCATCaa                                               >  1:2113385/1‑99 (MQ=255)
                                                  aaCATGAAAGGAACCTGAATTTTTTTTTTTTGTAATAACGTATCGCTATCATCAATCTACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACt                                             >  2:2335846/1‑100 (MQ=255)
                                                    cATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGc                                           >  1:2252452/1‑100 (MQ=255)
                                                    cATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGc                                           >  2:2176016/1‑100 (MQ=255)
                                                    cATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGc                                           <  1:2608444/100‑1 (MQ=255)
                                                     aTGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGc                                           <  1:1034234/99‑1 (MQ=255)
                                                      tGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCTTTATCATCAATCCACTCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCAc                                         <  1:2852105/100‑1 (MQ=255)
                                                          aGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATc                                     >  2:1866435/1‑100 (MQ=255)
                                                           ggAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCa                                    >  1:11957/1‑100 (MQ=255)
                                                               ccTGAATTTTTTCTGGTTGTAATGACGTATCGTTATCATTTATCCAAGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCActtc                                <  1:263214/100‑1 (MQ=255)
                                                                cTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCActtct                               <  1:2392423/100‑1 (MQ=255)
                                                                 tGACTTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACTGTTTCCCTTGTTCATCAACTGCACCATCActtctt                              <  2:2312/100‑1 (MQ=255)
                                                                 tGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCCCGCCTACAGCAGCGAAGAACGATTTCCCTGGTTCATAAACTGCACCAtaaactctt                              <  2:1196386/100‑9 (MQ=255)
                                                                 tGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCActtctt                              <  2:303982/100‑1 (MQ=255)
                                                                 tGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCActtctt                              >  1:520247/1‑100 (MQ=255)
                                                                 tGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACCCCTGCAGCAGCGAAGAACGATTTCCCTGGTTCAACAACTGCACCATCActtc                                >  1:2696319/1‑98 (MQ=255)
                                                                  gAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCttt                             >  1:2574366/1‑100 (MQ=255)
                                                                    aTTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTAt                           >  1:90684/1‑100 (MQ=255)
                                                                    aTTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTAt                           >  1:995492/1‑100 (MQ=255)
                                                                       ttttCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTCCACCATCACTTCTTTAtttt                        <  2:1846378/100‑1 (MQ=255)
                                                                        tttCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTc                       <  1:315046/100‑1 (MQ=255)
                                                                        tttCTGGTGGTAATGACGTATCGCCATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTAGTTTc                       <  2:535163/100‑1 (MQ=255)
                                                                         ttCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCt                      <  1:1204991/100‑1 (MQ=255)
                                                                          tCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTg                     <  2:2679696/100‑1 (MQ=255)
                                                                           cTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGTAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTgn                    <  2:1791596/100‑2 (MQ=255)
                                                                           cTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGc                    <  2:1202840/100‑1 (MQ=255)
                                                                            tggtagTACTTACGTATCGCTATCATCAATCCACGTCTGCAGCAGTGAAGATCGGTTTCCCTGTTTCATCAACTGTGTCATCACTTCTGTATTTTCTg                     <  2:373782/98‑1 (MQ=255)
                                                                             ggtggtAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTACCTGGTTCATCAACTGCACCAGCACTTCTTTATTTTCTGcc                   >  2:2831192/1‑99 (MQ=255)
                                                                              gtggtAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCgg                 <  1:1375961/100‑1 (MQ=255)
                                                                               tggtAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCggg                <  1:2034670/100‑1 (MQ=255)
                                                                                ggtAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGt               >  1:428883/1‑100 (MQ=255)
                                                                                   aaTTACGTATCGCTATCATCACTCCACGCCTCCAGCAGCGAAGAACGGTTTCCCTGGTTCATTAACTGCACCATCTCTTCTTTTTTTTCTGCCGGGTaaa            <  2:2043961/100‑1 (MQ=255)
                                                                                   aaTGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGCTCATCAACTGCACCATCACTTCTTTATTTTCGGCCTGGTaa             >  2:1505101/1‑99 (MQ=255)
                                                                                    aTGACGTATCGCCATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATAAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAaat           >  2:1331754/1‑100 (MQ=255)
                                                                                      gACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAaataa         <  2:2190607/100‑1 (MQ=255)
                                                                                      gACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAaataa         >  2:1566941/1‑100 (MQ=255)
                                                                                      gACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAaataa         >  2:2573551/1‑100 (MQ=255)
                                                                                      gACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAaataa         >  1:2023062/1‑100 (MQ=255)
                                                                                        cGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGATTCCCTGGCTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATg       >  2:2333482/1‑100 (MQ=255)
                                                                                          tATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATgcg     <  1:2685463/100‑1 (MQ=255)
                                                                                          tATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGAGTAAATAATgc      <  2:399136/99‑1 (MQ=255)
                                                                                            tCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGcgcg   >  2:1112112/1‑100 (MQ=255)
                                                                                            tCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGcgcg   >  2:1835858/1‑100 (MQ=255)
                                                                                            tCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGcgcg   <  1:1343596/100‑1 (MQ=255)
                                                                                            tCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGcgcg   <  2:3033646/100‑1 (MQ=255)
                                                                                            tCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGcgcg   <  1:490988/100‑1 (MQ=255)
                                                                                            tCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGcgcg   >  1:2708298/1‑100 (MQ=255)
                                                                                             cGCCATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGCGCGt  >  1:1333884/1‑100 (MQ=255)
                                                                                                |                                                                                                
TTTTTACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCCATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGCGCGT  >  NC_000913/2191336‑2191528

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