Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,636,747 T→C 8.8% T222A (ACC→GCC)  yhiN ← putative oxidoreductase YhiN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,636,7470TC 8.8% 249.7 / 12.5 91T222A (ACC→GCC) yhiNputative oxidoreductase YhiN
Reads supporting (aligned to +/- strand):  ref base T (40/43);  new base C (3/5);  total (43/48)
Fisher's exact test for biased strand distribution p-value = 7.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.60e-01

GCTGCTCATTCAGGAAGGTTTCGAGGTCCACATCCGGTAGCAGATTGATGCTGACAAATTCCCCCGGTTGCCAGTAGCTTGAAATCTGCAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCACCTGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGGCAGCACGTTGAGGCCAAATTGTTCGGCAATCTTATAACC  >  NC_000913/3636600‑3636894
                                                                                                                                                   |                                                                                                                                                   
gctgctCATTCAGGAAGGTTTCGAGGTCCACATCCGGTAGCAGATTGATGCTGACAAATTCCCCCGGTTGCCAGTAGCTTGAAATCTGCAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGc                                                                                                                                                   >  1:366199/1‑150 (MQ=255)
gctgctCATTCAGGAAGGTTTCGAGGTCCACATCCGGTAGCAGATTGATGCTGACAAATTCCCCCGGTTGCCAGTAGCTTGAAATCTGCAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGc                                                                                                                                                   >  2:321552/1‑150 (MQ=255)
  tgctCATTCAGGAAGGTTTCGAGGTCCACATCCGGTAGCAGATTGATGCTGACAAATTCCCCCGGTTGCCAGTAGCTTGAAATCTGCAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCg                                                                                                                                                 >  1:1545431/1‑150 (MQ=255)
   gctCATTCAGGAAGGTTTCGAGGTCCACATCCGGTAGCAGATTGATGCTGACAAATTCCCCCGGTTGCCAGTAGCTTGAAATCTGCAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGCGCCGt                                                                                                                                                >  1:1484840/1‑150 (MQ=255)
      cATTCAGGAAGGTTTCGAGGTCCACATCCGGTAGCAGATTGATGCTGACAAATTCCCCCGGTTGCCAGTAGCTTGAAATCTGCAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGtttt                                                                                                                                             >  1:684231/1‑150 (MQ=255)
      cATTCAGGAAGGTTTCGAGGTCCACATCCGGTAGCAGATTGATGCTGACAAATTCCCCCGGTTGCCAGTAGCTTGAAATCTGCAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGtttt                                                                                                                                             >  2:318097/1‑150 (MQ=255)
       aTTCAGGAAGGTTTCGAGGTCCACATCCGGTAGCAGATTGATGCTGACAAATTCCCCCGGTTGCCAGTAGCTTGAAATCTGCAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTc                                                                                                                                            <  2:684231/150‑1 (MQ=255)
          cAGGAAGGTTTCGAGGTCCACATCCGGTAGCAGATTGATGCTGACAAATTCCCCCGGTTGCCAGTAGCTTGAAATCTGCAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGc                                                                                                                                         >  2:223291/1‑150 (MQ=255)
          cAGGAAGGTTTCGAGGTCCACATCCGGTAGCAGATTGATGCTGACAAATTCCCCCGGTTGCCAGTAGCTTGAAATCTGCAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGc                                                                                                                                         >  1:1461396/1‑150 (MQ=255)
          cAGGAAGGTTTCGAGGTCCACATCCGGTAGCAGATTGATGCTGACAAATTCCCCCGGTTGCCAGTAGCTTGAAATCTGCAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGc                                                                                                                                         >  2:91737/1‑150 (MQ=255)
          cAGGAAGGTTTCGAGGTCCACATCCGGTAGCAGATTGATGCTGACAAATTCCCCCGGTTGCCAGTAGCTTGAAATCTGCAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGCACCGTTTTccgt                                                                                                                                         <  2:257851/150‑4 (MQ=255)
           aGGACGGTTTCGAGGTCCACATGCGGTAGCAGTTTTATGCTGACAAATTCCCCAGGTTGCCAGTAGCTTGAAATCTGCAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGCACCGTTTTccgtg                                                                                                                                        <  1:1366597/150‑5 (MQ=255)
           aGGAAGGTTTCGAGGTCCACTTCCGGTAGCAGATTGATGCTGACATATTCCCCGGGTTGCCAGTAGCTTGAAATCTGCAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTATGTTCTCACGGAAAACGGTGCCGTTTTCAGCg                                                                                                                                        <  1:42910/150‑1 (MQ=255)
           aGGAAGGTTTCGAGGTCCACATCCGGTAGCAGATTGATGCTGACAAATTCCCCCGGTTGCCAGTAGCTTGAAATCTGCAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGCACCGTTTTccgtg                                                                                                                                        <  2:88297/150‑5 (MQ=255)
            ggAAGGTTTCGAGGTCCACATCCGGTAGCAGATTGATGCTTACAAATTCCCCCGGTTGCCAGTAGCTTGAAATCTGCAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCgg                                                                                                                                       <  1:1407416/150‑1 (MQ=255)
            ggAAGGTTTCGAGGTCCACATCCGGTAGCAGATTGATGCTGACAAATTCCCCCGGTTGCCAGTAGCTTGAAATCTGCAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCgg                                                                                                                                       <  1:356665/150‑1 (MQ=255)
            ggAAGGTTTCGAGGTCCACATCCGGTAGCAGATTGATGCTGACAAATTCCCCCGGTTGCCAGTAGCTTGAAATCTGCAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGCACCGTTTTccgtga                                                                                                                                       <  1:897120/150‑6 (MQ=255)
             gAAGGTTTCGAGGTCCACATCCGGTAGCAGATTGATGCTGACAAATTCCCCCGGTTGCCAGTAGCTTGAAATCTGCAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGt                                                                                                                                      >  1:375074/1‑150 (MQ=255)
                   ttCGAGGTCCACATCCGGTAGCAGATTGATGCTGACAAATTCCCCCGGTTGCCAGTAGCTTGAAATCTGCAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCAc                                                                                                                                >  1:491457/1‑150 (MQ=255)
                     cGAGGTCCACATCCGGTAGCAGATTGATGCTGACAAATTCCCCCGGTTGCCAGTAGCTTGAAATCTGCAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACgg                                                                                                                              <  1:973405/150‑1 (MQ=255)
                           ccACATCCGGTAGCAGATTGATGCTGACAAATTCCCCCGGTTGCCAGTAGCTTGAAATCTGCAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATTACGGAAGGCa                                                                                                                        <  1:318097/150‑1 (MQ=255)
                           ccACATCCGGTAGCAGATTGATGCTGACAAATTCCCCCGGTTGCCAGTAGCTTGAAATCTGCAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCa                                                                                                                        >  2:648029/1‑150 (MQ=255)
                           ccACATCCGGTAGCAGATTGATGATGACAAATTCCCCCGGTTGCCAGTAGCTTGAAATCTGCAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCa                                                                                                                        >  1:1087307/1‑150 (MQ=255)
                            aagaTCCGGTAGCAGATTGATGCTGACAAATTCCCCCGGTTGCCAGTTGCTTGAAATCTGCAACACCGGCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGTAAAACGGTGCAGTTTTGAGCGGGAATCACGGAAGGCAc                                                                                                                       <  2:1307228/147‑1 (MQ=255)
                                    gTAGCAGATTGATGCTGACAAATTCCCCCGGTTGCCAGTAGCTTGAAATCTGCAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACAGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCAcgc                                                                                                               <  2:5089/150‑1 (MQ=255)
                                        cAGATTGATGCTGACAAATTCCCCCGGTTGCCAGTAGCTTGAAATCTGCAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCAcgcccgc                                                                                                           >  2:403856/1‑150 (MQ=255)
                                         aGATTGATGCTGACAAATTCCCCCGGTTGCCAGTAGCTTGAAATCTGCAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCAcgcccgcc                                                                                                          <  1:897389/150‑1 (MQ=255)
                                          gATTGATGCTGACAAATTCCCCCGGTTGCCAGTAGCTTGAAATCTGCAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCa                                                                                                         <  1:771920/150‑1 (MQ=255)
                                             tGATGCTGACAAATTCCCCCGGTTGCCAGTAGCTTGAAATCTGCAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCa                                                                                                      <  1:206903/150‑1 (MQ=255)
                                               atgatgACAAATTCCCCCGGTTGCCAGTAGCTTGAAATCTGCAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGACCGCCAGCAcc                                                                                                    <  2:1087307/150‑1 (MQ=255)
                                                  cTGACAAATTCCCCCGGTTGCCAGTAGCTTGAAATCTGCAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTc                                                                                                                                            >  2:810175/1‑107 (MQ=255)
                                                  cTGACAAATTCCCCCGGTTGCCAGTAGCTTGAAATCTGCAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTc                                                                                                                                            <  1:810175/107‑1 (MQ=255)
                                                  cTGACAAATTCCCCCGGTTGCCAGTAGCTTGAAATCTGCAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCACCTGt                                                                                                 >  1:1592875/1‑150 (MQ=255)
                                                   tGACAAATTCCCCCGGTTGCCAGTAGCTTGAAATCTGCAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCACCTGTa                                                                                                >  1:1495265/1‑150 (MQ=255)
                                                      cAAATTCCCCCGGTTGCCAGTAGCTTGAAATCTGCAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCACCTGTAAct                                                                                             <  2:715631/150‑1 (MQ=255)
                                                        aaTTCCCCCGGTTGCCAGTAGCTTGAAATCTGCAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCACCTGTAActct                                                                                           >  2:923843/1‑150 (MQ=255)
                                                          ttCCCCCGGTTGCCAGTAGCTTGAAATCTGCAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCCCGCCCGCCAGCACCTGTAACTCTTc                                                                                         >  2:802513/1‑150 (MQ=255)
                                                             ccccGGTTGCCAGTAGCTTGAAATCTGCAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCACCTGTAACTCTTCGAg                                                                                      >  2:1437216/1‑150 (MQ=255)
                                                               ccGGTTGCCAGTAGCTTGAAATCTGCAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCACCTGTAACTCTTCGAGCa                                                                                    <  2:1461396/150‑1 (MQ=255)
                                                                 ggTTGCCAGTAGCTTGAAATCTGCAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCACCTGTAACTCTTCGAGCAAc                                                                                  <  1:703215/150‑1 (MQ=255)
                                                                 ggTTGCCAGTAGCTTGAAATCTGCAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCACCTGTAACTCTTCGAGCAAc                                                                                  <  2:1545431/150‑1 (MQ=255)
                                                                      ccAGTAGCTTGAAATCTGCAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCACCTGTAACTCTTCGAGCAACGGttt                                                                             <  2:375074/150‑1 (MQ=255)
                                                                      ccAGTAGCTTGAAATCTGCAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCACCTGTAACTCTTCGAGCAACGGttt                                                                             <  2:565678/150‑1 (MQ=255)
                                                                       gtgagaacttactcttcacccaccgcggcttgtctggaccggcggtgGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCACCTGTAACTCTTCGAGCAACGGTTTa                                                                            >  2:481131‑M2/48‑150 (MQ=255)
                                                                           gaacttactcttcacccaccgcggcttgtctggaccggcggtgGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCCCCGCCCCGCCCGCCAGCACCTGTAACTCTTCGAGCACCGGTTTATGCa                                                                        >  1:1314745‑M2/44‑150 (MQ=255)
                                                                               tGAAATCTGCAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCACCTGTAACTCTTCGAGCAACGGTTTATGCAGAGt                                                                    <  1:403856/150‑1 (MQ=255)
                                                                                 aaaTCTGCAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGCGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCACCTGTAACTCTTCGAGCAACGGTTTATGCAGAGTGa                                                                  <  2:1484840/150‑1 (MQ=255)
                                                                                  ctcttcacccaccgcggcttgtctggaccggcggtgGGTGAAGAGTAAGTTCTCCCGGAAACCGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCACCTGTAACTCTTCGAGCAACGGTTTATGCAGAGTGaa                                                                 >  2:836576‑M2/37‑150 (MQ=255)
                                                                                   tcttcacccaccgcggcttgtctggaccggcggtgGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCACCTGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAAt                                                                >  2:1619546‑M2/36‑150 (MQ=255)
                                                                                   aTCTGCAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCACCTGTAACTCTTCGAGCAACGGTTTATGCAGAGTg                                                                   >  1:817200/1‑147 (MQ=255)
                                                                                   aTCTGCAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCACCTGTAACTCTTCGAGCAACGGTTTATGCAGAGTg                                                                   <  2:817200/147‑1 (MQ=255)
                                                                                    tCTGCAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCACCTGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATg                                                               <  1:923843/150‑1 (MQ=255)
                                                                                     cTGCAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCACCTGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATgg                                                              >  2:1264719/1‑150 (MQ=255)
                                                                                      tGCAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCACCTGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGc                                                             >  2:718135/1‑150 (MQ=255)
                                                                                      tGCAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCACCTGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGc                                                             >  2:752963/1‑150 (MQ=255)
                                                                                       gCAACACCGCCGGTCCAGACAAGCCGAGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCACCTCTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCa                                                            >  2:209664/1‑150 (MQ=255)
                                                                                       cacccaccgcggcttgtctggaccggcggtgGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCCCGGAAGGCACCGCCACGCCCGCCAGCACCTGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGgcc                                                            >  1:1433188‑M2/32‑149 (MQ=255)
                                                                                        cAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCACCTGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCAc                                                           <  1:1047060/150‑1 (MQ=255)
                                                                                        acccaccgcggcttgtctggaccggcggtgGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCACCTGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCAc                                                           >  2:308127‑M2/31‑150 (MQ=255)
                                                                                         aaCACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCACCCGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCAcc                                                          <  2:546571/150‑1 (MQ=255)
                                                                                          agACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGACACCGCCACGCCCGCCAGCACCTGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCa                                                         >  1:1538407‑M1/3‑150 (MQ=255)
                                                                                               GCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCACCTGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGAccc                                                    <  2:622273‑M1/150‑1 (MQ=255)
                                                                                                 cgGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCACCTGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCcgc                                                  <  1:1133653/150‑1 (MQ=255)
                                                                                                      cAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCACCTGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCggg                                             <  2:230046/150‑1 (MQ=255)
                                                                                                      cAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCACCTGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCggg                                             <  2:366199/150‑1 (MQ=255)
                                                                                                      cAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCACCTGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCggg                                             <  1:342798/150‑1 (MQ=255)
                                                                                                       aGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACg                                                                                                                               >  2:59908/1‑67 (MQ=255)
                                                                                                       aGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACg                                                                                                                               <  1:59908/67‑1 (MQ=255)
                                                                                                        gACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCACCTGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGCCCCGCGCGAGTc                                           >  1:208980/1‑150 (MQ=255)
                                                                                                           aaGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCACCTGTAACTCTTCGAGCAACGGTTTATGCAGAGTGGATGGCACCAGACCCGCGCGGGTCGgc                                        <  1:305022/150‑1 (MQ=255)
                                                                                                           aaGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCACCTGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGgc                                        >  2:641470/1‑150 (MQ=255)
                                                                                                           aaGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCACCTGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGgc                                        <  1:648029/150‑1 (MQ=255)
                                                                                                            aGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCACCTGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGgca                                       <  2:1592875/150‑1 (MQ=255)
                                                                                                             gCCGCGGTGGGGGAAGATTAAGTCCTCACGGAAACCGGTGCCGCGTTCAGCGGCACTCACTGAAGGCACCGCCACGCCCGCCACCACCTGTACCTCAGCGCGCACCGGTCATTGCAGACCCCCTGGCACCAGACCCGCGCGGGTCGgcag                                      <  2:730371/150‑1 (MQ=255)
                                                                                                                    tgGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCACCTGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGGCAGCACGTTg                               <  1:1591379‑M2/148‑1 (MQ=255)
                                                                                                                    tgGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGAACCTGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGGCAGCACGTTg                               <  2:1227410‑M2/148‑1 (MQ=255)
                                                                                                                     gGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCACCTGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGGCAGCACGTTGa                              >  2:642201‑M2/2‑150 (MQ=255)
                                                                                                                        tgAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGATGGCACCGCCACGCCCGCCAGCACCTGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGGCAGCACGTTTAGGc                           <  1:1626128/150‑1 (MQ=255)
                                                                                                                        tgAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCACCTGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGGCAGCACGTTGAGGc                           <  1:638276/150‑1 (MQ=255)
                                                                                                                        tgAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCACCTGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGGCAGCACGTTGAGGc                           <  2:975408/150‑1 (MQ=255)
                                                                                                                               tAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCACCTGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGGCAGCACGTTGAGGCCAAATTg                    >  2:1269878/1‑150 (MQ=255)
                                                                                                                                aaGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCACCTGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGGCAGCACGTTGAGGCCAAATTGt                   >  2:1355220/1‑150 (MQ=255)
                                                                                                                                  gTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCACCTGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGGCAGCACGTTGAGGCCAAATTGTTc                 <  1:832203/150‑1 (MQ=255)
                                                                                                                                      tcACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCACCTGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGGCAGCACGTTGAGGCCAAATTGTTCGGCa             >  2:535531/1‑150 (MQ=255)
                                                                                                                                       cgctgAAAACGGCACCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCACCTGTAACTCTTCGAGCAACGGATTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGGCAGCACGTTGAGGCCAAATTGTTCGGCaa            >  2:212046/5‑150 (MQ=255)
                                                                                                                                       cACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCACCTGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGGCAGCACGTTGAGGCCAAATTGTTCGGCaa            >  1:1303399/1‑150 (MQ=255)
                                                                                                                                        gctgAAAACGGCACCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCACCTGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGGCAGCACGTTGAGGCCAAATTGTTCGGCAAt           >  1:1108088/4‑150 (MQ=255)
                                                                                                                                         cGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCACCTGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGGCAGCACGTTGAGGCCAAATTGTTCGGCAATc          >  2:1575748/1‑150 (MQ=255)
                                                                                                                                          ggAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCACCTGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGGCAGCACGTTGAGGCCAAATTGTTCGGCAATCt         <  1:466677/150‑1 (MQ=255)
                                                                                                                                                 ggTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCACCTGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGGCAGCACGTTGAGGCCAAATTGTTCGGCAATCTTATAAcc  <  1:91737/150‑1 (MQ=255)
                                                                                                                                                 ggTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCACCTGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGGCAGCACGTTGAGGCCAAATTGTTCGGCAATCTTATAAcc  <  1:1264719/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCTGCTCATTCAGGAAGGTTTCGAGGTCCACATCCGGTAGCAGATTGATGCTGACAAATTCCCCCGGTTGCCAGTAGCTTGAAATCTGCAACACCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAAGGCACCGCCACGCCCGCCAGCACCTGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGGCAGCACGTTGAGGCCAAATTGTTCGGCAATCTTATAACC  >  NC_000913/3636600‑3636894

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: