Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,544,693 C→G 5.1% G51A (GGG→GCG) ‡ bioH ← pimeloyl‑ACP methyl ester carboxylesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,544,6930CG 5.1% 503.6 / 14.0 156G51A (GGG→GCG) ‡bioHpimeloyl‑ACP methyl ester carboxylesterase
Reads supporting (aligned to +/- strand):  ref base C (76/71);  new base G (5/3);  total (81/75)
Fisher's exact test for biased strand distribution p-value = 7.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.29e-01

Only 100 of 159 total aligned reads displayed.

CCTGAACACGCTCGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATT  >  NC_000913/3544544‑3544841
                                                                                                                                                     |                                                                                                                                                    
ccTGAACACGCTCGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGcc                                                                                                                                                      <  2:1183252/150‑1 (MQ=255)
  tGAACACGCTCGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCg                                                                                                                                                    <  1:52383/150‑1 (MQ=255)
   gAACACGCTCGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCCCCAAATCCCCGGCTACGCCCGa                                                                                                                                                   >  1:475079/1‑150 (MQ=255)
          cTCGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTTCGGGCGTAGCCggg                                                                                                                                            <  2:944834/150‑1 (MQ=255)
           tCGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGc                                                                                                                                           <  2:975654/150‑1 (MQ=255)
             gggATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAg                                                                                                                                         >  2:1544733/1‑150 (MQ=255)
             gggATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAg                                                                                                                                         >  2:1529479/1‑150 (MQ=255)
                aggggTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATAGCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTc                                                                                                                                      <  1:320312/148‑1 (MQ=255)
                   ggTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGTGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTcaac                                                                                                                                   >  2:1473832/1‑150 (MQ=255)
                   ggTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTcaac                                                                                                                                   >  2:1305433/1‑150 (MQ=255)
                    gTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTcaaca                                                                                                                                  <  2:1515973/150‑1 (MQ=255)
                      tAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAg                                                                                                                                >  2:99272/1‑150 (MQ=255)
                       aaCGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAgg                                                                                                                               >  1:55179/1‑150 (MQ=255)
                        aCGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGt                                                                                                                              <  1:899263/150‑1 (MQ=255)
                          ggtaTCTGGCTTGCCACCTGACCGCGCAGACTCCAGCCTAACCAATTTGCTTTATCAGGTGCCTGTTGCAGCAGGGATTCGGACATATCAGCAAGTGACAGCGCACCTAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACTAGGTGc                                                                                                                            <  2:1621697/147‑1 (MQ=255)
                           cAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCACATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCa                                                                                                                           >  2:1165530/1‑150 (MQ=255)
                            aaTCTGGCTTGCCACCAGCCCGACCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAg                                                                                                                          <  2:311674/150‑1 (MQ=255)
                            aaTCTGGCATGCCACCAGCCCGCCCAGACTCCAGCCTAAACAAATGGCTTTATCAGGTGCCTGTTGCAGCACTTCTTCGGCCAAATCAGCAAGTGACAGCGCCCCAAATCCCCGGCTACGCCCGCAGCCGGGCAGGTCAGCAAGGTGCac                                                                                                                          >  1:89343/1‑149 (MQ=255)
                                tGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTa                                                                                                                      >  2:1421628/1‑150 (MQ=255)
                                 ggCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTTCGGGCGTAGCCgggg                                                                                                                                           >  1:1554231/1‑127 (MQ=255)
                                 ggCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTTCGGGCGTAGCCgggg                                                                                                                                           <  2:1554231/128‑2 (MQ=255)
                                     tGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATg                                                                                                                 <  2:36839/150‑1 (MQ=255)
                                     tGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATg                                                                                                                 <  1:1313127/150‑1 (MQ=255)
                                      gCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGc                                                                                                                <  1:364810/150‑1 (MQ=255)
                                       ccaccaGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCg                                                                                                               >  2:164432/1‑150 (MQ=255)
                                       ccaccaGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCg                                                                                                               >  2:937176/1‑150 (MQ=255)
                                       ccaccaGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCg                                                                                                               <  2:740594/150‑1 (MQ=255)
                                       ccaccaGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCg                                                                                                               <  1:130355/150‑1 (MQ=255)
                                         accaGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAg                                                                                                             <  1:218729/150‑1 (MQ=255)
                                             gcccgcccAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTaa                                                                                                         >  1:1597866/1‑150 (MQ=255)
                                             gcccgcccAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCATCGGCCATATCAGCAAGTGACAGCGCACCAAATCGCCGGCTACGCCCGACGCCGGGCAGGTCCACAAGGTGCAGCGTACAATGCGAGCTaa                                                                                                         >  1:263159/1‑150 (MQ=255)
                                                cgcccAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGCCCGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGtt                                                                                                      >  2:1252083/1‑150 (MQ=255)
                                                  cccAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTcc                                                                                                    <  2:181545/150‑1 (MQ=255)
                                                        cTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCa                                                                                              >  2:1258233/1‑150 (MQ=255)
                                                            cgCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGc                                                                                          <  1:1413341/149‑1 (MQ=255)
                                                            aGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGc                                                                                          <  2:957657/150‑1 (MQ=255)
                                                            aGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGc                                                                                          >  1:1145195/1‑150 (MQ=255)
                                                             gCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCa                                                                                         >  1:1679079/1‑150 (MQ=255)
                                                             gCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCa                                                                                         <  1:1410951/150‑1 (MQ=255)
                                                              ccTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCaa                                                                                        <  1:1529479/150‑1 (MQ=255)
                                                                tAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACg                                                                                      <  1:164432/150‑1 (MQ=255)
                                                                     aaaTGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCcaca                                                                                 <  1:828559/150‑1 (MQ=255)
                                                                     aaaTGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTACAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCcaca                                                                                 >  1:514284/1‑150 (MQ=255)
                                                                           cTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTctg                                                                           <  2:1210495/150‑3 (MQ=255)
                                                                           cTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCgg                                                                           <  2:1597866/150‑1 (MQ=255)
                                                                            tttATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGACGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGc                                                                          >  1:104161/1‑150 (MQ=255)
                                                                            tttATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGc                                                                          <  1:194444/150‑1 (MQ=255)
                                                                              tATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCAt                                                                        >  1:165920/1‑150 (MQ=255)
                                                                              tATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCAt                                                                        >  1:952240/1‑150 (MQ=255)
                                                                              tATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCAt                                                                        >  1:1103947/1‑150 (MQ=255)
                                                                                tCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCATGGGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTc                                                                      <  1:437031/150‑1 (MQ=255)
                                                                                tCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTc                                                                      <  2:266932/150‑1 (MQ=255)
                                                                                tCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTc                                                                      <  1:1609071/150‑1 (MQ=255)
                                                                                   ggTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGt                                                                   <  2:570196/150‑1 (MQ=255)
                                                                                    gTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTc                                                                  <  1:660712/150‑1 (MQ=255)
                                                                                      gCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCCCACTTCGGCATTCAGTccc                                                                >  2:908165/1‑150 (MQ=255)
                                                                                      gCCTGTTGCAGCACGGCTTCGCCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCATCGCCACACTTCGGCATTCAGTccc                                                                >  2:1483980/1‑150 (MQ=255)
                                                                                        cTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGTGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCa                                                              <  1:940253/150‑1 (MQ=255)
                                                                                            tGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCg                                                          <  1:735020/150‑1 (MQ=255)
                                                                                              cagcaCGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTg                                                        <  1:1544733/150‑1 (MQ=255)
                                                                                               agcaCGGATTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCACGGATACGCCCGAAGCGGGGCAGGTCAACAAGGAGCAGAGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGCCATTCAGTCCCCATCCGTgc                                                       >  1:1071762/1‑150 (MQ=255)
                                                                                                     gCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCCCACTTCGGCATTCAGTCCCCATCCGTgcagcagc                                                 >  1:539382/1‑150 (MQ=255)
                                                                                                          ggCCATATCAGCAAGGGACAGCGCACCAAATCCCCGGCGACGCCCGAAGCCGTGCAGGTCAACAAGGTGCAGCGTAATATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAg                                            <  2:838874/150‑1 (MQ=255)
                                                                                                           gCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGa                                           >  2:1119426/1‑150 (MQ=255)
                                                                                                              atatCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGa                                        >  2:566376/1‑150 (MQ=255)
                                                                                                               tatCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCTGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGGAGCACAAGATgca                                       >  1:334218/1‑148 (MQ=255)
                                                                                                               tatCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAGAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATc                                                            >  2:785267/1‑129 (MQ=255)
                                                                                                               tatCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAGAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATc                                                            <  1:785267/129‑1 (MQ=255)
                                                                                                                atCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAAc                                      <  2:475079/150‑1 (MQ=255)
                                                                                                                  cAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGAGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACAt                                    >  1:106093/1‑150 (MQ=255)
                                                                                                                   aGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACAtt                                   <  1:99272/150‑1 (MQ=255)
                                                                                                                    gCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTc                                  >  1:1557366/1‑150 (MQ=255)
                                                                                                                    gCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGGGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTc                                  <  2:106093/150‑1 (MQ=255)
                                                                                                                      aaGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTccc                                <  2:419297/150‑1 (MQ=255)
                                                                                                                       agtgagtgCGCACCATATCCCCGGCTACGCACGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCGGTCCCCATCCGTGCAGCAGCACAAGATGAACATTcccc                               <  2:89343/150‑1 (MQ=255)
                                                                                                                       aGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTcccc                               <  1:224957/150‑1 (MQ=255)
                                                                                                                       aGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTcccc                               <  2:615271/150‑1 (MQ=255)
                                                                                                                       aGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTcccc                               <  1:1305433/150‑1 (MQ=255)
                                                                                                                       aGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTcccc                               >  2:549559/1‑150 (MQ=255)
                                                                                                                         tGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTg                             <  1:988011/150‑1 (MQ=255)
                                                                                                                          gACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGGTAACTGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGCACATTCCCCTgc                            >  2:299083/1‑149 (MQ=255)
                                                                                                                            cAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGAcc                          <  2:165920/150‑1 (MQ=255)
                                                                                                                            cAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGAcc                          >  1:1447969/1‑150 (MQ=255)
                                                                                                                            aaGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGAcc                          <  1:1037444/149‑1 (MQ=255)
                                                                                                                               cgcACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCTTCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCttt                       <  1:1115647/150‑1 (MQ=255)
                                                                                                                                gcACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTg                      <  1:937176/150‑1 (MQ=255)
                                                                                                                                 cACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTgg                     >  2:882591/1‑150 (MQ=255)
                                                                                                                                 cACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGACCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTgg                     >  2:653773/1‑150 (MQ=255)
                                                                                                                                   ccAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTc                   <  2:1393281/150‑1 (MQ=255)
                                                                                                                                     aaaTCCCCGGCTATGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAAAATTCCCCTGACCTTTGGTCTg                 >  2:1087469/1‑150 (MQ=255)
                                                                                                                                       ctgcccGGCTTCGGGCGTAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGcc               >  2:15656/4‑150 (MQ=255)
                                                                                                                                       aTCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGcc               <  1:1473832/150‑1 (MQ=255)
                                                                                                                                        tgcccGGCTTCGGGCGTAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGcca              >  1:1188484/3‑150 (MQ=255)
                                                                                                                                           ccGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGccacca           >  1:865616/1‑150 (MQ=255)
                                                                                                                                             ggCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGa         >  2:1616100/1‑150 (MQ=255)
                                                                                                                                             ggCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGa         >  2:1149478/1‑150 (MQ=255)
                                                                                                                                                  cGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTGAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTa    >  1:150974/1‑150 (MQ=255)
                                                                                                                                                  cGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTa    <  1:1258233/150‑1 (MQ=255)
                                                                                                                                                   gCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTAt   >  1:737262/1‑150 (MQ=255)
                                                                                                                                                    cccGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTAtt  <  1:503276/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCTGAACACGCTCGGGATGGGTTAACGCAATCTGGCTTGCCACCAGCCCGCCCAGACTCCAGCCTAACCAAATGGCTTTATCAGGTGCCTGTTGCAGCACGGCTTCGGCCATATCAGCAAGTGACAGCGCACCAAATCCCCGGCTACGCCCGAAGCCGGGCAGGTCAACAAGGTGCAGCGTAAAATGCGAGCTAAGTTCCTCGTCAATGCAACGCCACACTTCGGCATTCAGTCCCCATCCGTGCAGCAGCACAAGATGAACATTCCCCTGACCTTTGGTCTGCCACCAGATGTTATT  >  NC_000913/3544544‑3544841

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