Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,781,132 T→C D532G (GAC→GGC)  ypjA ← adhesin‑like autotransporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,781,1320TC98.7% 274.1 / NA 79D532G (GAC→GGC) ypjAadhesin‑like autotransporter
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (41/37);  total (41/38)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

TACATTTTTAGCTTCTGAATTCCTGATGTAGAAAGCATCACCATAACTGTTGGTTCCTTCGATAAGTGTTCCGGAGGTCGTGGAAGCAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGTCCTGAACCTCCAGAACGCCGCCATTATTAAGTCTGGTGCTATCTGTTTTAGCCTCCTTCTGGACAATCAGCTTACCGCCAGTATCAATGGTAGTATTTTTCGCCGAGGTTTTAGCCACTACCGTCAGTTCGC  >  NC_000913/2780996‑2781264
                                                                                                                                        |                                                                                                                                    
tACATTTTTAGCTTCTGAATTCCTGATGTAGAAAGCATCACCATAACTGTTGGTTCCTTCGATAAGTGTTCCGGAGGTCGTGGAAGCAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGccc                                                                                                                                    >  2:1598600/1‑139 (MQ=255)
       ttAGCTTCTGAATTCCTGATGTAGAAAGCATCACCATAACTGTTGGTTCCTTCGATAAGTGTTCCGGAGGTCGTGGAAGCAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGAACCt                                                                                                                             >  1:156498/1‑139 (MQ=255)
          gCTTCTGAATTCCTGATGTAGAAAGCATCACCATAACTGTTGGTTCCTTCGATAAGTGTTCCGGAGGTCGTGGAAGCAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGAACCTCCa                                                                                                                          >  2:809119/1‑139 (MQ=255)
          gCTTCTGAATTCCTGATGTAGAAAGCATCACCATAACTGTTGGTTCCTTCGATAAGTGTTCCGGAGGTCGTGAAAGCAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACAGCCCTGAACCTCCa                                                                                                                          >  1:616805/1‑139 (MQ=255)
                      cTGATGTAGAAAGCATCACCATAACTGTTGGTTCCTTCGATAAGTGTTCCGGAGGTCGTGGAAGCAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAACGCCGCCat                                                                                                              >  1:938222/1‑139 (MQ=255)
                        gATGTAGAAAGCATCACCATAACTGTTGGTTCCTTCGATAAGTGTTCCGGAGGTCGTGGAAGCAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAACGCCGCCatta                                                                                                            >  1:1324687/1‑139 (MQ=255)
                                 aGCATCACCATAACTGTTGGTTCCTTCGATAAGTGTTCCGGAGGTCGTGGAAGCAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAACGCCGCCattatt                                                                                                          <  1:1088744/132‑1 (MQ=255)
                                 aGCATCACCATAACTGTTGGTTCCTTCGATAAGTGTTCCGGAGGTCGTGGAAGCAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAACGCCGCCATTATTAAGTCTg                                                                                                   <  1:323259/139‑1 (MQ=255)
                                 aGCATCACCATAACTGTTGGTTCCTTCGATAAGTGTTCCGGAGGTCGTGGAACCAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAACGCCGCCattatt                                                                                                          >  2:1088744/1‑132 (MQ=255)
                                    aTCACCATAACTGTTGGTTCCTTCGATAAGTGTTCCGGAGGTCGTGGAAGCAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAACGCCGCCATTATTAAGTCTGGTg                                                                                                >  2:1082070/1‑139 (MQ=255)
                                    aTCACCATAACTGTTGGTTCCTTCGATAAGTGTTCCGGAGGTCGTGGAAGCAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAACGCCGCCATTATTAAGTCTGGTg                                                                                                <  1:1413764/139‑1 (MQ=255)
                                      cACCATAACTGTTGGTTCCTTCGATAAGTGTTCCGGAGGTCGTGGAAGCAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAAcgccgc                                                                                                                 <  2:898028/120‑1 (MQ=255)
                                      cACCATAACTGTTGGTTCCTTCGATAAGTGTTCCGGAGGTCGTGGAAGCAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAAcgccgc                                                                                                                 >  1:898028/1‑120 (MQ=255)
                                        ccATAACTGTTGGTTCCTTCGATAAGTGTTCCGGAGGTCGTGGAAGCAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAACGCCGCCATTATTAAGTCTGGTGCTAt                                                                                            <  1:281648/139‑1 (MQ=255)
                                             aCTGTTGGTTCCTTCGATAAGTGTTCCGGAGGTCGTGGAAGCAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAACGCCGCCa                                                                                                               >  1:1516644/1‑115 (MQ=255)
                                             aCTGTTGGTTCCTTCGATAAGTGTTCCGGAGGTCGTGGAAGCAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAACGCCGCCa                                                                                                               <  2:1516644/115‑1 (MQ=255)
                                                gTTGGTTCCTTCGATAAGTGTTCCGGAGGTCGTGGAAGCAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTACCCTCACCACCGCCCTGAACCTCCAGAACGCCGCCATTATTAAGTCTGGTGCTATCTGTTTTa                                                                                    >  2:206518/1‑139 (MQ=255)
                                                 ttGGTTCCTTCGATAAGTGTTCCGGAGGTCGTGGAAGCAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAACgccgcc                                                                                                                >  1:996162/1‑110 (MQ=255)
                                                 ttGGTTCCTTCGATAAGTGTTCCGGAGGTCGTGGAAGCAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAACgccgcc                                                                                                                <  2:996162/110‑1 (MQ=255)
                                                 ttGGTTCCTTCGATAAGGTTCCGGAGGTCGTGGAAGCAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAACGCCGCCATTATTAAGt                                                                                                      <  2:1086873/119‑1 (MQ=255)
                                                 ttGGTTCCTTCGATAAGGTTCCGGAGGTCGTGGAAGCAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAACGCCGCCATTATTAAGt                                                                                                      >  1:1086873/1‑119 (MQ=255)
                                                   ggTTCCTTCGATAAGTGTTCCGGAGGTCGTGGAAGCAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAACGCCGCCATTATTAAGTCTGGTGCTATCTGTTTTAgtc                                                                                 <  1:809119/139‑3 (MQ=255)
                                                   ggTTCCTTCGATAAGTGTTCCGGAGGTCGTGGAAGCAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAACGCCGCCATTATTAAGTCTGGTGCTATCTGTTTTAGcc                                                                                 <  1:94007/139‑1 (MQ=255)
                                                   ggTTCCTTCGATAAGTGTTCCGGAGGTCGTGGAAGCAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAACGCCGCCATTATTAAGTCTGGTGCTATCTGTTTTAGcc                                                                                 <  2:156498/139‑1 (MQ=255)
                                                        cTTCGATAAGTGTTCCGGAGGTCGTGGAAGCAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAACGCCGCCATTATTAAGTCTGGTGCTATCTGTTTTAGCCTCCtt                                                                            <  2:1324687/139‑1 (MQ=255)
                                                         ttCGATAAGTGTTCCGGAGGTCGTGGAAGCAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAACGCCGCCATTATTAAGTCTGGTGCTATCTGTTTTAGCCTCCTTc                                                                           >  1:1573266/1‑139 (MQ=255)
                                                         ttCGATAAGTGTTCCGGAGGTCGTGGAAGCAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAACGCCGCCATTATTAAGTCTGGTGCTATCTGTTTTAGCCTCCTTc                                                                           >  2:864116/1‑139 (MQ=255)
                                                         ttCGATAAGTATTCCGGAGGTCGTGGAAGCAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAACGCCGCCATTATTAAGTCTGGTGCTATCTGTTTTAGCCTCCTTc                                                                           >  1:140500/1‑139 (MQ=255)
                                                          tCGATAAGTGTTCCGGAGGTCGTGGAAGCAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTg                                                                                                                                  <  2:884419/83‑1 (MQ=255)
                                                          tCGATAAGTGTTCCGGAGGTCGTGGAAGCAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTg                                                                                                                                  >  1:884419/1‑83 (MQ=255)
                                                            gATAAGTGTTCCGGAGGTCGTGGAAGCAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAACGCCGCCATTATTAAGTCTGGTGCTATCTGTTTTAGCCTCCTTCTgg                                                                        <  2:938222/139‑1 (MQ=255)
                                                              tAAGTGTTCCGGAGGTCGTGGAAGCAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAACGCCGCCATTATTAAGTCTGGTGCTATCTGTTTTAGCCTCCTTCTGGAc                                                                      <  1:1082070/139‑1 (MQ=255)
                                                              tAAGTGTTCCGGAGGTCGTGGAAGCAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAACGCCGCCATTATTAAGTCTGGTGCTATCTGTTTTAGCCTCCTTCTGGAc                                                                      >  2:1124973/1‑139 (MQ=255)
                                                                 gtgtTCCGGAGGTCGTGGAAGCAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAACGCCGCCATTATTAAGTCTGGTGCTATCTGTTTTAGCCTCCTTCTGGACAAt                                                                   <  1:531954/139‑1 (MQ=255)
                                                                    ttCCGGAGGTCGTGGAAGCAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGa                                                                                                                                 >  2:683586/1‑74 (MQ=255)
                                                                    ttCCGGAGGTCGTGGAAGCAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGa                                                                                                                                 <  1:683586/74‑1 (MQ=255)
                                                                     tCCGGAGGTCGTGGAAGCAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAACGCCGCCATTATTAAGTCTGGTGCTATCTGTTTTAGCCTCCTTCTGGACAATCAGc                                                               >  2:562694/1‑139 (MQ=255)
                                                                      ccGGAGGTCGTGGAAGCAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAACGCCGCCa                                                                                                               <  2:263996/90‑1 (MQ=255)
                                                                      ccGGAGGTCGTGGAAGCAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAACGCCGCCa                                                                                                               >  1:263996/1‑90 (MQ=255)
                                                                       cGGAGGTCGTGGAAGCAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAACGCCGCCATTATTAAGTCTGGTGCTATCTGTTTTAGCCTCCTTCTGGACAATCAGCtt                                                             >  1:697537/1‑139 (MQ=255)
                                                                            gTCGTGGAAGCAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTg                                                                                                                                  <  1:1252185/65‑1 (MQ=255)
                                                                            gTCGTGGAAGCAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTg                                                                                                                                  >  2:1252185/1‑65 (MQ=255)
                                                                            gTCGTGGAAGCAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAACGCCGCCATTATTAAGTCTGGTGCTATCTGTTTTAGCCTCCTTCTGGACAATCAGCTTACCGc                                                        >  2:877610/1‑139 (MQ=255)
                                                                            gTCGTGGAAGCAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAACGCCGCCATTATTAAGTCTGGTGCTATCTGTTTTAGCCTCCTTCTGGACAATCAGCTTACCGc                                                        <  2:1573266/139‑1 (MQ=255)
                                                                              cGTGGAAGCAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAACGCCGCCATTATTAAGTCTGGTGCTATCTGTTTTAGCCTCCTTCTGGACAATCAGCTTACCGCCa                                                      <  1:1598600/139‑1 (MQ=255)
                                                                               gTGGAAGCAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAACGCCGCCATTATTAAGTCTGGTGCTAt                                                                                            <  2:622988/100‑1 (MQ=255)
                                                                               gTGGAAGCAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAACGCCGCCATTATTAAGTCTGGTGCTAt                                                                                            >  1:622988/1‑100 (MQ=255)
                                                                                     gCAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAACGCCGCCATTATTAAGTCTGGTGCTATCTGTTTTAGCCTCCTTCTGGACAATCAGCTTACCGCCAGTATCaa                                               <  2:391668/139‑1 (MQ=255)
                                                                                      cAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTg                                                                                                                                  >  1:1378088/1‑55 (MQ=255)
                                                                                      cAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTg                                                                                                                                  <  2:1378088/55‑1 (MQ=255)
                                                                                      cAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGa                                                                                                                                 >  1:912779/1‑56 (MQ=255)
                                                                                      cAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGa                                                                                                                                 <  2:912779/56‑1 (MQ=255)
                                                                                         ttAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAACGCCGCCATTATTAAGTCTGGTGCTATCTGTTTTAGCCTCCTTCTGGACa                                                                     <  2:931104/113‑1 (MQ=255)
                                                                                         ttAATGCGCCGCCGGATTGTTGCTCAACAGGCTTAGCCTCACCACTGCCCTGGACCTCCAGAACGCCGCCATTATTAAGTCTGGTGCTATCTGTTTTAGCCTCCTTCTGGACa                                                                     >  1:931104/1‑113 (MQ=255)
                                                                                             tGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGTCCTGAACCTCCAGAACGCCGCCATTATTAAGTCTGGTGCTATCTGTTTTAGCCTCCTTCTGGACAAt                                                                   <  2:1005674/111‑1 (MQ=255)
                                                                                             tGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAACGCCGCCATTATTAAGTCTGGTGCTATCTGTTTTAGCCTCCTTCTGGACAATCAGCTTACCGCCAGTATCAATGgtagta                                       >  1:572489/1‑139 (MQ=255)
                                                                                             ggcgcCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAACGCCGCCATTATTAAGTCTGGTGCTATCTGTTTTAGCCTCCTTCTGGACAAt                                                                   >  1:1005674/2‑111 (MQ=255)
                                                                                              gcgcCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGa                                                                                                                                 >  1:501359/1‑48 (MQ=255)
                                                                                              gcgcCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGa                                                                                                                                 <  2:501359/48‑1 (MQ=255)
                                                                                                 ccgccgGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAACGCCGCCATTATTAAGTCTGGTGCTATCTGTTTTAGcc                                                                                 >  2:1056678/1‑93 (MQ=255)
                                                                                                 ccgccgGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAACGCCGCCATTATTAAGTCTGGTGCTATCTGTTTTAGcc                                                                                 <  1:1056678/93‑1 (MQ=255)
                                                                                                   gccgGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAACGCCGCCATTATTAAGTCTGGTGCTATCTGTTTTAGCCTCCTTCTGGACAATCAGCTTACCGCCAGTATCAATGGTAGTATTTTTc                                 >  1:315189/1‑139 (MQ=255)
                                                                                                    ccgGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAAc                                                                                                                      <  2:1505713/53‑1 (MQ=255)
                                                                                                    ccgGATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAAc                                                                                                                      >  1:1505713/1‑53 (MQ=255)
                                                                                                       gATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAACGCCGCCATTATTAAGTCTGGTGCt                                                                                              >  2:85267/1‑74 (MQ=255)
                                                                                                       gATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAACGCCGCCATTATTAAGTCTGGTGCt                                                                                              <  1:85267/74‑1 (MQ=255)
                                                                                                       gATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAACGCCGCCATTATTAAGTCTGGTGCTATCTGTTTTAGCCTCCTTCTGGACAATCAGCTTACCGCCAGTATCAATGGTAGTATTTTTCGCCg                             <  2:315189/139‑1 (MQ=255)
                                                                                                       gATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAACGCCGCCATTATTAAGTCTGGTGCTATCTGTTTTAGCCTCCTTCTGGACAATCAGCTTACCGCCAGTATCAATGGTAGTATTTTTCGCCg                             <  1:877610/139‑1 (MQ=255)
                                                                                                       gATTGTTGCTCAACATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAACGCCGCCATTATTAAGTCTGGTACTATCTGTTTTAGCCTCCTTCTGGACAATCAGCTTACCGCCAGTATCAATGGTAGTATTTTTCGCCg                             <  1:206518/139‑1 (MQ=255)
                                                                                                            ttgCTCAACATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAACGCCGCCATTATTAAGTCTGGTGCTATCTGTTTTAGCCTCCTTCTGGACAATCAGCTTACCGCCAGTATCAATGGTAGTATTTTTCGCCGAGGtt                        >  1:502523/1‑139 (MQ=255)
                                                                                                                    cATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAAcgc                                                                                                                    >  1:169824/1‑39 (MQ=255)
                                                                                                                    cATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAAcgc                                                                                                                    <  2:169824/39‑1 (MQ=255)
                                                                                                                    cATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAACGCCGCCATTATTAAGTCTGGTGCTATCTGttt                                                                                      <  1:14889/69‑1 (MQ=255)
                                                                                                                    cATGCTTAGCCTCACCACCGCCCTGAACCTCCAGAACGCCGCCATTATTAAGTCTGGTGCTATCTGttt                                                                                      >  2:14889/1‑69 (MQ=255)
                                                                                                                      tGCTTAGCCTCACCACCGCCCTGAACCTCCAGAACGCCGCCATTATTAAGTCTGGTGCTATCTGTTTTAGCCTCCTTCTGGACAATCAGCTTACCGCCAGTATCAATGGTAGTATTTTTCGCCGAGGTTTTAGCCACTa              >  1:1371283/1‑139 (MQ=255)
                                                                                                                                caccacCGCCCTGAACCTCCAGAACGCCGCCATTATTAAGTCTGGTGCTAt                                                                                            <  1:1222187/51‑1 (MQ=255)
                                                                                                                                caccacCGCCCTGAACCTCCAGAACGCCGCCATTATTAAGTCTGGTGCTAt                                                                                            >  2:1222187/1‑51 (MQ=255)
                                                                                                                                 accaccGCCCTGAACCTCCAGAACGCCGCCATTATTAAGTCTGGTGCTATCTGTTTTAGCCTCCTTCTGGAc                                                                      <  2:527645/72‑1 (MQ=255)
                                                                                                                                 accaccGCCCTGAACCTCCAGAACGCCGCCATTATTAAGTCTGGTGCTATCTGTTTTAGCCTCCTTCTGGAc                                                                      >  1:527645/1‑72 (MQ=255)
                                                                                                                                  ccaccGCCCTGAACCTCCAGAACGCCGCCATTTTTAAGTCTGGTGCTATCTGTTTTAGCCTCCTTCTGGACAATCAGCTTACCGCCAGTATCAATGGTAGTATTTTTCGCCGAGGTTTTAGCCACTACCGTCAGTTCGc  >  1:959033/1‑139 (MQ=255)
                                                                                                                                        |                                                                                                                                    
TACATTTTTAGCTTCTGAATTCCTGATGTAGAAAGCATCACCATAACTGTTGGTTCCTTCGATAAGTGTTCCGGAGGTCGTGGAAGCAATTAATGCGCCGCCGGATTGTTGCTCAACATGCTTAGCCTCACCACCGTCCTGAACCTCCAGAACGCCGCCATTATTAAGTCTGGTGCTATCTGTTTTAGCCTCCTTCTGGACAATCAGCTTACCGCCAGTATCAATGGTAGTATTTTTCGCCGAGGTTTTAGCCACTACCGTCAGTTCGC  >  NC_000913/2780996‑2781264

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