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.98.5% 906.6 / 4.7 202intergenic (+89/‑175)yihK/yihLputative GTP‑binding factor/putative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base . (111/88);  total (113/89)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

Only 100 of 232 total aligned reads displayed.

CTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTA  >  NC_000913/4057759‑4057939
                                                                                            |                                                                                        
cTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGt                                                                                          >  1:3427763/1‑93 (MQ=255)
 tcaCGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGAAGATTAACTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGCTGGCAGGtt                                                                                         >  1:3151749/3‑93 (MQ=255)
    cGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttt                                                                                      <  2:2847636/93‑1 (MQ=255)
    cGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttt                                                                                      >  1:3787505/1‑93 (MQ=255)
    cGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttt                                                                                      >  1:139882/1‑93 (MQ=255)
      gAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttt                                                                                    >  1:1292148/1‑93 (MQ=255)
       aaaaCGTTCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTATTTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtttttttt                                                                                   <  2:4207328/93‑1 (MQ=255)
        aaaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttg                                                                                  >  1:3563050/1‑92 (MQ=255)
        aaaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttg                                                                                  >  2:1589682/1‑92 (MQ=255)
          aCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttgcc                                                                                >  1:4063923/1‑90 (MQ=255)
           cGATCGTCGCCGCGCCAACCGCGCGCCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttgcct                                                                               >  2:1977218/1‑89 (MQ=255)
            gATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTg                                                                             <  1:3959504/93‑1 (MQ=255)
            gATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTg                                                                             >  2:3973348/1‑93 (MQ=255)
             aTCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTCCTGt                                                                            >  2:3576423/1‑93 (MQ=255)
             aTCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGt                                                                            <  1:2377587/93‑1 (MQ=255)
             aTCGTCACCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGt                                                                            <  2:1647694/93‑1 (MQ=255)
               cgtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCg                                                                          >  1:2331886/1‑93 (MQ=255)
                gtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGt                                                                         >  1:731617/1‑93 (MQ=255)
                 tcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTg                                                                        <  1:1968106/93‑1 (MQ=255)
                      gcgcCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAgcgc                                                                   >  2:3408621/1‑93 (MQ=255)
                      gcgcCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAgcgc                                                                   <  1:2562175/93‑1 (MQ=255)
                      gcgcCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAgcgc                                                                   >  2:539961/1‑93 (MQ=255)
                        gcCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGcac                                                                 >  2:2938292/1‑93 (MQ=255)
                          cAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCacac                                                               >  2:3552062/1‑93 (MQ=255)
                           aaCCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACAcc                                                              >  1:3157261/1‑93 (MQ=255)
                              cgcgCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAAt                                                           <  1:2968238/93‑1 (MQ=255)
                                  cACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGa                                                       >  2:1226313/1‑93 (MQ=255)
                                   aCCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                      >  1:15067/1‑93 (MQ=255)
                                     cGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAATCCTGCCAGCGATGGAAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTa                                                    <  1:3319217/93‑1 (MQ=255)
                                     cGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTa                                                    >  2:3483377/1‑93 (MQ=255)
                                      gAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAt                                                   >  2:1071668/1‑93 (MQ=255)
                                      gAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAt                                                   >  1:1074786/1‑93 (MQ=255)
                                      gAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAt                                                   >  2:4628326/1‑93 (MQ=255)
                                      gAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAt                                                   >  1:16847/1‑93 (MQ=255)
                                        aaGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtct                                                 >  1:4735394/1‑93 (MQ=255)
                                        aaGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtct                                                 <  1:1195400/93‑1 (MQ=255)
                                        aaGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtct                                                 >  2:2449305/1‑93 (MQ=255)
                                         aGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtctc                                                >  2:1729591/1‑93 (MQ=255)
                                         aGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtctc                                                <  1:135391/93‑1 (MQ=255)
                                          gACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTcc                                               >  1:482074/1‑93 (MQ=255)
                                           aCGATTATTTTCTTTTTTTAGTCATAAAAACCCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTccc                                              <  2:3524134/93‑1 (MQ=255)
                                           aCGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTccc                                              <  1:4155375/93‑1 (MQ=255)
                                           aCGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTccc                                              >  2:3760144/1‑93 (MQ=255)
                                            cGATTAATTTCTTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCt                                             <  2:4404264/93‑1 (MQ=255)
                                            cGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCt                                             >  1:935045/1‑93 (MQ=255)
                                             gATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCtt                                            >  2:1310657/1‑93 (MQ=255)
                                              gttAATTCCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTta                                           >  2:2882224/2‑93 (MQ=255)
                                              aTTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTta                                           >  1:3182948/1‑93 (MQ=255)
                                              aTTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTta                                           <  1:3767727/93‑1 (MQ=255)
                                              aTTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTta                                           <  1:2266646/93‑1 (MQ=255)
                                               ttAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTtat                                          >  1:4017904/1‑93 (MQ=255)
                                               ttAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTGTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTtat                                          >  1:4215157/1‑93 (MQ=255)
                                                tAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTtata                                         >  2:3153819/1‑93 (MQ=255)
                                                tAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTtata                                         <  1:3333045/93‑1 (MQ=255)
                                                 aaTTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATaa                                        >  2:2750910/1‑93 (MQ=255)
                                                 aaTTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATaa                                        <  2:4564884/93‑1 (MQ=255)
                                                   tttCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGt                                      <  1:455772/93‑1 (MQ=255)
                                                     tCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCt                                    <  2:4257130/93‑1 (MQ=255)
                                                     tCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCt                                    >  1:1623435/1‑93 (MQ=255)
                                                       gTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgt                                  <  2:1181030/93‑1 (MQ=255)
                                                       gTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgt                                  >  1:2201245/1‑93 (MQ=255)
                                                         tttttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                                <  1:3847287/93‑1 (MQ=255)
                                                           tttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGtt                              <  2:4617718/93‑1 (MQ=255)
                                                           tttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGtt                              >  1:809185/1‑93 (MQ=255)
                                                            ttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTa                             <  1:4186072/93‑1 (MQ=255)
                                                             ttgTCAAAAAAAACCCGCCAGGGATGGCAGGTTTGTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTa                             <  2:2640533/91‑1 (MQ=255)
                                                             tAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTAc                            <  2:3361805/93‑1 (MQ=255)
                                                               gTCATAAAAAACCTGCCAGCGGTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGt                          <  1:3726979/93‑1 (MQ=255)
                                                               gTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGt                          >  1:1885204/1‑93 (MQ=255)
                                                                 cATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCa                        >  1:2699567/1‑93 (MQ=255)
                                                                   tAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAAt                      >  1:4600838/1‑93 (MQ=255)
                                                                   tAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAAt                      <  2:1429331/93‑1 (MQ=255)
                                                                    aaaaaaCCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATa                     <  1:2786596/93‑1 (MQ=255)
                                                                    aaaaaaCCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATa                     <  1:2000719/93‑1 (MQ=255)
                                                                     aaaaaCCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATag                    >  2:49370/1‑93 (MQ=255)
                                                                     aaaaaCCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATag                    <  2:789483/93‑1 (MQ=255)
                                                                     aaaaaCCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATag                    >  1:2792824/1‑93 (MQ=255)
                                                                      aaaaCCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGATTATCTCCCTTATAAGGCCGAGGATACGTCAATaga                   >  2:4221314/1‑93 (MQ=255)
                                                                      aaaaCCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATaga                   <  1:1223971/93‑1 (MQ=255)
                                                                       aaaCCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGGGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagag                  <  1:879721/93‑1 (MQ=255)
                                                                       aaaCCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagag                  <  1:1141750/93‑1 (MQ=255)
                                                                        aaCCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                 >  2:2706101/1‑93 (MQ=255)
                                                                        aaCCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                 <  1:4585477/93‑1 (MQ=255)
                                                                        aaCCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                 >  2:1049985/1‑93 (MQ=255)
                                                                         aCCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagagag                >  2:1860105/1‑93 (MQ=255)
                                                                          ccTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGc               <  2:2287596/93‑1 (MQ=255)
                                                                           cTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCt              <  2:2278392/93‑1 (MQ=255)
                                                                           cTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCt              >  1:1746344/1‑93 (MQ=255)
                                                                            tGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCtt             >  2:2217328/1‑93 (MQ=255)
                                                                             gCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCttt            >  1:1716503/1‑93 (MQ=255)
                                                                             gCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCttt            >  2:255755/1‑93 (MQ=255)
                                                                             gCCAGCGATGGCAGG‑‑TTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTg           <  1:3659858/93‑1 (MQ=255)
                                                                               cAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGt          >  2:2679214/1‑93 (MQ=255)
                                                                                aGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTc         >  1:713482/1‑93 (MQ=255)
                                                                                 gCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTca        >  1:2612818/1‑93 (MQ=255)
                                                                                   gATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTcaca      >  2:2091529/1‑93 (MQ=255)
                                                                                    aTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATTAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCCCaa     >  1:2830712/1‑93 (MQ=255)
                                                                                     tGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAAt    >  2:3269023/1‑93 (MQ=255)
                                                                                      ggCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAAtt   <  2:58125/93‑1 (MQ=255)
                                                                                      ggCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAAtt   >  2:259488/1‑93 (MQ=255)
                                                                                            |                                                                                        
CTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTA  >  NC_000913/4057759‑4057939

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