Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,057,859 (T)9→8 intergenic (+97/‑167) yihK → / → yihL putative GTP‑binding factor/putative transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,057,8510T.99.3% 2558.8 / 7.7 572intergenic (+89/‑175)yihK/yihLputative GTP‑binding factor/putative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (2/2);  new base . (320/246);  total (322/250)
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 670 total aligned reads displayed.

CTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTA  >  NC_000913/4057759‑4057939
                                                                                            |                                                                                        
cTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGt                                                                                          >  1:8237829/1‑93 (MQ=255)
cTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGt                                                                                          <  2:6455685/93‑1 (MQ=255)
 tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtt                                                                                         >  1:8045649/1‑93 (MQ=255)
 tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtt                                                                                         >  2:9098143/1‑93 (MQ=255)
  gACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttt                                                                                        <  2:736818/93‑1 (MQ=255)
  gACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttt                                                                                        >  2:9395590/1‑93 (MQ=255)
  gACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttt                                                                                        >  1:1918468/1‑93 (MQ=255)
    cGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttt                                                                                      <  1:4415615/93‑1 (MQ=255)
     ggAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtttttt                                                                                     >  1:3129920/1‑93 (MQ=255)
     ggAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtttttt                                                                                     >  2:6280870/1‑93 (MQ=255)
      gAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttt                                                                                    >  1:1377935/1‑93 (MQ=255)
        aaaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttg                                                                                  <  2:9223541/93‑2 (MQ=255)
          aCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttgcc                                                                                <  2:1302669/93‑4 (MQ=255)
          aCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttgcc                                                                                >  1:8489526/1‑90 (MQ=255)
           cGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttgcct                                                                               <  2:7695347/93‑5 (MQ=255)
           cGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttgcct                                                                               <  1:4235320/93‑5 (MQ=255)
             aTCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGt                                                                            >  1:1540906/1‑93 (MQ=255)
              tcgtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTc                                                                           <  1:2381668/93‑1 (MQ=255)
              tcgccgCCCCCCCAACCGCGCCCCGAATGACGATTAATTTCTTTTTTTAGTCATAAAAATCCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTc                                                                           <  2:7215648/93‑1 (MQ=255)
                 tcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTg                                                                        <  1:865221/93‑1 (MQ=255)
                  cgccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGa                                                                       <  2:4657136/93‑1 (MQ=255)
                      gcgcCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAgcgc                                                                   >  1:7710697/1‑93 (MQ=255)
                        gcCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGcac                                                                 >  1:10015544/1‑93 (MQ=255)
                        gcCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGcac                                                                 <  1:3372316/93‑1 (MQ=255)
                         ccAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGcacc                                                                >  1:8342611/1‑92 (MQ=255)
                         ccAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGcaca                                                                <  2:8218310/93‑1 (MQ=255)
                          cAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCacac                                                               <  1:5245696/93‑1 (MQ=255)
                             ccGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCaa                                                            >  1:6037318/1‑93 (MQ=255)
                             ccGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCaa                                                            <  2:6762148/93‑1 (MQ=255)
                              cgcgCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAAt                                                           >  2:8745387/1‑93 (MQ=255)
                              cgcgCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAAt                                                           <  1:5032668/93‑1 (MQ=255)
                              cgcgCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAAt                                                           <  1:4875983/93‑1 (MQ=255)
                               gcgcACCGAAAGACAATTAATTTCGTTTTTTATTCATAAAAACCCTGCCAGCGATGGCAGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATg                                                          <  1:2253307/93‑1 (MQ=255)
                               gccaaCCGAAAGAGGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATg                                                          <  1:2102858/89‑1 (MQ=255)
                                  cACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGa                                                       >  2:342199/1‑93 (MQ=255)
                                   aCCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                      >  1:6566418/1‑93 (MQ=255)
                                   aCCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                      >  1:2657595/1‑93 (MQ=255)
                                   aCCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCCGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                      >  2:9861057/1‑93 (MQ=255)
                                    ccGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGt                                                     >  2:3084142/1‑93 (MQ=255)
                                    ccGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGt                                                     >  2:3953221/1‑93 (MQ=255)
                                    ccGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGt                                                     >  1:9976022/1‑93 (MQ=255)
                                    ccGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGt                                                     <  1:4004750/93‑1 (MQ=255)
                                     cGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTa                                                    >  2:4332395/1‑93 (MQ=255)
                                     cGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTa                                                    <  2:3213991/93‑1 (MQ=255)
                                     cGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTa                                                    >  2:1683286/1‑93 (MQ=255)
                                      gAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAt                                                   <  2:375631/93‑1 (MQ=255)
                                       aaaGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtc                                                  <  2:6967485/93‑1 (MQ=255)
                                        aaGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtct                                                 <  1:3017490/93‑1 (MQ=255)
                                         aGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtctc                                                >  1:7284817/1‑93 (MQ=255)
                                           aCGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTccc                                              <  1:1839442/93‑1 (MQ=255)
                                            cGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCt                                             >  2:9647294/1‑93 (MQ=255)
                                            cGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCt                                             >  2:4137930/1‑93 (MQ=255)
                                            cGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCt                                             >  2:8553561/1‑93 (MQ=255)
                                             gATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCtt                                            >  2:5056837/1‑93 (MQ=255)
                                             gATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCtt                                            >  1:6160291/1‑93 (MQ=255)
                                             gATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCtt                                            >  1:2918401/1‑93 (MQ=255)
                                              aTTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTta                                           <  2:1073342/93‑1 (MQ=255)
                                                 aaTTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATaa                                        <  2:886737/93‑1 (MQ=255)
                                                   tttCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGt                                      >  1:7333224/1‑93 (MQ=255)
                                                      cGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtg                                   <  1:2778163/93‑1 (MQ=255)
                                                       gTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgt                                  <  1:4665220/93‑1 (MQ=255)
                                                          ttttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGt                               <  1:6437264/93‑1 (MQ=255)
                                                           tttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGtt                              >  2:9596937/1‑93 (MQ=255)
                                                            ttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTa                             >  2:2903611/1‑93 (MQ=255)
                                                            ttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTa                             <  1:6257573/93‑1 (MQ=255)
                                                            ttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTa                             <  1:9645005/93‑1 (MQ=255)
                                                            ttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTa                             >  2:569521/1‑93 (MQ=255)
                                                              aGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACg                           <  2:8334213/93‑1 (MQ=255)
                                                               gTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGt                          >  2:2311450/1‑93 (MQ=255)
                                                               gTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGt                          <  2:2759351/93‑1 (MQ=255)
                                                                tCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTc                         >  2:2098846/1‑93 (MQ=255)
                                                                 cATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCa                        >  1:740164/1‑93 (MQ=255)
                                                                  aTAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCaa                       <  1:5457468/93‑1 (MQ=255)
                                                                  aTAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCaa                       >  2:4531929/1‑93 (MQ=255)
                                                                  aTAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCaa                       <  1:7017817/93‑1 (MQ=255)
                                                                  aTAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCaa                       <  1:9888137/93‑1 (MQ=255)
                                                                   tAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCCCACCCATGCGAGTATCTCCCCTATAAGTCTGTGGGTACGTCCAt                      >  2:1606398/1‑93 (MQ=255)
                                                                   tAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAAt                      >  1:5744850/1‑93 (MQ=255)
                                                                    aaaaaaCCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATa                     >  2:1805163/1‑93 (MQ=255)
                                                                    aaaaaaCCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATa                     >  1:7834169/1‑93 (MQ=255)
                                                                     aaaaaCCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATag                    >  2:5622388/1‑93 (MQ=255)
                                                                     aaaaaCCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATag                    >  2:2559150/1‑93 (MQ=255)
                                                                      aaaaCCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATaga                   >  2:1196340/1‑93 (MQ=255)
                                                                      aaaaCCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATaga                   <  2:5938301/93‑1 (MQ=255)
                                                                       aaaCCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagag                  <  1:8926152/93‑1 (MQ=255)
                                                                        aaCCTGCCAGCGATGGCAGCTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                 <  1:3685897/93‑1 (MQ=255)
                                                                         aCCTGGCAGCGATGGCAGGTTTTTTGTGCCTGTCGTGAGTGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagagag                <  1:6804068/93‑1 (MQ=255)
                                                                         aCCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagagcg                >  2:4820393/1‑91 (MQ=255)
                                                                         aCCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagagag                >  1:9425017/1‑93 (MQ=255)
                                                                         aCCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagagag                >  2:7527430/1‑93 (MQ=255)
                                                                             tccAGCGATGTCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCttt            <  2:6094540/92‑1 (MQ=255)
                                                                               cAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGt          >  1:3846130/1‑93 (MQ=255)
                                                                                 gCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTca        >  1:511714/1‑93 (MQ=255)
                                                                                    aTGGCAGGCTGTTTTGTCCTGGCGTGGCCGAACGCCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACaa     <  1:1291557/93‑1 (MQ=255)
                                                                                    aTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACaa     <  2:1431335/93‑1 (MQ=255)
                                                                                    aTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACaa     >  1:5592171/1‑93 (MQ=255)
                                                                                     tGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAAt    <  2:517616/93‑1 (MQ=255)
                                                                                      ggCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAAtt   >  2:4911999/1‑93 (MQ=255)
                                                                                       gCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTa  <  2:8011401/93‑1 (MQ=255)
                                                                                       gCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTa  >  1:9158450/1‑93 (MQ=255)
                                                                                            |                                                                                        
CTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTA  >  NC_000913/4057759‑4057939

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: