Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 645,034 A→G V257A (GTG→GCG)  ybdS ← YbdS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913645,0340AG99.4% 499.4 / ‑4.7 161V257A (GTG→GCG) ybdSYbdS
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (66/94);  total (66/95)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.75e-01

Only 100 of 166 total aligned reads displayed.

TCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCACGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAGCTAATCTGGATACCG  >  NC_000913/644944‑645123
                                                                                          |                                                                                         
tctctGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGcgcg                                                                                         <  2:4375569/93‑1 (MQ=255)
    tGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAgg                                                                                     >  2:2284505/1‑93 (MQ=255)
    tGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAgg                                                                                     >  1:1815021/1‑93 (MQ=255)
    tGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAgg                                                                                     >  1:3215558/1‑93 (MQ=255)
      cGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGaa                                                                                   <  1:1292932/93‑1 (MQ=255)
      cGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCACGTAGGaa                                                                                   <  1:3343691/93‑1 (MQ=255)
       gtgtCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGaaa                                                                                  <  2:4432966/93‑1 (MQ=255)
       gtgtCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGaaa                                                                                  <  1:3149153/93‑1 (MQ=255)
        tgtCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAg                                                                                 <  2:1734984/93‑1 (MQ=255)
         gtCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGc                                                                                <  1:2188120/93‑1 (MQ=255)
          tCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGcc                                                                               <  1:3622330/93‑1 (MQ=255)
          tCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGcc                                                                               <  1:2608655/93‑1 (MQ=255)
           cAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCa                                                                              <  1:4432001/93‑1 (MQ=255)
            aGCGCACGCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCAc                                                                             <  2:1064644/93‑1 (MQ=255)
               gcACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGgc                                                                          <  2:2684814/93‑1 (MQ=255)
                 aCCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGgcgc                                                                        <  1:214235/93‑1 (MQ=255)
                    cATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAAt                                                                     <  1:4239296/93‑1 (MQ=255)
                    cATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAAt                                                                     >  2:1160489/1‑93 (MQ=255)
                    cATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAAt                                                                     <  1:689642/93‑1 (MQ=255)
                     aTGGTTTTCAGTTCATCCCCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATg                                                                    <  2:393623/93‑1 (MQ=255)
                     aTGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATg                                                                    >  2:2618216/1‑93 (MQ=255)
                      tGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCTCTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGa                                                                   <  1:4342248/93‑1 (MQ=255)
                      tGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGa                                                                   >  1:2995024/1‑93 (MQ=255)
                      tGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGa                                                                   <  1:1220462/93‑1 (MQ=255)
                        gTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATa                                                                 <  1:4103848/93‑1 (MQ=255)
                        gTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATa                                                                 >  1:3971799/1‑93 (MQ=255)
                         ttttCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATaa                                                                <  1:4099441/93‑1 (MQ=255)
                          tttCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTTATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAg                                                               <  1:1010097/93‑1 (MQ=255)
                          tttCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAg                                                               <  1:4631676/93‑1 (MQ=255)
                           ttCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGc                                                              <  2:2433438/93‑1 (MQ=255)
                             cAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCaa                                                            >  1:2770399/1‑93 (MQ=255)
                              aGTTCATCCCCCTCCCAGGTTGCCACTTCTTCACTGTGTGTTATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAg                                                           <  2:3938741/93‑1 (MQ=255)
                               gTTCATCCCCCGCCCTGGTTGCCACTTCTTCACTGTGTGTAATTTCCGGTTTGTACAGCCCGTAGGAAAGCCACGGCGCAATGATAAGCAAGa                                                          <  2:3205105/93‑1 (MQ=255)
                               gTTCATCACGCGGCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGa                                                          <  1:531480/93‑1 (MQ=255)
                               gTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGa                                                          <  2:2216638/93‑1 (MQ=255)
                               gTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGa                                                          >  1:963330/1‑93 (MQ=255)
                               gTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGa                                                          >  1:1606595/1‑93 (MQ=255)
                               gTTCAGCGCCCGCCCAGGTTGCCACTTCTTCACTGTGTGTTATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGa                                                          <  2:2396489/93‑1 (MQ=255)
                                ttCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGAt                                                         >  1:2523844/1‑93 (MQ=255)
                                ttCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGAt                                                         >  2:3965067/1‑93 (MQ=255)
                                ttCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGAt                                                         >  1:2405112/1‑93 (MQ=255)
                                ttCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGAt                                                         >  2:1183413/1‑93 (MQ=255)
                                  catcaCCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATaa                                                       >  2:190632/1‑93 (MQ=255)
                                  catcaCCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATaa                                                       <  1:3495783/93‑1 (MQ=255)
                                      ccccGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAAcccc                                                   <  2:3619955/92‑1 (MQ=255)
                                      aCCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAAcccc                                                   >  1:2017641/1‑93 (MQ=255)
                                        ccGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCaa                                                 <  1:934558/93‑1 (MQ=255)
                                         cGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAAc                                                <  1:2607214/93‑1 (MQ=255)
                                         cGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGAGCGTAGGAAAGCCACGGGGCAATGATAAGCAAGATAGCCCCAAc                                                >  2:4140699/1‑93 (MQ=255)
                                            ccAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGAGCGTAGGAAAGACACGGCGCAATGATAAGCAAGATACCCCCAACCgg                                             >  2:134678/1‑93 (MQ=255)
                                              aGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCa                                           <  2:589714/93‑1 (MQ=255)
                                              aGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCa                                           >  1:3574755/1‑93 (MQ=255)
                                               ggTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAg                                          >  2:1692667/1‑93 (MQ=255)
                                                gTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAgg                                         >  1:4577387/1‑93 (MQ=255)
                                                gTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAgg                                         >  1:363237/1‑93 (MQ=255)
                                                gTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAgg                                         >  2:35511/1‑93 (MQ=255)
                                                 ttGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGa                                        >  2:18963/1‑93 (MQ=255)
                                                  tGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGaa                                       >  2:3357770/1‑93 (MQ=255)
                                                  tGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGaa                                       <  1:1925692/93‑1 (MQ=255)
                                                     cACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCa                                    <  1:653097/93‑1 (MQ=255)
                                                     cACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCa                                    <  2:873296/93‑1 (MQ=255)
                                                     cACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCa                                    >  1:2198433/1‑93 (MQ=255)
                                                         tcttcACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAgg                                <  1:102381/93‑1 (MQ=255)
                                                           ttcACTGTGGGGGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGaa                              <  2:145384/93‑1 (MQ=255)
                                                           ttcACTATGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCATTGATAAGCAAGATATCCCCAACCGGCAGGAAGCAGAGGaa                              <  2:2627347/93‑1 (MQ=255)
                                                            tcACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAAc                             <  2:1762663/93‑1 (MQ=255)
                                                            tcACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAAc                             >  1:1883510/1‑93 (MQ=255)
                                                             cACTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAAcc                            >  1:3731416/1‑93 (MQ=255)
                                                              aCTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCa                           >  1:1460239/1‑93 (MQ=255)
                                                               cTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCAc                          >  2:4209300/1‑93 (MQ=255)
                                                               cTGTGTGTGATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCAc                          >  1:1536321/1‑93 (MQ=255)
                                                               cTGTGTGTGATTTCCGGTTTGTACAGCGCGGAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCAc                          >  2:2072769/1‑93 (MQ=255)
                                                                tgtgtgtgATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACt                         >  2:859359/1‑93 (MQ=255)
                                                                tgtgtgtgATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACt                         >  2:1215924/1‑93 (MQ=255)
                                                                 gtgtgtgATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTg                        <  2:260541/93‑1 (MQ=255)
                                                                    tgtgATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCaa                     <  2:2793374/93‑1 (MQ=255)
                                                                    tgtgATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCaa                     >  2:3503687/1‑93 (MQ=255)
                                                                      tgATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAAcc                   <  1:2789590/93‑1 (MQ=255)
                                                                      tgATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGAAAcc                   >  2:940182/1‑93 (MQ=255)
                                                                       gATTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCa                  <  1:248372/93‑1 (MQ=255)
                                                                        aTTTCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAg                 <  1:4467070/93‑1 (MQ=255)
                                                                         tttCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAGc                <  2:1867075/93‑1 (MQ=255)
                                                                          ttCCGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAGCt               >  2:3671108/1‑93 (MQ=255)
                                                                            ccGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAGCTaa             <  2:378736/93‑1 (MQ=255)
                                                                             cGGTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAGCTAAt            <  2:1655918/93‑1 (MQ=255)
                                                                              ggTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAGCTAATc           <  2:4623774/93‑1 (MQ=255)
                                                                              ggTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAGCTAATc           <  1:434312/93‑1 (MQ=255)
                                                                               gTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAGCTAATCt          >  2:199423/1‑93 (MQ=255)
                                                                               gTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAGCTAATCt          <  2:653668/93‑1 (MQ=255)
                                                                               gTTTGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAGCTAATCt          <  2:781835/93‑1 (MQ=255)
                                                                                tttGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGCGGAACCACTGCAACCACCTAATCTg         >  2:3754857/1‑93 (MQ=255)
                                                                                tttGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAGCTAATCTg         <  2:1332410/93‑1 (MQ=255)
                                                                                 ttGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAGCTAATCTgg        >  2:1958213/1‑93 (MQ=255)
                                                                                 ttGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAGCTAATCTgg        <  2:2547496/93‑1 (MQ=255)
                                                                                 ttGTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAGCTAATCTgg        <  1:2675849/93‑1 (MQ=255)
                                                                                   gTACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAGCTAATCTGGAt      <  2:1504367/93‑1 (MQ=255)
                                                                                    tACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAGCTAATCTGGATa     >  2:213699/1‑93 (MQ=255)
                                                                                    tACAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAGCTAATCTGGATa     <  2:2883255/93‑1 (MQ=255)
                                                                                      cAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAGCTAATCTGGATAcc   <  2:3132631/93‑1 (MQ=255)
                                                                                      cAGCGCGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAGCTAATCTGGATAcc   <  2:4529656/93‑1 (MQ=255)
                                                                                          |                                                                                         
TCTCTGCGTGTCAGCGCACCCATGGTTTTCAGTTCATCACCCGCCCAGGTTGCCACTTCTTCACTGTGTGTGATTTCCGGTTTGTACAGCACGTAGGAAAGCCACGGCGCAATGATAAGCAAGATAACCCCAACCGGCAGGAAGCAGAGGAACCACTGCAACCAGCTAATCTGGATACCG  >  NC_000913/644944‑645123

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