Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,180,036 A→T 6.3% Q98L (CAG→CTG)  ygiB → DUF1190 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,180,0360AT 6.3% 223.1 / 4.6 63Q98L (CAG→CTG) ygiBDUF1190 family protein
Reads supporting (aligned to +/- strand):  ref base A (29/30);  new base T (2/2);  total (31/32)
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

GCTGACGACTGTTCAGCTGCAAACCCAGGCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGCACAGCAGCCGCTGTTC  >  NC_000913/3179891‑3180184
                                                                                                                                                 |                                                                                                                                                    
gCTGACGACTGTTCAGCTGCAAACCCAGGCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCcagca                                                                                                                                                   <  4:344302/149‑1 (MQ=255)
      gacTGTTCAGCTGCAAACCCAGGCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCAcc                                                                                                                                             >  2:652546/1‑149 (MQ=255)
      gacTGTTCAGCTGCAAACCCAGGCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCAcc                                                                                                                                             >  2:638275/1‑149 (MQ=255)
      gacTGTTCAGCTGCAAACCCAGGCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCAcc                                                                                                                                             >  2:440306/1‑149 (MQ=255)
              aGCTGCAAACCCAGGCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAgg                                                                                                                                     <  1:300982/149‑1 (MQ=255)
              aGCTGCAAACCCAGGCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAgg                                                                                                                                     <  1:608843/149‑1 (MQ=255)
               gCTGCAAACCCAGGCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGc                                                                                                                                    >  2:36926/1‑149 (MQ=255)
                 tGCAAACCCAGGCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCtg                                                                                                                                  >  4:166598/1‑149 (MQ=255)
                  gCAAACCCAGGCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTCCTGCGCCGAAATACGCGACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCtgg                                                                                                                                 <  1:47050/149‑1 (MQ=255)
                  gCAAACCCAGGCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCtgg                                                                                                                                 >  2:599750/1‑149 (MQ=255)
                  gCAAACCCAGGCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCtgg                                                                                                                                 >  4:384909/1‑149 (MQ=255)
                    aaaCCCAGGCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCtggca                                                                                                                               <  3:384909/149‑1 (MQ=255)
                    aaaCCCAGGCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCtggca                                                                                                                               <  4:615558/149‑1 (MQ=255)
                      aCCCAGGCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCtggcatg                                                                                                                             >  4:70671/1‑149 (MQ=255)
                        ccAGGCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCtggcatggc                                                                                                                           >  4:24082/1‑149 (MQ=255)
                        ccAGGCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCtggcatggc                                                                                                                           >  3:332568/1‑149 (MQ=255)
                            gCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCa                                                                                                                       <  2:267516/149‑1 (MQ=255)
                            gCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCa                                                                                                                       >  1:543691/1‑149 (MQ=255)
                                 aGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGaaaa                                                                                                                  <  1:36926/149‑1 (MQ=255)
                                    gcCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCa                                                                                                               >  3:63805/1‑149 (MQ=255)
                                     ccGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGTATGGCACCAGAAAACCAg                                                                                                              <  4:393501/149‑1 (MQ=255)
                                     ccGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAg                                                                                                              <  1:440306/149‑1 (MQ=255)
                                     ccGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAg                                                                                                              <  2:546542/149‑1 (MQ=255)
                                       gAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGgc                                                                                                            <  3:166598/149‑1 (MQ=255)
                                            tACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAgg                                                                                                       <  4:332568/149‑1 (MQ=255)
                                               caccGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGccc                                                                                                    >  1:81537/1‑149 (MQ=255)
                                                accGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCa                                                                                                   <  3:643847/149‑1 (MQ=255)
                                                accGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCa                                                                                                   <  1:338819/149‑1 (MQ=255)
                                                         aacaaTGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAg                                                                                          >  1:247164/1‑149 (MQ=255)
                                                            aaTGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCgg                                                                                       <  1:691420/149‑1 (MQ=255)
                                                              tGCGCTGAAAGAAGCCGAACGTACCGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGa                                                                                     >  3:440175/1‑149 (MQ=255)
                                                                   tGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTc                                                                                <  4:160115/149‑1 (MQ=255)
                                                                    gAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCt                                                                               <  2:358272/149‑1 (MQ=255)
                                                                     aaaGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTg                                                                              <  2:81537/149‑1 (MQ=255)
                                                                      aagaagCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGccc                                                                                                                                        <  1:235440/90‑1 (MQ=255)
                                                                      aagaagCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGccc                                                                                                                                        >  2:235440/1‑90 (MQ=255)
                                                                        gaagCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGAt                                                                           <  2:651545/149‑1 (MQ=255)
                                                                          agCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGc                                                                         >  4:301760/1‑149 (MQ=255)
                                                                                           cGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGtt                                                          >  2:619755/1‑147 (MQ=255)
                                                                                           cGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGtt                                                          <  1:619755/147‑1 (MQ=255)
                                                                                            gAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACa                                                       <  1:652546/149‑1 (MQ=255)
                                                                                                     cACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGc                                              <  4:543045/149‑1 (MQ=255)
                                                                                                               gACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATggg                                    >  3:384166/1‑149 (MQ=255)
                                                                                                                aCTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGggc                                   <  1:599750/149‑1 (MQ=255)
                                                                                                                     gtTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGggcggcgg                              >  4:590986/1‑149 (MQ=255)
                                                                                                                      ttGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGggcggcggc                             >  3:488977/1‑149 (MQ=255)
                                                                                                                       tGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGa                                                                            <  4:102927/101‑1 (MQ=255)
                                                                                                                       tGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGa                                                                            >  3:102927/1‑101 (MQ=255)
                                                                                                                         cTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACa                                                       <  1:626146/120‑1 (MQ=255)
                                                                                                                         cTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACa                                                       >  2:626146/1‑120 (MQ=255)
                                                                                                                               ttGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGAttt                    <  3:69/149‑1 (MQ=255)
                                                                                                                                 ggTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGc                  <  1:638275/149‑1 (MQ=255)
                                                                                                                                 ggTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGc                  <  4:440175/149‑1 (MQ=255)
                                                                                                                                  gTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGca                 >  3:304711/1‑149 (MQ=255)
                                                                                                                                      aGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGCAcagc             >  2:692112/1‑149 (MQ=255)
                                                                                                                                        gTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGCAcagcag           >  4:250404/1‑149 (MQ=255)
                                                                                                                                        gTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGCAcagcag           >  3:305452/1‑149 (MQ=255)
                                                                                                                                         ctggTGCCTGCTGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGca                 <  4:12945/139‑1 (MQ=255)
                                                                                                                                         ctggTGCCTGCTGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGca                 >  3:12945/4‑142 (MQ=255)
                                                                                                                                           ggTGCCTGCTGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGa                                        <  3:554169/116‑1 (MQ=255)
                                                                                                                                           ggTGCCTGCTGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGa                                        >  4:554169/2‑117 (MQ=255)
                                                                                                                                            gTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGCACAGCAGCCGc       >  3:409358/1‑149 (MQ=255)
                                                                                                                                            gTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGCACAGCAGCCGc       >  2:464403/1‑149 (MQ=255)
                                                                                                                                             tGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGCACAGCAGCCGCt      <  1:464403/149‑1 (MQ=255)
                                                                                                                                                 agcagGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGCACAGCAGCCGCTGTtc  >  3:538940/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GCTGACGACTGTTCAGCTGCAAACCCAGGCAAAAGCGCCGAATGTACCACCGCGTACAACAATGCGCTGAAAGAAGCCGAACGTACTGCGCCGAAATACGCCACCCGTGAAGACTGTGTTGCTGAATTTGGTGAAGGTCAGTGCCAGCAGGCACCAGCCCAGGCTGGCATGGCACCAGAAAACCAGGCGCAGGCCCAGCAATCCAGCGGGAGTTTCTGGATGCCGCTGATGGCCGGTTACATGATGGGGCGTCTGATGGGCGGCGGCGCGGGATTTGCACAGCAGCCGCTGTTC  >  NC_000913/3179891‑3180184

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