Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,167,075 T→C 19.1% D215G (GAT→GGT)  ygiS ← putative ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,167,0750TC19.1% 199.5 / 39.3 99D215G (GAT→GGT) ygiSputative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (37/43);  new base C (10/9);  total (47/52)
Fisher's exact test for biased strand distribution p-value = 7.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTAGCGCCAGATACTCAACCTGTTGCAATACTGTATGTTGCGCATCGCGGTACTTTGGATTTTTGCGTGCAGTAATCTTTTCGTTAACTACCCACTGATCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCTTACTCCAGCTATCGCCATGTTTAGCGATGACATGATGAGGAACCGGGAACAGCGTCGGC  >  NC_000913/3166976‑3167173
                                                                                                   |                                                                                                  
tCTAGCGCCAGATACTCAACCTGTTGCAATACTGTATGTTGCGCATCGCGGTACTTTGGATTTTTGCGTGCAGTAATCTTTTCGTTAACTACCCACTGAt                                                                                                    >  2:758989/1‑100 (MQ=255)
 cTAGCGCCAGATACTCAACCTGTTGCAATACTGTATGTTGCGCATCGCGGTACTTTGGATTTTTGCGTGCAGTAATCTTTTCGTTAACTACCCACTGATc                                                                                                   >  1:69992/1‑100 (MQ=255)
 cTAGCGCCAGATACTCAACCTGTTGCAATACTGTATGTTGCGCATCGCGGTACTTTGGATTTTTGCGTGCAGTAATCTTTTCGTTAACTACCCACTGATc                                                                                                   >  2:2202788/1‑100 (MQ=255)
 cTAGCGCCAGATACTCAACCTGTTGCAATACTGTATGTTGCGCATCGCGGTACTTTGGATTTTTGCGTGCAGTAATCTTTTCGTTAACTACCCACTGATc                                                                                                   <  2:2563481/100‑1 (MQ=255)
  tAGCGCCAGATACTCAACCTGTTGCAATACTGTATGTTGCGCATCGCGGTACTTTGGATTTTTGCGTGCAGTAATCTTTTCGTTAACTACCCACTGATCa                                                                                                  <  2:2402598/100‑1 (MQ=255)
  tAGCGCCAGATACTCAACCTGTTGCAATACTGTATGTTGCGCATCGCGGTACTTTGGATTTTTGCGTGCAGTAATCTTTTCGTTAACTACCCACTGATCa                                                                                                  >  2:2542719/1‑100 (MQ=255)
     cgcCAGATACTCAACCTGTTGCAATACTGTATGTTGCGCATCGCGGTACTTTGGATTTTTGCGTGCAGTAATCTTTTCGTTAACTACCCACTGATCAAGc                                                                                               >  1:685098/1‑100 (MQ=255)
       ccAGATACTCAACCTGTTGCAATACTGTATGTTGCGCATCGCGGTACTTTGGATTTTTGCGTGCAGTAATCTTTTCGTTAACTCCCCACTGACCAAGcac                                                                                             >  2:2258720/1‑100 (MQ=255)
       ccAGATACTCAACCTGTTGCAATACTGTATGTTGCGCATCGCGGTACTTTGGATTTTTGCGTGCAGTAATCTTTTCGTTAACTACCCACTGATCAAGcac                                                                                             <  2:2086674/100‑1 (MQ=255)
          gATACTCAACCTGTTGCAATACTGTATGTTGCGCATCGCGGTACTTTGGATTTTTGCGTGCAGTAATCTTTTCGTTAACTACCCACTGATCAAGCACaaa                                                                                          >  2:501042/1‑100 (MQ=255)
            tACTCAACCTGTTGCAATACTGTATGTTGCGCATCGCGGTACTTTGGATTTTTGCGTGCAGTAATCTTTTCGTTAACTACCCACTGATCAAGCACAAAgg                                                                                        >  1:410326/1‑100 (MQ=255)
            tACTCAACCTGTTGCAATACTGTATGTTGCGCATCGCGGTACTTTGGATTTTTGCGTGCAGTAATCTTTTCGTTAACTACCCACTGATCAAGCACAAAgg                                                                                        >  1:623092/1‑100 (MQ=255)
            tACTCAACCTGTTGCAATACTGTATGTTGCGCATCGCGGTACTTTGGATTTTTGCGTGCAGTAATCTTTTCGTTAACTACCCACTGACCAAGCACAAAgg                                                                                        <  1:692265/100‑1 (MQ=255)
             aCTCGACCTGTTGCAATACTGTATGTTGCGCATCGCGGTACTTTGGATCTTTGCGTGCAGTAATCTTTTCGTTAACTACCCACTGATCAAGCACAAAGGc                                                                                       <  1:1269934/100‑1 (MQ=255)
              cTCAACCTGTTGCAATACTGTATGTTGCGCATCGCGGTACTTTGGATTTTTGCGTGCAGTAATCTTTTCGTTAACTACCCACTGATCAAGCACAAAGGCa                                                                                      >  2:2374206/1‑100 (MQ=255)
              cTCAACCTGTTGCAATACTGTATGTTGCGCATCGCGGTACTTTGGATATTTGCGTGCAGTAATCTTTTCATTAACTACCCACTGCTCACGCACCAAgg                                                                                        <  1:1025099/98‑1 (MQ=255)
                cAACCTGTTGCAATACTGTATGTTGCGCATCGCGGTACTTTGGATTTTTGCGTGCAGTAATCTTTTCGTTAACTACCCACTGATCAAGCACAAAGGCAcc                                                                                    >  1:558751/1‑100 (MQ=255)
                cAACCTGTTGCAATACTGTATGTTGCGCATCGCGGTACTTTGGATTTTTGCGTGCAGTAATCTTTTCGTTAACTACCCACTGACCAAGCACAAAGGCAcc                                                                                    >  1:1291330/1‑100 (MQ=255)
                  aCCTGTTGCAATACTGTATGTTGCGCATCGCGGTACTTTGGATTTTTGCGTGCAGTAATCTTTTCGTTAACTACCCACTGATCAAGCACAAAGGCACCGt                                                                                  <  1:758986/100‑1 (MQ=255)
                  aCCTGTTGCAATACTGTATGTTGCGCATCGCGGTACCTTGGATTTTTGCGTGCAGTAATCTTTTCGTTAACTACCCCCTGACCAAGCACAAAGGCACCGt                                                                                  >  1:2472403/1‑100 (MQ=255)
                   ccTGTTGCAATACTGTATGTTGCGCATCGCGGTACTTTGGATTTTTGCGTGCAGTAATCTTTTCGTTAACTACCCACTGATCAAGCACAAAGGCACCGt                                                                                  <  1:220823/99‑1 (MQ=255)
                     tGTTGCAATACTGTATGTTGCGCATCGCGGTACTTTGGATTTTTGCGTGCAGTAATCTTTTCGTTAACTACCCACTGACCAAGCACAAAGGCACCGTTGt                                                                               <  2:905317/100‑1 (MQ=255)
                     tGTTGCAATACTGTATGTTGCGCATCGCGGTACTTTGGATTTTTGCGTGCAGTAATCTTTTCGTTAACTACCCACTGACCAAGCACAAAGGCACCGTTGt                                                                               <  1:2501124/100‑1 (MQ=255)
                      gTTGCAATACTGTATGGTGCGCATCGCGGTACTTGGGATTTTTGCGTGCAGTAATCTTTGCGTTAACTACCCACGGATCAAGCGCAAAGGCACCGTTGTa                                                                              <  1:502377/100‑1 (MQ=255)
                         gCAATACTGTATGTTGCGCATCGCGGTACTTTGGATTTTTGCGTGCAGTAATCTTTTCGTTAACTACCCACTGATCAAGCACAAAGGCACCGTTGTAAAc                                                                           >  2:652403/1‑100 (MQ=255)
                         gCAATACTGTATGTTGCGCATCGCGGTACTTTGGATTTTTGCGTGCAGTAATCTTTTCGTTAACTACCCACTGACCAAGCACAAAGGCACCGTTGTAAAc                                                                           >  2:1882161/1‑100 (MQ=255)
                          cAATACTGTATGTTGCGCATCGCGGTACTTTGGATTTTTGCGTGCAGTAATGTTTTCGTTAACTACCCACTGACCAAGCACAAAGGCACCGTTGTAAAcc                                                                          >  1:1445788/1‑100 (MQ=255)
                          cAATACTGTATGTTGCGCATCGCGGTACTTTGGATTTTTGCGTGCAGTAATCTTTTCGTTAACTACCCACTGATCAAGCACAAAGGCACCGTTGTAAAcc                                                                          >  1:106233/1‑100 (MQ=255)
                          cAATACTGTATGTTGCGCATCGCGGTACTTTGGATTTTTGCGTGCAGTAATCTTTTCGTTAACTACCCACTGATCAAGCACAAAGGCACCGTTGTAAAcc                                                                          <  2:619704/100‑1 (MQ=255)
                          cAATACTGTATGTTGCGCATCGCGGTACTTTGGATTTTTGCGTGCAGTAATCTTTTCGTTAACTACCCACTGATCAAGCACAAAGGCACCGTTGTAAAcc                                                                          >  1:553098/1‑100 (MQ=255)
                          cAATACTGTATGTTGCGCATCGCGGTACTTTGGATTTTTGCGTGCAGTAATCTTTTCGTTAACTACCCACTGATCAAGCACAAAGGCACCGTTGTAAAcc                                                                          >  1:771919/1‑100 (MQ=255)
                          cAATACTGTATGTTGCGCATCGCGGTACTTTGGATTTTTGCGTGCAGTAATCTTTTCGTTAACTACCCACTGACCAAGCACAAAGGCACCGTTGTAAAcc                                                                          >  2:980676/1‑100 (MQ=255)
                           aaTACTGTATGTTGCGCATCGCGGTACTTTGGATTTTTGCGTGCAGTAATCTTTTCGTTAACTACCCACTGATCAAGCACAAAGGCACCGTTGTAAAcc                                                                          <  2:679318/99‑1 (MQ=255)
                                tGTATGTTGCGCATCGCGGTACTTTGGATTTTTGCGTGCAGTAATCTTTTCGTTAACTACCCACTGATCAAGCACAAAGGCACCGTTGTAAACCATGTtc                                                                    <  1:2537668/100‑1 (MQ=255)
                                 gTATGTTGCGCATCGCGGTACTTTGGATTTTTGCGTGCAGTAATCTTTTCGTTAACTACCCACTGATCAAGCACAAAGGCACCGTTGTAAACCATGTtct                                                                   >  1:2220012/1‑100 (MQ=255)
                                  tATGTTGCGCATCGCGGTACTTTGGATTTTTGCGTGCAGTAATCTTTTCGTTAACTACCCACTGATCAAGCACAAAGGCACCGTTGTAAACCATGTtctc                                                                  >  1:26912/1‑100 (MQ=255)
                                     gTTGCGCATCGCGGTACTTTGGATTTTTGCGTGCAGTAATCTTTTCGTTAACTACCCACTGATCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTgg                                                               <  2:2115185/100‑1 (MQ=255)
                                     gTTGCGCATCGCGGTACTTTGGATTTTTGCGTGCAGTAATCTTTTCGTTAACTACCCACTGATCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTgg                                                               >  2:1237631/1‑100 (MQ=255)
                                     gTTGCGCATCGCGGTACTTTGGATTTTTGCGTGCAGTAATCTTTTCGTTAACTACCCACTGATCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTgg                                                               <  2:1318462/100‑1 (MQ=255)
                                       tGCGCATCGCGGTACTTTGGATTTTTGCGTGCAGTAATCTTTTCGTTAACTACCCACTGGTCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCt                                                             <  1:208010/100‑1 (MQ=255)
                                       tGCGCATCGCGGTACTTTGGATTTTTGCGTGCAGTAATCTTTTCGTTAACTACCCACTGATCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCt                                                             <  2:1060002/100‑1 (MQ=255)
                                       tGCGCATCGCGGTACTTTGGATTTTTGCGTGCAGTAATCTTTTCGTTAACTACCCACTGATCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCt                                                             <  2:349154/100‑1 (MQ=255)
                                       tGCGCATCGCGGTACTTTGGATTTTTGCGTGCAGTAATCTTTTCGTTAACTACCCACTGATCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCt                                                             <  1:623095/100‑1 (MQ=255)
                                       tGCGCATCGCGGTACTTTGGATTTTTGCGTGCAGTAATCTTTTCGTTAACTACCCACTGATCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCt                                                             <  2:1544529/100‑1 (MQ=255)
                                       tGCGCATCGCGGTACTTTGGATTTTTGCGTGCAGTAAGCTTTTCGTTAACTACCCACTGACCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCt                                                             >  1:1386998/1‑100 (MQ=255)
                                        gcgcATCGCGGTACTTTGGATTTTTGCGTGCAGTAATCTTTTCGTTAACTACCCACTGATCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCtt                                                            <  1:400894/100‑1 (MQ=255)
                                           cATCGCGGTACTTTGGATTTTTGCGTGCAGTAATCTTTTCGTTAACTACCCACTGATCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCTTACt                                                         <  1:73704/100‑1 (MQ=255)
                                            aTCGCGGTACTTTGGATTTTTGCGTGCAGTAATCTTTTCGTTAACTACCCACTGATCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCTTACTc                                                        >  1:2456551/1‑100 (MQ=255)
                                             tCGCGGTACTTTGGATTTTTGCGTGCAGTAATCTTTTCGTTAACTACCCACTGATCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCTTACTcc                                                       <  1:378623/100‑1 (MQ=255)
                                               gcgGTACTTTGGATTTTTGCGTGCAGTAATCTTTTCGTTAACTACCCACTGATCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCTTACTCCAg                                                     >  2:814098/1‑100 (MQ=255)
                                               gcgGTACTTTGGATTTTTGCGTGCAGTAATCTTTTCGTTAACTACCCACTGATCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCTTACTCCAg                                                     >  2:686348/1‑100 (MQ=255)
                                                      tttGGCTTTTTGCGTGCAGTAATCTTTTCGTTAACTACCCACTGATCAAGCACAAAGGCACCGTTATAAACCATGTTCTCTGGCTTACTCCAGCTATCGc                                              >  1:1294032/1‑100 (MQ=255)
                                                       ttGGATTTTTGCGTGCAGTAATCTTTTCGTTAACTACCCACTGATCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCTTACTCCAGCTATCGcc                                             <  2:1621206/100‑1 (MQ=255)
                                                       ttGGATTTTTGCGTGCAGTAATCTTTTCGTTAACTACCCACTGATCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCTTACTCCAGCTATCGcc                                             >  1:1725871/1‑100 (MQ=255)
                                                        tGGATTTTTGCGTGCAGTAATCTTTTCGTTAATTACCCACTGATCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCTTACTCCAGCTATCGCCa                                            <  2:464782/100‑1 (MQ=255)
                                                          gATTTTTGCGTGCAGTAATCTTTTCGTTAACTACCCACTGATCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCTTACTCCAGCTATCGCCATg                                          >  2:676693/1‑100 (MQ=255)
                                                            tttttGCGTGCAGTAATCTTTTCGTTAACTACCCACTGATCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCTTACTCCAGCTATCGCCATGtt                                        >  2:1344081/1‑100 (MQ=255)
                                                               ttGCGTGCAGTAATCTTTTCGTTAACTACCCACTGATCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCTTACTCCAGCTATCGCCATGTTTAg                                     <  2:647780/100‑1 (MQ=255)
                                                                tGCGTGCAGTAATCTTTTCGTTAACTACCCACTGCTCAAGCTCAAAGGCACCGTTGTAAACCATGTTCTCTGGCTTACTCCAGCTATCGCCATGTTTAGc                                    <  2:567630/100‑1 (MQ=255)
                                                                tGCGTGCAGTAATCTTTTCGTTAACTACCCACTGATCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCTTACTCCAGCTATCGCCATGTTTAGc                                    <  2:1256790/100‑1 (MQ=255)
                                                                    tGCAGTAATCTTTTCGTTAACTACCCACTGACCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCTTACTCCAGCTATCGCCATGTTTAGCGATg                                <  2:185760/100‑1 (MQ=255)
                                                                      cAGTAATCTTTTCGTTAACTACCCACTGATCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCTTACTCCAGCTATCGCCATGTTTAGCGATGAc                              <  1:721669/100‑1 (MQ=255)
                                                                      cAGTAATCTTTTCGTTAACTACCCACTGATCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCGTAATTCAGCTATCGCCATGTTTAGCGATGAc                              >  1:985725/1‑100 (MQ=255)
                                                                      cAGTAATCTCTTCGTTAACTACCCACTGATCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCTTACTCCAGCTATCGCCATGTTTAGCGATGAc                              <  1:852259/100‑1 (MQ=255)
                                                                       aGTAATCTTTTCGTTAACTACCCACTGGCCAAGCACAAGTGCACCGTTGTAAACCATGTGCTCTGGCTTAGTCCAGCTATCGCCATGTTTAGCGATGAc                              <  1:141215/99‑1 (MQ=255)
                                                                        gTAATCTTTTCGTTAACTACCCACTGATCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCTTACTCCAGCTATCGCCATGTTTAGCGATGACat                            <  2:1621522/100‑1 (MQ=255)
                                                                          aaTCTTTTCGTTAACTACCCACTGATCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCTTACTCCAGCTATCGCCATGTTTAGCGATGACatga                          <  1:1921485/100‑1 (MQ=255)
                                                                          aaTCTTTTCGTTAACTACCCACTGACCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCTTACTCCAGCTATCGCCATGTTTAGCGATGACatga                          <  1:852033/100‑1 (MQ=255)
                                                                            tCTTTTCGTTAACTACCCACTGATCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCTTACTCCAGCTATCGCCATGTTTAGCGATGACatgatg                        <  2:1874382/100‑1 (MQ=255)
                                                                            tCTTTTCGTTAACTACCCACTGATCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCTTACTCCAGCTATCGCCATGTTTAGCGATGACatgatg                        <  1:52770/100‑1 (MQ=255)
                                                                            tCTTTTCGTTAACTACCCACTGACCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCTTACTCCAGCTATCGCCATGTTTAGCGATGACatgatg                        >  2:704180/1‑100 (MQ=255)
                                                                             cTTTTCGTTAACTACCCACTGATCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCTTACTCCAGCTATCGCCATGTTTAGCGATGACAtgatga                       <  1:1796559/100‑1 (MQ=255)
                                                                             cTTTTCGTTAACTACCCACTGATCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCTTACTCCAGCTATCGCCATGTTTAGCGATGACAtgatga                       <  2:819554/100‑1 (MQ=255)
                                                                               tttCGTTAACTACCCACTGATCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCTTACTCCAGCTATCGCCATGTTTAGCGATGACATGATGAgg                     <  2:1230257/100‑1 (MQ=255)
                                                                               tttCGTTAACTACCCACTGACCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCTTACTCCAGCTATCGCCATGTTTAGCGATGACATGATGAgg                     <  1:1297117/100‑1 (MQ=255)
                                                                                ttCGTTAACTAACCACTGATCAAGCACAAAGGAACCGTTGTAAACCATGTTCTCTGGCTTAATCCAGCTATCGCCATGTTTAGCGATGACATGATGAGGa                    <  1:926622/100‑1 (MQ=255)
                                                                                 tCGTTAACTCCCCACTGATCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCTTACTCCAGCTATCGCCATGTTTAGCGATGACATGATGAGGa                    >  1:712551/1‑99 (MQ=255)
                                                                                 tCGTTAACTACCCACTGACCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCTTACTCCAGCTATCGCCATGTTTAGCGATGACATGATGAGGaa                   >  2:2577775/1‑100 (MQ=255)
                                                                                  cGTTAACTACCCACTGATCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCTTACTCCAGCTATCGCCATGTTTAGCGATGACATGATGAGGAAc                  >  1:2056377/1‑100 (MQ=255)
                                                                                   gTTAACTACCCACTGATCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCTTACTCCAGCTATCGCCATGTTTAGCGATGACATGATGAGGAAcc                 >  2:2348611/1‑100 (MQ=255)
                                                                                     tAACTCCCCACTGATCAAGCACCAAGGCACCGTTGTAAACCATGTTCTCTGGCTTACTCCAGCTGTGGCCATGTTTAGCGATGACATGATGAGGAAGCgg               <  1:1969683/100‑1 (MQ=255)
                                                                                      aaCTACCCACTGATCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCTTACTCCAGCTATCGCCATGTTTAGCGATGACATGATGAGGAACCggg              >  2:1561540/1‑100 (MQ=255)
                                                                                      aaCTACCCACTGATCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCTTACTCCAGCTATCGCCATGTTTAGCGATGACATGATGAGGAACCggg              >  1:2537796/1‑100 (MQ=255)
                                                                                       aCTACCCACTGATCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCTTACTCCAGCTATCGCCATGTTTAGCGATGACATGATGAGGAACCGGGa             <  2:2537815/100‑1 (MQ=255)
                                                                                       aCTACCCACTGATCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCTTACTCCAGCTATCGCCATGTTTAGCGATGACATGATGAGGAACCGGGa             <  2:1064532/100‑1 (MQ=255)
                                                                                       aCTACCCACTGATCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCTTACTCCAGCTATCGCCATGTTTAGCGATGACATGATGAGGAACCGGGa             >  2:1044422/1‑100 (MQ=255)
                                                                                       aCTACCCACTGACCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCTTACTCCAGCTATCGCCATGTTTAGCGATGACATGATGAGGAACCGGGa             <  1:1201805/100‑1 (MQ=255)
                                                                                         tACCCACTGATCAAGCACAAAGCCACCGGTGTAAACCATGTTCTCTGGCTTACTACAGCCATCGCCATGTTTAGCGATGACATGATGAGGAACCGGGAAc           <  2:673912/100‑1 (MQ=255)
                                                                                         tACCCACTGACCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCTTACTCCAGCTAGCGCCATGTTTAGCGATGACATGATGAGGAACCGGGAAc           <  1:36245/100‑1 (MQ=255)
                                                                                         atctcACTGATCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCTTACTCCAGCTATCGCCATGTTTAGCGATGACATGATGAGGAACCGGGAAc           <  2:1047426/96‑1 (MQ=255)
                                                                                          aCCCACTGATCAAGCACAACGGCACCGTTGTAAACCATGTTCTCTGCCTTACTCCAGCTATCGCCATATTTAGCGATGACATGCTGAGGAACCGGGAAc           >  1:755628/1‑99 (MQ=255)
                                                                                          aCCCACTGATCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCTTACTCCAGCTATCGCCATGTTTAGCGATGACATGATGAGGAACCGGGAACa          >  1:2519997/1‑100 (MQ=255)
                                                                                           cccACTGATCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCTTACTCCAGCTATCGCCATGTTTAGCGATGACATGATGAGGAACCGGGAACAg         <  1:2542700/100‑1 (MQ=255)
                                                                                             cACTGATCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCTTACTCCAGCTATCGCCATGTTTAGCGATGACATGATGAGGAACCGGGAACAgat       >  1:1047423/1‑98 (MQ=255)
                                                                                             cACTGACCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCTTACTCCAGCTATCGCCATGTTTAGCGATGACATGATGAGGAACCGGGAACAGCg       >  2:2360730/1‑100 (MQ=255)
                                                                                                tGATCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCTTACTCCAGCTATCGCCATGTTTAGCGATGACATGATGAGGAACCGGGAACAGCGTCg    <  2:2056390/100‑1 (MQ=255)
                                                                                                 gATCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCTTACTCCAGCTATCGCCATGTTTAGCGATGACATGATGAGGAACCGGGAACATCGTCgg   >  1:968114/1‑100 (MQ=255)
                                                                                                 gATCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCTTACTCCAGCTATCGCCATGTTTAGCGATGACATGATGAGGAACCGGGAACAGCGTCgg   >  1:1324803/1‑100 (MQ=255)
                                                                                                  aTCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCTTACTCCAGCTATCGCCATGTTTAGCGATGACATGATGAGGAACCGGGAACAGCGTCGGc  >  2:1700651/1‑100 (MQ=255)
                                                                                                  aTCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCTTACTCCAGCTATCGCCATGTTTAGCGATGACATGATGAGGAACCGGGAACAGCGTCGGc  >  1:1242941/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
TCTAGCGCCAGATACTCAACCTGTTGCAATACTGTATGTTGCGCATCGCGGTACTTTGGATTTTTGCGTGCAGTAATCTTTTCGTTAACTACCCACTGATCAAGCACAAAGGCACCGTTGTAAACCATGTTCTCTGGCTTACTCCAGCTATCGCCATGTTTAGCGATGACATGATGAGGAACCGGGAACAGCGTCGGC  >  NC_000913/3166976‑3167173

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